- Timestamp:
- Aug 2, 2024, 9:58:25 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/nxgrad_gam_loc.f90
r5136 r5159 1 1 SUBROUTINE nxgrad_gam_loc( klevel, rot, x, y ) 2 ! 2 3 3 ! P. Le Van 4 ! 4 5 5 ! ******************************************************************** 6 6 ! calcul du gradient tourne de pi/2 du rotationnel du vect.v … … 8 8 ! rot est un argument d'entree pour le s-prog 9 9 ! x et y sont des arguments de sortie pour le s-prog 10 ! 10 11 11 USE parallel_lmdz 12 12 USE lmdz_comgeom 13 13 14 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 15 USE lmdz_paramet 14 16 IMPLICIT NONE 15 17 ! 16 INCLUDE "dimensions.h" 17 INCLUDE "paramet.h" 18 19 18 20 INTEGER :: klevel 19 21 REAL :: rot( ijb_v:ije_v,klevel ) … … 23 25 external ismin,ismax 24 26 INTEGER :: ijb,ije 25 ! 27 26 28 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 27 29 DO l = 1,klevel 28 ! 30 29 31 ijb=ij_begin 30 32 ije=ij_end … … 34 36 y( ij,l ) = (rot( ij,l ) - rot( ij-1,l )) * cvscuvgam( ij ) 35 37 END DO 36 ! 38 37 39 ! ..... correction pour y ( 1,j,l ) ...... 38 ! 40 39 41 ! .... y(1,j,l)= y(iip1,j,l) .... 40 42 !DIR$ IVDEP … … 42 44 y( ij,l ) = y( ij +iim,l ) 43 45 END DO 44 ! 46 45 47 ijb=ij_begin 46 48 ije=ij_end+iip1 … … 63 65 ENDDO 64 66 ENDIF 65 ! 67 66 68 END DO 67 69 !$OMP END DO NOWAIT
Note: See TracChangeset
for help on using the changeset viewer.