Changeset 4428


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

Add mbkmax to dynspg_ts and trazdf_imp

Location:
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90

    r3432 r4428  
    209209      DO jj = 1, jpj 
    210210         DO ji = 1, jpi 
    211             DO jk = 1, jpkm1 
     211            DO jk = 1, mbkmax(ji,jj)-1 ! jpkm1 
    212212#else 
    213213      DO jk = 1, jpkm1 
     
    242242      DO jj = 2, jpjm1 
    243243         DO ji = 2, jpim1 
    244             DO jk = 1, jpkm1 
     244            DO jk = 1, mbkmax(ji,jj)-1 ! jpkm1 
    245245#else 
    246246      DO jk = 1, jpkm1                    ! -------------------------- 
     
    691691            zu_asp = ( zu_sum(ji,jj) - ub_b(ji,jj) ) / z2dt_b 
    692692            zv_asp = ( zv_sum(ji,jj) - vb_b(ji,jj) ) / z2dt_b 
    693             DO jk=1,jpkm1 
     693            DO jk=1,mbkmax(ji,jj)-1 ! jpkm1 
    694694               ua(ji,jj,jk) = ua(ji,jj,jk) + zu_asp 
    695695               va(ji,jj,jk) = va(ji,jj,jk) + zv_asp 
  • 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.