Changeset 9546


Ignore:
Timestamp:
2018-05-04T07:57:05+02:00 (2 years ago)
Author:
gm
Message:

dev_merge_2017: add a TKE parameter to control the activation of nn_etau below sea-ice

Location:
branches/2017/dev_merge_2017/NEMOGCM
Files:
7 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/AGRIF_NORDIC/EXP00/namelist_cfg

    r9531 r9546  
    357357&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
    358358!----------------------------------------------------------------------- 
     359      rn_eice     =   0       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4    
    359360/ 
    360361 
  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/EXP00/namelist_cfg

    r9533 r9546  
    350350/ 
    351351!----------------------------------------------------------------------- 
     352&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     353!----------------------------------------------------------------------- 
     354      rn_eice     =   0       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4    
     355/ 
     356!----------------------------------------------------------------------- 
    352357&namzdf_iwm    !    internal wave-driven mixing parameterization        (ln_zdfiwm =T) 
    353358!----------------------------------------------------------------------- 
  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/cpp_ORCA2_LIM3_PISCES.fcm

    r9303 r9546  
    1 bld::tool::fppkeys   key_lim3 key_top key_iomput key_mpp_mpi key_nosignedzero 
     1bld::tool::fppkeys   key_lim3 key_top key_iomput key_mpp_mpi 
  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/SHARED/namelist_ref

    r9532 r9546  
    783783/ 
    784784!----------------------------------------------------------------------- 
    785 &namtra_mle    !   mixed layer eddy parametrisation (Fox-Kemper)        (default: OFF) 
     785&namtra_mle    !   mixed layer eddy parametrisation (Fox-Kemper)       (default: OFF) 
    786786!----------------------------------------------------------------------- 
    787787   ln_mle      = .false.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     
    10401040                              !        = 0  constant 10 m length scale 
    10411041                              !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
     1042      rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4    
    10421043/ 
    10431044!----------------------------------------------------------------------- 
  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/cfg.txt

    r9526 r9546  
    1010ORCA2_LIM3 OPA_SRC LIM_SRC_3 NST_SRC 
    1111AGRIF_NORDIC OPA_SRC LIM_SRC_3 NST_SRC 
     12GYRE_PISCES OPA_SRC TOP_SRC 
    1213ORCA2_LIM3_PISCES OPA_SRC LIM_SRC_3 TOP_SRC NST_SRC 
    13 GYRE_PISCES OPA_SRC TOP_SRC 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90

    r9490 r9546  
    145145   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   ht_0            ,    ht_n             !: t-depth              [m] 
    146146   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   hu_0  ,    hu_b ,    hu_n ,    hu_a   !: u-depth              [m] 
    147    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   hv_0  ,    hv_b ,    hv_n ,    hv_a   !: u-depth              [m] 
     147   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   hv_0  ,    hv_b ,    hv_n ,    hv_a   !: v-depth              [m] 
    148148   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::           r1_hu_b , r1_hu_n , r1_hu_a   !: inverse of u-depth [1/m] 
    149149   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::           r1_hv_b , r1_hv_n , r1_hv_a   !: inverse of v-depth [1/m] 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

    r9440 r9546  
    7575   INTEGER  ::      nn_htau   ! type of tke profile of penetration (=0/1) 
    7676   REAL(wp) ::      rn_efr    ! fraction of TKE surface value which penetrates in the ocean 
     77   REAL(wp) ::      rn_eice   ! =0 ON below sea-ice, =4 OFF when ice fraction > 1/4    
    7778   LOGICAL  ::   ln_lc     ! Langmuir cells (LC) as a source term of TKE or not 
    7879   REAL(wp) ::      rn_lc     ! coef to compute vertical velocity of Langmuir cells 
     
    298299                  !                                           ! TKE Langmuir circulation source term 
    299300                  en(ji,jj,jk) = en(ji,jj,jk) + rdt * MAX(0.,1._wp - 4.*fr_i(ji,jj) ) * ( zwlc * zwlc * zwlc )   & 
    300                      &                              / zhlc(ji,jj) * wmask(ji,jj,jk) 
    301 !!gm                     &                          / zhlc(ji,jj) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     301                     &                              / zhlc(ji,jj) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    302302               END DO 
    303303            END DO 
     
    397397       
    398398      IF( nn_etau == 1 ) THEN           !* penetration below the mixed layer (rn_efr fraction) 
    399          DO jk = 2, jpkm1 
     399         DO jk = 2, jpkm1                       ! rn_eice =0 ON below sea-ice, =4 OFF when ice fraction > 0.25 
    400400            DO jj = 2, jpjm1 
    401401               DO ji = fs_2, fs_jpim1   ! vector opt. 
    402402                  en(ji,jj,jk) = en(ji,jj,jk) + rn_efr * en(ji,jj,1) * EXP( -pdepw(ji,jj,jk) / htau(ji,jj) )   & 
    403                      &                                 * MAX(0.,1._wp - 4.*fr_i(ji,jj) )  * wmask(ji,jj,jk) 
    404 !!gm                     &                             * MAX(0.,1._wp - 4.*fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    405 !!gm    multiplication by surface tmask useless (already includes in MAX( 0, 1-4*fr_i ) 
     403                     &                                 * MAX(0.,1._wp - rn_eice *fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    406404               END DO 
    407405            END DO 
     
    412410               jk = nmln(ji,jj) 
    413411               en(ji,jj,jk) = en(ji,jj,jk) + rn_efr * en(ji,jj,1) * EXP( -pdepw(ji,jj,jk) / htau(ji,jj) )   & 
    414                   &                                 * MAX(0.,1._wp - 4.*fr_i(ji,jj) )  * wmask(ji,jj,jk) 
    415 !!gm                  &                             * MAX(0.,1._wp - 4.*fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     412                  &                                 * MAX(0.,1._wp - rn_eice *fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    416413            END DO 
    417414         END DO 
     
    426423                  zdif = rhftau_scl * MAX( 0._wp, zdif + rhftau_add )  ! apply some modifications... 
    427424                  en(ji,jj,jk) = en(ji,jj,jk) + zbbrau * zdif * EXP( -pdepw(ji,jj,jk) / htau(ji,jj) )   & 
    428                      &                        * MAX(0.,1._wp - 4.*fr_i(ji,jj) ) * wmask(ji,jj,jk) 
    429 !!gm                     &                    * MAX(0.,1._wp - 4.*fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     425                     &                        * MAX(0.,1._wp - rn_eice *fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    430426               END DO 
    431427            END DO 
     
    683679         WRITE(numout,*) '          type of tke penetration profile            nn_htau   = ', nn_htau 
    684680         WRITE(numout,*) '          fraction of TKE that penetrates            rn_efr    = ', rn_efr 
     681         WRITE(numout,*) '          below sea-ice:  =0 ON                      rn_eice   = ', rn_eice 
     682         WRITE(numout,*) '          =4 OFF when ice fraction > 1/4   ' 
    685683         IF( ln_drg ) THEN 
    686684            WRITE(numout,*) 
Note: See TracChangeset for help on using the changeset viewer.