LMDZ
divgrad.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE divgrad (klevel,h, lh, divgra )
5  IMPLICIT NONE
6 c
7 c=======================================================================
8 c
9 c Auteur : P. Le Van
10 c ----------
11 c
12 c lh
13 c calcul de (div( grad )) de h .....
14 c h et lh sont des arguments d'entree pour le s-prog
15 c divgra est un argument de sortie pour le s-prog
16 c
17 c=======================================================================
18 c
19 c declarations:
20 c -------------
21 c
22 #include "dimensions.h"
23 #include "paramet.h"
24 #include "comgeom.h"
25 #include "comdissipn.h"
26 #include "logic.h"
27 c
28  INTEGER klevel
29  REAL h( ip1jmp1,klevel ), divgra( ip1jmp1,klevel )
30 c
31  REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)
32 
33  INTEGER l,ij,iter,lh
34 c
35 c
36 c
37  CALL scopy ( ip1jmp1*klevel,h,1,divgra,1 )
38 c
39  DO 10 iter = 1,lh
40 
41  CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1 )
42 
43  CALL grad (klevel,divgra, ghx , ghy )
44  CALL diverg (klevel, ghx , ghy , divgra )
45 
46  CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1)
47 
48  DO 5 l = 1,klevel
49  DO 4 ij = 1, ip1jmp1
50  divgra( ij,l ) = - cdivh * divgra( ij,l )
51  4 CONTINUE
52  5 CONTINUE
53 c
54  10 CONTINUE
55  RETURN
56  END
grad
subroutine grad(klevel, pg, pgx, pgy)
Definition: grad.F:5
ip1jmp1
!$Header llmm1 INTEGER ip1jmp1
Definition: paramet.h:14
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
divgrad
subroutine divgrad(klevel, h, lh, divgra)
Definition: divgrad.F:5
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
filtreg
subroutine filtreg(champ, nlat, nbniv, ifiltre, iaire, griscal, iter)
Definition: filtreg.F:6