Changeset 5086 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergst.F
- Timestamp:
- Jul 19, 2024, 7:54:50 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergst.F
r4593 r5086 30 30 c 31 31 c 32 DO 10l = 1,klevel32 DO l = 1,klevel 33 33 c 34 DO 1ij = iip2, ip1jm - 134 DO ij = iip2, ip1jm - 1 35 35 div( ij + 1, l ) = x(ij+1,l) - x(ij,l)+ y(ij-iim,l)-y(ij+1,l) 36 1 CONTINUE36 END DO 37 37 c 38 38 c .... correction pour div( 1,j,l) ...... … … 40 40 c 41 41 CDIR$ IVDEP 42 DO 3ij = iip2,ip1jm,iip142 DO ij = iip2,ip1jm,iip1 43 43 div( ij,l ) = div( ij + iim,l ) 44 3 CONTINUE44 END DO 45 45 c 46 46 c .... calcul aux poles ..... 47 47 c 48 48 c 49 DO 5i = 1,iim49 DO i = 1,iim 50 50 aiy1(i)= y(i,l) 51 51 aiy2(i)= y(i+ip1jmi1,l) 52 5 CONTINUE52 END DO 53 53 sumypn = SSUM ( iim,aiy1,1 ) 54 54 sumyps = SSUM ( iim,aiy2,1 ) 55 DO 7i = 1,iip155 DO i = 1,iip1 56 56 div( i , l ) = - sumypn/iim 57 57 div( i + ip1jm, l ) = sumyps/iim 58 7 CONTINUE58 END DO 59 59 c 60 10 CONTINUE60 END DO 61 61 RETURN 62 62 END
Note: See TracChangeset
for help on using the changeset viewer.