Changeset 5048 for branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref
- Timestamp:
- 2015-02-02T11:28:50+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref
r5046 r5048 2 2 !! NEMO/LIM3 : 1 - dynamics/advection/thermo (namicerun) 3 3 !! namelists 2 - ice intialisation (namiceini) 4 !! 3 - ice dynamic 4 !! 3 - ice dynamics (namicedyn) 5 5 !! 4 - ice advection (namicetrp) 6 !! 5 - thermodynamic 6 !! 5 - thermodynamics (namicethd) 7 7 !! 6 - ice salinity (namicesal) 8 8 !! 7 - mechanical redistribution of ice (namiceitdme) … … 17 17 cn_icerst_out = "restart_ice" ! suffix of ice restart name (output) 18 18 ln_limdyn = .true. ! ice dynamics (T) or thermodynamics only (F) 19 amax = 0.999 ! maximum ice concentration20 ln_nicep = .false. ! Ice points output for debug (yes or no)21 ln_limdiahsb = .false. 19 amax = 0.999 ! maximum tolerated ice concentration 20 ln_nicep = .false. ! ice points output for debug (yes or no) 21 ln_limdiahsb = .false. ! check the heat and salt budgets (T) or not (F) 22 22 ln_limdiaout = .true. ! output the heat and salt budgets (T) or not (F) 23 23 / 24 24 !----------------------------------------------------------------------- 25 &namiceini ! ice initiali sation25 &namiceini ! ice initialization 26 26 !----------------------------------------------------------------------- 27 27 ln_limini = .false. ! activate ice initialization (T) or not (F) 28 28 thres_sst = 0.0 ! threshold water temperature for initial sea ice 29 hts_ini_n = 0.3 ! initial snow thickness in the north30 hts_ini_s = 0.3 ! " " south31 hti_ini_n = 3.0 ! initial ice thickness in the north32 hti_ini_s = 1.0 ! " " south33 ati_ini_n = 0.9 ! initial ice concentration in the north34 ati_ini_s = 0.9 ! " " south35 smi_ini_n = 6.3 ! initial ice salinity in the north36 smi_ini_s = 6.3 ! " " south37 tmi_ini_n = 270. ! initial ice/snw temp in the north38 tmi_ini_s = 270. ! initial ice/snw temp in the south29 hts_ini_n = 0.3 ! initial snow thickness (North) 30 hts_ini_s = 0.3 ! " " (South) 31 hti_ini_n = 3.0 ! initial ice thickness (North) 32 hti_ini_s = 1.0 ! " " (South) 33 ati_ini_n = 0.9 ! initial ice concentration (North) 34 ati_ini_s = 0.9 ! " " (South) 35 smi_ini_n = 6.3 ! initial ice salinity (North) 36 smi_ini_s = 6.3 ! " " (South) 37 tmi_ini_n = 270. ! initial ice/snw temp (North) 38 tmi_ini_s = 270. ! initial ice/snw temp (South) 39 39 / 40 40 !----------------------------------------------------------------------- 41 &namicedyn ! ice dynamic 41 &namicedyn ! ice dynamics 42 42 !----------------------------------------------------------------------- 43 cw = 5.0e-03 ! drag coefficient for oceanic stress44 pstar = 2.0e+04 ! 1st bulk-rheology parameter 43 cw = 5.0e-03 ! ice-ocean drag coefficient 44 pstar = 2.0e+04 ! 1st bulk-rheology parameter (N/m) 45 45 c_rhg = 20.0 ! 2nd bulk-rhelogy parameter 46 creepl = 1.0e-12 ! creep limit 46 creepl = 1.0e-12 ! creep limit (s-1) 47 47 ecc = 2.0 ! eccentricity of the elliptical yield curve 48 ahi0 = 350.e0 ! horizontal eddy diffusivity coefficient for sea-ice [m2/s] !depends on resolution49 nevp = 120 ! number of iterations for subcycling in EVP50 relast = 0.333 ! ratio of elastic timescale over ice time step (1/3 if nevp=120 ;1/9 if nevp=300)51 hminrhg = 0.001 ! ice volume (a*h in m) below which ice velocity equal ocean velocity48 ahi0 = 350.e0 ! horizontal eddy diffusivity coefficient for sea-ice (m2/s), depends on resolution 49 nevp = 120 ! number of EVP subcycles 50 relast = 0.333 ! ratio of elastic timescale to ice time step (1/3 if nevp=120, 1/9 if nevp=300) 51 hminrhg = 0.001 ! ice volume (a*h in m) below which ice velocity equals ocean velocity 52 52 / 53 53 !----------------------------------------------------------------------- 54 &namicethd ! ice thermodynamic 54 &namicethd ! ice thermodynamics 55 55 !----------------------------------------------------------------------- 56 hiccrit = 0.1 ! ice thickness for lateral accretion 57 ! caution 1.0, 1.0 best value to be used!!! (gilles G.) ???? 58 fraz_swi = 0 ! use of frazil ice collection thickness in function of wind (1.0) or not (0.0) 56 hiccrit = 0.1 ! thickness for new ice formation in open water 57 fraz_swi = 0 ! use frazil ice collection thickness as a function of wind (1.0) or not (0.0) 59 58 maxfrazb = 0.0 ! maximum portion of frazil ice collecting at the ice bottom 60 vfrazb = 0.41 66667! thresold drift speed for frazil ice collecting at the ice bottom59 vfrazb = 0.417 ! thresold drift speed for frazil ice collecting at the ice bottom 61 60 Cfrazb = 5.0 ! squeezing coefficient for frazil ice collecting at the ice bottom 62 61 hiclim = 0.10 ! minimum thickness of the first thickness category, must be smaller than hiccrit 63 62 betas = 0.6 ! exponent in lead-ice fractionation of snow precipitation 0.66 64 63 ! betas = 1 -> equipartition, betas < 1 -> more on leads 65 kappa_i = 1.0 ! extinction radiation parameter in sea ice ( 1.0)64 kappa_i = 1.0 ! extinction radiation parameter in sea ice (m-1) 66 65 nconv_i_thd = 50 ! maximal number of iterations for heat diffusion computation 67 maxer_i_thd = 0.0001 ! maximal error in temperature for heat diffusion computation 68 thcon_i_swi = 1 ! switch for computation of thermal conductivity in the ice 69 ! (0) Untersteiner (1964), (1) Pringle et al. (2007) 66 maxer_i_thd = 0.0001 ! temperature error tolerance after heat diffusion 67 thcon_i_swi = 1 ! sea ice thermal conductivity, Untersteiner_64 (0) or Pringle_et_al_07 (1) 68 nn_monocat = 0 ! virtual ITD mono-category parameterizations (1) or not (0); 1 ok for jpl=1 only 69 ! 2=replace ridging by piling, 3=activate G(he) only, 4=activate lateral melting only) --- those are temporary test options 70 70 / 71 71 !----------------------------------------------------------------------- 72 72 &namicesal ! ice salinity 73 73 !----------------------------------------------------------------------- 74 num_sal = 2 ! salinity option: 1 -> S = bulk_sal 75 ! 2 -> S = S(z,t) with a simple parameterization 76 ! 3 -> S = S(z) profile of Scwharzacher [1959] 77 ! 4 -> S = S(h) Cox and Weeks [1974] 74 num_sal = 2 ! S=bulk_sal (1), S(z,t) - (2), S(z) Schwarzacher_59 (3) 78 75 bulk_sal = 4.0 ! if 1 is used, it represents the ice salinity 79 sal_G = 5.00 ! restoring salinity for GD80 time_G = 1.728e+6 ! restoring time for GD81 sal_F = 2.00 ! restoring salinity for flushing82 time_F = 8.640e+5 ! restoring time for flushing83 s_i_max = 20.0 ! Maximum salinity84 s_i_min = 0.1 ! Minimum tolerated ice salinity76 sal_G = 5.00 ! restoring salinity, gravity drainage (g/kg) 77 time_G = 1.728e+6 ! restoring time scale, gravity drainage (s) 78 sal_F = 2.00 ! restoring salinity, flushing (g/kg) 79 time_F = 8.640e+5 ! restoring time, flushing (s) 80 s_i_max = 20.0 ! maximum tolerated ice salinity (g/kg) 81 s_i_min = 0.1 ! minimum tolerated ice salinity (g/kg) 85 82 / 86 83 !----------------------------------------------------------------------- 87 84 &namiceitdme ! parameters for mechanical redistribution of ice 88 85 !----------------------------------------------------------------------- 89 ridge_scheme_swi = 0 ! which ridging scheme using (1=Rothrock,else=Hibler79)90 Cs = 0.50 ! shearing energy contributionto ridging91 Cf = 17.0 ! ratio of ridging work to PEchange in ridging92 fsnowrdg = 0.5 ! snow fraction that survives in ridging93 fsnowrft = 0.5 ! snow fraction that survives in rafting94 Gstar = 0.15 ! fractional area of thin ice being ridged 95 astar = 0.05 ! equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice)96 Hstar = 100.0 ! parameter determining the maximum thickness of ridged ice97 raft_swi = 1 ! rafting or not98 hparmeter = 0.75 ! threshold thickness for rafting or not86 ridge_scheme_swi = 0 ! ice strength parameteriztaion, Hibler_79 (0), Rothrock_75 (1) 87 Cs = 0.50 ! ratio of shearing energy contributing to ridging 88 Cf = 17.0 ! ratio of ridging work to potential energy change in ridging 89 fsnowrdg = 0.5 ! snow volume fraction that survives in ridging 90 fsnowrft = 0.5 ! snow volume fraction that survives in rafting 91 Gstar = 0.15 ! fractional area of thin ice being ridged (partfun_swi = 0) 92 astar = 0.05 ! measures fractional area of thin ice being ridged (partfun_swi = 1) 93 Hstar = 100.0 ! determines the maximum thickness of ridged ice (m) 94 raft_swi = 1 ! rafting activated (1) or not (0) 95 hparmeter = 0.75 ! threshold thickness for rafting (m) 99 96 Craft = 5.0 ! coefficient used in the rafting function 100 ridge_por = 0.3 ! initial porosity of the ridged ice (typically 0.30)101 partfun_swi = 1 ! participation function linear, TH75 (0) or exponential Letal07 (1)102 brinstren_swi = 0 ! (1) use brine volume to diminish ice strength97 ridge_por = 0.3 ! porosity of newly ridged ice 98 partfun_swi = 1 ! participation function: linear, Thorndike et al 75 (0) or exponential, Lipscomb 07 (1) 99 brinstren_swi = 0 ! ice strength is a function of brine volume fraction (1) or not (0)
Note: See TracChangeset
for help on using the changeset viewer.