- Timestamp:
- 2017-01-04T17:47:47+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.F90
r7508 r7525 78 78 CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv ) 79 79 !$OMP PARALLEL DO schedule(static) private(jk, jj, ji) 80 DO jk = 1, jpk81 DO jj = 1, jpj82 DO ji = 1, jpi83 ztrdu(ji,jj,jk) = ua(ji,jj,jk)84 ztrdv(ji,jj,jk) = va(ji,jj,jk)85 END DO86 END DO87 END DO80 DO jk = 1, jpk 81 DO jj = 1, jpj 82 DO ji = 1, jpi 83 ztrdu(ji,jj,jk) = ua(ji,jj,jk) 84 ztrdv(ji,jj,jk) = va(ji,jj,jk) 85 END DO 86 END DO 87 END DO 88 88 ENDIF 89 89 … … 145 145 IF( l_trddyn ) THEN ! save the vertical advection trends for diagnostic 146 146 !$OMP PARALLEL DO schedule(static) private(jk, jj, ji) 147 148 149 150 151 152 153 154 147 DO jk = 1, jpk 148 DO jj = 1, jpj 149 DO ji = 1, jpi 150 ztrdu(ji,jj,jk) = ua(ji,jj,jk) - ztrdu(ji,jj,jk) 151 ztrdv(ji,jj,jk) = va(ji,jj,jk) - ztrdv(ji,jj,jk) 152 END DO 153 END DO 154 END DO 155 155 CALL trd_dyn( ztrdu, ztrdv, jpdyn_zad, kt ) 156 156 CALL wrk_dealloc( jpi, jpj, jpk, ztrdu, ztrdv )
Note: See TracChangeset
for help on using the changeset viewer.