Changeset 11896


Ignore:
Timestamp:
2019-11-13T12:40:52+01:00 (7 months ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo: update ISOMIP+ test case configuration

Location:
NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/EXPREF/namelist_cfg

    r11889 r11896  
    6262   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    6363   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    64    sn_tem = 'TS_init',  -12.    , 'Tinit',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    65    sn_sal = 'TS_init',  -12.    , 'Sinit',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     64   sn_tem  = 'nemo_base_WARM.nc'         ,  -12.    , 'Tinit',   .false.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     65   sn_sal  = 'nemo_base_WARM.nc'         ,  -12.    , 'Sinit',   .false.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     66 
     67   sn_dmpt = 'nemo_base_WARM.nc'         ,  -12.    , 'Tinit',   .false.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     68   sn_dmps = 'nemo_base_WARM.nc'         ,  -12.    , 'Sinit',   .false.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    6669/ 
    6770!----------------------------------------------------------------------- 
     
    111114 
    112115   ln_usr      = .true.   !  user defined formulation                  (T => check usrdef_sbc) 
     116   nn_fwb      = 1 
    113117/ 
    114118!----------------------------------------------------------------------- 
     
    181185         !                       ! ad15 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    182186         !                       ! hj99 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    183          rn_gammat0  = 0.045     ! gammat coefficient used in blk formula 
    184          rn_gammas0  = 1.29e-3   ! gammas coefficient used in blk formula 
     187         rn_gammat0  = 0.0215    ! gammat coefficient used in blk formula 
     188         rn_gammas0  = 0.614e-3  ! gammas coefficient used in blk formula 
    185189         rn_vtide    = 0.01      ! tidal velocity [m/s] 
    186190         ! 
     
    227231   !                       !  free slip  !   partial slip  !   no slip   ! strong slip 
    228232   rn_shlat    =  2.       !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    229    ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs. 
    230233/ 
    231234!----------------------------------------------------------------------- 
     
    270273!----------------------------------------------------------------------- 
    271274   rn_Cd0      =  2.5e-3   !  drag coefficient [-] 
     275   rn_ke0      =  0.0e-3   !  background kinetic energy  [m2/s2] (non-linear cases) 
    272276/ 
    273277!----------------------------------------------------------------------- 
     
    275279!----------------------------------------------------------------------- 
    276280   rn_Cd0      =  2.5e-3    !  drag coefficient [-] 
     281   rn_ke0      =  0.0e-3    !  background kinetic energy  [m2/s2] (non-linear cases) 
    277282/ 
    278283!----------------------------------------------------------------------- 
     
    300305&nameos        !   ocean Equation Of Seawater                           (default: NO selection) 
    301306!----------------------------------------------------------------------- 
    302    ln_eos80    = .true.          !  = Use EOS80 
     307   ln_teos10   = .false.         !  = Use TEOS-10 
     308   ln_eos80    = .false.         !  = Use EOS80 
     309   ln_leos     = .true.          !  = Use S-EOS (simplified Eq.) 
     310                                 ! 
     311   !                     ! S-EOS coefficients (ln_seos=T): 
     312   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     313   !                     ! L-EOS coefficients (ln_seos=T): 
     314   !                             !  rd(T,S,Z)*rau0 = rau0*(-a0*dT+b0*dS) 
     315   rn_a0       =  3.7330e-5      !  thermal expension coefficient 
     316   rn_b0       =  7.8430e-4      !  saline  expension coefficient 
    303317/ 
    304318!----------------------------------------------------------------------- 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/MY_SRC/sbcfwb.F90

    r11889 r11896  
    105105         ! 
    106106         IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 
     107            y_fwfnow(1) = local_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) - snwice_fmass(:,:) ) ) 
     108            CALL mpp_delay_sum( 'sbcfwb', 'fwb', y_fwfnow(:), z_fwfprv(:), kt == nitend - nn_fsbc + 1 ) 
     109            z_fwfprv(1) = z_fwfprv(1) / area 
     110            zcoef = z_fwfprv(1) * rcp 
     111            emp(:,:) = emp(:,:) - z_fwfprv(1)        * tmask(:,:,1) 
     112            qns(:,:) = qns(:,:) + zcoef * sst_m(:,:) * tmask(:,:,1) ! account for change to the heat budget due to fw correction 
     113         ENDIF 
     114         ! 
     115      CASE ( 4 )                             !==  global mean fwf set to zero (ISOMIP case) ==! 
     116         ! 
     117         IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 
    107118            z_fwf = glob_sum( 'sbcfwb',  e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) - snwice_fmass(:,:) ) ) 
    108119            z_fwf = z_fwf / area 
    109120            zcoef = z_fwf * rcp 
    110             emp(:,:) = emp(:,:) - z_fwf              * tmask(:,:,1) 
    111             qns(:,:) = qns(:,:) + zcoef * sst_m(:,:) * tmask(:,:,1) ! account for change to the heat budget due to fw correction 
    112             sfx(:,:) = sfx(:,:) + z_fwf * sss_m(:,:) * tmask(:,:,1) 
     121            emp(:,:) = emp(:,:) - z_fwf            * tmask(:,:,1) ! (Eq. 34 AD2015) 
     122            qns(:,:) = qns(:,:) + zcoef * ( -1.9 ) * tmask(:,:,1) ! (Eq. 35 AD2015) ! could be sst_m if we don't want any bouyancy fluxes 
     123            sfx(:,:) = sfx(:,:) + z_fwf * ( 33.8 ) * tmask(:,:,1) ! (Eq. 36 AD2015) ! could be sss_m if we don't want any bouyancy fluxes 
    113124         ENDIF 
    114125         ! 
Note: See TracChangeset for help on using the changeset viewer.