Context Navigation

hemin40_param_list.dat @ 146

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

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

Download in other formats:

Original Format