Ignore:
Timestamp:
2020-07-09T12:39:43+02:00 (3 months ago)
Author:
clem
Message:

merge with r4.0-HEAD at r13278

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/OCE/ZDF/zdftke.F90

    r13249 r13279  
    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 
     
    7979   REAL(wp) ::   rn_emin0  ! surface minimum value of tke   [m2/s2] 
    8080   REAL(wp) ::   rn_bshear ! background shear (>0) currently a numerical threshold (do not change it) 
    81    LOGICAL  ::   ln_drg    ! top/bottom friction forcing flag  
    8281   INTEGER  ::   nn_etau   ! type of depth penetration of surface tke (=0/1/2/3) 
    8382   INTEGER  ::      nn_htau   ! type of tke profile of penetration (=0/1) 
     
    247246      ! Note that stress averaged is done using an wet-only calculation of u and v at t-point like in zdfsh2 
    248247      ! 
    249       IF( ln_drg ) THEN       !== friction used as top/bottom boundary condition on TKE 
    250          ! 
    251          DO jj = 2, jpjm1           ! bottom friction 
     248      IF( .NOT.ln_drg_OFF ) THEN    !== friction used as top/bottom boundary condition on TKE 
     249         ! 
     250         DO jj = 2, jpjm1              ! bottom friction 
    252251            DO ji = fs_2, fs_jpim1     ! vector opt. 
    253252               zmsku = ( 2. - umask(ji-1,jj,mbkt(ji,jj)) * umask(ji,jj,mbkt(ji,jj)) ) 
     
    717716         &                 rn_emin0, rn_bshear, nn_mxl   , ln_mxl0  ,  & 
    718717         &                 rn_mxl0 , nn_mxlice, rn_mxlice,             & 
    719          &                 nn_pdl  , ln_drg   , ln_lc    , rn_lc,      & 
     718         &                 nn_pdl  , ln_lc    , rn_lc,                 & 
    720719         &                 nn_etau , nn_htau  , rn_efr   , nn_eice   
    721720      !!---------------------------------------------------------------------- 
     
    760759            END SELECT 
    761760         ENDIF 
    762          WRITE(numout,*) '      top/bottom friction forcing flag            ln_drg    = ', ln_drg 
    763761         WRITE(numout,*) '      Langmuir cells parametrization              ln_lc     = ', ln_lc 
    764762         WRITE(numout,*) '         coef to compute vertical velocity of LC     rn_lc  = ', rn_lc 
     
    775773            CALL ctl_stop( 'zdf_tke_init: wrong value for nn_eice, should be 0,1,2, or 3') 
    776774         END SELECT       
    777          IF( ln_drg ) THEN 
     775         IF( .NOT.ln_drg_OFF ) THEN 
    778776            WRITE(numout,*) 
    779777            WRITE(numout,*) '   Namelist namdrg_top/_bot:   used values:' 
Note: See TracChangeset for help on using the changeset viewer.