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