Context Navigation

paramci15ref1 @ 4

Last change on this file since 4 was 4, checked in by dumas, 10 years ago
initial import GRISLI trunk
File size: 11.3 KB

Line
1	! Parametres du run :ci15ref1 greeneem15
2	! sources dans SOURCES
3
4	!___________________________________________________________
5	&runpar ! nom du bloc parametres du run
6
7	runname = "ci15ref1"
8	icompteur = 0
9	iout = 1
10	reprcptr = "../Fichier-CPTR/ci15int1.CPTR"
11	num_tracebug = 166
12	comment_run ="neem v43. statio. nv mk str. rapp 0.05. tob 1000. clim CISM"
13	/
14	! runname : nom de l experience (8 caracteres)
15	! icompteur : reprise dans un fichier 0 -> non, 1 -> oui, 2 -> T et Hw
16	! icompteur : 3-> T seulement
17	! iout : 1-> sortie cptr pour reprise,2 -> sortie nc pour reprise
18	! reprcptr : nom du fichier
19	! itracebug: 1-> ecriture de traces au debut des routines
20	! num_tracebug numero d unite ecriture traces
21	! comment_run: commentaire court sur le run
22
23	!___________________________________________________________
24	&grdline ! bloc grounding line
25
26	igrdline = 0
27	/
28	! igrdline : 1 ligne d echouage fixée, sinon 0
29
30	!___________________________________________________________
31	&timesteps ! bloc time steps
32
33	tend = 0.0000000E+00
34	tbegin = -40000.00
35	dtmin = 2.0000001E-03
36	dtmax = 1.000000
37	dtt = 5.000000
38	testdiag = 2.5000000E-02
39	/
40	! tous les temps en annees. tbegin et tend : debut et fin du run
41	! pour equation masse, pas de temps mini -> dtmin, maxi -> dtmax
42	! dtt : pas de temps long
43	! testdiag, pour gérer le pas de temps dynamique dt
44	! ordres de grandeur (a moduler selon dx) :
45	! 40 km dtmin=2.e-2, dtmax=1., dtt=5., tesdiag=0.02
46
47	!___________________________________________________________
48	&toponeen ! module lect_topo_greeneem
49	topo_ref =groen_grisli_15km.dat
50	topo_dep =groen_grisli_15km.dat
51	grid_topo =groen_grisli_grid_xylonlat_15km.dat
52	ghf0 =
53	ghf_cism_15km.dat
54
55	/
56	! topo_ref= topo actuelle
57	! ghf0 flux geothermique en mW/m2
58
59	!___________________________________________________________
60	&clim_perturb_mois ! module climat_perturb_mois_mod
61	filtr_topo =surface_grisli_15km.dat
62	filtr_t =CNRM_3.3_tt_80-99_15km.dat
63	filtr_p =precipes_cism_v1.1_15km.dat
64	/
65	! fichiers temperature et precip : 12 mois
66	! faire un bloc namelist par snapshot
67
68	!___________________________________________________________
69	&lapse_month ! module climat_perturb_mois_mod
70	filtr_lr =climato_actu_grad-temp_15km.dat
71	/
72	! lapse rates janvier -> decembre en deg/m
73
74	!___________________________________________________________
75	&clim_pert ! module climat_perturb_mod
76
77	coefT = 1.000000
78	rappact = 5.0000001E-02
79	rretroac = 1
80	rapbmshelf = 5.000000
81	mincoefbmelt = 0.0000000E+00
82	maxcoefbmelt = 2.000000
83	filforc = forc_traceurs-stationnaire.dat
84	/
85
86	!___________________________________________________________
87	! ablation_month PDD base Tann et Tjuly
88	&ABLATION_MONTH
89	CICE = 8.0000004E-03,
90	CSNOW = 4.9999999E-03,
91	CSI = 0.6000000 ,
92	SIGMA = 5.000000
93	/
94
95
96	! Cice and Csnow, melting factors for ice and snow
97	! sigma variabilite Tday
98	! csi proportion of melted water that can refreeze
99	!___________________________________________________________
100	&tracer ! module tracer_mod
101	tracer .dat file : ../INPUT/tracer.dat
102	tracer .out file : ../RESULTATS/tracer.out
103	tracer dep. file : ../RESULTATS/deparchive.out
104	/
105	declaration des repertoires traceurs
106
107	!___________________________________________________________
108	&eaubasale1 ! nom du premier bloc eau basale
109
110	ecoulement_eau = T
111	hwatermax = 5000.000
112	infiltr = 1.0000000E-03
113	/
114	! ecoulement eau : .false. -> modele bucket, sinon equ. diffusion
115	! hwatermax : hauteur d eau basale maximum dans le sediment (m)
116	! infiltr est la quantite d eau qui peut s infiltrer dans le sol (m/an)
117
118	!___________________________________________________________
119	&param_hydr ! nom du bloc parametres hydrauliques
120
121	hmax_till = 20.00000
122	poro_till = 0.5000000
123	kond0 = 1.0000000E-06
124	/
125	! hmax_till (m) : epaisseur max du sediment
126	! poro_till : porosite du sediment
127	! conductivite du sediment : kond0 (m/s)
128
129	Calcul Des Temperatures
130	Nfracq= 1
131	Chaleur Demi Maille, Iq= 1
132
133	Si Base Froide, La Chaleur De Glissement Est En Exp(Delta T * Ecart_phid)
134	Ecart_phid= 0.5000000
135	---------------------------------------------------------------------------
136	&BMELT_SEUIL
137	BM_GRZ = 3.000000 ,
138	BMSHELF_PLATEAU = 1.000000 ,
139	BMSHELF_ABYSSES = 10.00000 ,
140	DEPTH_TALUS = -2000.000
141	/
142	!___________________________________________________________
143	&bmelt_seuil ! module bmelt_seuil_prof
144	bm_grz = 3.000000
145	bmshelf_plateau = 1.000000
146	bmshelf_abysses = 10.00000
147	depth_talus = -2000.000
148	/
149	! Pour l actuel : bm_grz a la grounding line
150	! bmshelf_plateau sur le plateau continental
151	! bmshelf_abysses pour les grandes profondeurs
152	! depth_talus, negative, separation entre les 2 domaines
153
154	!___________________________________________________________
155	!___________________________________________________________
156	! loi de deformation 1 module deformation_mod_2lois
157
158	&LOIDEF_1
159	EXPOSANT_1 = 3.000000 ,
160	TEMP_TRANS_1 = -6.500000 ,
161	ENHANC_FACT_1 = 3.000000 ,
162	COEF_COLD_1 = 1.6599999E-16,
163	Q_COLD_1 = 78200.00 ,
164	COEF_WARM_1 = 2.0000000E-16,
165	Q_WARM_1 = 95450.00
166	/
167	! exposant (glen), temperature de transition (ttrans)
168	! enhancement factor (sf)
169	! pour les temperatures inf. a Temp_trans :
170	! coef_cold (Bat1) et Q_cold (Q1)
171	! pour les temperatures sup. a Temp_trans :
172	! coef_warm (Bat2) et Q_warm (Q2)
173	!___________________________________________________________
174	!___________________________________________________________
175	! loi de deformation 2 module deformation_mod_2lois
176
177	&LOIDEF_2
178	EXPOSANT_2 = 1.000000 ,
179	TEMP_TRANS_2 = -10.00000 ,
180	ENHANC_FACT_2 = 1.000000 ,
181	COEF_COLD_2 = 8.3129997E-08,
182	Q_COLD_2 = 40000.00 ,
183	COEF_WARM_2 = 8.3129997E-08,
184	Q_WARM_2 = 60000.00
185	/
186	! exposant (glen), temperature de transition (ttrans)
187	! enhancement factor (sf)
188	! pour les temperatures inf. a Temp_trans :
189	! coef_cold (Bat1) et Q_cold (Q1)
190	! pour les temperatures sup. a Temp_trans :
191	! coef_warm (Bat2) et Q_warm (Q2)
192	!___________________________________________________________
193	!___________________________________________________________
194	& meca_SIA_L1 ! bloc resol_meca
195
196	i_resolmeca = 2
197	/
198	! i_resolmeca type d association entre SIA et L1
199	! i_resolmeca=0 chacun dans sa zone
200	! i_resolmeca=1 dans les zones stream, addition si uxdef > uxL1 (MIS11 Cairns)
201	! i_resolmeca=2 addition systematique dans les zones stream
202
203	!___________________________________________________________
204	&diagno_rheol ! nom du bloc diagno_rheol
205
206	sf01 = 0.1250000
207	sf03 = 0.1250000
208	pvimin = 1000.000
209	/
210	! coefficients par rapport a la loi glace posee
211	! sf01 : coefficient viscosite loi lineaire
212	! sf03 : coefficient viscosite loi n=3
213	! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3
214	! tres petit par rapport aux valeurs standards ~ 1.e10
215
216	!___________________________________________________________
217	! glissement module sliding_Bindschadler
218
219	&SLID_BINDSH
220	KWEERT = 5.0000001E-11,
221	LOIGLISS = 2,
222	COEFBMAX = 10.00000
223	/
224	! kweert : coefficent, loigliss le type de loi
225	! coefbmax : facteur de normalisation pour influence eau
226	!___________________________________________________________
227	!___________________________________________________________
228	! spinup module no_spinup
229	&SPINUP
230	ISPINUP = 1
231	/
232
233
234	! ispinup = 0 run standard
235	! ispinup = 1 equilibre temperature avec vitesses grisli
236	! ispinup = 2 conservation de la masse avec vitesses bilan
237	! ispinup = 3 equilibre temperature avec vitesses bilan
238
239	!___________________________________________________________
240	! module dragging_neem
241	&DRAG_NEEM
242	HWATSTREAM = 50.00000 ,
243	CF = 9.9999997E-06,
244	TOBMAX = 500.0000 ,
245	TOBLIM = 70000.00 ,
246	MASQUE_STREAM = mstream_15k.dat
247	/
248	! hwatstream (m) : critere de passage en stream vitesses de bilan
249	! si hwater > hwatstream
250	! cf coefficient de la loi de frottement fonction Neff
251	! seulement pour les points cotiers
252	! tobmax : (Pa) frottement maxi sous les streams
253	! toblim : (Pa) pour les iles
254	! masque_stream : nom du fichier qui contient le masque de streams autorises
255
256	Conservation de la masse avec equation advection-diffusion
257	la repartition depend de adv_frac
258	>1 -> advection seule
259	0 -> diffusion seule
260	0<*<1 -> fraction de l advection
261	-1 -> zones diffusion + zones advecttion
262	adv_frac= 2.000000
263	partie diffusive
264	le type de schema temporel diffusif depend de omega
265	0 -> explicite, 0.5 -> Crank-Nicolson
266	1 -> semi implicite, >1 -> over-implicite
267	omega = 2.500000
268
269	partie advective
270	le schéma temporel advectif dépend de mu
271	0 -> explicite, 0.5 -> Crank-Nicolson
272	1 -> semi implicite, >1 -> over-implicite
273	mu_adv = 1.000000
274	le schéma spatial dépend de upwind
275	1 -> schema upwind, 0.5 -> schema centre
276	upwind = 1.000000
277	------------------------------------------------------
278	!___________________________________________________________
279	&calving ! nom du bloc calving méthode Vincent
280
281	Hcoup = 250.0000
282	ifrange = 4
283	meth_hcoup = 0
284	/
285	! Hcoup epaisseur de coupure, le max si attache a coefbmshelf
286	! ifrange=0 -> pas de traitement particulier sur les bords
287	! ifrange=1 -> traitement de Vincent avec ice shelves frangeants partout
288	! ifrange=2 -> ice shelves frangeants seulement si bm-bmelt positif
289	! meth_hcoup pour faire eventuellement varier Hcoup avec le climat
290	! Hcoup_min= 100.0000
291

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

Download in other formats:

Original Format