LMDZ
nxgraro2.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE nxgraro2 (klevel,xcov, ycov, lr, grx, gry )
5 c
6 c p.le van .
7 c ***********************************************************
8 c lr
9 c calcul de( nxgrad(rot) ) du vect. v ....
10 c
11 c xcov et ycov etant les compos. covariantes de v
12 c ***********************************************************
13 c xcov , ycov et lr sont des arguments d'entree pour le s-prog
14 c grx et gry sont des arguments de sortie pour le s-prog
15 c
16 c
17  IMPLICIT NONE
18 c
19 #include "dimensions.h"
20 #include "paramet.h"
21 #include "comdissipn.h"
22 c
23 c ...... variables en arguments .......
24 c
25  INTEGER klevel
26  REAL xcov( ip1jmp1,klevel ), ycov( ip1jm,klevel )
27  REAL grx( ip1jmp1,klevel ), gry( ip1jm,klevel )
28 c
29 c ...... variables locales ........
30 c
31  REAL rot(ip1jm,llm) , signe, nugradrs
32  INTEGER l,ij,iter,lr
33 c ........................................................
34 c
35 c
36 c
37  signe = (-1.)**lr
38  nugradrs = signe * crot
39 c
40  CALL SCOPY ( ip1jmp1* klevel, xcov, 1, grx, 1 )
41  CALL SCOPY ( ip1jm * klevel, ycov, 1, gry, 1 )
42 c
43  CALL rotatf ( klevel, grx, gry, rot )
44 c
45  CALL laplacien_rot ( klevel, rot, rot,grx,gry )
46 
47 c
48 c ..... Iteration de l'operateur laplacien_rotgam .....
49 c
50  DO iter = 1, lr -2
51  CALL laplacien_rotgam ( klevel, rot, rot )
52  ENDDO
53 c
54 c
55  CALL filtreg( rot, jjm, klevel, 2,1, .false.,1)
56  CALL nxgrad ( klevel, rot, grx, gry )
57 c
58  DO l = 1, klevel
59  DO ij = 1, ip1jm
60  gry( ij,l ) = gry( ij,l ) * nugradrs
61  ENDDO
62  DO ij = 1, ip1jmp1
63  grx( ij,l ) = grx( ij,l ) * nugradrs
64  ENDDO
65  ENDDO
66 c
67  RETURN
68  END
ip1jmp1
!$Header llmm1 INTEGER ip1jmp1
Definition: paramet.h:14
l
c c $Id c c calculs statistiques distribution nuage ftion du regime dynamique c c Ce calcul doit etre fait a partir de valeurs mensuelles CALL nbregdyn DO kmaxm1 DO l
Definition: calcul_REGDYN.h:13
nxgrad
subroutine nxgrad(klevel, rot, x, y)
Definition: nxgrad.F:5
laplacien_rotgam
subroutine laplacien_rotgam(klevel, rotin, rotout)
Definition: laplacien_rotgam.F:5
nxgraro2
subroutine nxgraro2(klevel, xcov, ycov, lr, grx, gry)
Definition: nxgraro2.F:5
ip1jm
!$Header llmm1 INTEGER ip1jm
Definition: paramet.h:14
false
!$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 false
Definition: calcul_STDlev.h:26
ij
!$Id klon klev DO klon!IM klev DO klon klon nbp_lat DO nbp_lon ij
Definition: write_histday_seri.h:96
filtreg
subroutine filtreg(champ, nlat, nbniv, ifiltre, iaire, griscal, iter)
Definition: filtreg.F:6
de
do llm!au dessus de
Definition: 1D_nudge_sandu_astex.h:1