Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/nxgrad_loc.f90
- 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_loc.f90
r5136 r5159 1 1 SUBROUTINE nxgrad_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 ),x( ijb_u:ije_u,klevel ) … … 21 23 INTEGER :: l,ij 22 24 INTEGER :: ijb,ije 23 ! 24 ! 25 26 25 27 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 26 28 DO l = 1,klevel 27 ! 29 28 30 ijb=ij_begin 29 31 ije=ij_end … … 33 35 y( ij,l ) = ( rot( ij,l ) - rot( ij-1,l ) ) * cvsurcuv( ij ) 34 36 END DO 35 ! 37 36 38 ! ..... correction pour y ( 1,j,l ) ...... 37 ! 39 38 40 ! .... y(1,j,l)= y(iip1,j,l) .... 39 41 !DIR$ IVDEP … … 41 43 y( ij,l ) = y( ij +iim,l ) 42 44 END DO 43 ! 45 44 46 ijb=ij_begin 45 47 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.