source: LMDZ.3.3/trunk/libf/dyn3d/rotatst.F @ 346

Last change on this file since 346 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.1 KB
Line 
1      SUBROUTINE rotatst (klevel,x, y, rot )
2c
3c  P. Le Van
4c
5c    *****************************************************************
6c     .. calcule le rotationnel a tous les niveaux d'1 vecteur de comp. x et y ..
7c         x  et  y etant des composantes  covariantes  .....
8c    *****************************************************************
9c        x  et y     sont des arguments d'entree pour le s-prog
10c        rot          est  un argument  de sortie pour le s-prog
11c
12      IMPLICIT NONE
13c
14      INTEGER klevel
15#include "dimensions.h"
16#include "paramet.h"
17
18      REAL rot( ip1jm,klevel )
19      REAL x( ip1jmp1,klevel ), y( ip1jm,klevel )
20      INTEGER  l, ij
21c
22c
23      DO 5 l = 1,klevel
24c
25      DO 1 ij = 1, ip1jm - 1
26      rot( ij,l )  =  (  y( ij+1 , l )  -  y( ij,l )   +
27     *                 x(ij +iip1, l )  -  x( ij,l )  )
28   1  CONTINUE
29c
30c    .... correction pour rot( iip1,j,l)  ....
31c
32c    ....   rot(iip1,j,l)= rot(1,j,l) ...
33CDIR$ IVDEP
34      DO 2 ij = iip1, ip1jm, iip1
35      rot( ij,l ) = rot( ij -iim,l )
36   2  CONTINUE
37c
38   5  CONTINUE
39      RETURN
40      END
Note: See TracBrowser for help on using the repository browser.