source: trunk/LMDZ.COMMON/libf/dyn3d_common/nxgrad_gam.F @ 3553

Last change on this file since 3553 was 1300, checked in by emillour, 11 years ago

Common dynamics:
Some updates to keep up with LMDZ5 Earth model evolution (up to LMDZ5 rev 1955).
Main change is the introduction of a "dyn3d_common" directory
to store files common to dyn3d and dyn3dpar.
See file "DOC/chantiers/commit_importants.log" for detailed list
of changes. These changes do not change results on test cases.
EM

File size: 1.2 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE nxgrad_gam( klevel, rot, x, y )
5c
6c  P. Le Van
7c
8c   ********************************************************************
9c      calcul du gradient tourne de pi/2 du rotationnel du vect.v
10c   ********************************************************************
11c       rot          est un argument  d'entree pour le s-prog
12c       x  et y    sont des arguments de sortie pour le s-prog
13c
14      IMPLICIT NONE
15c
16#include "dimensions.h"
17#include "paramet.h"
18#include "comgeom.h"
19      INTEGER klevel
20      REAL rot( ip1jm,klevel ),x( ip1jmp1,klevel ),y(ip1jm,klevel )
21      INTEGER   l,ij
22c
23      DO 10 l = 1,klevel
24c
25      DO 1  ij = 2, ip1jm
26      y( ij,l ) = (rot( ij,l ) - rot( ij-1,l )) * cvscuvgam( ij )
27   1  CONTINUE
28c
29c    ..... correction pour  y ( 1,j,l )  ......
30c
31c    ....    y(1,j,l)= y(iip1,j,l) ....
32CDIR$ IVDEP
33      DO 2  ij = 1, ip1jm, iip1
34      y( ij,l ) = y( ij +iim,l )
35   2  CONTINUE
36c
37      DO 4  ij = iip2,ip1jm
38      x( ij,l ) = (rot( ij,l ) - rot( ij -iip1,l )) * cuscvugam( ij )
39   4  CONTINUE
40      DO 6 ij = 1,iip1
41      x(    ij    ,l ) = 0.
42      x( ij +ip1jm,l ) = 0.
43   6  CONTINUE
44c
45  10  CONTINUE
46      RETURN
47      END
Note: See TracBrowser for help on using the repository browser.