LMDZ
gradiv.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE gradiv(klevel, xcov, ycov, ld, gdx, gdy )
5 c
6 c Auteur : P. Le Van
7 c
8 c ***************************************************************
9 c
10 c ld
11 c calcul de (grad (div) ) du vect. v ....
12 c
13 c xcov et ycov etant les composant.covariantes de v
14 c ****************************************************************
15 c xcov , ycov et ld sont des arguments d'entree pour le s-prog
16 c gdx et gdy sont des arguments de sortie pour le s-prog
17 c
18 c
19  IMPLICIT NONE
20 c
21 #include "dimensions.h"
22 #include "paramet.h"
23 #include "comdissipn.h"
24 #include "logic.h"
25 
26  INTEGER klevel
27 c
28  REAL xcov( ip1jmp1,klevel ), ycov( ip1jm,klevel )
29  REAL gdx( ip1jmp1,klevel ), gdy( ip1jm,klevel )
30 
31  REAL div(ip1jmp1,llm)
32 
33  INTEGER l,ij,iter,ld
34 c
35 c
36 c
37  CALL scopy( ip1jmp1*klevel,xcov,1,gdx,1 )
38  CALL scopy( ip1jm*klevel, ycov,1,gdy,1 )
39 c
40  DO 10 iter = 1,ld
41 c
42  CALL diverg( klevel, gdx , gdy, div )
43  CALL filtreg( div, jjp1, klevel, 2,1, .true.,2 )
44  CALL grad( klevel, div, gdx, gdy )
45 c
46  DO 5 l = 1, klevel
47  DO 3 ij = 1, ip1jmp1
48  gdx( ij,l ) = - gdx( ij,l ) * cdivu
49  3 CONTINUE
50  DO 4 ij = 1, ip1jm
51  gdy( ij,l ) = - gdy( ij,l ) * cdivu
52  4 CONTINUE
53  5 CONTINUE
54 c
55  10 CONTINUE
56  RETURN
57  END
grad
subroutine grad(klevel, pg, pgx, pgy)
Definition: grad.F:5
ip1jmp1
!$Header llmm1 INTEGER ip1jmp1
Definition: paramet.h:14
cdivu
!$Header cdivu
Definition: comdissipn.h:11
llm
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd w_profd omega_profd!do llm!print llm l llm
Definition: 1D_interp_cases.h:103
diverg
subroutine diverg(klevel, x, y, div)
Definition: diverg.F:5
scopy
subroutine scopy(n, sx, incx, sy, incy)
Definition: cray.F:9
ip1jm
!$Header llmm1 INTEGER ip1jm
Definition: paramet.h:14
jjp1
!$Header jjp1
Definition: paramet.h:14
true
!$Id itapm1 ENDIF!IM on interpole les champs sur les niveaux STD de pression!IM a chaque pas de temps de la physique c!positionnement de l argument logique a false c!pour ne pas recalculer deux fois la meme chose!c!a cet effet un appel a plevel_new a ete deplace c!a la fin de la serie d appels c!la boucle DO nlevSTD a ete internalisee c!dans d ou la creation de cette routine c c!CALL ulevSTD CALL &zphi philevSTD CALL &zx_rh rhlevSTD!DO klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon du jour ou toutes les read_climoz CALL true
Definition: calcul_STDlev.h:139
gradiv
subroutine gradiv(klevel, xcov, ycov, ld, gdx, gdy)
Definition: gradiv.F:5
filtreg
subroutine filtreg(champ, nlat, nbniv, ifiltre, iaire, griscal, iter)
Definition: filtreg.F:6