Changeset 5246 for LMDZ6/trunk/libf/dyn3d_common/rotat_nfil.f90
- Timestamp:
- Oct 21, 2024, 2:58:45 PM (23 hours ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d_common/rotat_nfil.f90
r5245 r5246 2 2 ! $Header$ 3 3 ! 4 5 c 6 c Auteur : P.Le Van 7 c**************************************************************8 c. Calcule le rotationnel non filtre ,9 ca tous les niveaux d'1 vecteur de comp. x et y ..10 cx et y etant des composantes covariantes ...11 c********************************************************************12 cklevel, x et y sont des arguments d'entree pour le s-prog13 crot est un argument de sortie pour le s-prog14 c 15 16 c 17 18 19 20 c 21 c..... variables en arguments ......22 c 23 INTEGERklevel24 REALrot( ip1jm,klevel )25 REALx( ip1jmp1,klevel ), y( ip1jm,klevel )26 c 27 c... variables locales ...28 c 29 INTEGERl, ij30 c 31 c 32 DO 10l = 1,klevel33 c 34 35 rot( ij,l ) = y( ij+1 , l ) - y( ij,l ) +36 * x(ij +iip1, l ) - x( ij,l )37 38 c 39 c.... correction pour rot( iip1,j,l) ....40 c.... rot(iip1,j,l)= rot(1,j,l) ...41 CDIR$ IVDEP42 43 44 45 c 46 10 CONTINUE4 SUBROUTINE rotat_nfil (klevel, x, y, rot ) 5 ! 6 ! Auteur : P.Le Van 7 !************************************************************** 8 !. Calcule le rotationnel non filtre , 9 ! a tous les niveaux d'1 vecteur de comp. x et y .. 10 ! x et y etant des composantes covariantes ... 11 !******************************************************************** 12 ! klevel, x et y sont des arguments d'entree pour le s-prog 13 ! rot est un argument de sortie pour le s-prog 14 ! 15 IMPLICIT NONE 16 ! 17 INCLUDE "dimensions.h" 18 INCLUDE "paramet.h" 19 INCLUDE "comgeom.h" 20 ! 21 ! ..... variables en arguments ...... 22 ! 23 INTEGER :: klevel 24 REAL :: rot( ip1jm,klevel ) 25 REAL :: x( ip1jmp1,klevel ), y( ip1jm,klevel ) 26 ! 27 ! ... variables locales ... 28 ! 29 INTEGER :: l, ij 30 ! 31 ! 32 DO l = 1,klevel 33 ! 34 DO ij = 1, ip1jm - 1 35 rot( ij,l ) = y( ij+1 , l ) - y( ij,l ) + & 36 x(ij +iip1, l ) - x( ij,l ) 37 ENDDO 38 ! 39 ! .... correction pour rot( iip1,j,l) .... 40 ! .... rot(iip1,j,l)= rot(1,j,l) ... 41 !DIR$ IVDEP 42 DO ij = iip1, ip1jm, iip1 43 rot( ij,l ) = rot( ij -iim,l ) 44 ENDDO 45 ! 46 END DO 47 47 48 49 END 48 RETURN 49 END SUBROUTINE rotat_nfil
Note: See TracChangeset
for help on using the changeset viewer.