Context Navigation

hemin40_param_list.dat @ 83

Last change on this file since 83 was 83, checked in by aquiquet, 8 years ago
iLOVECLIM branch: updated standard module choix + standard param_list for Hemin-40 geometry
File size: 16.9 KB

Line
1	! Parametres du run : geoplace GI2S_5T Groenland Ice2sea 5 km cut Tamsin
2
3	!___________________________________________________________
4	&runpar ! nom du bloc parametres du run
5	runname = "THEBEAUTIFULRUNNAME" ! 8 caracteres
6	icompteur = THERESTARTCOMPTEUR
7	iout = THERESTARTIOUT
8	reprcptr = "THERESTARTFILEGRISLI"
9	itracebug = 0
10	num_tracebug = 166
11	comment_run = "Test couple"
12	/
13	! runname : nom de l experience (8 caracteres)
14	! icompteur : reprise dans un fichier 0 -> non, 1 -> oui, 2 -> T et Hwat
15	! 3-> T seulement
16	! iout : 1-> sortie cptr pour reprise,2 -> sortie nc pour reprise
17	! reprcptr : nom du fichier
18	!___________________________________________________________
19	&grdline ! bloc grounding line
20
21	igrdline = 0
22	/
23	! igrdline : 1 ligne d echouage fixée, sinon 0
24	! paleo grounding line -> 2
25	!___________________________________________________________
26	&timesteps ! bloc timestep
27
28	tend = 100.
29	tbegin = 0. !1.e10 ! si tbegin > 1.e9 on prend le temps du fichier cptr
30	dtmin = 2.e-3
31	dtmax = 1.
32	dtt = 1.
33	testdiag = 0.025 !0.016
34	/
35	! tous les temps en annees. tbegin et tend : debut et fin du run
36	! pour equation masse, pas de temps mini -> dtmin, maxi -> dtmax
37	! dtt : pas de temps long
38	! testdiag, pour gerer le pas de temps dynamique dt
39	! ordres de grandeur (a moduler selon dx) :
40	! 40 km dtmin=2.e-3, dtmax=1., dtt=5., tesdiag=0.02
41
42	!___________________________________________________________
43	&topo_file
44	topo_ref = "Greenland_Bamber-et-al_composite_hemin40.nc" !"etopo1_grisliN_GMT.g40"
45	topo_dep = "Greenland_Bamber-et-al_composite_hemin40.nc" !"etopo1_grisliN_GMT.g40"
46	grid_topo = "coord-nord-40km.dat"
47	ghf_fich = "ijphi_hemin40.nc"
48	/
49	! file1 : topo de depart
50	! file2 : topo de reference
51	!___________________________________________________________
52	&mass_conserv ! conservation de la masse
53	adv_frac = 2. ! 2-> advection seule
54	V_limit = 3000. ! maximum velocity for mass conservation : depend de la calotte
55	/
56	! Conservation de la masse avec equation advection-diffusion '
57	! la repartition depend de adv_frac'
58	! >1 -> advection seule'
59	! 0 -> diffusion seule'
60	! 0<*<1 -> fraction de l advection'
61	! -1 -> zones diffusion + zones advecttion'
62	! V_limit depend de la calotte : typiquement 3000 en Antarctique, 10000 au Groenland
63
64	!--------------------------------------------------------------------
65	&eaubasale1 ! nom du premier bloc eau basale
66
67	ecoulement_eau = T
68	hwatermax = 5000.000
69	infiltr = 1.e-3 !1.0000001E-03
70	/
71	! ecoulement eau : .false. -> modele bucket, sinon equ. diffusion
72	! hwatermax : hauteur d eau basale maximum dans le sediment (m)
73	! infiltr est la quantite d eau qui peut s infiltrer dans le sol (m/an)
74
75	!___________________________________________________________
76	&param_hydr ! nom du bloc parametres hydrauliques
77
78	hmax_till = 20.00000
79	poro_till = 0.5000000
80	kond0 = 1.000000E-06
81	/
82	! hmax_till (m) : epaisseur max du sediment
83	! poro_till : porosite du sediment
84	! conductivite du sediment : kond0 (m/s)
85
86	!____________________________________________________________
87	&drag_hwat_cont ! nom du bloc dragging hwater contigu
88
89	hwatstream = 50.
90	cf = 1.e-5 ! 1.e-4
91	betamax = 1000.
92	toblim = 0.7e5 ! 0.25e5
93	/
94	! hwatstream (m) : critere de passage en stream en partant de la cote'
95	! si hwater > hwatstream '
96	! cf coefficient de la loi de frottement fonction Neff'
97	! seulement pour les points cotiers
98	! betamax : (Pa) frottement maxi sous les streams '
99	! toblim : (Pa) pour les iles '
100
101	!____________________________________________________________
102	&drag_hwat_sedim ! nom du bloc dragging hwater sediment
103
104	file_sedim = 'sediment_ij_hemin40.dat'
105	hwatstream = 50.
106	cf = 2.e-5 ! 1.e-4
107	betamax = 1000.
108	toblim = 0.7e5 ! 0.25e5
109	seuil_hwater = 250.
110	seuil_sedim = 150.
111	seuil_neff = 5.e5 !350.e5
112	coef_sedim = 1.
113	coef_gz = 10.
114	coef_ile = 0.1
115	/
116	! hwatstream (m) : critere de passage en stream en partant de la cote'
117	! si hwater > hwatstream '
118	! cf coefficient de la loi de frottement fonction Neff'
119	! tobmax : (Pa) frottement maxi '
120	! toblim : (Pa) tes pour les iles '
121	! seuil_hwater (m) : seuil hwater pour avoir glissement sur zone sediment'
122	! seuil_sedim (m) : seuil epaisseur sediment pour avoir glissement sur zone sediment'
123	! coef_sedim : coef frottement zones stream sediments
124	! coef_gz : coef frottement zones stream std
125	! coef_ile : coef frottement zones iles
126	!____________________________________________________________
127	&drag_neff_slope ! nom du bloc dragging neff slope
128
129	cf = 2.e-5 ! 1.e-4
130	betamax = 1000.
131	betamin = 10.
132	toblim = 0.7e5 ! 0.25e5
133	tostick = 1.e5
134	seuil_neff = 5.e5 !350.e5
135	coef_gz = 10.
136	coef_ile = 0.1
137	slope_fich = "ETOPO1_slope_R2min_hemin40.nc"
138	expo_slope = 1.
139	pente_min = 500.
140	pente_max = 2000.
141	/
142	! cf coefficient de la loi de frottement fonction Neff'
143	! betamax : (Pa) frottement maxi '
144	! betamin : (Pa) frottement mini '
145	! toblim : (Pa) tes pour les iles '
146	! tostick : (Pa) pour les points non flgzmx '
147	! seuil_neff (Pa) : seuil sur la pression effective pour avoir glissement'
148	! coef_gz : coef frottement zones stream std
149	! coef_ile : coef frottement zones iles
150
151	!___________________________________________________________
152	&drag_jorge ! nom du bloc dragging_jorge
153
154	hwatstream = 400.
155	cf = 2.e-5 ! 1.e-4
156	tobmax = 1000.
157	toblim = 0.7e5 ! 0.25e5
158	seuil_sedim = 1. ! 0.0001
159	/
160	! hwatstream (m) : critere de passage en stream en partant de la cote'
161	! si hwater > hwatstream '
162	! cf coefficient de la loi de frottement fonction Neff'
163	! tobmax : (Pa) frottement maxi '
164	! toblim : (Pa) tes pour les iles '
165	!___________________________________________________________
166
167	&calving ! nom du bloc calving méthode Vincent
168
169	Hcoup = 200 ! tres petit quand shelves fixes sinon 250
170	ifrange = 4
171	meth_Hcoup = 0
172	/
173
174	! Hcoup epaisseur de coupure
175	! ifrange=0 -> pas de traitement particulier sur les bords'
176	! ifrange=1 -> traitement de Vincent avec ice shelves frangeants'
177	! ifrange=2 -> ice shelves frangeant seulement si bm-bmelt positif
178	!____________________________________________________________
179	! loi de deformation 1 module deformation_mod_2lois
180	&loidef_1
181
182	exposant_1 = 3.
183	temp_trans_1 = -6.5
184	enhanc_fact_1 = 3. ! Grice2sea 3.
185	coef_cold_1 = 1.660E-16
186	Q_cold_1 = 7.820E+04
187	coef_warm_1 = 2.000E-16
188	Q_warm_1 = 9.545E+04
189	/
190	! exposant (glen), temperature de transition (ttrans)
191	! enhancement factor (sf)
192	! pour les temperatures inf. a Temp_trans :
193	! coef_cold (Bat1) et Q_cold (Q1)
194	! pour les temperatures sup. a Temp_trans :
195	! coef_warm (Bat2) et Q_warm (Q2)
196	!________________________________________________________
197	! loi de deformation 2 module deformation_mod_2lois
198	&loidef_2
199
200	exposant_2 = 1.
201	temp_trans_2 = -10.
202	enhanc_fact_2 = 1. ! Grice2sea 3.
203	coef_cold_2 = 8.313E-08
204	Q_cold_2 = 4.000E+04
205	coef_warm_2 = 8.313E-08
206	Q_warm_2 = 6.000E+04
207	/
208	! exposant (glen), temperature de transition (ttrans)
209	! enhancement factor (sf)
210	! pour les temperatures inf. a Temp_trans :
211	! coef_cold (Bat1) et Q_cold (Q1)
212	! pour les temperatures sup. a Temp_trans :
213	! coef_warm (Bat2) et Q_warm (Q2)
214	!___________________________________________________________
215
216	&diagno_rheol ! nom du bloc diagno_rheol
217
218	sf01 = 0.125 ! 0.125
219	sf03 = 0.125 ! 0.125
220	pvimin = 1.5e3 ! 1.e3
221	/
222	! coefficients par rapport a la loi glace posee
223	! sf01 : coefficient viscosite loi lineaire
224	! sf03 : coefficient viscosite loi n=3
225	! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3
226	! tres petit par rapport aux valeurs standards ~ 1.e10
227	!___________________________________________________________
228
229	! glissement module sliding_Bindschadler
230	&slid_bindsh ! nom du bloc
231
232	kweert = 0. ! 5.e-11
233	loigliss = 2
234	coefbmax = 10.
235	/
236
237	! kweert : coefficent, loigliss le type de loi
238	! coefbmax : facteur de normalisation pour influence eau
239
240	!___________________________________________________________
241	&snap_forcage_mois ! module climat_forcage_mois_mod'
242
243	filtr_t1 = 'PnGrmin360_SE_2431_2440_1M_histmthij-t2m.dat'
244	filtr_p1 = 'PnGrmin360_SE_2431_2440_1M_histmthij-precip.dat'
245	/
246	! fichiers temperature et precip : 12 mois
247
248	!___________________________________________________________
249	&lapse_month ! module climat_forcage_mois_mod'
250
251	lect_lapse = 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0 ! 12 mois
252	/
253	! laspe rates janvier -> decembre en deg/km
254	!___________________________________________________________
255	&clim_coupl_atm ! nom du bloc
256
257	mincoefbmelt = 0.
258	maxcoefbmelt = 2.
259	filforc = "signal-cycle-hemin.dat"
260	snow = .true.
261	! mincoefbmelt : borne mini de coefbmshelf
262	! maxcoefbmelt : borne maxi de coefbmshelf
263	! filforc : fichier de forcage nivx des mers
264	! snow : .true. si neige, .false. si precip
265	/
266	!___________________________________________________________
267	&ablation ! module ablation_mod'
268
269	pdd_type = 1
270	annual = .false.
271	Cice = 0.008
272	Csnow = 0.003
273	csi = 0.6
274	sigma_ice = 5.0
275	/
276	! pdd_type : 0 reeh, 1 Fausto, 2 Tarasov
277	! annual : T = annuel, F = mensuel
278	! Cice and Csnow, melting factors for ice and snow
279	! sigma variabilite Tday
280	! csi proportion of melted water that can refreeze
281
282	!___________________________________________________________
283	&clim_pert ! nom du bloc
284
285	coefT = 1.
286	rappact = 0.07 ! 0.05->0.75 0.07 -> 0.5
287	retroac = 1
288	rapbmshelf = 5.
289	mincoefbmelt = 0.
290	maxcoefbmelt = 2.
291	filforc = 'forcage-900k-2007-Parrenin-Bassinot.dat' !'LGM_permanent.dat'
292	/
293	!___________________________________________________________
294	&clim_pert_massb ! nom du bloc
295
296	coefT = 1.
297	filforc = 'quiquet-et-al_CP2013_forcage.dat' !'forcage-900k-2007-Parrenin-Bassinot.dat' !'LGM_permanent.dat'
298	pertsmb = 1 !1
299	rapsmb = 0.07 !7
300	/
301	! pertsmb = 1, smb pert via tpert / pertsmb = 0, smb = cst
302	!----------------------------------------------------------
303	&meca_SIA_L1 ! bloc resol_meca
304
305	i_resolmeca = 2
306	/
307	! i_resolmeca type d association entre SIA et L1'
308	! i_resolmeca=0 chacun dans sa zone'
309	! i_resolmeca=1 dans les zones stream, addition si uxdef > uxL1 (MIS11 Cairns)
310	! i_resolmeca=2 addition systematique dans les zones stream
311
312	!----------------------------------------------------------
313	! PDD base Tann et Tjuly module ablation_ann
314
315	&ablation_ann
316	Cice = 0.008
317	Csnow = 0.005
318	Sigma = 5.
319	csi = 0.6
320	/
321	! Cice and Csnow, melting factors for ice and snow
322	! sigma variabilite Tday
323	! csi proportion of melted water that can refreeze
324	!----------------------------------------------------------
325	! pour le module tracer (suivi d'isotopes)
326	&tracer
327	file_tr_dat = "../INPUT/tracer.dat"
328	file_tr_out = "../RESULTATS/tracer.out"
329	file_tr_dep = "../RESULTATS/deparchive.out"
330	lambdab = 0.0001
331	agemax = 1000000
332	coefT_tra = 1.
333	rappact_tra = 0.07
334	filforc = "quiquet-et-al_CP2013_forcage.dat" !"forcage-900k-2007-Parrenin-Bassinot.dat"
335	pert_type = 0
336	/
337	!----------------------------------------------------------
338	! module lect_topo_green_gen
339	&topo_groen_gen
340	topo_surf = 'Surface_mean_5km_Tamsin_ieq-param.grd' ! surface
341	correc_surf = 'no' ! ice-real correction
342	topo_bed = 'Bedrock_mean_5km_Tamsin.grd' ! bedrock
343	mask_grounded = 'Mask_5km_checked_Tamsin-Fjords.grd' ! mask
344	longitude = 'long_ZBL_CISM_5km.grd' ! longitude ! a garder en .dat
345	latitude = 'lat_ZBL_CISM_5km.grd' ! latitude ! car donne xmin ymin
346	heatflux = 'ghf-cism.grd' ! geothermal heat flux
347	/
348	! remarques : l'epaisseur n'est pas lue mais calculee a partir du masque
349	!
350	!----------------------------------------------------------
351	! module lect_climref_Ice2sea
352	&clim_smb_T_gen
353	smb_file = 'smb_MARv3.5.2-25km_1958_1989mean_hemin40.nc' !'smb_ERAI_MAR_highres_1989-2008mean_hemin40.nc' ! smb (en kg/m2/an
354	coef_smb_unit = 1.0893e-3 ! 1.0989e-3 ! 1.0989e-3 pour1/910 1.0893e-3 pour 1/918 (valeur habituelle)
355	temp_annual_file ='Tann-cism-hemin40.nc' ! annual surface temperature
356	/
357	!----------------------------------------------------------
358	! module climat_Grice2sea_mod
359	&lapse_rates
360	T_lapse_rate = 0. !-6.309e-3 ! celui pour la T annuelle de Fausto
361	/
362	!----------------------------------------------------------
363	&clim_snap
364	nb_snap = 200
365	time_depart_snaps = 2000000
366	ecart_snap = 1.
367	file_smb_snap = 'smb_ECHAM5_A1B_MAR_highres_all_time.nc'
368	massb_time = 0 ! 0=fixe, 1:interpolation snapshots
369	/
370	!----------------------------------------------------------
371	&grad_smb !routine init_grad_smb
372	file_grad_smb = 'GRISLI_gradients_MAR_00.txt'
373	/
374	!---------------------------------------------------------
375	&bmelt_seuil ! module bmelt_seuil_prof
376	bm_grz = 3. !3
377	bmshelf_plateau = 1. !0.4
378	bmshelf_abysses = 10.
379	depth_talus = -1000.
380	/
381	! Pour l actuel : bm_grz a la grounding line
382	! bmshelf_plateau sur le plateau continental
383	! bmshelf_abysses pour les grandes profondeurs
384	! depth_talus, negative, separation entre les 2 domaines
385	!___________________________________________________________
386	&beta_prescr ! dragging_prescr_beta'
387	beta_c_file = 'beta-GISCAM23_20ans.grd' ! 'beta_Fab_10avril12_bord.grd' ! 'beta_fab_comb_alpha_v4.grd'
388	beta_limgz = .5e6 ! -1.e12 pour plastic ? ! .5e6 en visqueux
389	beta_min = 1 !10. ! for grounded ice (Pa)
390	beta_mult = 1. ! coefficient multiplicateur
391	/
392	! read beta on centered grid
393	! beta_file : nom des fichiers qui contiennent les betamx et betamy
394	! beta-estime-run07.dat calcul direct beta sans modif
395	! beta-staggered-15km-opt1.dat premier test en utilisant beta calcule avec iterations
396	! above beta_limgz, gzmx is false if negative : no sliding
397	! betamx_file = not read anymore
398	! betamy_file = not read anymore
399	!___________________________________________________________
400	&beta_iter_vitbil ! beta_iter_vitbil
401
402	time_iter = 5 ! temps de demarrage des iterations
403	nb_iter_vitbil = 2 ! nombre d'iterations par pas de temps (dtt)
404	coef_iter_vitbil = 1. ! coefficient pour rapport des vitesses <=1
405	Umag_bil_file = 'ebauche_vitbil_GrIce2sea.grd' ! fichier des vitesses de bilan centrees
406	/
407	!
408	!___________________________________________________________
409	&spinup ! warning : 2 different modules
410	ispinup = 0
411	/
412	! with module no_spinup
413	! ispinup = 0 run standard ou calcul du beta
414	! ispinup = 1 temperature equilibrium with grisli velocities'
415	!
416	! with module spinup_vitbil
417	! ispinup = 2 conservation de la masse avec vitesses bilan '
418	! ispinup = 3 equilibre temperature avec vitesses bilan'
419	!
420	!___________________________________________________________
421	&vitbil_upwind ! nom du bloc vitbil calcule sur le mailles staggered
422	balance_Ux_file = 'Uy_moy_stag_5km.nc'
423	balance_Uy_file = 'Ux_moy_stag_5km.nc'
424	/
425	!balance_Ux_file = 'Ux_moy_stag_5km.nc'
426	!balance_Uy_file = 'Uy_moy_stag_5km.nc'
427	!balance_vel_file = 'Umoy-for-calc-beta_LBq15-06_ZBL_15km.dat'
428	! balance_vel_file = 'Uslid-for-calc-beta_LBq15-06_ZBL_15km.dat'
429	! balance_vel_file = 'Umoy-for-calc-beta_LBq15-06_ZBL_15km.dat'
430	! balance_vel_file = 'Vitbil-stag-pour-spinup_merge_runs_Lbq_ZBL_15km.dat'
431	! balance velocities on staggered grid Ux, Uy
432	!___________________________________________________________
433	&output_regions
434	nbregions = 7
435	region_file = 'Basins_Greenland_aurel.grd'
436	/
437	! nbregions = nombre de regions
438	! region_file = carte des regions
439	!_______________________________________________________________________________

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

Download in other formats:

Original Format