My Project
 All Classes Files Functions Variables Macros
nxgrad_gam.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE nxgrad_gam( klevel, rot, x, y )
5 c
6 c P. Le Van
7 c
8 c ********************************************************************
9 c calcul du gradient tourne de pi/2 du rotationnel du vect.v
10 c ********************************************************************
11 c rot est un argument d'entree pour le s-prog
12 c x et y sont des arguments de sortie pour le s-prog
13 c
14  IMPLICIT NONE
15 c
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
22 c
23  DO 10 l = 1,klevel
24 c
25  DO 1 ij = 2, ip1jm
26  y( ij,l ) = (rot( ij,l ) - rot( ij-1,l )) * cvscuvgam( ij )
27  1 CONTINUE
28 c
29 c ..... correction pour y ( 1,j,l ) ......
30 c
31 c .... y(1,j,l)= y(iip1,j,l) ....
32 CDIR$ IVDEP
33  DO 2 ij = 1, ip1jm, iip1
34  y( ij,l ) = y( ij +iim,l )
35  2 CONTINUE
36 c
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
44 c
45  10 CONTINUE
46  RETURN
47  END