Changeset 4446


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

Use of mbkmax in ldf_slp→ldf_slp_mxl

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/LDF/ldfslp.F90

    r4427 r4446  
    172172      DO jj = 1, jpjm1           !==   i- & j-gradient of density   ==! 
    173173         DO ji = 1, jpim1 
    174             DO jk = 1, jpkf 
     174            DO jk = 1, mbkmax(ji,jj) ! jpkf 
    175175#else 
    176176      DO jk = 1, jpkf            !==   i- & j-gradient of density   ==! 
     
    202202         DO ji = 1, jpi 
    203203            zdzr(ji,jj,1) = 0._wp   !==   Local vertical density gradient at T-point   == !   (evaluated from N^2) 
    204             DO jk = 2, jpkfm1 
     204            DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    205205               zdzr(ji,jj,jk) = zm1_g * ( prd(ji,jj,jk) + 1._wp )              & 
    206206                  &                   * ( pn2(ji,jj,jk) + pn2(ji,jj,jk+1) ) * ( 1._wp - 0.5_wp * tmask(ji,jj,jk+1) ) 
     
    231231      DO jj = 2, jpjm1               !* Slopes at u and v points 
    232232         DO ji = 2, jpim1 
    233             DO jk = 2, jpkfm1 
     233            DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    234234#else 
    235235      DO jk = 2, jpkfm1                            !* Slopes at u and v points 
     
    275275      DO jj = 2, jpjm1, MAX(1, jpj-3)                        ! rows jj=2 and =jpjm1 only 
    276276         DO ji = 2, jpim1   
    277             DO jk = 2, jpkfm1 
     277            DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    278278#else 
    279279      DO jk = 2, jpkfm1 
     
    297297      DO jj = 3, jpj-2                               ! other rows 
    298298         DO ji = 2, jpim1 
    299             DO jk = 2, jpkfm1 
     299            DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    300300#else 
    301301         DO jj = 3, jpj-2                               ! other rows 
     
    319319      DO jj = 2, jpjm1 
    320320         DO ji = 2, jpim1 
    321             DO jk = 2, jpkfm1 
     321            DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    322322#else 
    323323         !                                        !* decrease along coastal boundaries 
     
    340340      DO jj = 2, jpjm1 
    341341         DO ji = 2, jpim1 
    342             DO jk = 2, jpkfm1 
     342            DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    343343#else 
    344344      DO jk = 2, jpkfm1 
     
    384384      DO jj = 2, jpjm1, MAX(1, jpj-3)                           ! rows jj=2 and =jpjm1 only 
    385385         DO ji = 2, jpim1 
    386             DO jk = 2, jpkfm1 
     386            DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    387387#else 
    388388      DO jk = 2, jpkfm1 
     
    407407      DO jj = 3, jpj-2                                  ! other rows 
    408408         DO ji = 2, jpim1 
    409             DO jk = 2, jpkfm1 
     409            DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    410410#else 
    411411         DO jj = 3, jpj-2                               ! other rows 
     
    430430      DO jj = 2, jpjm1 
    431431         DO ji = 2, jpim1 
    432             DO jk = 2, jpkfm1 
     432            DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 
    433433#else 
    434434         !                                        !* decrease along coastal boundaries 
     
    818818      DO jj = 1, jpj 
    819819         DO ji = 1, jpi 
    820             DO jk = 1, jpkf          ! =1 inside the mixed layer, =0 otherwise 
     820            DO jk = 1, mbkmax(ji,jj) ! jpkf ! =1 inside the mixed layer, =0 otherwise 
    821821#else 
    822822      DO jk = 1, jpkf                ! =1 inside the mixed layer, =0 otherwise 
     
    858858            ! 
    859859            !                          !- vertical density gradient for u- and v-slopes (from dzr at T-point) 
     860#if defined key_z_first 
     861            iku = MIN(  MAX( 1, nmln(ji,jj) , nmln(ji+1,jj) ) , mbkmax(ji,jj)-1  ) ! ML (MAX of T-pts, bound by mbkmax(ji,jj)) 
     862            ikv = MIN(  MAX( 1, nmln(ji,jj) , nmln(ji,jj+1) ) , mbkmax(ji,jj)-1  ) !  
     863#else 
    860864            iku = MIN(  MAX( 1, nmln(ji,jj) , nmln(ji+1,jj) ) , jpkm1  )   ! ML (MAX of T-pts, bound by jpkm1) 
    861865            ikv = MIN(  MAX( 1, nmln(ji,jj) , nmln(ji,jj+1) ) , jpkm1  )   !  
     866#endif 
    862867            zbu = 0.5_wp * ( p_dzr(ji,jj,iku) + p_dzr(ji+1,jj  ,iku) ) 
    863868            zbv = 0.5_wp * ( p_dzr(ji,jj,ikv) + p_dzr(ji  ,jj+1,ikv) ) 
     
    875880            !                    !==   i- & j-slopes at w-points just below the Mixed Layer   ==! 
    876881            ! 
     882#if defined key_z_first 
     883            ik   = MIN( nmln(ji,jj) + 1, mbkmax(ji,jj) ) 
     884#else 
    877885            ik   = MIN( nmln(ji,jj) + 1, jpkf ) 
     886#endif 
    878887            ikm1 = MAX( 1, ik-1 ) 
    879888            !                          !- vertical density gradient for w-slope (from N^2) 
Note: See TracChangeset for help on using the changeset viewer.