Changeset 5086 for LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv2.F
- Timestamp:
- Jul 19, 2024, 7:54:50 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv2.F
r2600 r5086 34 34 c 35 35 c 36 DO 5l = 1,llm36 DO l = 1,llm 37 37 c 38 DO 2ij = iip2, ip1jm - 138 DO ij = iip2, ip1jm - 1 39 39 du(ij,l) = du(ij,l) + 0.5* ( teta( ij,l ) + teta( ij+1,l ) ) * 40 40 * ( pkf( ij,l ) - pkf(ij+1,l) ) + bern(ij,l) - bern(ij+1,l) 41 2 CONTINUE41 END DO 42 42 c 43 43 c … … 46 46 c 47 47 CDIR$ IVDEP 48 DO 3ij = iip1+ iip1, ip1jm, iip148 DO ij = iip1+ iip1, ip1jm, iip1 49 49 du( ij,l ) = du( ij - iim,l ) 50 3 CONTINUE50 END DO 51 51 c 52 52 c 53 DO 4ij = 1,ip1jm53 DO ij = 1,ip1jm 54 54 dv( ij,l) = dv(ij,l) + 0.5 * ( teta(ij,l) + teta( ij+iip1,l ) ) * 55 55 * ( pkf(ij+iip1,l) - pkf( ij,l ) ) 56 56 * + bern( ij+iip1,l ) - bern( ij ,l ) 57 4 CONTINUE57 END DO 58 58 c 59 5 CONTINUE59 END DO 60 60 c 61 61 RETURN
Note: See TracChangeset
for help on using the changeset viewer.