!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! SI3 configuration namelist: Overwrites SHARED/namelist_ice_ref
!!              1 - Generic parameters                 (nampar)
!!              2 - Ice thickness discretization       (namitd)
!!              3 - Ice dynamics                       (namdyn)
!!              4 - Ice ridging/rafting                (namdyn_rdgrft)
!!              5 - Ice rheology                       (namdyn_rhg)
!!              6 - Ice advection                      (namdyn_adv)
!!              7 - Ice surface boundary conditions    (namsbc)
!!              8 - Ice thermodynamics                 (namthd)
!!              9 - Ice heat diffusion                 (namthd_zdf)
!!             10 - Ice lateral melting                (namthd_da)
!!             11 - Ice growth in open water           (namthd_do)
!!             12 - Ice salinity                       (namthd_sal)
!!             13 - Ice melt ponds                     (namthd_pnd)
!!             14 - Ice initialization                 (namini)
!!             15 - Ice/snow albedos                   (namalb)
!!             16 - Ice diagnostics                    (namdia)
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!------------------------------------------------------------------------------
&nampar         !   Generic parameters
!------------------------------------------------------------------------------
   jpl              =   5 ! 10         !  number of ice  categories
   nlay_i           =   3 ! 5          !  number of ice  layers
/
!------------------------------------------------------------------------------
&namitd         !   Ice discretization
!------------------------------------------------------------------------------
   ln_cat_hfn       = .true.          !  ice categories are defined by a function following rn_himean**(-0.05)
      rn_himean     =   1.5 !2.0           !  expected domain-average ice thickness (m)
   rn_himin         =   0.01           !  minimum ice thickness (m) allowed
/
!------------------------------------------------------------------------------
&namdyn         !   Ice dynamics
!------------------------------------------------------------------------------
   ln_landfast_L16  = .true.          !  landfast: parameterization from Lemieux 2016
      rn_lf_tensile =   0.2 !0.05          !        isotropic tensile strength [0-0.5??]
/
!------------------------------------------------------------------------------
&namdyn_rdgrft  !   Ice ridging/rafting
!------------------------------------------------------------------------------
          ! -- ice_rdgrft_strength -- !
   ln_str_H79       = .true.          !  ice strength param.: Hibler_79   => P = pstar*<h>*exp(-c_rhg*A)                      
      rn_pstar      =   2.7e+04 !3.5e+04 !2.0e+04       !     ice strength thickness parameter [N/m2]
      rn_crhg       =   20.0 !30.0 !20.0          !     ice strength conc. parameter (-)
/
!------------------------------------------------------------------------------
&namdyn_rhg     !   Ice rheology
!------------------------------------------------------------------------------
   ln_rhg_EVP       = .true.          !  EVP rheology
      ln_aEVP       = .false.          !     adaptive rheology (Kimmritz et al. 2016 & 2017)
      rn_creepl     =   2.0e-9        !     creep limit [1/s]
      rn_ecc        =   2.0           !     eccentricity of the elliptical yield curve          
      nn_nevp       = 120             !     number of EVP subcycles                             
      rn_relast     =   0.333         !     ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast 
                                      !        advised value: 1/3 (nn_nevp=100) or 1/9 (nn_nevp=300)
   nn_rhg_chkcvg    =   0             !  check convergence of rheology
                                      !     = 0  no check
                                      !     = 1  check at the main time step (output xml: uice_cvg)
                                      !     = 2  check at both main and rheology time steps (additional output: ice_cvg.nc)
                                      !          this option 2 asks a lot of communications between cpu
/
!------------------------------------------------------------------------------
&namdyn_adv     !   Ice advection
!------------------------------------------------------------------------------
   ln_adv_Pra       = .false.         !  Advection scheme (Prather)
   ln_adv_UMx       = .true.          !  Advection scheme (Ultimate-Macho)
      nn_UMx        =   5             !     order of the scheme for UMx (1-5 ; 20=centered 2nd order)
/
!------------------------------------------------------------------------------
&namsbc         !   Ice surface boundary conditions
!------------------------------------------------------------------------------
   nn_snwfra        =   1 !2             !  calculate the fraction of ice covered by snow (for zdf and albedo)
   nn_qtrice        =   0 !1             !  Solar flux transmitted thru the surface scattering layer:
                                      !     = 0  Grenfell and Maykut 1977 (depends on cloudiness and is 0 when there is snow) 
                                      !     = 1  Lebrun 2019 (equals 0.3 anytime with different melting/dry snw conductivities)
/
!------------------------------------------------------------------------------
&namthd         !   Ice thermodynamics
!------------------------------------------------------------------------------
   ln_icedA         = .false.          !  activate lateral melting param. (T) or not (F)
   ln_leadhfx       = .false.          !  heat in the leads is used to melt sea-ice before warming the ocean
/
!------------------------------------------------------------------------------
&namthd_zdf     !   Ice heat diffusion
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namthd_da      !   Ice lateral melting
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namthd_do      !   Ice growth in open water
!------------------------------------------------------------------------------
   rn_hinew         =   0.1 !0.02          !  thickness for new ice formation in open water (m), must be larger than rn_himin
   ln_frazil        = .false. !.true.          !  Frazil ice parameterization (ice collection as a function of wind)
/
!------------------------------------------------------------------------------
&namthd_sal     !   Ice salinity
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namthd_pnd     !   Melt ponds
!------------------------------------------------------------------------------
   ln_pnd           = .false. !.true.          !  activate melt ponds or not
      ln_pnd_LEV     = .false.         !  level ice melt ponds (from Flocco et al 2007,2010 & Holland et al 2012)
/
!------------------------------------------------------------------------------
&namini         !   Ice initialization
!------------------------------------------------------------------------------
   ln_iceini        = .true.         !  activate ice initialization (T) or not (F)
   nn_iceini_file   =   0             !     0 = Initialise sea ice based on SSTs
                                      !     1 = Initialise sea ice from single category netcdf file
                                      !     2 = Initialise sea ice from multi category restart file
   rn_thres_sst     =   0.5           !  max delta temp. above Tfreeze with initial ice = (sst - tfreeze)
   rn_hts_ini_n     =   0.1           !  initial real snow thickness (m), North
   rn_hti_ini_n     =   1.0           !  initial real ice thickness  (m), North
   rn_ati_ini_n     =   0.9           !  initial ice concentration   (-), North
   rn_smi_ini_n     =   6.0           !  initial ice salinity     (g/kg), North
   rn_tmi_ini_n     = 267.            !  initial ice/snw temperature (K), North
   ! -- for nn_iceini_file = 1
!           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
!           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
   sn_hti = 'iicethic_ORAS5'        , -1 ,'iicethic' ,  .false.  , .false., 'monthly'  , '' , '', ''
   sn_hts = 'isnowinvent_ORAS5'     , -1  ,'isnow'   ,  .false.  , .false., 'monthly'  , '' , '', ''
   sn_ati = 'ileadfra_ORAS5'        , -1  ,'ileadfra',  .false.  , .false., 'monthly'  , '' , '', ''
   sn_smi = 'NOT USED'              , -12 ,'smi'     ,  .false.  , .true., 'yearly'  , '' , '', ''
   sn_tmi = 'NOT USED'              , -12 ,'tmi'     ,  .false.  , .true., 'yearly'  , '' , '', ''
   sn_tsu = 'NOT USED'              , -12 ,'tsu'     ,  .false.  , .true., 'yearly'  , '' , '', ''
   sn_tms = 'NOT USED'              , -12 ,'tms'     ,  .false.  , .true., 'yearly'  , '' , '', ''
   !      melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd)
   sn_apd = 'NOT USED'              , -12 ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', ''
   sn_hpd = 'NOT USED'              , -12 ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', ''
   sn_hld = 'NOT USED'              , -12 ,'hld'   ,  .false.  , .true., 'yearly'  , '' , '', ''
   cn_dir='./'
/
!------------------------------------------------------------------------------
&namalb         !   albedo parameters
!------------------------------------------------------------------------------
   !                                  !                          !  obs range (cloud-sky)
   rn_alb_sdry      =  0.87 !0.87 !0.87 !0.85          !  dry snow albedo         :  0.85 -- 0.87
   rn_alb_smlt      =  0.82 !0.82 !0.82 !0.75          !  melting snow albedo     :  0.72 -- 0.82
   rn_alb_idry      =  0.65 !0.72 !0.65 !0.60          !  dry ice albedo          :  0.54 -- 0.65
   rn_alb_imlt      =  0.61 !0.65 !0.58 !0.50          !  bare puddled ice albedo :  0.49 -- 0.58
   rn_alb_dpnd      =  0.61 !0.65 !0.58 !0.27          !  ponded ice albedo       :  0.10 -- 0.30 
/
!------------------------------------------------------------------------------
&namdia         !   Diagnostics
!------------------------------------------------------------------------------
/