- Timestamp:
- Aug 2, 2024, 9:58:25 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotat_nfil.f90
r5136 r5159 3 3 4 4 SUBROUTINE rotat_nfil(klevel, x, y, rot ) 5 ! 5 6 6 ! Auteur : P.Le Van 7 7 !************************************************************** … … 12 12 ! klevel, x et y sont des arguments d'entree pour le s-prog 13 13 ! rot est un argument de sortie pour le s-prog 14 ! 14 15 15 USE lmdz_comgeom 16 16 17 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 18 USE lmdz_paramet 17 19 IMPLICIT NONE 18 20 ! 19 INCLUDE "dimensions.h" 20 INCLUDE "paramet.h" 21 ! 21 22 23 22 24 ! ..... variables en arguments ...... 23 ! 25 24 26 INTEGER :: klevel 25 27 REAL :: rot( ip1jm,klevel ) 26 28 REAL :: x( ip1jmp1,klevel ), y( ip1jm,klevel ) 27 ! 29 28 30 ! ... variables locales ... 29 ! 31 30 32 INTEGER :: l, ij 31 ! 32 ! 33 34 33 35 DO l = 1,klevel 34 ! 36 35 37 DO ij = 1, ip1jm - 1 36 38 rot( ij,l ) = y( ij+1 , l ) - y( ij,l ) + & 37 39 x(ij +iip1, l ) - x( ij,l ) 38 40 ENDDO 39 ! 41 40 42 ! .... correction pour rot( iip1,j,l) .... 41 43 ! .... rot(iip1,j,l)= rot(1,j,l) ... … … 44 46 rot( ij,l ) = rot( ij -iim,l ) 45 47 ENDDO 46 ! 48 47 49 END DO 48 50
Note: See TracChangeset
for help on using the changeset viewer.