Changeset 4431 for branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90

Timestamp:

2014-02-04T13:14:23+01:00 (10 years ago)

Author:

trackstand2

Message:

Use of mbkmax in eosbn2 and trazdf_imp

File:

• branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90 (modified) (1 diff)

branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90

@@ -307 +307 @@
          ! third recurrence:    Xk = (Zk - Sk Xk+1 ) / Tk   (result is the after tracer)
 #if defined key_z_first
-               pta(ji,jj,jpkm1,jn) = pta(ji,jj,jpkm1,jn) / zwt(ji,jj,jpkm1) * tmask(ji,jj,jpkm1)
-               DO jk = mbkmax(ji,jj)-2, 1, -1 ! jpk-2
-                  pta(ji,jj,jk,jn) = ( pta(ji,jj,jk,jn) - zws(ji,jj,jk) * pta(ji,jj,jk+1,jn) )   &
+!               pta(ji,jj,jpkm1,jn) = pta(ji,jj,jpkm1,jn) / zwt(ji,jj,jpkm1) * tmask(ji,jj,jpkm1)
+               jk = mbkmax(ji,jj)-1
+               IF(jk>0)THEN
+                  pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) / zwt(ji,jj,jk) * tmask(ji,jj,jk)
+                  DO jk = mbkmax(ji,jj)-2, 1, -1 ! jpk-2
+                     pta(ji,jj,jk,jn) = ( pta(ji,jj,jk,jn) - zws(ji,jj,jk) * pta(ji,jj,jk+1,jn) )   &
                      &             / zwt(ji,jj,jk) * tmask(ji,jj,jk)
-               END DO
+                  END DO
+               END IF
             END DO
          END DO