New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4428 for branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90 – NEMO

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

Add mbkmax to dynspg_ts and trazdf_imp

File:
1 edited

Legend:

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

    r4409 r4428  
    144144                  DO ji = 1, jpi 
    145145                     zwt(ji,jj,1) = 0._wp 
    146                      DO jk = 2, jpk 
     146                     DO jk = 2, mbkmax(ji,jj) ! jpk 
    147147                        zwt(ji,jj,jk) = avt  (ji,jj,jk) 
    148148                     END DO 
     
    153153                  DO ji = 1, jpi 
    154154                     zwt(ji,jj,1) = 0._wp 
    155                      DO jk = 2, jpk 
     155                     DO jk = 2, mbkmax(ji,jj) ! jpk 
    156156                        zwt(ji,jj,jk) = fsavs(ji,jj,jk) 
    157157                     END DO 
     
    172172               DO jj = 2, jpjm1 
    173173                  DO ji = 2, jpim1 
    174                      DO jk = 2, jpkm1 
     174                     DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1 
    175175#else 
    176176               DO jk = 2, jpkm1 
     
    186186               DO jj = 2, jpjm1 
    187187                  DO ji = 2, jpim1 
    188                      DO jk = 2, jpkm1 
     188                     DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1 
    189189#else 
    190190               DO jk = 2, jpkm1 
     
    204204            DO jj = 2, jpjm1 
    205205               DO ji = 2, jpim1 
    206                   DO jk = 1, jpkm1 
     206                  DO jk = 1, mbkmax(ji,jj) ! jpkm1 
    207207                     ! after scale factor at T-point 
    208208                     ze3ta =  ( 1. - r_vvl ) +        r_vvl   * fse3t_a(ji,jj,jk)  
     
    250250#if defined key_z_first 
    251251                  zwt(ji,jj,1) = zwd(ji,jj,1) 
    252                   DO jk = 2, jpkm1 
     252                  DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1 
    253253                    zwt(ji,jj,jk) = zwd(ji,jj,jk) - zwi(ji,jj,jk) * zws(ji,jj,jk-1) / zwt(ji,jj,jk-1) 
    254254                  END DO 
     
    279279               ze3tn = ( 1. - r_vvl ) + r_vvl * fse3t(ji,jj,1) 
    280280               pta(ji,jj,1,jn) = ze3tb * ptb(ji,jj,1,jn) + p2dt(1) * ze3tn * pta(ji,jj,1,jn) 
    281                DO jk = 2, jpkm1 
     281               DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1 
    282282                  ze3tb = ( 1. - r_vvl ) + r_vvl * fse3t_b(ji,jj,jk) 
    283283                  ze3tn = ( 1. - r_vvl ) + r_vvl * fse3t  (ji,jj,jk) 
     
    308308#if defined key_z_first 
    309309               pta(ji,jj,jpkm1,jn) = pta(ji,jj,jpkm1,jn) / zwt(ji,jj,jpkm1) * tmask(ji,jj,jpkm1) 
    310                DO jk = jpk-2, 1, -1 
     310               DO jk = mbkmax(ji,jj)-2, 1, -1 ! jpk-2 
    311311                  pta(ji,jj,jk,jn) = ( pta(ji,jj,jk,jn) - zws(ji,jj,jk) * pta(ji,jj,jk+1,jn) )   & 
    312312                     &             / zwt(ji,jj,jk) * tmask(ji,jj,jk) 
Note: See TracChangeset for help on using the changeset viewer.