New URL for NEMO forge! http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.

namelist_ice_lim3_ref in branches/2017/dev_r8183_ICEMODEL/NEMOGCM/CONFIG/SHARED – NEMO

Context Navigation

source: branches/2017/dev_r8183_ICEMODEL/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref @ 8514

Last change on this file since 8514 was 8514, checked in by clem, 7 years ago
changes in style - part5 - almost done
File size: 17.1 KB

Rev	Line
[3875]	1	!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[7646]	2	!! LIM3 namelist:
[8512]	3	!! 1 - Generic parameters (namice_run)
	4	!! 2 - Ice thickness discretization (namice_itd)
	5	!! 3 - Ice ridging/rafting (namice_rdgrft)
	6	!! 4 - Ice rheology (namice_rhg)
	7	!! 5 - Ice advection (namice_adv)
	8	!! 6 - Ice thermodynamics (namice_thd)
	9	!! 7 - Ice salinity (namice_sal)
	10	!! 8 - Ice melt ponds (namice_mp)
	11	!! 9 - Ice initialization (namice_ini)
	12	!! 10 - Ice/snow albedos (namice_alb)
	13	!! 11 - Ice diagnostics (namice_dia)
[7646]	14	!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5124]	15	!
	16	!------------------------------------------------------------------------------
[8512]	17	&namice_run ! Generic parameters
[5124]	18	!------------------------------------------------------------------------------
[8512]	19	jpl = 5 ! number of ice categories
	20	nlay_i = 2 ! number of ice layers
	21	nlay_s = 1 ! number of snow layers (only 1 is working)
	22	nn_monocat = 0 ! virtual ITD mono-category parameterizations (1-4 => jpl = 1 only) or not (0)
	23	! 2: simple piling instead of ridging --- temporary option
	24	! 3: activate G(he) only --- temporary option
	25	! 4: activate extra lateral melting only --- temporary option
	26	rn_amax_n = 0.997 ! maximum tolerated ice concentration NH
	27	rn_amax_s = 0.997 ! maximum tolerated ice concentration SH
	28	cn_icerst_in = "restart_ice" ! suffix of ice restart name (input)
	29	cn_icerst_out = "restart_ice" ! suffix of ice restart name (output)
	30	cn_icerst_indir = "." ! directory to read input ice restarts
	31	cn_icerst_outdir = "." ! directory to write output ice restarts
[3875]	32	/
[5124]	33	!------------------------------------------------------------------------------
[8512]	34	&namice_itd ! Ice discretization
[7646]	35	!------------------------------------------------------------------------------
[8512]	36	rn_himean = 2.0 ! expected domain-average ice thickness (m)
[7646]	37	/
	38	!------------------------------------------------------------------------------
[8512]	39	&namice_rdgrft ! Ice ridging/rafting
[5124]	40	!------------------------------------------------------------------------------
[8512]	41	! -- ice_rdgrft_strength -- !
[8514]	42	ln_str_H79 = .true. ! ice strength param.: Hibler_79 => P = pstar<h>exp(-c_rhg*A)
[8512]	43	rn_pstar = 2.0e+04 ! ice strength thickness parameter (N/m2)
	44	rn_crhg = 20.0 ! ice strength conc. parameter (-)
[8514]	45	ln_str_R75 = .false. ! ice strength param.: Rothrock_75 => P = Cfcoeffintegral(wr.h^2)
[8512]	46	rn_perdg = 17.0 ! ridging work divided by pot. energy change in ridging
	47	! -- ice_rdgrft -- !
	48	rn_cs = 0.5 ! fraction of shearing energy contributing to ridging
	49	! -- ice_rdgrft_prep -- !
	50	ln_partf_lin = .false. ! Linear ridging participation function (Thorndike et al, 1975)
	51	rn_gstar = 0.15 ! fractional area of thin ice being ridged
	52	ln_partf_exp = .true. ! Exponential ridging participation function (Lipscomb, 2007)
	53	rn_astar = 0.03 ! exponential measure of ridging ice fraction [set to 0.05 if hstar=100]
	54	ln_ridging = .true. ! ridging activated (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO
	55	rn_hstar = 25.0 ! determines the maximum thickness of ridged ice (m) (Hibler, 1980)
	56	rn_porordg = 0.3 ! porosity of newly ridged ice (Lepparanta et al., 1995)
	57	rn_fsnwrdg = 0.5 ! snow volume fraction that survives in ridging
	58	rn_fpndrdg = 1.0 ! pond fraction that survives in ridging (small a priori)
	59	ln_rafting = .true. ! rafting activated (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO
	60	rn_hraft = 0.75 ! threshold thickness for rafting (m)
	61	rn_craft = 5.0 ! squeezing coefficient used in the rafting function
	62	rn_fsnwrft = 0.5 ! snow volume fraction that survives in rafting
	63	rn_fpndrft = 1.0 ! pond fraction that survives in rafting (0.5 a priori)
[3875]	64	/
[5124]	65	!------------------------------------------------------------------------------
[8512]	66	&namice_rhg ! Ice rheology
[5124]	67	!------------------------------------------------------------------------------
[8512]	68	rn_ishlat = 2. ! free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2)
	69	rn_cio = 5.0e-03 ! ice-ocean drag coefficient (-)
	70	rn_creepl = 1.0e-12 ! creep limit (s-1)
	71	rn_ecc = 2.0 ! eccentricity of the elliptical yield curve
	72	nn_nevp = 120 ! number of EVP subcycles
	73	rn_relast = 0.333 ! ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast
	74	! advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300)
	75	ln_landfast = .false. ! landfast ice parameterization (T or F)
	76	rn_gamma = 0.15 ! fraction of ocean depth that ice must reach to initiate landfast
	77	! recommended range: [0.1 ; 0.25]
	78	rn_icebfr = 10. ! maximum bottom stress per unit area of contact (N/m2)
	79	! a very large value ensures ice velocity=0 even with a small contact area
	80	! recommended range: ?? (should be greater than atm-ice stress => >0.1 N/m2)
	81	rn_lfrelax = 1.e-5 ! relaxation time scale to reach static friction (s-1)
[3875]	82	/
[5124]	83	!------------------------------------------------------------------------------
[8512]	84	&namice_adv ! Ice advection
[5124]	85	!------------------------------------------------------------------------------
[8512]	86	ln_icedyn = .true. ! ice dynamics (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO
	87	nn_icedyn = 2 ! switch for ice dynamics
	88	! 2: total
	89	! 1: advection only (no diffusion, no ridging/rafting)
	90	! 0: advection only (as 1 but with prescribed velocity, bypass rheology)
	91	rn_uice = 0.00001 ! prescribed ice u-velocity
	92	rn_vice = -0.00001 ! prescribed ice v-velocity
	93	! -- ice_adv -- !
	94	ln_adv_Pra = .false. ! Advection scheme (Prather)
	95	ln_adv_UMx = .true. ! Advection scheme (Ultimate-Macho)
[8514]	96	nn_UMx = 5 ! order of the scheme for UMx (1-5 ; 20=centered 2nd order)
[3875]	97	/
[5124]	98	!------------------------------------------------------------------------------
[8512]	99	&namice_thd ! Ice thermodynamics
[5124]	100	!------------------------------------------------------------------------------
[8514]	101	ln_icethd = .true. ! ice thermo (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO
	102	! -- limthd_dif -- !
	103	rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice (m-1)
	104	ln_cndi_U64 = .false. ! sea ice thermal conductivity: k = k0 + beta.S/T (Untersteiner, 1964)
	105	ln_cndi_P07 = .true. ! sea ice thermal conductivity: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007)
	106	ln_dqns_i = .true. ! change the surface non-solar flux with surface temperature (T) or not (F)
	107	rn_cnd_s = 0.31 ! thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971)
	108	! Obs: 0.1-0.5 (Lecomte et al, JAMES 2013)
	109	! -- limthd_dh -- !
	110	ln_icedH = .true. ! activate ice thickness change from growing/melting (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO
	111	rn_blow_s = 0.66 ! mesure of snow blowing into the leads
	112	! = 1 => no snow blowing, < 1 => some snow blowing
	113	! -- limthd_da -- !
	114	ln_icedA = .true. ! activate lateral melting param. (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO
	115	rn_beta = 1.0 ! coef. beta for lateral melting param. Recommended range=[0.8-1.2]
	116	! => decrease = more melt and melt peaks toward higher concentration (A~0.5 for beta=1 ; A~0.8 for beta=0.2)
	117	! 0.3 = best fit for western Fram Strait and Antarctica
	118	! 1.4 = best fit for eastern Fram Strait
	119	rn_dmin = 8. ! minimum floe diameter for lateral melting param. Recommended range=[6-10]
	120	! => 6 vs 8m = +40% melting at the peak (A~0.5)
	121	! 10 vs 8m = -20% melting
	122	! -- limthd_lac -- !
	123	ln_icedO = .true. ! activate ice growth in open-water (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO
	124	rn_hinew = 0.1 ! thickness for new ice formation in open water (m)
	125	ln_frazil = .false. ! Frazil ice parameterization (ice collection as a function of wind)
	126	rn_maxfraz = 1.0 ! maximum fraction of frazil ice collecting at the ice base
	127	rn_vfraz = 0.417 ! thresold drift speed for frazil ice collecting at the ice bottom (m/s)
	128	rn_Cfraz = 5.0 ! squeezing coefficient for frazil ice collecting at the ice bottom
[8512]	129	! -- limitd_th -- !
[8514]	130	rn_himin = 0.1 ! minimum ice thickness (m) used in remapping, must be smaller than rn_hnewice
	131	! -- icestp -- !
	132	nn_iceflx = -1 ! Redistribute heat flux over ice categories
	133	! =-1 Do nothing (needs N(cat) fluxes)
	134	! ==> forced mode only
	135	! = 0 Average N(cat) fluxes then apply the average over the N(cat) ice
	136	! ==> forced and coupled modes
	137	! = 1 Average N(cat) fluxes then redistribute over the N(cat) ice using T-ice and albedo sensitivity
	138	! ==> forced mode only
	139	! = 2 Redistribute a single flux over categories
	140	! ==> coupled mode only
[3875]	141	/
[5124]	142	!------------------------------------------------------------------------------
[8512]	143	&namice_sal ! Ice salinity
[5124]	144	!------------------------------------------------------------------------------
[8514]	145	! -- limthd_sal -- !
	146	ln_icedS = .true. ! activate gravity drainage and flushing (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO
	147	nn_icesal = 2 ! ice salinity option
	148	! 1: constant ice salinity (S=rn_icesal)
	149	! 2: varying salinity parameterization S(z,t)
	150	! 3: prescribed salinity profile S(z), Schwarzacher, 1959
	151	rn_icesal = 4. ! (nn_icesal=1) ice salinity (g/kg)
	152	rn_sal_gd = 5. ! restoring ice salinity, gravity drainage (g/kg)
	153	rn_time_gd = 1.73e+6 ! restoring time scale, gravity drainage (s)
	154	rn_sal_fl = 2. ! restoring ice salinity, flushing (g/kg)
	155	rn_time_fl = 8.64e+5 ! restoring time scale, flushing (s)
	156	rn_simax = 20. ! maximum tolerated ice salinity (g/kg)
	157	rn_simin = 0.1 ! minimum tolerated ice salinity (g/kg)
[3875]	158	/
[5124]	159	!------------------------------------------------------------------------------
[8233]	160	&namicemp ! Melt ponds
	161	!------------------------------------------------------------------------------
[8514]	162	! -- limmp -- !
	163	ln_pnd = .false. ! active melt ponds
	164	ln_pnd_rad = .false. ! active melt ponds radiative coupling
	165	ln_pnd_fw = .false. ! active melt ponds freshwater coupling
	166	nn_pnd_scheme = 0 ! type of melt pond scheme : =0 prescribed ( Tsu=0 ), =1 empirical, =2 topographic
	167	rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 : (0<rn_apnd<1, nn_pnd_scheme = 0)
	168	rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 : (0<rn_apnd<1, nn_pnd_scheme = 0)
[8233]	169	/
	170	!------------------------------------------------------------------------------
[8512]	171	&namice_ini ! Ice initialization
[5124]	172	!------------------------------------------------------------------------------
[8514]	173	! -- limistate -- !
	174	ln_iceini = .true. ! activate ice initialization (T) or not (F)
	175	ln_iceini_file = .false. ! netcdf file provided for initialization (T) or not (F)
	176	rn_thres_sst = 2.0 ! maximum water temperature with initial ice (degC)
	177	rn_hts_ini_n = 0.3 ! initial real snow thickness (m), North
	178	rn_hts_ini_s = 0.3 ! " " South
	179	rn_hti_ini_n = 3.0 ! initial real ice thickness (m), North
	180	rn_hti_ini_s = 1.0 ! " " South
	181	rn_ati_ini_n = 0.9 ! initial ice concentration (-), North
	182	rn_ati_ini_s = 0.9 ! " " South
	183	rn_smi_ini_n = 6.3 ! initial ice salinity (g/kg), North
	184	rn_smi_ini_s = 6.3 ! " " South
	185	rn_tmi_ini_n = 270. ! initial ice/snw temperature (K), North
	186	rn_tmi_ini_s = 270. ! " " South
[8512]	187
	188	sn_hti = 'Ice_initialization' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', ''
	189	sn_hts = 'Ice_initialization' , -12 ,'hts' , .false. , .true., 'yearly' , '' , '', ''
	190	sn_ati = 'Ice_initialization' , -12 ,'ati' , .false. , .true., 'yearly' , '' , '', ''
	191	sn_tsu = 'Ice_initialization' , -12 ,'tsu' , .false. , .true., 'yearly' , '' , '', ''
	192	sn_tmi = 'Ice_initialization' , -12 ,'tmi' , .false. , .true., 'yearly' , '' , '', ''
	193	sn_smi = 'Ice_initialization' , -12 ,'smi' , .false. , .true., 'yearly' , '' , '', ''
	194	cn_dir='./'
[3875]	195	/
[8512]	196	!------------------------------------------------------------------------------
	197	&namice_alb ! albedo parameters
	198	!------------------------------------------------------------------------------
[8514]	199	nn_ice_alb = 1 ! parameterization of ice/snow albedo
	200	! 0: Shine & Henderson-Sellers (JGR 1985), giving clear-sky albedo
	201	! 1: "home made" based on Brandt et al. (JClim 2005)
	202	! and Grenfell & Perovich (JGR 2004), giving cloud-sky albedo
	203	! 2: as 1 with melt ponds
	204	rn_alb_sdry = 0.85 ! dry snow albedo : 0.80 (nn_ice_alb = 0); 0.85 (nn_ice_alb = 1); obs 0.85-0.87 (cloud-sky)
	205	rn_alb_smlt = 0.75 ! melting snow albedo : 0.65 ( '' ) ; 0.75 ( '' ) ; obs 0.72-0.82 ( '' )
	206	rn_alb_idry = 0.60 ! dry ice albedo : 0.72 ( '' ) ; 0.60 ( '' ) ; obs 0.54-0.65 ( '' )
	207	rn_alb_imlt = 0.50 ! bare puddled ice albedo : 0.53 ( '' ) ; 0.50 ( '' ) ; obs 0.49-0.58 ( '' )
	208	rn_alb_dpnd = 0.27 ! ponded ice albedo : 0.25 ( '' ) ; 0.27 ( '' ) ; obs 0.10-0.30 ( '' )
[8321]	209	/
[8512]	210	!------------------------------------------------------------------------------
	211	&namice_dia ! Diagnostics
	212	!------------------------------------------------------------------------------
[8514]	213	ln_icediachk = .false. ! check online the heat, mass & salt budgets (T) or not (F)
	214	ln_icediahsb = .false. ! output the heat, mass & salt budgets (T) or not (F)
	215	ln_icectl = .false. ! ice points output for debug (T or F)
	216	iiceprt = 10 ! i-index for debug
	217	jiceprt = 10 ! j-index for debug
[8512]	218	/

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format