Ignore:
Timestamp:
2020-07-08T16:27:49+02:00 (3 months ago)
Author:
clem
Message:

should solve ticket #2435

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/r4.0/r4.0-HEAD/src/OCE/ZDF/zdftke.F90

    r12703 r13268  
    2828   !!            3.6  !  2014-11  (P. Mathiot) add ice shelf capability 
    2929   !!            4.0  !  2017-04  (G. Madec)  remove CPP ddm key & avm at t-point only  
    30    !!             -   !  2017-05  (G. Madec)  add top/bottom friction as boundary condition (ln_drg) 
     30   !!             -   !  2017-05  (G. Madec)  add top/bottom friction as boundary condition 
    3131   !!---------------------------------------------------------------------- 
    3232 
     
    7171   REAL(wp) ::   rn_emin0  ! surface minimum value of tke   [m2/s2] 
    7272   REAL(wp) ::   rn_bshear ! background shear (>0) currently a numerical threshold (do not change it) 
    73    LOGICAL  ::   ln_drg    ! top/bottom friction forcing flag  
    7473   INTEGER  ::   nn_etau   ! type of depth penetration of surface tke (=0/1/2/3) 
    7574   INTEGER  ::      nn_htau   ! type of tke profile of penetration (=0/1) 
     
    229228      ! Note that stress averaged is done using an wet-only calculation of u and v at t-point like in zdfsh2 
    230229      ! 
    231       IF( ln_drg ) THEN       !== friction used as top/bottom boundary condition on TKE 
    232          ! 
    233          DO jj = 2, jpjm1           ! bottom friction 
     230      IF( .NOT.ln_drg_OFF ) THEN    !== friction used as top/bottom boundary condition on TKE 
     231         ! 
     232         DO jj = 2, jpjm1              ! bottom friction 
    234233            DO ji = fs_2, fs_jpim1     ! vector opt. 
    235234               zmsku = ( 2. - umask(ji-1,jj,mbkt(ji,jj)) * umask(ji,jj,mbkt(ji,jj)) ) 
     
    644643      INTEGER ::   ios 
    645644      !! 
    646       NAMELIST/namzdf_tke/ rn_ediff, rn_ediss , rn_ebb , rn_emin  ,          & 
    647          &                 rn_emin0, rn_bshear, nn_mxl , ln_mxl0  ,          & 
    648          &                 rn_mxl0 , nn_pdl   , ln_drg , ln_lc    , rn_lc,   & 
     645      NAMELIST/namzdf_tke/ rn_ediff, rn_ediss , rn_ebb , rn_emin  ,  & 
     646         &                 rn_emin0, rn_bshear, nn_mxl , ln_mxl0  ,  & 
     647         &                 rn_mxl0 , nn_pdl   , ln_lc  , rn_lc,      & 
    649648         &                 nn_etau , nn_htau  , rn_efr , rn_eice   
    650649      !!---------------------------------------------------------------------- 
     
    676675         WRITE(numout,*) '         surface mixing length = F(stress) or not    ln_mxl0   = ', ln_mxl0 
    677676         WRITE(numout,*) '         surface  mixing length minimum value        rn_mxl0   = ', rn_mxl0 
    678          WRITE(numout,*) '      top/bottom friction forcing flag            ln_drg    = ', ln_drg 
    679677         WRITE(numout,*) '      Langmuir cells parametrization              ln_lc     = ', ln_lc 
    680678         WRITE(numout,*) '         coef to compute vertical velocity of LC     rn_lc  = ', rn_lc 
     
    684682         WRITE(numout,*) '          below sea-ice:  =0 ON                      rn_eice   = ', rn_eice 
    685683         WRITE(numout,*) '          =4 OFF when ice fraction > 1/4   ' 
    686          IF( ln_drg ) THEN 
     684         IF( .NOT.ln_drg_OFF ) THEN 
    687685            WRITE(numout,*) 
    688686            WRITE(numout,*) '   Namelist namdrg_top/_bot:   used values:' 
Note: See TracChangeset for help on using the changeset viewer.