Changeset 4447


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

Added mbkmax to dyn_ldf_iso

File:
1 edited

Legend:

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

    r3432 r4447  
    169169         DO jj = 2, jpjm1       ! set the slopes of iso-level 
    170170            DO ji = fs_2, fs_jpim1     
    171                DO jk = 1, jpk        
     171               DO jk = 1, mbkmax(ji,jj) ! jpk        
    172172#else 
    173173         DO jk = 1, jpk         ! set the slopes of iso-level 
     
    368368               ! when jk is 1 but that doesn't matter. 
    369369!DIR$ SAFE_ADDRESS 
    370                DO jk = 1, jpkm1, 1 
     370               DO jk = 1, mbkmax(ji,jj)-1 ! jpkm1, 1 
    371371 
    372372                  ! Vertical u- and v-shears at level jk and jk+1 
     
    465465! when jk is 1 but that doesn't matter. 
    466466!DIR$ SAFE_ADDRESS 
    467                DO jk = 1, jpkm1, 1 
     467               DO jk = 1, mbkmax(ji,jj)-1, 1 ! jpkm1, 1 
    468468 
    469469                  zdk1v    = ( vb(ji  ,jj  ,jk) -vb(ji  ,jj  ,jk+1) ) * vmask(ji  ,jj,jk+1) 
     
    566566               ! when jk is 1 but that doesn't matter. 
    567567!DIR$ SAFE_ADDRESS 
    568                DO jk = 1, jpkm1, 1 
     568               DO jk = 1, mbkmax(ji,jj)-1, 1 ! jpkm1, 1 
    569569 
    570570                  ! Vertical u- and v-shears at level jk and jk+1 
     
    661661! when jk is 1 but that doesn't matter. 
    662662!DIR$ SAFE_ADDRESS 
    663                DO jk = 1, jpkm1, 1 
     663               DO jk = 1, mbkmax(ji,jj)-1, 1 ! jpkm1, 1 
    664664 
    665665                  zdk1v    = ( vb(ji  ,jj  ,jk) -vb(ji  ,jj  ,jk+1) ) * vmask(ji  ,jj,jk+1) 
     
    918918 
    919919         DO ji = 2, jpi 
    920             DO jk = 1, jpk 
     920            DO jk = 1, mbkmax(ji,jj) ! jpk 
    921921               ! i-gradient of u at jj 
    922922               zdiu (ji,jk) = tmask(ji,jj  ,jk) * ( ub(ji,jj  ,jk) - ub(ji-1,jj  ,jk) ) 
     
    931931 
    932932         DO ji = 1, jpim1 
    933             DO jk = 1, jpk 
     933            DO jk = 1, mbkmax(ji,jj) ! jpk 
    934934               ! i-gradient of v at jj 
    935935               zdiv (ji,jk) = fmask(ji,jj  ,jk) * ( vb(ji+1,jj,jk) - vb(ji  ,jj  ,jk) ) 
     
    951951         ! interior (2=<jk=<jpk-1) on U field 
    952952         DO ji = 2, jpim1 
    953             DO jk = 2, jpkm1 
     953            DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1 
    954954               zcoef0= 0.5 * aht0 * umask(ji,jj,jk) 
    955955 
     
    977977         ! interior (2=<jk=<jpk-1) on V field 
    978978         DO ji = 2, jpim1 
    979             DO jk = 2, jpkm1 
     979            DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1 
    980980               zcoef0= 0.5 * aht0 * vmask(ji,jj,jk) 
    981981 
     
    10051005         ! ------------------------------------------------------------------- 
    10061006!         DO ji = 2, jpim1 
    1007             zfuw(ji,1)   = 0.0_wp 
    1008             zfuw(ji,jpk) = 0.0_wp 
    1009             zfvw(ji,1)   = 0.0_wp 
    1010             zfvw(ji,jpk) = 0.0_wp 
    1011  
    1012             DO jk = 1, jpkm1 
     1007            zfuw(ji,1)             = 0.0_wp 
     1008            zfuw(ji,mbkmax(ji,jj)) = 0.0_wp 
     1009            zfvw(ji,1)             = 0.0_wp 
     1010            zfvw(ji,mbkmax(ji,jj)) = 0.0_wp 
     1011 
     1012            DO jk = 1, mbkmax(ji,jj)-1 ! jpkm1 
    10131013               ! volume elements 
    10141014               zbu = e1u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,jk) 
Note: See TracChangeset for help on using the changeset viewer.