My Project
 All Classes Files Functions Variables Macros
grad.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE grad(klevel, pg,pgx,pgy )
5 c
6 c P. Le Van
7 c
8 c ******************************************************************
9 c .. calcul des composantes covariantes en x et y du gradient de g
10 c
11 c ******************************************************************
12 c pg est un argument d'entree pour le s-prog
13 c pgx et pgy sont des arguments de sortie pour le s-prog
14 c
15  IMPLICIT NONE
16 c
17 #include "dimensions.h"
18 #include "paramet.h"
19  INTEGER klevel
20  REAL pg( ip1jmp1,klevel )
21  REAL pgx( ip1jmp1,klevel ) , pgy( ip1jm,klevel )
22  INTEGER l,ij
23 c
24 c
25  DO 6 l = 1,klevel
26 c
27  DO 2 ij = 1, ip1jmp1 - 1
28  pgx( ij,l ) = pg( ij +1,l ) - pg( ij,l )
29  2 CONTINUE
30 c
31 c .... correction pour pgx(ip1,j,l) ....
32 c ... pgx(iip1,j,l)= pgx(1,j,l) ....
33 CDIR$ IVDEP
34  DO 3 ij = iip1, ip1jmp1, iip1
35  pgx( ij,l ) = pgx( ij -iim,l )
36  3 CONTINUE
37 c
38  DO 4 ij = 1,ip1jm
39  pgy( ij,l ) = pg( ij,l ) - pg( ij +iip1,l )
40  4 CONTINUE
41 c
42  6 CONTINUE
43  RETURN
44  END