- Timestamp:
- 2017-09-12T20:46:13+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref
r8516 r8517 18 18 &namice_run ! Generic parameters 19 19 !------------------------------------------------------------------------------ 20 jpl = 5! number of ice categories21 nlay_i = 2! number of ice layers22 nlay_s = 1! number of snow layers (only 1 is working)23 nn_monocat = 0! virtual ITD mono-category parameterizations (1-4 => jpl = 1 only) or not (0)20 jpl = 5 ! number of ice categories 21 nlay_i = 2 ! number of ice layers 22 nlay_s = 1 ! number of snow layers (only 1 is working) 23 nn_monocat = 0 ! virtual ITD mono-category parameterizations (1-4 => jpl = 1 only) or not (0) 24 24 ! 2: simple piling instead of ridging --- temporary option 25 25 ! 3: activate G(he) only --- temporary option 26 26 ! 4: activate extra lateral melting only --- temporary option 27 ln_icedyn = .true. ! ice dynamics (T) or not (F) 28 ln_icethd = .true. ! ice thermo (T) or not (F) 27 29 rn_amax_n = 0.997 ! maximum tolerated ice concentration NH 28 30 rn_amax_s = 0.997 ! maximum tolerated ice concentration SH … … 41 43 &namice_dyn ! Ice dynamics 42 44 !------------------------------------------------------------------------------ 43 ln_icedyn = .true. ! ice dynamics (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 44 nn_icedyn = 2 ! switch for ice dynamics 45 ! 2: total 46 ! 1: advection only (no diffusion, no ridging/rafting) 47 ! 0: advection only (as 1 but with prescribed velocity, bypass rheology) 48 rn_uice = 0.00001 ! prescribed ice u-velocity 49 rn_vice = -0.00001 ! prescribed ice v-velocity 45 ln_dynFULL = .true. ! dyn.: full ice dynamics (rheology + advection + ridging/rafting + correction) 46 ln_dynRHGADV = .false. ! dyn.: no ridge/raft & no corrections (rheology + advection) 47 ln_dynADV = .false. ! dyn.: only advection w prescribed vel.(rn_uvice + advection) 48 rn_uice = 0.00001 ! prescribed ice u-velocity 49 rn_vice = 0. ! prescribed ice v-velocity 50 50 rn_ishlat = 2. ! free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2) 51 51 rn_cio = 5.0e-03 ! ice-ocean drag coefficient (-) … … 53 53 rn_gamma = 0.15 ! fraction of ocean depth that ice must reach to initiate landfast 54 54 ! recommended range: [0.1 ; 0.25] 55 rn_icebfr = 10. ! maximum bottom stress per unit area of contact (N/m2)55 rn_icebfr = 10. ! maximum bottom stress per unit area of contact [N/m2] 56 56 ! a very large value ensures ice velocity=0 even with a small contact area 57 57 ! recommended range: ?? (should be greater than atm-ice stress => >0.1 N/m2) 58 rn_lfrelax = 1.e-5 ! relaxation time scale to reach static friction (s-1)58 rn_lfrelax = 1.e-5 ! relaxation time scale to reach static friction [s-1] 59 59 / 60 60 !------------------------------------------------------------------------------ … … 63 63 ! -- ice_rdgrft_strength -- ! 64 64 ln_str_H79 = .true. ! ice strength param.: Hibler_79 => P = pstar*<h>*exp(-c_rhg*A) 65 rn_pstar = 2.0e+04 ! ice strength thickness parameter (N/m2)65 rn_pstar = 2.0e+04 ! ice strength thickness parameter [N/m2] 66 66 rn_crhg = 20.0 ! ice strength conc. parameter (-) 67 67 ln_str_R75 = .false. ! ice strength param.: Rothrock_75 => P = Cf*coeff*integral(wr.h^2) … … 74 74 ln_partf_exp = .true. ! Exponential ridging participation function (Lipscomb, 2007) 75 75 rn_astar = 0.03 ! exponential measure of ridging ice fraction [set to 0.05 if hstar=100] 76 ln_ridging = .true. ! ridging activated (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO77 rn_hstar = 25.0 ! determines the maximum thickness of ridged ice (m)(Hibler, 1980)76 ln_ridging = .true. ! ridging activated (T) or not (F) 77 rn_hstar = 25.0 ! determines the maximum thickness of ridged ice [m] (Hibler, 1980) 78 78 rn_porordg = 0.3 ! porosity of newly ridged ice (Lepparanta et al., 1995) 79 79 rn_fsnwrdg = 0.5 ! snow volume fraction that survives in ridging 80 80 rn_fpndrdg = 1.0 ! pond fraction that survives in ridging (small a priori) 81 ln_rafting = .true. ! rafting activated (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO82 rn_hraft = 0.75 ! threshold thickness for rafting (m)81 ln_rafting = .true. ! rafting activated (T) or not (F) 82 rn_hraft = 0.75 ! threshold thickness for rafting [m] 83 83 rn_craft = 5.0 ! squeezing coefficient used in the rafting function 84 84 rn_fsnwrft = 0.5 ! snow volume fraction that survives in rafting … … 89 89 !------------------------------------------------------------------------------ 90 90 ln_rhg_EVP = .true. ! EVP rheology 91 rn_creepl = 1.0e-12 ! creep limit (s-1)91 rn_creepl = 1.0e-12 ! creep limit [1/s] 92 92 rn_ecc = 2.0 ! eccentricity of the elliptical yield curve 93 93 nn_nevp = 120 ! number of EVP subcycles … … 105 105 &namice_thd ! Ice thermodynamics 106 106 !------------------------------------------------------------------------------ 107 ln_icethd = .true. ! ice thermo (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO108 107 ! -- icethd_dif -- ! 109 rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice (m-1)108 rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice [1/m] 110 109 ln_cndi_U64 = .false. ! sea ice thermal conductivity: k = k0 + beta.S/T (Untersteiner, 1964) 111 110 ln_cndi_P07 = .true. ! sea ice thermal conductivity: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007) … … 114 113 ! Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) 115 114 ! -- icethd_dh -- ! 116 ln_icedH = .true. ! activate ice thickness change from growing/melting (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO115 ln_icedH = .true. ! activate ice thickness change from growing/melting (T) or not (F) 117 116 rn_blow_s = 0.66 ! mesure of snow blowing into the leads 118 117 ! = 1 => no snow blowing, < 1 => some snow blowing 119 118 ! -- icethd_da -- ! 120 ln_icedA = .true. ! activate lateral melting param. (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO119 ln_icedA = .true. ! activate lateral melting param. (T) or not (F) 121 120 rn_beta = 1.0 ! coef. beta for lateral melting param. Recommended range=[0.8-1.2] 122 121 ! => decrease = more melt and melt peaks toward higher concentration (A~0.5 for beta=1 ; A~0.8 for beta=0.2) … … 127 126 ! 10 vs 8m = -20% melting 128 127 ! -- icethd_lac -- ! 129 ln_icedO = .true. ! activate ice growth in open-water (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO128 ln_icedO = .true. ! activate ice growth in open-water (T) or not (F) 130 129 rn_hinew = 0.1 ! thickness for new ice formation in open water (m), must be larger than rn_hnewice 131 130 ln_frazil = .false. ! Frazil ice parameterization (ice collection as a function of wind) … … 148 147 !------------------------------------------------------------------------------ 149 148 ! -- icethd_sal -- ! 150 ln_icedS = .true. ! activate gravity drainage and flushing (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO149 ln_icedS = .true. ! activate gravity drainage and flushing (T) or not (F) 151 150 nn_icesal = 2 ! ice salinity option 152 151 ! 1: constant ice salinity (S=rn_icesal)
Note: See TracChangeset
for help on using the changeset viewer.