!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !! NEMO/LIM3 : 1 - dynamics/advection/thermo (namicerun) !! namelists 2 - ice intialisation (namiceini) !! 3 - ice dynamic (namicedyn) !! 4 - ice advection (namicetrp) !! 5 - thermodynamic (namicethd) !! 6 - ice salinity (namicesal) !! 7 - mechanical redistribution of ice (namiceitdme) !! 8 - ice diagnostics (namicedia) !! 9 - ice outputs (namiceout) !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !----------------------------------------------------------------------- &namicerun ! Share parameters for dynamics/advection/thermo !----------------------------------------------------------------------- cn_icerst_in = "restart_ice" ! suffix of ice restart name (input) cn_icerst_out = "restart_ice" ! suffix of ice restart name (output) ln_limdyn = .true. ! ice dynamics (T) or thermodynamics only (F) amax = 0.999 ! maximum ice concentration ln_nicep = .false. ! Ice points output for debug (yes or no) ln_limdiahsb = .false. ! check the heat and salt budgets (T) or not (F) ln_limdiaout = .true. ! output the heat and salt budgets (T) or not (F) / !----------------------------------------------------------------------- &namiceini ! ice initialisation !----------------------------------------------------------------------- ln_limini = .false. ! activate ice initialization (T) or not (F) thres_sst = 0.0 ! threshold water temperature for initial sea ice hts_ini_n = 0.3 ! initial snow thickness in the north hts_ini_s = 0.3 ! " " south hti_ini_n = 3.0 ! initial ice thickness in the north hti_ini_s = 1.0 ! " " south ati_ini_n = 0.9 ! initial ice concentration in the north ati_ini_s = 0.9 ! " " south smi_ini_n = 6.3 ! initial ice salinity in the north smi_ini_s = 6.3 ! " " south tmi_ini_n = 270. ! initial ice/snw temp in the north tmi_ini_s = 270. ! initial ice/snw temp in the south / !----------------------------------------------------------------------- &namicedyn ! ice dynamic !----------------------------------------------------------------------- cw = 5.0e-03 ! drag coefficient for oceanic stress pstar = 2.0e+04 ! 1st bulk-rheology parameter c_rhg = 20.0 ! 2nd bulk-rhelogy parameter creepl = 1.0e-12 ! creep limit ecc = 2.0 ! eccentricity of the elliptical yield curve ahi0 = 350.e0 ! horizontal eddy diffusivity coefficient for sea-ice [m2/s] ! depends on resolution nevp = 120 ! number of iterations for subcycling in EVP relast = 0.333 ! ratio of elastic timescale over ice time step (1/3 if nevp=120 ; 1/9 if nevp=300) hminrhg = 0.001 ! ice volume (a*h in m) below which ice velocity equal ocean velocity / !----------------------------------------------------------------------- &namicethd ! ice thermodynamic !----------------------------------------------------------------------- hiccrit = 0.1 ! ice thickness for lateral accretion ! caution 1.0, 1.0 best value to be used!!! (gilles G.) ???? fraz_swi = 0 ! use of frazil ice collection thickness in function of wind (1.0) or not (0.0) maxfrazb = 0.0 ! maximum portion of frazil ice collecting at the ice bottom vfrazb = 0.4166667 ! thresold drift speed for frazil ice collecting at the ice bottom Cfrazb = 5.0 ! squeezing coefficient for frazil ice collecting at the ice bottom hiclim = 0.10 ! minimum thickness of the first thickness category, must be smaller than hiccrit betas = 0.6 ! exponent in lead-ice fractionation of snow precipitation 0.66 ! betas = 1 -> equipartition, betas < 1 -> more on leads kappa_i = 1.0 ! extinction radiation parameter in sea ice (1.0) nconv_i_thd = 50 ! maximal number of iterations for heat diffusion computation maxer_i_thd = 0.0001 ! maximal error in temperature for heat diffusion computation thcon_i_swi = 1 ! switch for computation of thermal conductivity in the ice ! (0) Untersteiner (1964), (1) Pringle et al. (2007) / !----------------------------------------------------------------------- &namicesal ! ice salinity !----------------------------------------------------------------------- num_sal = 2 ! salinity option: 1 -> S = bulk_sal ! 2 -> S = S(z,t) with a simple parameterization ! 3 -> S = S(z) profile of Scwharzacher [1959] ! 4 -> S = S(h) Cox and Weeks [1974] bulk_sal = 4.0 ! if 1 is used, it represents the ice salinity sal_G = 5.00 ! restoring salinity for GD time_G = 1.728e+6 ! restoring time for GD sal_F = 2.00 ! restoring salinity for flushing time_F = 8.640e+5 ! restoring time for flushing s_i_max = 20.0 ! Maximum salinity s_i_min = 0.1 ! Minimum tolerated ice salinity / !----------------------------------------------------------------------- &namiceitdme ! parameters for mechanical redistribution of ice !----------------------------------------------------------------------- ridge_scheme_swi = 0 ! which ridging scheme using (1=Rothrock,else=Hibler79) Cs = 0.50 ! shearing energy contribution to ridging Cf = 17.0 ! ratio of ridging work to PE change in ridging fsnowrdg = 0.5 ! snow fraction that survives in ridging fsnowrft = 0.5 ! snow fraction that survives in rafting Gstar = 0.15 ! fractional area of thin ice being ridged astar = 0.05 ! equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice) Hstar = 100.0 ! parameter determining the maximum thickness of ridged ice raft_swi = 1 ! rafting or not hparmeter = 0.75 ! threshold thickness for rafting or not Craft = 5.0 ! coefficient used in the rafting function ridge_por = 0.3 ! initial porosity of the ridged ice (typically 0.30) partfun_swi = 1 ! participation function linear, TH75 (0) or exponential Letal07 (1) brinstren_swi = 0 ! (1) use brine volume to diminish ice strength