- Timestamp:
- 2014-02-04T13:14:23+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90
r4429 r4431 307 307 ! third recurrence: Xk = (Zk - Sk Xk+1 ) / Tk (result is the after tracer) 308 308 #if defined key_z_first 309 pta(ji,jj,jpkm1,jn) = pta(ji,jj,jpkm1,jn) / zwt(ji,jj,jpkm1) * tmask(ji,jj,jpkm1) 310 DO jk = mbkmax(ji,jj)-2, 1, -1 ! jpk-2 311 pta(ji,jj,jk,jn) = ( pta(ji,jj,jk,jn) - zws(ji,jj,jk) * pta(ji,jj,jk+1,jn) ) & 309 ! pta(ji,jj,jpkm1,jn) = pta(ji,jj,jpkm1,jn) / zwt(ji,jj,jpkm1) * tmask(ji,jj,jpkm1) 310 jk = mbkmax(ji,jj)-1 311 IF(jk>0)THEN 312 pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) / zwt(ji,jj,jk) * tmask(ji,jj,jk) 313 DO jk = mbkmax(ji,jj)-2, 1, -1 ! jpk-2 314 pta(ji,jj,jk,jn) = ( pta(ji,jj,jk,jn) - zws(ji,jj,jk) * pta(ji,jj,jk+1,jn) ) & 312 315 & / zwt(ji,jj,jk) * tmask(ji,jj,jk) 313 END DO 316 END DO 317 END IF 314 318 END DO 315 319 END DO
Note: See TracChangeset
for help on using the changeset viewer.