Ignore:
Timestamp:
2014-02-04T13:14:23+01:00 (7 years ago)
Author:
trackstand2
Message:

Use of mbkmax in eosbn2 and trazdf_imp

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90

    r4429 r4431  
    307307         ! third recurrence:    Xk = (Zk - Sk Xk+1 ) / Tk   (result is the after tracer) 
    308308#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) )   & 
    312315                     &             / zwt(ji,jj,jk) * tmask(ji,jj,jk) 
    313                END DO 
     316                  END DO 
     317               END IF 
    314318            END DO 
    315319         END DO 
Note: See TracChangeset for help on using the changeset viewer.