Changeset 4444


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

Use of mbkmax in tra_adv_tvd

File:
1 edited

Legend:

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

    r4409 r4444  
    130130      END IF 
    131131      ! 
    132       zwi(:,:,1:jpkf) = 0.e0 
     132#if defined key_z_first 
     133      DO jj = 1, jpj, 1 
     134         DO ji = 1, jpi, 1 
     135            zwi(ji,jj,1:mbkmax(ji,jj)) = 0.e0_wp 
     136         END DO 
     137      END DO 
     138#else 
     139      zwi(:,:,1:jpk) = 0.e0_wp 
     140#endif 
    133141      ! 
    134142      !                                                          ! =========== 
     
    139147         ! ARPDBG: FINISS was using jpk correct here given that it can be  
    140148         ! below the ocean floor?? 
     149#if defined key_z_first 
     150         DO jj = 1, jpj, 1 
     151            DO ji = 1, jpi, 1 
     152               zwx(ji,jj,mbkmax(ji,jj)) = 0.0_wp 
     153               zwy(ji,jj,mbkmax(ji,jj)) = 0.0_wp 
     154               zwz(ji,jj,mbkmax(ji,jj)) = 0.0_wp 
     155               zwi(ji,jj,mbkmax(ji,jj)) = 0.0_wp 
     156            END DO             
     157         END DO 
     158#else 
    141159         zwx(:,:,jpkf) = 0.e0    ;    zwz(:,:,jpkf) = 0.e0 
    142160         zwy(:,:,jpkf) = 0.e0    ;    zwi(:,:,jpkf) = 0.e0 
    143  
     161#endif 
    144162         ! 2. upstream advection with initial mass fluxes & intermediate update 
    145163         ! -------------------------------------------------------------------- 
     
    149167         DO jj = 1, jpjm1 
    150168            DO ji = 1, jpim1 
    151                DO jk = 1, jpkfm1 
     169               DO jk = 1, mbkmax(ji,jj)-1 ! jpkfm1 
    152170#else 
    153171         DO jk = 1, jpkfm1 
     
    177195         DO jj = 1, jpj 
    178196            DO ji = 1, jpi 
    179                DO jk = 2, jpkfm1 
     197               DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    180198#else 
    181199         DO jk = 2, jpkfm1 
     
    196214         DO jj = 2, jpjm1 
    197215            DO ji = 2, jpim1 
    198                DO jk = 1, jpkfm1 
     216               DO jk = 1, mbkmax(ji,jj)-1 ! jpkfm1 
    199217                  z2dtt = p2dt(jk) 
    200218#else 
     
    240258         DO jj = 1, jpjm1 
    241259            DO ji = 1, jpim1 
    242                DO jk = 1, jpkfm1 
     260               DO jk = 1, mbkmax(ji,jj)-1 ! jpkfm1 
    243261#else 
    244262         DO jk = 1, jpkfm1 
     
    259277            DO ji = 1, jpi 
    260278               zwz(ji,jj,1) = 0.e0   ! Surface value 
    261                DO jk = 2, jpkfm1 
     279               DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    262280#else 
    263281         zwz(:,:,1) = 0.e0           ! Surface value 
     
    291309         DO jj = 2, jpjm1 
    292310            DO ji = 2, jpim1 
    293                DO jk = 1, jpkfm1 
     311               DO jk = 1, mbkmax(ji,jj)-1 ! jpkfm1 
    294312#else 
    295313         DO jk = 1, jpkfm1 
Note: See TracChangeset for help on using the changeset viewer.