!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !! ESIM namelist: !! 1 - Generic parameters (namice_run) !! 2 - Ice thickness discretization (namice_itd) !! 3 - Ice dynamics (namice_dyn) !! 4 - Ice ridging/rafting (namice_rdgrft) !! 5 - Ice rheology (namice_rhg) !! 6 - Ice advection (namice_adv) !! 7 - Ice thermodynamics (namice_thd) !! 8 - Ice salinity (namice_sal) !! 9 - Ice melt ponds (namice_mp) !! 10 - Ice initialization (namice_ini) !! 11 - Ice/snow albedos (namice_alb) !! 12 - Ice diagnostics (namice_dia) !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ! !------------------------------------------------------------------------------ &namice_run ! Generic parameters !------------------------------------------------------------------------------ jpl = 5 ! number of ice categories nlay_i = 2 ! number of ice layers nlay_s = 1 ! number of snow layers (only 1 is working) nn_monocat = 0 ! virtual ITD mono-category parameterizations (1-4 => jpl = 1 only) or not (0) ! 2: simple piling instead of ridging --- temporary option ! 3: activate G(he) only --- temporary option ! 4: activate extra lateral melting only --- temporary option ln_icedyn = .true. ! ice dynamics (T) or not (F) ln_icethd = .true. ! ice thermo (T) or not (F) rn_amax_n = 0.997 ! maximum tolerated ice concentration NH rn_amax_s = 0.997 ! maximum tolerated ice concentration SH cn_icerst_in = "restart_ice" ! suffix of ice restart name (input) cn_icerst_out = "restart_ice" ! suffix of ice restart name (output) cn_icerst_indir = "." ! directory to read input ice restarts cn_icerst_outdir = "." ! directory to write output ice restarts / !------------------------------------------------------------------------------ &namice_itd ! Ice discretization !------------------------------------------------------------------------------ rn_himean = 2.0 ! expected domain-average ice thickness (m) rn_himin = 0.1 ! minimum ice thickness (m) used in remapping / !------------------------------------------------------------------------------ &namice_dyn ! Ice dynamics !------------------------------------------------------------------------------ ln_dynFULL = .true. ! dyn.: full ice dynamics (rheology + advection + ridging/rafting + correction) ln_dynRHGADV = .false. ! dyn.: no ridge/raft & no corrections (rheology + advection) ln_dynADV = .false. ! dyn.: only advection w prescribed vel.(rn_uvice + advection) rn_uice = 0.00001 ! prescribed ice u-velocity rn_vice = 0. ! prescribed ice v-velocity rn_ishlat = 2. ! free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2) rn_cio = 5.0e-03 ! ice-ocean drag coefficient (-) ln_landfast = .false. ! landfast ice parameterization (T or F) rn_gamma = 0.15 ! fraction of ocean depth that ice must reach to initiate landfast ! recommended range: [0.1 ; 0.25] rn_icebfr = 10. ! maximum bottom stress per unit area of contact [N/m2] ! a very large value ensures ice velocity=0 even with a small contact area ! recommended range: ?? (should be greater than atm-ice stress => >0.1 N/m2) rn_lfrelax = 1.e-5 ! relaxation time scale to reach static friction [s-1] / !------------------------------------------------------------------------------ &namice_rdgrft ! Ice ridging/rafting !------------------------------------------------------------------------------ ! -- ice_rdgrft_strength -- ! ln_str_H79 = .true. ! ice strength param.: Hibler_79 => P = pstar**exp(-c_rhg*A) rn_pstar = 2.0e+04 ! ice strength thickness parameter [N/m2] rn_crhg = 20.0 ! ice strength conc. parameter (-) ln_str_R75 = .false. ! ice strength param.: Rothrock_75 => P = Cf*coeff*integral(wr.h^2) rn_perdg = 17.0 ! ridging work divided by pot. energy change in ridging ! -- ice_rdgrft -- ! rn_csrdg = 0.5 ! fraction of shearing energy contributing to ridging ! -- ice_rdgrft_prep -- ! ln_partf_lin = .false. ! Linear ridging participation function (Thorndike et al, 1975) rn_gstar = 0.15 ! fractional area of thin ice being ridged ln_partf_exp = .true. ! Exponential ridging participation function (Lipscomb, 2007) rn_astar = 0.03 ! exponential measure of ridging ice fraction [set to 0.05 if hstar=100] ln_ridging = .true. ! ridging activated (T) or not (F) rn_hstar = 25.0 ! determines the maximum thickness of ridged ice [m] (Hibler, 1980) rn_porordg = 0.3 ! porosity of newly ridged ice (Lepparanta et al., 1995) rn_fsnwrdg = 0.5 ! snow volume fraction that survives in ridging rn_fpndrdg = 1.0 ! pond fraction that survives in ridging (small a priori) ln_rafting = .true. ! rafting activated (T) or not (F) rn_hraft = 0.75 ! threshold thickness for rafting [m] rn_craft = 5.0 ! squeezing coefficient used in the rafting function rn_fsnwrft = 0.5 ! snow volume fraction that survives in rafting rn_fpndrft = 1.0 ! pond fraction that survives in rafting (0.5 a priori) / !------------------------------------------------------------------------------ &namice_rhg ! Ice rheology !------------------------------------------------------------------------------ ln_rhg_EVP = .true. ! EVP rheology rn_creepl = 1.0e-12 ! 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 (rn_nevp=120) or 1/9 (rn_nevp=300) / !------------------------------------------------------------------------------ &namice_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) / !------------------------------------------------------------------------------ &namice_thd ! Ice thermodynamics !------------------------------------------------------------------------------ ! -- icethd_dif -- ! rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice [1/m] ln_cndi_U64 = .false. ! sea ice thermal conductivity: k = k0 + beta.S/T (Untersteiner, 1964) ln_cndi_P07 = .true. ! sea ice thermal conductivity: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007) ln_dqns_i = .true. ! change the surface non-solar flux with surface temperature (T) or not (F) rn_cnd_s = 0.31 ! thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971) ! Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) ! -- icethd_dh -- ! ln_icedH = .true. ! activate ice thickness change from growing/melting (T) or not (F) rn_blow_s = 0.66 ! mesure of snow blowing into the leads ! = 1 => no snow blowing, < 1 => some snow blowing ! -- icethd_da -- ! ln_icedA = .true. ! activate lateral melting param. (T) or not (F) rn_beta = 1.0 ! coef. beta for lateral melting param. Recommended range=[0.8-1.2] ! => decrease = more melt and melt peaks toward higher concentration (A~0.5 for beta=1 ; A~0.8 for beta=0.2) ! 0.3 = best fit for western Fram Strait and Antarctica ! 1.4 = best fit for eastern Fram Strait rn_dmin = 8. ! minimum floe diameter for lateral melting param. Recommended range=[6-10] ! => 6 vs 8m = +40% melting at the peak (A~0.5) ! 10 vs 8m = -20% melting ! -- icethd_lac -- ! ln_icedO = .true. ! activate ice growth in open-water (T) or not (F) rn_hinew = 0.1 ! thickness for new ice formation in open water (m), must be larger than rn_hnewice ln_frazil = .false. ! Frazil ice parameterization (ice collection as a function of wind) rn_maxfraz = 1.0 ! maximum fraction of frazil ice collecting at the ice base rn_vfraz = 0.417 ! thresold drift speed for frazil ice collecting at the ice bottom (m/s) rn_Cfraz = 5.0 ! squeezing coefficient for frazil ice collecting at the ice bottom ! -- icestp -- ! nn_iceflx = -1 ! Redistribute heat flux over ice categories ! =-1 Do nothing (needs N(cat) fluxes) ! ==> forced mode only ! = 0 Average N(cat) fluxes then apply the average over the N(cat) ice ! ==> forced and coupled modes ! = 1 Average N(cat) fluxes then redistribute over the N(cat) ice using T-ice and albedo sensitivity ! ==> forced mode only ! = 2 Redistribute a single flux over categories ! ==> coupled mode only / !------------------------------------------------------------------------------ &namice_sal ! Ice salinity !------------------------------------------------------------------------------ ! -- icethd_sal -- ! ln_icedS = .true. ! activate gravity drainage and flushing (T) or not (F) nn_icesal = 2 ! ice salinity option ! 1: constant ice salinity (S=rn_icesal) ! 2: varying salinity parameterization S(z,t) ! 3: prescribed salinity profile S(z), Schwarzacher, 1959 rn_icesal = 4. ! (nn_icesal=1) ice salinity (g/kg) rn_sal_gd = 5. ! restoring ice salinity, gravity drainage (g/kg) rn_time_gd = 1.73e+6 ! restoring time scale, gravity drainage (s) rn_sal_fl = 2. ! restoring ice salinity, flushing (g/kg) rn_time_fl = 8.64e+5 ! restoring time scale, flushing (s) rn_simax = 20. ! maximum tolerated ice salinity (g/kg) rn_simin = 0.1 ! minimum tolerated ice salinity (g/kg) / !------------------------------------------------------------------------------ &namicemp ! Melt ponds !------------------------------------------------------------------------------ ! -- limmp -- ! ln_pnd = .false. ! active melt ponds ln_pnd_rad = .false. ! active melt ponds radiative coupling ln_pnd_fw = .false. ! active melt ponds freshwater coupling nn_pnd_scheme = 0 ! type of melt pond scheme : =0 prescribed ( Tsu=0 ), =1 empirical, =2 topographic rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 : (0