Changeset 5086 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dudv1_loc.F
- Timestamp:
- Jul 19, 2024, 7:54:50 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dudv1_loc.F
r4593 r5086 28 28 29 29 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 30 DO 10l = 1,llm30 DO l = 1,llm 31 31 c 32 32 ijb=ij_begin … … 36 36 if (pole_sud) ije=ij_end-iip1 37 37 38 DO 2 ij = ijb, ije-138 DO ij = ijb, ije-1 39 39 du( ij,l ) = 0.125 *( vorpot(ij-iip1, l) + vorpot( ij, l) ) * 40 40 * ( pbarv(ij-iip1, l) + pbarv(ij-iim, l) + 41 41 * pbarv( ij , l) + pbarv(ij+ 1 , l) ) 42 2 CONTINUE42 END DO 43 43 44 44 … … 46 46 if (pole_nord) ijb=ij_begin 47 47 48 DO 3 ij = ijb, ije-148 DO ij = ijb, ije-1 49 49 dv( ij+1,l ) = - 0.125 *( vorpot(ij, l) + vorpot(ij+1, l) ) * 50 50 * ( pbaru(ij, l) + pbaru(ij+1 , l) + 51 51 * pbaru(ij+iip1, l) + pbaru(ij+iip2, l) ) 52 3 CONTINUE52 END DO 53 53 c 54 54 c .... correction pour dv( 1,j,l ) ..... … … 56 56 c 57 57 CDIR$ IVDEP 58 DO 4ij = ijb, ije, iip158 DO ij = ijb, ije, iip1 59 59 dv( ij,l ) = dv( ij + iim, l ) 60 4 CONTINUE60 END DO 61 61 c 62 10 CONTINUE62 END DO 63 63 c$OMP END DO NOWAIT 64 64 RETURN
Note: See TracChangeset
for help on using the changeset viewer.