source: LMDZ.3.3/branches/rel-LF/libf/dyn3d/grad.F @ 363

Last change on this file since 363 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.0 KB
Line 
1      SUBROUTINE  grad(klevel, pg,pgx,pgy )
2c
3c      P. Le Van
4c
5c    ******************************************************************
6c     .. calcul des composantes covariantes en x et y du gradient de g
7c
8c    ******************************************************************
9c             pg        est un   argument  d'entree pour le s-prog
10c       pgx  et  pgy    sont des arguments de sortie pour le s-prog
11c
12      IMPLICIT NONE
13c
14#include "dimensions.h"
15#include "paramet.h"
16      INTEGER klevel
17      REAL  pg( ip1jmp1,klevel )
18      REAL pgx( ip1jmp1,klevel ) , pgy( ip1jm,klevel )
19      INTEGER  l,ij
20c
21c
22      DO 6 l = 1,klevel
23c
24      DO 2  ij = 1, ip1jmp1 - 1
25      pgx( ij,l ) = pg( ij +1,l ) - pg( ij,l )
26   2  CONTINUE
27c
28c    .... correction pour  pgx(ip1,j,l)  ....
29c    ...    pgx(iip1,j,l)= pgx(1,j,l)  ....
30CDIR$ IVDEP
31      DO 3  ij = iip1, ip1jmp1, iip1
32      pgx( ij,l ) = pgx( ij -iim,l )
33   3  CONTINUE
34c
35      DO 4 ij = 1,ip1jm
36      pgy( ij,l ) = pg( ij,l ) - pg( ij +iip1,l )
37   4  CONTINUE
38c
39   6  CONTINUE
40      RETURN
41      END
Note: See TracBrowser for help on using the repository browser.