Context Navigation

orchidee.default @ 7485

Last change on this file since 7485 was 7485, checked in by josefine.ghattas, 2 years ago
Changed keyword PARAM_FILE into PARAM_SPMIP_FILE to more explicit.
File size: 73.0 KB

Line
1	#
2	#
3	# WARNING !!!
4	# DO NOT MODIFY THIS FILE.
5	# THIS FILE IS ONLY PROVIDING INFORMATION ABOUT DEFAULT PARAMETER SETTINGS IN ORCHIDEE.
6	#
7	#*******************************************************************************************
8	# Namelist for ORCHIDEE
9	#*******************************************************************************************
10	#
11	# For more details, see : http://forge.ipsl.jussieu.fr/orchidee/wiki/Documentation/OrchideeParameters
12	#
13	# Note : [m] : meters; [K] : Kelvin degrees; [C] : Celsius degrees
14	#
15
16	#*******************************************************************************************
17	# ORCHIDEE driver parameters (read in Off-line mode only)
18	#*******************************************************************************************
19
20	# LWDOWN_CONS ([FLAG]) : Conserve longwave downwelling radiation in the forcing {}
21	LWDOWN_CONS = n
22
23	# SWDOWN_CONS ([FLAG]) : Conserve shortwave downwelling radiation in the forcing {}
24	SWDOWN_CONS = LWDOWN_CONS
25
26	# FORCING_FILE ([FILE] ) : Name of file containing the forcing data {[-]}
27	FORCING_FILE = forcing_file.nc
28
29	# DT_SECHIBA ([seconds]) : Time-step of the SECHIBA component {NOT(WEATHERGEN)}
30	DT_SECHIBA = 1800.
31
32	# RESTART_FILEIN ([FILE]) : Name of restart to READ for initial conditions {[-]}
33	RESTART_FILEIN = NONE
34
35	# RESTART_FILEOUT ([FILE]) : Name of restart files to be created by the driver {[-]}
36	RESTART_FILEOUT = driver_rest_out.nc
37
38	# DRIVER_reset_time ([FLAG]) : Overwrite time values from the driver restart file {[-]}
39	DRIVER_reset_time = n
40
41	# TIME_SKIP ([seconds, days, months, years]) : Time in the forcing file at which the model is started. {[-]}
42	TIME_SKIP = 0
43
44	# TIME_LENGTH ([seconds, days, months, years]) : Length of the integration in time. {[-]}
45	TIME_LENGTH = Full length of the forcing file
46
47	# RELAXATION ([FLAG]) : method of forcing {[-]}
48	RELAXATION = n
49
50	# RELAX_A ([days?]) : Time constant of the relaxation layer {RELAXATION}
51	RELAX_A = 1.0
52
53	# SPREAD_PREC ([nb of dt_sechiba]) : On how long we spread the precipitation {[orchidee_ol]}
54	SPREAD_PREC = Half of the forcing time step or uniform, depending on dt_force and dt_sechiba
55
56	# ATM_CO2 ([ppm]) : Value to precribe atmosoheric CO2 {[FORCE_CO2_VEG=y or Offline mode]}
57	ATM_CO2 = 350.
58
59	# SPREAD_PREC ([FLAG]) : Spread the precipitation. {orchideedriver}
60	SPREAD_PREC = Spread the precipitation.
61
62	# SPREAD_PREC_SEC ([-]) : On how long we spread the precipitation {orchideedriver and nb_spread < 0}
63	SPREAD_PREC_SEC = 3600
64
65	# SPREAD_PREC_CONT ([s]) : Take into account precipitation on the next forcing step for spreading it. {orchideedriver}
66	SPREAD_PREC_CONT = FALSE
67
68	# ALLOW_WEATHERGEN ([FLAG]) : Allow weather generator to create data {[-]}
69	ALLOW_WEATHERGEN = n
70
71	# DT_WEATHGEN ([seconds]) : Calling frequency of weather generator {ALLOW_WEATHERGEN}
72	DT_WEATHGEN = 1800.
73
74	# LIMIT_WEST ([Degrees] ) : Western limit of region {[-]}
75	LIMIT_WEST = -180.
76
77	# LIMIT_EAST ([Degrees] ) : Eastern limit of region {[-]}
78	LIMIT_EAST = 180.
79
80	# LIMIT_NORTH ([Degrees]) : Northern limit of region {[-]}
81	LIMIT_NORTH = 90.
82
83	# LIMIT_SOUTH ([Degrees]) : Southern limit of region {[-]}
84	LIMIT_SOUTH = -90.
85
86	# MERID_RES ([Degrees]) : North-South Resolution {ALLOW_WEATHERGEN}
87	MERID_RES = 2.
88
89	# ZONAL_RES ([Degrees] ) : East-West Resolution {ALLOW_WEATHERGEN}
90	ZONAL_RES = 2.
91
92	# HEIGHT_LEV1 ([m]) : Height at which T and Q are given {offline mode}
93	HEIGHT_LEV1 = 2.0
94
95	# HEIGHT_LEVW ([m]) : Height at which the wind is given {offline mode}
96	HEIGHT_LEVW = 10.0
97
98	# NBUFF (-) : Number of time steps of data to buffer between each reading of the forcing file {OFF_LINE}
99	NBUFF = 1
100
101	# IPPREC ([-] ) : Use prescribed values {ALLOW_WEATHERGEN}
102	IPPREC = 0
103
104	# WEATHGEN_PRECIP_EXACT ([FLAG]) : Exact monthly precipitation {ALLOW_WEATHERGEN}
105	WEATHGEN_PRECIP_EXACT = n
106
107	# DUMP_WEATHER ([FLAG]) : Write weather from generator into a forcing file {ALLOW_WEATHERGEN }
108	DUMP_WEATHER = n
109
110	# DUMP_WEATHER_FILE ([FILE]) : Name of the file that contains the weather from generator {DUMP_WEATHER}
111	DUMP_WEATHER_FILE = weather_dump.nc
112
113	# DUMP_WEATHER_GATHERED ([FLAG]) : Dump weather data on gathered grid {DUMP_WEATHER}
114	DUMP_WEATHER_GATHERED = y
115
116	# HEIGHT_LEV1_DUMP ([m]) : {DUMP_WEATHER}
117	HEIGHT_LEV1_DUMP = 10.
118
119	#*******************************************************************************************
120	# ORCHIDEE parameters
121	#*******************************************************************************************
122
123	# SOILTYPE_CLASSIF ([-]) : Type of soil texture classification (for hydraulic and thermic properties) {!IMPOSE_VEG}
124	SOILTYPE_CLASSIF = zobler
125
126	# RIVER_ROUTING ([FLAG]) : Decides if we route the water or not {OK_SECHIBA}
127	RIVER_ROUTING = y
128
129	# DO_IRRIGATION ([FLAG]) : Should we compute an irrigation flux {RIVER_ROUTING }
130	DO_IRRIGATION = n
131
132	# DO_FLOODPLAINS ([FLAG] ) : Should we include floodplains {RIVER_ROUTING }
133	DO_FLOODPLAINS = n
134
135	# STOMATE_OK_STOMATE ([FLAG]) : Activate STOMATE? {OK_SECHIBA}
136	STOMATE_OK_STOMATE = y
137
138	# DO_WOOD_HARVEST ([FLAG]) : Activate Wood Harvest ? {OK_STOMATE}
139	DO_WOOD_HARVEST = y
140
141	# STOMATE_OK_DGVM ([FLAG]) : Activate DGVM? {OK_STOMATE}
142	STOMATE_OK_DGVM = n
143
144	# CHEMISTRY_BVOC ([FLAG]) : Activate calculations for BVOC {OK_SECHIBA}
145	CHEMISTRY_BVOC = n
146
147	# CHEMISTRY_LEAFAGE ([FLAG]) : Activate LEAFAGE? {CHEMISTRY_BVOC}
148	CHEMISTRY_LEAFAGE = n
149
150	# CANOPY_EXTINCTION ([FLAG]) : Use canopy radiative transfer model? {CHEMISTRY_BVOC }
151	CANOPY_EXTINCTION = n
152
153	# CANOPY_MULTILAYER ([FLAG]) : Use canopy radiative transfer model with multi-layers {CANOPY_EXTINCTION }
154	CANOPY_MULTILAYER = n
155
156	# NOx_RAIN_PULSE ([FLAG]) : Calculate NOx emissions with pulse? {CHEMISTRY_BVOC }
157	NOx_RAIN_PULSE = n
158
159	# NOx_BBG_FERTIL ([FLAG]) : Calculate NOx emissions with bbg fertilizing effect? {CHEMISTRY_BVOC }
160	NOx_BBG_FERTIL = n
161
162	# NOx_FERTILIZERS_USE ([FLAG] ) : Calculate NOx emissions with fertilizers use? {CHEMISTRY_BVOC }
163	NOx_FERTILIZERS_USE = n
164
165	# NVM ([-]) : number of PFTs {OK_SECHIBA or OK_STOMATE}
166	NVM = 13
167
168	# IMPOSE_PARAM ([FLAG]) : Do you impose the values of the parameters? {OK_SECHIBA or OK_STOMATE}
169	IMPOSE_PARAM = y
170
171	# DEPTH_MAX_T (m) : Maximum depth of the soil thermodynamics {}
172	DEPTH_MAX_T = 90.0
173
174	# DEPTH_MAX_H (m) : Maximum depth of soil moisture {}
175	DEPTH_MAX_H = 2.0
176
177	# DEPTH_TOPTHICK (m) : Thickness of upper most Layer {}
178	DEPTH_TOPTHICK = 9.77517107e-04
179
180	# DEPTH_CSTTHICK (m) : Depth at which constant layer thickness start {}
181	DEPTH_CSTTHICK = DEPTH_MAX_H
182
183	# REFINEBOTTOM (-) : Depth at which the hydrology layers will be refined towards the bottom. {}
184	REFINEBOTTOM = .FALSE.
185
186	# DEPTH_GEOM (m) : Depth at which we resume geometrical increases for temperature {}
187	DEPTH_GEOM = DEPTH_MAX_H
188
189	# RATIO_GEOM_BELOW (-) : Ratio of the geometrical series defining the thickness below DEPTH_GEOM {}
190	RATIO_GEOM_BELOW = 2
191
192	# ALMA_OUTPUT ([FLAG]) : Should the output follow the ALMA convention {OK_SECHIBA}
193	ALMA_OUTPUT = n
194
195	# OUTPUT_FILE ([FILE]) : Name of file in which the output is going to be written {OK_SECHIBA}
196	OUTPUT_FILE = sechiba_history.nc
197
198	# WRITE_STEP ([seconds]) : Frequency in seconds for sechiba_history.nc file with IOIPSL {OK_SECHIBA, NOT XIOS_ORCHIDEE_OK}
199	WRITE_STEP = 86400.
200
201	# SECHIBA_HISTLEVEL ([-]) : SECHIBA history output level (0..10) {OK_SECHIBA and HF}
202	SECHIBA_HISTLEVEL = 5
203
204	# SECHIBA_HISTFILE2 ([FLAG]) : Flag to switch on histfile 2 for SECHIBA (hi-frequency ?) {OK_SECHIBA}
205	SECHIBA_HISTFILE2 = n
206
207	# WRITE_STEP2 ([seconds]) : Frequency in seconds at which to WRITE output {SECHIBA_HISTFILE2}
208	WRITE_STEP2 = 1800.0
209
210	# SECHIBA_OUTPUT_FILE2 ([FILE]) : Name of file in which the output number 2 is going to be written {SECHIBA_HISTFILE2}
211	SECHIBA_OUTPUT_FILE2 = sechiba_out_2.nc
212
213	# SECHIBA_HISTLEVEL2 ([-] ) : SECHIBA history 2 output level (0..10) {SECHIBA_HISTFILE2}
214	SECHIBA_HISTLEVEL2 = 1
215
216	# STOMATE_OUTPUT_FILE ([FILE]) : Name of file in which STOMATE's output is going to be written {OK_STOMATE}
217	STOMATE_OUTPUT_FILE = stomate_history.nc
218
219	# STOMATE_HIST_DT ([days]) : STOMATE history time step {OK_STOMATE}
220	STOMATE_HIST_DT = 10.
221
222	# STOMATE_IPCC_OUTPUT_FILE ([FILE]) : Name of file in which STOMATE's output is going to be written {OK_STOMATE}
223	STOMATE_IPCC_OUTPUT_FILE = stomate_ipcc_history.nc
224
225	# STOMATE_IPCC_HIST_DT ([days]) : STOMATE IPCC history time step {OK_STOMATE}
226	STOMATE_IPCC_HIST_DT = 0.
227
228	# OK_HISTSYNC ([FLAG]) : Syncronize and write IOIPSL output files at each time step {}
229	OK_HISTSYNC = FALSE
230
231	# STOMATE_HISTLEVEL ([-]) : STOMATE history output level (0..10) {OK_STOMATE}
232	STOMATE_HISTLEVEL = 10
233
234	# SECHIBA_restart_in ([FILE]) : Name of restart to READ for initial conditions {OK_SECHIBA }
235	SECHIBA_restart_in = NONE
236
237	# SECHIBA_rest_out ([FILE]) : Name of restart files to be created by SECHIBA {OK_SECHIBA}
238	SECHIBA_rest_out = sechiba_rest_out.nc
239
240	# STOMATE_RESTART_FILEIN ([FILE]) : Name of restart to READ for initial conditions of STOMATE {STOMATE_OK_STOMATE}
241	STOMATE_RESTART_FILEIN = NONE
242
243	# STOMATE_RESTART_FILEOUT ([FILE]) : Name of restart files to be created by STOMATE {STOMATE_OK_STOMATE}
244	STOMATE_RESTART_FILEOUT = stomate_rest_out.nc
245
246	# FORCE_CO2_VEG ([FLAG]) : Flag to force the value of atmospheric CO2 for vegetation. {Only in coupled mode}
247	FORCE_CO2_VEG = FALSE
248
249	# TAU_OUTFLOW ([days]) : Number of days over which the coastal- and riverflow will be distributed {Only in coupled mode}
250	TAU_OUTFLOW = 0
251
252	# ECCENTRICITY ([-]) : Use prescribed values {ALLOW_WEATHERGEN}
253	ECCENTRICITY = 0.016724
254
255	# PERIHELIE ([-]) : Use prescribed values {ALLOW_WEATHERGEN}
256	PERIHELIE = 102.04
257
258	# OBLIQUITY ([Degrees]) : Use prescribed values {ALLOW_WEATHERGEN}
259	OBLIQUITY = 23.446
260
261	# PFT_TO_MTC ([-]) : correspondance array linking a PFT to MTC {OK_SECHIBA or OK_STOMATE}
262	PFT_TO_MTC = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13
263
264	# PFT_NAME ([-]) : Name of a PFT {OK_SECHIBA or OK_STOMATE}
265	PFT_NAME = bare ground, tropical broad-leaved evergreen, tropical broad-leaved raingreen, temperate needleleaf evergreen, temperate broad-leaved evergreen, temperate broad-leaved summergreen, boreal needleleaf evergreen, boreal broad-leaved summergreen, boreal needleleaf summergreen, C3 grass, C4 grass, C3 agriculture, C4 agriculture
266
267	# LEAF_TAB ([-] ) : leaf type : 1 {OK_STOMATE}
268	LEAF_TAB = 4, 1, 1, 2, 1, 1, 2, 1, 2, 3, 3, 3, 3
269
270	# PHENO_MODEL ([-] ) : which phenology model is used? (tabulated) {OK_STOMATE}
271	PHENO_MODEL = none, none, moi, none, none, ncdgdd, none, ncdgdd, ngd, moigdd, moigdd, moigdd, moigdd
272
273	# SECHIBA_LAI ([m^2/m^2]) : laimax for maximum lai(see also type of lai interpolation) {OK_SECHIBA or IMPOSE_VEG}
274	SECHIBA_LAI = 0., 8., 8., 4., 4.5, 4.5, 4., 4.5, 4., 2., 2., 2., 2.
275
276	# LLAIMIN ([m^2/m^2]) : laimin for minimum lai(see also type of lai interpolation) {OK_SECHIBA or IMPOSE_VEG}
277	LLAIMIN = 0., 8., 0., 4., 4.5, 0., 4., 0., 0., 0., 0., 0., 0.
278
279	# SLOWPROC_HEIGHT ([m] ) : prescribed height of vegetation {OK_SECHIBA}
280	SLOWPROC_HEIGHT = 0., 30., 30., 20., 20., 20., 15., 15., 15., .5, .6, 1., 1.
281
282	# Z0_OVER_HEIGHT ([-] ) : factor to calculate roughness height from height of canopy {OK_SECHIBA}
283	Z0_OVER_HEIGHT = 0., 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625
284
285	# RATIO_Z0M_Z0H ([-]) : Ratio between z0m and z0h {OK_SECHIBA}
286	RATIO_Z0M_Z0H = 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0
287
288	# TYPE_OF_LAI ([-]) : Type of behaviour of the LAI evolution algorithm {OK_SECHIBA}
289	TYPE_OF_LAI = inter, inter, inter, inter, inter, inter, inter, inter, inter, inter, inter, inter, inter
290
291	# NATURAL ([BOOLEAN]) : natural? {OK_SECHIBA, OK_STOMATE}
292	NATURAL = y, y, y, y, y, y, y, y, y, y, y, n, n
293
294	# IS_C4 ([BOOLEAN]) : flag for C4 vegetation types {OK_SECHIBA or OK_STOMATE}
295	IS_C4 = n, n, n, n, n, n, n, n, n, n, n, y, n, y
296
297	# VCMAX_FIX ([micromol/m^2/s] ) : values used for vcmax when STOMATE is not activated {OK_SECHIBA and NOT(OK_STOMATE)}
298	VCMAX_FIX = 0., 40., 50., 30., 35., 40.,30., 40., 35., 60., 60., 70., 70.
299
300	# DOWNREG_CO2 ([-]) : coefficient for CO2 downregulation (unitless) {OK_CO2 and DOWNREGULATION_CO2}
301	DOWNREG_CO2 = 0., 0.38, 0.38, 0.28, 0.28, 0.28, 0.22, 0.22, 0.22, 0.26, 0.26, 0.26, 0.26
302
303	# DOWNREG_CO2_NEW ([-]) : coefficient for CO2 downregulation (unitless) {OK_CO2 and DOWNREGULATION_CO2_NEW}
304	DOWNREG_CO2_NEW = 0., 0.35, 0.35, 0.26, 0.26, 0.26, 0.20, 0.20, 0.20, 0.24, 0.03, 0.24, 0.03
305
306	# E_KmC ([J mol-1]) : Energy of activation for KmC {}
307	E_KmC = -9999., 79430., 79430., 79430., 79430., 79430., 79430., 79430., 79430., 79430., 79430., 79430., 79430.
308
309	# E_KmO ([J mol-1]) : Energy of activation for KmO {}
310	E_KmO = -9999., 36380., 36380., 36380., 36380., 36380., 36380., 36380., 36380., 36380., 36380., 36380., 36380.
311
312	# E_Sco ([J mol-1]) : Energy of activation for Sco {}
313	E_Sco = -9999., -24460., -24460., -24460., -24460., -24460., -24460., -24460., -24460., -24460., -24460., -24460., -24460.
314
315	# E_gamma_star ([J mol-1]) : Energy of activation for gamma_star {}
316	E_gamma_star = -9999., 37830., 37830., 37830., 37830., 37830., 37830., 37830., 37830., 37830., 37830., 37830., 37830.
317
318	# E_Vcmax ([J mol-1]) : Energy of activation for Vcmax {}
319	E_Vcmax = -9999., 71513., 71513., 71513., 71513., 71513., 71513., 71513., 71513., 71513., 67300., 71513., 67300.
320
321	# E_Jmax ([J mol-1]) : Energy of activation for Jmax {}
322	E_Jmax = -9999., 49884., 49884., 49884., 49884., 49884., 49884., 49884., 49884., 49884., 77900., 49884., 77900.
323
324	# aSV ([J K-1 mol-1]) : a coefficient of the linear regression (a+bT) defining the Entropy term for Vcmax {}
325	aSV = -9999., 668.39, 668.39, 668.39, 668.39, 668.39, 668.39, 668.39, 668.39, 668.39, 641.64, 668.39, 641.64
326
327	# bSV ([J K-1 mol-1 Â°C-1]) : b coefficient of the linear regression (a+bT) defining the Entropy term for Vcmax {}
328	bSV = -9999., -1.07, -1.07, -1.07, -1.07, -1.07, -1.07, -1.07, -1.07, -1.07, 0., -1.07, 0.
329
330	# TPHOTO_MIN ([-]) : minimum photosynthesis temperature (deg C) {OK_STOMATE}
331	TPHOTO_MIN = -9999., -4., -4., -4., -4.,-4.,-4., -4., -4., -4., -4., -4., -4.
332
333	# TPHOTO_MAX ([-]) : maximum photosynthesis temperature (deg C) {OK_STOMATE}
334	TPHOTO_MAX = -9999., 55., 55., 55., 55., 55., 55., 55., 55., 55., 55., 55., 55.
335
336	# aSJ ([J K-1 mol-1]) : a coefficient of the linear regression (a+bT) defining the Entropy term for Jmax {}
337	aSJ = -9999., 659.70, 659.70, 659.70, 659.70, 659.70, 659.70, 659.70, 659.70, 659.70, 630., 659.70, 630.
338
339	# bSJ ([J K-1 mol-1 Â°C-1]) : b coefficient of the linear regression (a+bT) defining the Entropy term for Jmax {}
340	bSJ = -9999., -0.75, -0.75, -0.75, -0.75, -0.75, -0.75, -0.75, -0.75, -0.75, 0., -0.75, 0.
341
342	# D_Vcmax ([J mol-1]) : Energy of deactivation for Vcmax {}
343	D_Vcmax = -9999., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 192000., 200000., 192000.
344
345	# D_Jmax ([J mol-1]) : Energy of deactivation for Jmax {}
346	D_Jmax = -9999., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 192000., 200000., 192000.
347
348	# E_gm ([J mol-1] ) : Energy of activation for gm { }
349	E_gm = -9999., 49600., 49600., 49600., 49600., 49600., 49600., 49600., 49600., 49600., -9999., 49600., -9999.
350
351	# S_gm ([J K-1 mol-1] ) : Entropy term for gm { }
352	S_gm = -9999., 1400., 1400., 1400., 1400., 1400., 1400., 1400., 1400., 1400., -9999., 1400., -9999.
353
354	# D_gm ([J mol-1] ) : Energy of deactivation for gm { }
355	D_gm = -9999., 437400., 437400., 437400., 437400., 437400., 437400., 437400., 437400., 437400., -9999., 437400., -9999.
356
357	# E_Rd ([J mol-1]) : Energy of activation for Rd {}
358	E_Rd = -9999., 46390., 46390., 46390., 46390., 46390., 46390., 46390., 46390., 46390., 46390., 46390., 46390.
359
360	# VCMAX25 ([micromol/m^2/s]) : Maximum rate of Rubisco activity-limited carboxylation at 25Â°C {OK_STOMATE}
361	VCMAX25 = -9999., 45.0, 45.0, 35.0, 40.0, 50.0, 45.0, 35.0, 35.0, 50.0, 50.0, 60.0, 60.0
362
363	# ARJV ([mu mol e- (mu mol CO2)-1]) : a coefficient of the linear regression (a+bT) defining the Jmax25/Vcmax25 ratio {OK_STOMATE}
364	ARJV = -9999., 2.59, 2.59, 2.59, 2.59, 2.59, 2.59, 2.59, 2.59, 2.59, 1.715, 2.59, 1.715
365
366	# BRJV ([(mu mol e- (mu mol CO2)-1) (Â°C)-1]) : b coefficient of the linear regression (a+bT) defining the Jmax25/Vcmax25 ratio {OK_STOMATE}
367	BRJV = -9999., -0.035, -0.035, -0.035, -0.035, -0.035, -0.035, -0.035, -0.035, -0.035, 0., -0.035, 0.
368
369	# KmC25 ([ubar]) : MichaelisâMenten constant of Rubisco for CO2 at 25Â°C {}
370	KmC25 = -9999., 404.9, 404.9, 404.9, 404.9, 404.9, 404.9, 404.9, 404.9, 404.9, 650., 404.9, 650.
371
372	# KmO25 ([ubar]) : MichaelisâMenten constant of Rubisco for O2 at 25Â°C {}
373	KmO25 = -9999., 278400., 278400., 278400., 278400., 278400., 278400., 278400., 278400., 278400., 450000., 278400., 450000.
374
375	# Sco25 ([bar bar-1]) : Relative CO2 /O2 specificity factor for Rubisco at 25ÃÂ°C {}
376	Sco25 = -9999., 2800., 2800., 2800., 2800., 2800., 2800., 2800., 2800., 2800., 2590., 2800., 2590.
377
378	# gm25 ([mol m-2 s-1 bar-1] ) : Mesophyll diffusion conductance at 25ÃÂÃÂ°C { }
379	gm25 = -9999., 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, -9999., 0.4, -9999.
380
381	# gamma_star25 ([ubar]) : Ci-based CO2 compensation point in the absence of Rd at 25Â°C (ubar) {}
382	gamma_star25 = -9999., 42.75, 42.75, 42.75, 42.75, 42.75, 42.75, 42.75, 42.75, 42.75, 42.75, 42.75, 42.75
383
384	# a1 ([-]) : Empirical factor involved in the calculation of fvpd {}
385	a1 = -9999., 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.72, 0.85, 0.72
386
387	# b1 ([-]) : Empirical factor involved in the calculation of fvpd {}
388	b1 = -9999., 0.14, 0.14, 0.14, 0.14, 0.14, 0.14, 0.14, 0.14, 0.14, 0.20, 0.14, 0.20
389
390	# g0 ([mol mâ2 sâ1 barâ1]) : Residual stomatal conductance when irradiance approaches zero {}
391	g0 = -9999., 0.00625, 0.00625, 0.00625, 0.00625, 0.00625, 0.00625, 0.00625, 0.00625, 0.00625, 0.01875, 0.00625, 0.01875
392
393	# h_protons ([mol mol-1]) : Number of protons required to produce one ATP {}
394	h_protons = -9999., 4., 4., 4., 4., 4., 4., 4., 4., 4., 4., 4., 4.
395
396	# fpsir ([-]) : Fraction of PSII eâ transport rate partitioned to the C4 cycle {}
397	fpsir = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.4, -9999., 0.4
398
399	# fQ ([-]) : Fraction of electrons at reduced plastoquinone that follow the Q-cycle {}
400	fQ = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 1., -9999., 1.
401
402	# fpseudo ([-]) : Fraction of electrons at PSI that follow pseudocyclic transport {}
403	fpseudo = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.1, -9999., 0.1
404
405	# kp ([mol mâ2 sâ1 barâ1]) : Initial carboxylation efficiency of the PEP carboxylase {}
406	kp = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.7, -9999., 0.7
407
408	# alpha ([-]) : Fraction of PSII activity in the bundle sheath {}
409	alpha = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.1, -9999., 0.1
410
411	# gbs ([mol mâ2 sâ1 barâ1]) : Bundle-sheath conductance {}
412	gbs = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.003, -9999., 0.003
413
414	# theta ([â]) : Convexity factor for response of J to irradiance {}
415	theta = -9999., 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7
416
417	# alpha_LL ([mol eâ (mol photon)â1]) : Conversion efficiency of absorbed light into J at strictly limiting light {}
418	alpha_LL = -9999., 0.372, 0.372, 0.372, 0.372, 0.372, 0.372, 0.372, 0.372, 0.372, 0.372, 0.372, 0.372
419
420	# STRESS_VCMAX ([-]) : Stress on vcmax {OK_SECHIBA or OK_STOMATE}
421	STRESS_VCMAX = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.
422
423	# STRESS_GS ([-]) : Stress on gs {OK_SECHIBA or OK_STOMATE}
424	STRESS_GS = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.
425
426	# STRESS_GM ([-]) : Stress on gm {OK_SECHIBA or OK_STOMATE}
427	STRESS_GM = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.
428
429	# EXT_COEFF ([-]) : extinction coefficient of the Monsi&Seaki relationship (1953) {OK_SECHIBA or OK_STOMATE}
430	EXT_COEFF = .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5
431
432	# EXT_COEFF_VEGETFRAC ([-]) : extinction coefficient used for the calculation of the bare soil fraction {OK_SECHIBA or OK_STOMATE}
433	EXT_COEFF_VEGETFRAC = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.
434
435	# HYDROL_HUMCSTE ([m]) : Root profile {OK_SECHIBA}
436	HYDROL_HUMCSTE = humcste_ref2m or humcste_ref4m depending on zmaxh
437
438	# PREF_SOIL_VEG ([-] ) : The soil tile number for each vegetation {OK_SECHIBA or OK_STOMATE}
439	PREF_SOIL_VEG = 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3
440
441	# RSTRUCT_CONST ([s/m]) : Structural resistance {OK_SECHIBA}
442	RSTRUCT_CONST = 0.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, 2.5, 2.0, 2.0, 2.0
443
444	# KZERO ([kg/m^2/s]) : A vegetation dependent constant used in the calculation of the surface resistance. {OK_SECHIBA}
445	KZERO = 0.0, 12.E-5, 12.E-5, 12.e-5, 12.e-5, 25.e-5, 12.e-5,25.e-5, 25.e-5, 30.e-5, 30.e-5, 30.e-5, 30.e-5
446
447	# RVEG_PFT ([-]) : Artificial parameter to increase or decrease canopy resistance. {OK_SECHIBA}
448	RVEG_PFT = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.
449
450	# WMAX_VEG ([kg/m^3]) : Maximum field capacity for each of the vegetations (Temporary): max quantity of water {OK_SECHIBA}
451	WMAX_VEG = 150., 150., 150., 150., 150., 150., 150.,150., 150., 150., 150., 150., 150.
452
453	# PERCENT_THROUGHFALL_PFT ([%]) : Percent by PFT of precip that is not intercepted by the canopy. Default value depend on run mode. {OK_SECHIBA}
454	PERCENT_THROUGHFALL_PFT = Case offline [0. 0. 0....] else [30. 30. 30.....]
455
456	# SNOWA_AGED_VIS ([-]) : Minimum snow albedo value for each vegetation type after aging (dirty old snow), visible albedo {OK_SECHIBA}
457	SNOWA_AGED_VIS = 0.74, 0.0, 0.0, 0.08, 0.24, 0.07, 0.18, 0.18, 0.33, 0.57, 0.57, 0.57, 0.57
458
459	# SNOWA_AGED_NIR ([-]) : Minimum snow albedo value for each vegetation type after aging (dirty old snow), near infrared albedo {OK_SECHIBA}
460	SNOWA_AGED_NIR = 0.50, 0.0, 0.0, 0.10, 0.37, 0.08, 0.16, 0.17, 0.27, 0.44, 0.44, 0.44, 0.44
461
462	# SNOWA_DEC_VIS ([-]) : Decay rate of snow albedo value for each vegetation type as it will be used in condveg_snow, visible albedo {OK_SECHIBA}
463	SNOWA_DEC_VIS = 0.21, 0.0, 0.0, 0.14, 0.08, 0.17, 0.05, 0.06, 0.09, 0.15, 0.15, 0.15, 0.15
464
465	# SNOWA_DEC_NIR ([-]) : Decay rate of snow albedo value for each vegetation type as it will be used in condveg_snow, near infrared albedo {OK_SECHIBA}
466	SNOWA_DEC_NIR = 0.13, 0.0, 0.0, 0.10, 0.10, 0.16, 0.04, 0.07, 0.08, 0.12, 0.12, 0.12, 0.12
467
468	# ALB_LEAF_VIS ([-]) : leaf albedo of vegetation type, visible albedo {OK_SECHIBA}
469	ALB_LEAF_VIS = 0.00, 0.04, 0.04, 0.04, 0.04, 0.03, 0.03, 0.03, 0.03, 0.06, 0.06, 0.06, 0.06
470
471	# ALB_LEAF_NIR ([-]) : leaf albedo of vegetation type, near infrared albedo {OK_SECHIBA}
472	ALB_LEAF_NIR = 0.00, 0.23, 0.18, 0.18, 0.20, 0.24, 0.15, 0.26, 0.20, 0.24, 0.27, 0.28, 0.26
473
474	# ISO_ACTIVITY ([-]) : Biogenic activity for each age class : isoprene {CHEMISTRY_BVOC}
475	ISO_ACTIVITY = 0.5, 1.5, 1.5, 0.5
476
477	# METHANOL_ACTIVITY ([-]) : Isoprene emission factor for each age class : methanol {CHEMISTRY_BVOC}
478	METHANOL_ACTIVITY = 1., 1., 0.5, 0.5
479
480	# EM_FACTOR_ISOPRENE ([ugC/g/h] ) : Isoprene emission factor {CHEMISTRY_BVOC}
481	EM_FACTOR_ISOPRENE = 0., 24., 24., 8., 16., 45., 8., 18., 0.5, 12., 18., 5., 5.
482
483	# EM_FACTOR_MONOTERPENE ([ugC/g/h] ) : Monoterpene emission factor {CHEMISTRY_BVOC }
484	EM_FACTOR_MONOTERPENE = 0., 2.0, 2.0, 1.8, 1.4, 1.6, 1.8, 1.4, 1.8, 0.8, 0.8, 0.22, 0.22
485
486	# C_LDF_MONO ([]) : Monoterpenes fraction dependancy to light {CHEMISTRY_BVOC}
487	C_LDF_MONO = 0.6
488
489	# C_LDF_SESQ ([]) : Sesquiterpenes fraction dependancy to light {CHEMISTRY_BVOC}
490	C_LDF_SESQ = 0.5
491
492	# C_LDF_METH ([]) : Methanol fraction dependancy to light {CHEMISTRY_BVOC}
493	C_LDF_METH = 0.8
494
495	# C_LDF_ACET ([]) : Acetone fraction dependancy to light {CHEMISTRY_BVOC}
496	C_LDF_ACET = 0.2
497
498	# EM_FACTOR_APINENE ([ugC/g/h] ) : Alfa pinene emission factor {CHEMISTRY_BVOC }
499	EM_FACTOR_APINENE = 0., 1.35, 1.35, 0.85, 0.95, 0.75, 0.85, 0.60, 1.98, 0.30, 0.30, 0.09, 0.09
500
501	# EM_FACTOR_BPINENE ([ugC/g/h] ) : Beta pinene emission factor {CHEMISTRY_BVOC }
502	EM_FACTOR_BPINENE = 0., 0.30, 0.30, 0.35, 0.25, 0.20, 0.35, 0.12, 0.45, 0.16, 0.12, 0.05, 0.05
503
504	# EM_FACTOR_LIMONENE ([ugC/g/h] ) : Limonene emission factor {CHEMISTRY_BVOC}
505	EM_FACTOR_LIMONENE = 0., 0.25, 0.25, 0.20, 0.25, 0.14, 0.20, 0.135, 0.11, 0.19, 0.42, 0.03, 0.03
506
507	# EM_FACTOR_MYRCENE ([ugC/g/h] ) : Myrcene emission factor {CHEMISTRY_BVOC}
508	EM_FACTOR_MYRCENE = 0., 0.20, 0.20, 0.12, 0.11, 0.065, 0.12, 0.036, 0.075, 0.08, 0.085, 0.015, 0.015
509
510	# EM_FACTOR_SABINENE ([ugC/g/h] ) : Sabinene emission factor {CHEMISTRY_BVOC}
511	EM_FACTOR_SABINENE = 0., 0.20, 0.20, 0.12, 0.17, 0.70, 0.12, 0.50, 0.09, 0.085, 0.075, 0.02, 0.02
512
513	# EM_FACTOR_CAMPHENE ([ugC/g/h] ) : Camphene emission factor {CHEMISTRY_BVOC}
514	EM_FACTOR_CAMPHENE = 0., 0.15, 0.15, 0.10, 0.10, 0.01, 0.10, 0.01, 0.07, 0.07, 0.08, 0.01, 0.01
515
516	# EM_FACTOR_3CARENE ([ugC/g/h] ) : 3-Carene emission factor {CHEMISTRY_BVOC}
517	EM_FACTOR_3CARENE = 0., 0.13, 0.13, 0.42, 0.02, 0.055, 0.42,0.025, 0.125, 0.085, 0.085, 0.065, 0.065
518
519	# EM_FACTOR_TBOCIMENE ([ugC/g/h] ) : T-beta-ocimene emission factor {CHEMISTRY_BVOC}
520	EM_FACTOR_TBOCIMENE = 0., 0.25, 0.25, 0.13, 0.09, 0.26, 0.13, 0.20, 0.085, 0.18, 0.18, 0.01, 0.01
521
522	# EM_FACTOR_OTHERMONOT ([ugC/g/h] ) : Other monoterpenes emission factor {CHEMISTRY_BVOC}
523	EM_FACTOR_OTHERMONOT = 0., 0.17, 0.17, 0.11, 0.11, 0.125, 0.11, 0.274, 0.01, 0.15, 0.155, 0.035, 0.035
524
525	# EM_FACTOR_SESQUITERP ([ugC/g/h] ) : Sesquiterpenes emission factor {CHEMISTRY_BVOC}
526	EM_FACTOR_SESQUITERP = 0., 0.45, 0.45, 0.13, 0.3, 0.36, 0.15, 0.3, 0.25, 0.6, 0.6, 0.08, 0.08
527
528	# C_BETA_MONO ([]) : Monoterpenes temperature dependency coefficient {CHEMISTRY_BVOC}
529	C_BETA_MONO = 0.1
530
531	# C_BETA_SESQ ([]) : Sesquiterpenes temperature dependency coefficient {CHEMISTRY_BVOC}
532	C_BETA_SESQ = 0.17
533
534	# C_BETA_METH ([]) : Methanol temperature dependency coefficient {CHEMISTRY_BVOC}
535	C_BETA_METH = 0.08
536
537	# C_BETA_ACET ([]) : Acetone temperature dependency coefficient {CHEMISTRY_BVOC}
538	C_BETA_ACET = 0.1
539
540	# C_BETA_OXYVOC ([]) : Other oxygenated BVOC temperature dependency coefficient {CHEMISTRY_BVOC}
541	C_BETA_OXYVOC = 0.13
542
543	# EM_FACTOR_ORVOC ([ugC/g/h] ) : ORVOC emissions factor {CHEMISTRY_BVOC }
544	EM_FACTOR_ORVOC = 0., 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5
545
546	# EM_FACTOR_OVOC ([ugC/g/h] ) : OVOC emissions factor {CHEMISTRY_BVOC}
547	EM_FACTOR_OVOC = 0., 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5
548
549	# EM_FACTOR_MBO ([ugC/g/h] ) : MBO emissions factor {CHEMISTRY_BVOC }
550	EM_FACTOR_MBO = 0., 2.e-5, 2.e-5, 1.4, 2.e-5, 2.e-5, 0.14, 2.e-5, 2.e-5, 2.e-5, 2.e-5, 2.e-5, 2.e-5
551
552	# EM_FACTOR_METHANOL ([ugC/g/h] ) : Methanol emissions factor {CHEMISTRY_BVOC }
553	EM_FACTOR_METHANOL = 0., 0.8, 0.8, 1.8, 0.9, 1.9, 1.8, 1.8, 1.8, 0.7, 0.9, 2., 2.
554
555	# EM_FACTOR_ACETONE ([ugC/g/h] ) : Acetone emissions factor {CHEMISTRY_BVOC }
556	EM_FACTOR_ACETONE = 0., 0.25, 0.25, 0.3, 0.2, 0.33, 0.3, 0.25, 0.25, 0.2, 0.2, 0.08, 0.08
557
558	# EM_FACTOR_ACETAL ([ugC/g/h] ) : Acetaldehyde emissions factor {CHEMISTRY_BVOC}
559	EM_FACTOR_ACETAL = 0., 0.2, 0.2, 0.2, 0.2, 0.25, 0.25, 0.16, 0.16, 0.12, 0.12, 0.035, 0.02
560
561	# EM_FACTOR_FORMAL ([ugC/g/h] ) : Formaldehyde emissions factor {CHEMISTRY_BVOC }
562	EM_FACTOR_FORMAL = 0., 0.04, 0.04, 0.08, 0.04, 0.04, 0.04, 0.04, 0.04, 0.025, 0.025, 0.013, 0.013
563
564	# EM_FACTOR_ACETIC ([ugC/g/h] ) : Acetic Acid emissions factor {CHEMISTRY_BVOC }
565	EM_FACTOR_ACETIC = 0., 0.025, 0.025,0.025,0.022,0.08,0.025,0.022,0.013,0.012,0.012,0.008,0.008
566
567	# EM_FACTOR_FORMIC ([ugC/g/h] ) : Formic Acid emissions factor {CHEMISTRY_BVOC}
568	EM_FACTOR_FORMIC = 0., 0.015, 0.015, 0.02, 0.02, 0.025, 0.025, 0.015, 0.015,0.010,0.010,0.008,0.008
569
570	# EM_FACTOR_NO_WET ([ngN/m^2/s]) : NOx emissions factor wet soil emissions and exponential dependancy factor {CHEMISTRY_BVOC}
571	EM_FACTOR_NO_WET = 0., 2.6, 0.06, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.36, 0.36, 0.36, 0.36
572
573	# EM_FACTOR_NO_DRY ([ngN/m^2/s] ) : NOx emissions factor dry soil emissions and exponential dependancy factor {CHEMISTRY_BVOC}
574	EM_FACTOR_NO_DRY = 0., 8.60, 0.40, 0.22, 0.22, 0.22, 0.22, 0.22, 0.22, 2.65, 2.65, 2.65, 2.65
575
576	# LARCH ([-] ) : Larcher 1991 SAI/LAI ratio {CHEMISTRY_BVOC }
577	LARCH = 0., 0.015, 0.015, 0.003, 0.005, 0.005, 0.003, 0.005, 0.003, 0.005, 0.005, 0.008, 0.008
578
579	# SLA ([m^2/gC]) : specif leaf area {OK_STOMATE}
580	SLA = 1.5E-2, 1.53E-2, 2.6E-2, 9.26E-3, 2E-2, 2.6E-2, 9.26E-3, 2.6E-2, 1.9E-2, 2.6E-2, 2.6E-2, 2.6E-2, 2.6E-2
581
582	# AVAILABILITY_FACT ([-] ) : Calculate dynamic mortality in lpj_gap, pft dependent parameter {OK_STOMATE }
583	AVAILABILITY_FACT = -9999., 0.14, 0.14, 0.10, 0.10, 0.10, 0.05, 0.05, 0.05, -9999., -9999., -9999., -9999.
584
585	# R0 ([-] ) : Standard root allocation {OK_STOMATE }
586	R0 = -9999., .30, .30, .30, .30, .30, .30, .30, .30, .30, .30, .30, .30
587
588	# S0 ([-] ) : Standard sapwood allocation {OK_STOMATE }
589	S0 = -9999., .25, .25, .30, .30, .30, .30, .30, .30, .30, .30, .30, .30
590
591	# FRAC_GROWTHRESP ([-]) : fraction of GPP which is lost as growth respiration {OK_STOMATE}
592	FRAC_GROWTHRESP = -9999., 0.35, 0.35, 0.28, 0.28, 0.28, 0.35, 0.35, 0.35, 0.28, 0.28, 0.28, 0.28
593
594	# MAINT_RESP_SLOPE_C ([-]) : slope of maintenance respiration coefficient (1/K), constant c of aT^2+bT+c , tabulated {OK_STOMATE}
595	MAINT_RESP_SLOPE_C = -9999., 0.12, 0.12, 0.16, 0.16, 0.16, 0.25, 0.25, 0.25, 0.16, 0.12, 0.16, 0.12
596
597	# MAINT_RESP_SLOPE_B ([-]) : slope of maintenance respiration coefficient (1/K), constant b of aT^2+bT+c , tabulated {OK_STOMATE}
598	MAINT_RESP_SLOPE_B = -9999., .0, .0, .0, .0, .0, .0, .0, .0, -.00133, .0, -.00133, .0
599
600	# MAINT_RESP_SLOPE_A ([-]) : slope of maintenance respiration coefficient (1/K), constant a of aT^2+bT+c , tabulated {OK_STOMATE}
601	MAINT_RESP_SLOPE_A = -9999., .0, .0, .0, .0, .0, .0, .0, .0, .0, .0, .0, .0
602
603	# CM_ZERO_LEAF ([g/g/day]) : maintenance respiration coefficient at 0 deg C, for leaves, tabulated {OK_STOMATE}
604	CM_ZERO_LEAF = -9999., 2.35E-3, 2.62E-3, 1.01E-3, 2.35E-3, 2.62E-3, 1.01E-3,2.62E-3, 2.05E-3, 2.62E-3, 2.62E-3, 2.62E-3, 2.62E-3
605
606	# CM_ZERO_SAPABOVE ([g/g/day]) : maintenance respiration coefficient at 0 deg C,for sapwood above, tabulated {OK_STOMATE}
607	CM_ZERO_SAPABOVE = -9999., 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4
608
609	# CM_ZERO_SAPBELOW ([g/g/day]) : maintenance respiration coefficient at 0 deg C, for sapwood below, tabulated {OK_STOMATE}
610	CM_ZERO_SAPBELOW = -9999., 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4
611
612	# CM_ZERO_HEARTABOVE ([g/g/day]) : maintenance respiration coefficient at 0 deg C, for heartwood above, tabulated {OK_STOMATE }
613	CM_ZERO_HEARTABOVE = -9999., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.
614
615	# CM_ZERO_HEARTBELOW ([g/g/day] ) : maintenance respiration coefficient at 0 deg C,for heartwood below, tabulated {OK_STOMATE }
616	CM_ZERO_HEARTBELOW = -9999., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.
617
618	# CM_ZERO_ROOT ([g/g/day] ) : maintenance respiration coefficient at 0 deg C, for roots, tabulated {OK_STOMATE}
619	CM_ZERO_ROOT = -9999.,1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3,1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3
620
621	# CM_ZERO_FRUIT ([g/g/day] ) : maintenance respiration coefficient at 0 deg C, for fruits, tabulated {OK_STOMATE}
622	CM_ZERO_FRUIT = -9999., 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4,1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4
623
624	# CM_ZERO_CARBRES ([g/g/day] ) : maintenance respiration coefficient at 0 deg C, for carbohydrate reserve, tabulated {OK_STOMATE}
625	CM_ZERO_CARBRES = -9999., 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4,1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4
626
627	# FLAM ([-]) : flamability: critical fraction of water holding capacity {OK_STOMATE}
628	FLAM = -9999., .15, .25, .25, .25, .25, .25, .25, .25, .25, .25, .35, .35
629
630	# RESIST ([-]) : fire resistance {OK_STOMATE}
631	RESIST = -9999., .95, .90, .12, .50, .12, .12, .12, .12, .0, .0, .0, .0
632
633	# COEFF_LCCHANGE_1 ([-]) : Coeff of biomass export for the year {OK_STOMATE}
634	COEFF_LCCHANGE_1 = -9999., 0.897, 0.897, 0.597, 0.597, 0.597, 0.597, 0.597, 0.597, 0.597, 0.597, 0.597, 0.597
635
636	# COEFF_LCCHANGE_10 ([-]) : Coeff of biomass export for the decade {OK_STOMATE}
637	COEFF_LCCHANGE_10 = -9999., 0.103, 0.103, 0.299, 0.299, 0.299, 0.299, 0.299, 0.299, 0.299, 0.403, 0.299, 0.403
638
639	# COEFF_LCCHANGE_100 ([-]) : Coeff of biomass export for the century {OK_STOMATE}
640	COEFF_LCCHANGE_100 = -9999., 0., 0., 0.104, 0.104, 0.104, 0.104, 0.104, 0.104, 0.104, 0., 0.104, 0.
641
642	# LAI_MAX_TO_HAPPY ([-]) : threshold of LAI below which plant uses carbohydrate reserves {OK_STOMATE}
643	LAI_MAX_TO_HAPPY = -9999., .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5
644
645	# LAI_MAX ([m^2/m^2]) : maximum LAI, PFT-specific {OK_STOMATE}
646	LAI_MAX = -9999., 7.0, 5.0, 5.0, 4.0, 5.0, 3.5, 4.0, 3.0, 2.5, 2.0, 5.0, 5.0
647
648	# PHENO_TYPE ([-]) : type of phenology, 0 {OK_STOMATE}
649	PHENO_TYPE = 0, 1, 3, 1, 1, 2, 1, 2, 2, 4, 4, 2, 3
650
651	# PHENO_GDD_CRIT_C ([-]) : critical gdd, tabulated (C), constant c of aT^2+bT+c {OK_STOMATE}
652	PHENO_GDD_CRIT_C = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 270., 400., 125., 400.
653
654	# PHENO_GDD_CRIT_B ([-]) : critical gdd, tabulated (C), constant b of aT^2+bT+c {OK_STOMATE}
655	PHENO_GDD_CRIT_B = -9999., -9999., -9999., -9999., -9999., -9999., -9999.,-9999., -9999., 6.25, 0., 0., 0.
656
657	# PHENO_GDD_CRIT_A ([-]) : critical gdd, tabulated (C), constant a of aT^2+bT+c {OK_STOMATE}
658	PHENO_GDD_CRIT_A = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.03125, 0., 0., 0.
659
660	# PHENO_MOIGDD_T_CRIT ([C]) : Average temperature threashold for C4 grass used in pheno_moigdd {OK_STOMATE}
661	PHENO_MOIGDD_T_CRIT = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 22.0, -9999., -9999.
662
663	# NGD_CRIT ([days]) : critical ngd, tabulated. Threshold -5 degrees {OK_STOMATE}
664	NGD_CRIT = -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0., -9999., -9999., -9999., -9999., -9999.
665
666	# NCDGDD_TEMP ([C] ) : critical temperature for the ncd vs. gdd function in phenology {OK_STOMATE}
667	NCDGDD_TEMP = -9999., -9999., -9999., -9999., -9999., 5., -9999., 0., -9999., -9999., -9999., -9999., -9999.
668
669	# HUM_FRAC ([%]) : critical humidity (relative to min/max) for phenology {OK_STOMATE}
670	HUM_FRAC = -9999., -9999., .5, -9999., -9999., -9999., -9999., -9999., -9999., .5, .5, .5,.5
671
672	# HUM_MIN_TIME ([days]) : minimum time elapsed since moisture minimum {OK_STOMATE}
673	HUM_MIN_TIME = -9999., -9999., 50., -9999., -9999., -9999., -9999., -9999., -9999., 35., 35., 75., 75.
674
675	# TAU_SAP ([days]) : sapwood -> heartwood conversion time {OK_STOMATE}
676	TAU_SAP = -9999., 730., 730., 730., 730., 730., 730., 730., 730., -9999., -9999., -9999., -9999.
677
678	# TAU_LEAFINIT ([days]) : time to attain the initial foliage using the carbohydrate reserve {OK_STOMATE}
679	TAU_LEAFINIT = -9999., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10.
680
681	# TAU_FRUIT ([days]) : fruit lifetime {OK_STOMATE}
682	TAU_FRUIT = -9999., 90., 90., 90., 90., 90., 90., 90., 90., -9999., -9999., -9999., -9999.
683
684	# ECUREUIL ([-]) : fraction of primary leaf and root allocation put into reserve {OK_STOMATE}
685	ECUREUIL = -9999., .0, 1., .0, .0, 1., .0, 1., 1., 1., 1., 1., 1.
686
687	# ALLOC_MIN ([-]) : minimum allocation above/below {OK_STOMATE}
688	ALLOC_MIN = -9999., 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, -9999., -9999., -9999., -9999.
689
690	# ALLOC_MAX ([-]) : maximum allocation above/below {OK_STOMATE}
691	ALLOC_MAX = -9999., 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, -9999., -9999., -9999., -9999.
692
693	# DEMI_ALLOC ([-]) : mean allocation above/below {OK_STOMATE}
694	DEMI_ALLOC = -9999., 5., 5., 5., 5., 5., 5., 5., 5., -9999., -9999., -9999., -9999.
695
696	# LEAFLIFE_TAB ([years]) : leaf longevity {OK_STOMATE}
697	LEAFLIFE_TAB = -9999., .5, 2., .33, 1., 2., .33, 2., 2., 2., 2., 2., 2.
698
699	# LEAFFALL ([days]) : length of death of leaves, tabulated {OK_STOMATE}
700	LEAFFALL = -9999., -9999., 10., -9999., -9999., 10., -9999., 10., 10., 10., 10., 10., 10.
701
702	# LEAFAGECRIT ([days]) : critical leaf age, tabulated {OK_STOMATE}
703	LEAFAGECRIT = -9999., 730., 180., 910., 730., 180., 910., 180., 180., 120., 120., 90., 90.
704
705	# SENESCENCE_TYPE ([-]) : type of senescence, tabulated {OK_STOMATE}
706	SENESCENCE_TYPE = none, none, dry, none, none, cold, none, cold, cold, mixed, mixed, mixed, mixed
707
708	# SENESCENCE_HUM ([-] ) : critical relative moisture availability for senescence {OK_STOMATE}
709	SENESCENCE_HUM = -9999., -9999., .3, -9999., -9999., -9999., -9999., -9999., -9999., .2, .2, .3, .2
710
711	# NOSENESCENCE_HUM ([-]) : relative moisture availability above which there is no humidity-related senescence {OK_STOMATE}
712	NOSENESCENCE_HUM = -9999., -9999., .8, -9999., -9999., -9999., -9999., -9999., -9999., .3, .3, .3, .3
713
714	# MAX_TURNOVER_TIME ([days]) : maximum turnover time for grasse {OK_STOMATE}
715	MAX_TURNOVER_TIME = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 80., 80., 80., 80.
716
717	# MIN_TURNOVER_TIME ([days]) : minimum turnover time for grasse {OK_STOMATE}
718	MIN_TURNOVER_TIME = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 10., 10., 10., 10.
719
720	# MIN_LEAF_AGE_FOR_SENESCENCE ([days] ) : minimum leaf age to allow senescence g {OK_STOMATE}
721	MIN_LEAF_AGE_FOR_SENESCENCE = -9999., -9999., 90., -9999., -9999., 90., -9999., 60., 60., 30., 30., 30., 30.
722
723	# SENESCENCE_TEMP_C ([-]) : critical temperature for senescence (C), constant c of aT^2+bT+c, tabulated {OK_STOMATE}
724	SENESCENCE_TEMP_C = -9999., -9999., -9999., -9999., -9999., 12., -9999., 7., 2., -1.375, 5., 5., 10.
725
726	# SENESCENCE_TEMP_B ([-]) : critical temperature for senescence (C), constant b of aT^2+bT+c ,tabulated {OK_STOMATE }
727	SENESCENCE_TEMP_B = -9999., -9999., -9999., -9999., -9999., 0., -9999., 0., 0., .1, 0., 0., 0.
728
729	# SENESCENCE_TEMP_A ([-] ) : critical temperature for senescence (C), constant a of aT^2+bT+c , tabulated {OK_STOMATE}
730	SENESCENCE_TEMP_A = -9999., -9999., -9999., -9999., -9999., 0., -9999., 0., 0.,.00375, 0., 0., 0.
731
732	# GDD_SENESCENCE ([days] ) : minimum gdd to allow senescence of crops {OK_STOMATE}
733	GDD_SENESCENCE = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 950., 4000.
734
735	# ALWAYS_INIT ([BOOLEAN]) : Take carbon from atmosphere if carbohydrate reserve too small {OK_STOMATE}
736	ALWAYS_INIT = y, y, y, y, y, y, y, y, y, y, n, y, y
737
738	# RESIDENCE_TIME ([years]) : residence time of trees {OK_DGVM and NOT(LPJ_GAP_CONST_MORT)}
739	RESIDENCE_TIME = -9999., 30.0, 30.0, 40.0, 40.0, 40.0, 80.0, 80.0, 80.0, 0.0, 0.0, 0.0, 0.0
740
741	# TMIN_CRIT ([C]) : critical tmin, tabulated {OK_STOMATE}
742	TMIN_CRIT = -9999., 0.0, 0.0, -30.0, -14.0, -30.0, -45.0, -45.0, -9999., -9999., -9999., -9999., -9999.
743
744	# TCM_CRIT ([C]) : critical tcm, tabulated {OK_STOMATE}
745	TCM_CRIT = -9999., -9999., -9999., 5.0, 15.5, 15.5, -8.0, -8.0, -8.0, -9999., -9999., -9999., -9999.
746
747	# HERBIVORES ([FLAG]) : herbivores allowed? {OK_STOMATE }
748	HERBIVORES = n
749
750	# TREAT_EXPANSION ([FLAG]) : treat expansion of PFTs across a grid cell? {OK_STOMATE }
751	TREAT_EXPANSION = n
752
753	# LPJ_GAP_CONST_MORT ([FLAG]) : Constant mortality {OK_STOMATE AND NOT OK_DGVM}
754	LPJ_GAP_CONST_MORT = y/n depending on OK_DGVM
755
756	# HARVEST_AGRI ([FLAG]) : Harvest model for agricultural PFTs. {OK_STOMATE }
757	HARVEST_AGRI = y
758
759	# FIRE_DISABLE ([FLAG]) : no fire allowed {OK_STOMATE }
760	FIRE_DISABLE = y
761
762	# SPINUP_ANALYTIC (BOOLEAN ) : Activation of the analytic resolution of the spinup. {OK_STOMATE}
763	SPINUP_ANALYTIC = n
764
765	# AGRICULTURE ([FLAG]) : agriculture allowed? {OK_SECHIBA or OK_STOMATE}
766	AGRICULTURE = y
767
768	# IMPOSE_VEG ([FLAG]) : Should the vegetation be prescribed ? {OK_SECHIBA or OK_STOMATE}
769	IMPOSE_VEG = n
770
771	# IMPOSE_SOILT ([FLAG]) : Should the soil type be prescribed ? {}
772	IMPOSE_SOILT = n
773
774	# LAI_MAP ([FLAG]) : Read the LAI map {OK_SECHIBA or OK_STOMATE}
775	LAI_MAP = n
776
777	# VEGET_REINIT ([FLAG] ) : Reset veget_year counter (obsolet) {VEGET_UPDATE > 0Y}
778	VEGET_REINIT = y
779
780	# VEGETMAP_RESET ([FLAG] ) : Flag to change vegetation map without activating LAND USE change for carbon fluxes and reset carbon related variables to zero {}
781	VEGETMAP_RESET = n
782
783	# VEGET_YEAR ([FLAG] ) : Year of the vegetation map to be read {}
784	VEGET_YEAR = 1
785
786	# MAXMASS_SNOW ([kg/m^2] ) : The maximum mass of a snow {OK_SECHIBA}
787	MAXMASS_SNOW = 3000.
788
789	# SNOWCRI ([kg/m^2] ) : Sets the amount above which only sublimation occures {OK_SECHIBA}
790	SNOWCRI = 1.5
791
792	# MIN_WIND ([m/s]) : Minimum wind speed {OK_SECHIBA}
793	MIN_WIND = 0.1
794
795	# MAX_SNOW_AGE ([days?]) : Maximum period of snow aging {OK_SECHIBA}
796	MAX_SNOW_AGE = 50.
797
798	# SNOW_TRANS ([m] ) : Transformation time constant for snow {OK_SECHIBA}
799	SNOW_TRANS = 0.2
800
801	# OK_NUDGE_MC ([FLAG]) : Activate nudging of soil moisture {}
802	OK_NUDGE_MC = n
803
804	# NUDGE_TAU_MC ([-]) : Relaxation time for nudging of soil moisture expressed in fraction of the day {OK_NUDGE_MC}
805	NUDGE_TAU_MC = 1
806
807	# OK_NUDGE_SNOW ([FLAG]) : Activate nudging of snow variables {}
808	OK_NUDGE_SNOW = n
809
810	# NUDGE_TAU_SNOW ([-]) : Relaxation time for nudging of snow variables {OK_NUDGE_SNOW}
811	NUDGE_TAU_SNOW = 1
812
813	# NUDGE_INTERPOL_WITH_XIOS ([FLAG]) : Activate reading and interpolation with XIOS for nudging fields {OK_NUDGE_MC or OK_NUDGE_SNOW}
814	NUDGE_INTERPOL_WITH_XIOS = n
815
816	# HEIGHT_DISPLACEMENT ([m] ) : Magic number which relates the height to the displacement height. {OK_SECHIBA }
817	HEIGHT_DISPLACEMENT = 0.75
818
819	# Z0_BARE ([m] ) : bare soil roughness length {OK_SECHIBA }
820	Z0_BARE = 0.01
821
822	# Z0_ICE ([m] ) : ice roughness length {OK_SECHIBA }
823	Z0_ICE = 0.001
824
825	# TCST_SNOWA ([days]) : Time constant of the albedo decay of snow {OK_SECHIBA }
826	TCST_SNOWA = 10.0
827
828	# SNOWCRI_ALB ([cm] ) : Critical value for computation of snow albedo {OK_SECHIBA}
829	SNOWCRI_ALB = 10.
830
831	# VIS_DRY ([-] ) : The correspondance table for the soil color numbers and their albedo {OK_SECHIBA }
832	VIS_DRY = 0.24, 0.22, 0.20, 0.18, 0.16, 0.14, 0.12, 0.10, 0.27
833
834	# NIR_DRY ([-] ) : The correspondance table for the soil color numbers and their albedo {OK_SECHIBA }
835	NIR_DRY = 0.48, 0.44, 0.40, 0.36, 0.32, 0.28, 0.24, 0.20, 0.55
836
837	# VIS_WET ([-] ) : The correspondance table for the soil color numbers and their albedo {OK_SECHIBA }
838	VIS_WET = 0.12, 0.11, 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.15
839
840	# NIR_WET ([-] ) : The correspondance table for the soil color numbers and their albedo {OK_SECHIBA }
841	NIR_WET = 0.24, 0.22, 0.20, 0.18, 0.16, 0.14, 0.12, 0.10, 0.31
842
843	# ALBSOIL_VIS ([-] ) : {OK_SECHIBA }
844	ALBSOIL_VIS = 0.18, 0.16, 0.16, 0.15, 0.12, 0.105, 0.09, 0.075, 0.25
845
846	# ALBSOIL_NIR ([-] ) : {OK_SECHIBA }
847	ALBSOIL_NIR = 0.36, 0.34, 0.34, 0.33, 0.30, 0.25, 0.20, 0.15, 0.45
848
849	# ALB_DEADLEAF ([-] ) : albedo of dead leaves, VIS+NIR {OK_SECHIBA }
850	ALB_DEADLEAF = 0.12, 0.35
851
852	# ALB_ICE ([-] ) : albedo of ice, VIS+NIR {OK_SECHIBA}
853	ALB_ICE = 0.60, 0.20
854
855	# CONDVEG_SNOWA ([-]) : The snow albedo used by SECHIBA {OK_SECHIBA}
856	CONDVEG_SNOWA = 1.E+20
857
858	# ALB_BARE_MODEL ([FLAG]) : Switch bare soil albedo dependent (if TRUE) on soil wetness {OK_SECHIBA}
859	ALB_BARE_MODEL = n
860
861	# ALB_BG_MODIS ([FLAG]) : Read bare soil albedo from file with background MODIS data {OK_SECHIBA}
862	ALB_BG_MODIS = y
863
864	# IMPOSE_AZE ([FLAG]) : Should the surface parameters be prescribed {OK_SECHIBA}
865	IMPOSE_AZE = n
866
867	# CONDVEG_Z0 ([m]) : Surface roughness {IMPOSE_AZE}
868	CONDVEG_Z0 = 0.15
869
870	# ROUGHHEIGHT ([m] ) : Height to be added to the height of the first level {IMPOSE_AZE}
871	ROUGHHEIGHT = 0.0
872
873	# CONDVEG_ALBVIS ([-]) : SW visible albedo for the surface {IMPOSE_AZE}
874	CONDVEG_ALBVIS = 0.25
875
876	# CONDVEG_ALBNIR ([-] ) : SW near infrared albedo for the surface {IMPOSE_AZE}
877	CONDVEG_ALBNIR = 0.25
878
879	# CONDVEG_EMIS ([-] ) : Emissivity of the surface for LW radiation {IMPOSE_AZE}
880	CONDVEG_EMIS = 1.0
881
882	# ROUGH_DYN ([FLAG]) : Account for a dynamic roughness height {OK_SECHIBA}
883	ROUGH_DYN = y
884
885	# C1 ([-] ) : Constant used in the formulation of the ratio of {ROUGH_DYN}
886	C1 = 0.32
887
888	# C2 ([-] ) : Constant used in the formulation of the ratio of {ROUGH_DYN}
889	C2 = 0.264
890
891	# C3 ([-] ) : Constant used in the formulation of the ratio of {ROUGH_DYN}
892	C3 = 15.1
893
894	# Cdrag_foliage ([-] ) : Drag coefficient of the foliage {ROUGH_DYN}
895	Cdrag_foliage = 0.2
896
897	# Ct ([-] ) : Heat transfer coefficient of the leaf {ROUGH_DYN}
898	Ct = 0.01
899
900	# Prandtl ([-] ) : Prandtl number used in the calculation of Ct* {ROUGH_DYN}
901	Prandtl = 0.71
902
903	# xansmax ([-] ) : maximum snow albedo {OK_SECHIBA}
904	xansmax = 0.85
905
906	# xansmin ([-] ) : minimum snow albedo {OK_SECHIBA}
907	xansmin = 0.50
908
909	# xans_todry ([S-1] ) : albedo decay rate for the dry snow {OK_SECHIBA}
910	xans_todry = 0.008
911
912	# xans_t ([S-1] ) : albedo decay rate for the wet snow {OK_SECHIBA}
913	xans_t = 0.24
914
915	# xrhosmax ([-] ) : maximum snow density {OK_SECHIBA}
916	xrhosmax = 750
917
918	# xwsnowholdmax1 ([-] ) : snow holding capacity 1 {OK_SECHIBA}
919	xwsnowholdmax1 = 0.03
920
921	# xwsnowholdmax2 ([-] ) : snow holding capacity 2 {OK_SECHIBA}
922	xwsnowholdmax2 = 0.10
923
924	# xsnowrhohold ([kg/m3] ) : snow density {OK_SECHIBA}
925	xsnowrhohold = 200.0
926
927	# ZSNOWTHRMCOND1 ([W/m/K] ) : Thermal conductivity Coef 1 {OK_SECHIBA}
928	ZSNOWTHRMCOND1 = 0.02
929
930	# ZSNOWTHRMCOND2 ([W m5/(kg2 K)] ) : Thermal conductivity Coef 2 {OK_SECHIBA}
931	ZSNOWTHRMCOND2 = 2.5E-6
932
933	# ZSNOWTHRMCOND_AVAP ([W/m/K] ) : Thermal conductivity Coef 1 water vapor {OK_SECHIBA}
934	ZSNOWTHRMCOND_AVAP = -0.06023
935
936	# ZSNOWTHRMCOND_BVAP ([W/m] ) : Thermal conductivity Coef 2 water vapor {OK_SECHIBA}
937	ZSNOWTHRMCOND_BVAP = -2.5425
938
939	# ZSNOWTHRMCOND_CVAP ([K] ) : Thermal conductivity Coef 3 water vapor {OK_SECHIBA}
940	ZSNOWTHRMCOND_CVAP = -289.99
941
942	# ZSNOWCMPCT_RHOD ([kg/m3]) : Snow compaction coefficent {OK_SECHIBA}
943	ZSNOWCMPCT_RHOD = 150.0
944
945	# ZSNOWCMPCT_ACM ([1/s]) : Coefficent for the thermal conductivity {OK_SECHIBA}
946	ZSNOWCMPCT_ACM = 2.8e-6
947
948	# ZSNOWCMPCT_BCM ([1/K]) : Coefficent for the thermal conductivity {OK_SECHIBA}
949	ZSNOWCMPCT_BCM = 0.04
950
951	# ZSNOWCMPCT_CCM ([m3/kg] ) : Coefficent for the thermal conductivity {OK_SECHIBA}
952	ZSNOWCMPCT_CCM = 460.
953
954	# ZSNOWCMPCT_V0 ([Pa/s]) : Vapor coefficent for the thermal conductivity {OK_SECHIBA}
955	ZSNOWCMPCT_V0 = 3.7e7
956
957	# ZSNOWCMPCT_VT ([1/K]) : Vapor coefficent for the thermal conductivity {OK_SECHIBA}
958	ZSNOWCMPCT_VT = 0.081
959
960	# ZSNOWCMPCT_VR ([m3/kg]) : Vapor coefficent for the thermal conductivity {OK_SECHIBA}
961	ZSNOWCMPCT_VR = 0.018
962
963	# CB ([-] ) : Constant of the Louis scheme {OK_SECHIBA}
964	CB = 5.0
965
966	# CC ([-] ) : Constant of the Louis scheme {OK_SECHIBA}
967	CC = 5.0
968
969	# CD ([-] ) : Constant of the Louis scheme {OK_SECHIBA}
970	CD = 5.0
971
972	# RAYT_CSTE ([W.m^{-2}] ) : Constant in the computation of surface resistance {OK_SECHIBA}
973	RAYT_CSTE = 125
974
975	# DEFC_PLUS ([K.W^{-1}] ) : Constant in the computation of surface resistance {OK_SECHIBA}
976	DEFC_PLUS = 23.E-3
977
978	# DEFC_MULT ([K.W^{-1}] ) : Constant in the computation of surface resistance {OK_SECHIBA}
979	DEFC_MULT = 1.5
980
981	# NLAI ([-] ) : Number of LAI levels {OK_SECHIBA}
982	NLAI = 20
983
984	# LAIMAX ([m^2/m^2] ) : Maximum LAI {OK_SECHIBA}
985	LAIMAX =
986
987	# DEW_VEG_POLY_COEFF ([-] ) : coefficients of the polynome of degree 5 for the dew {OK_SECHIBA}
988	DEW_VEG_POLY_COEFF = 0.887773, 0.205673, 0.110112, 0.014843, 0.000824, 0.000017
989
990	# DOWNREGULATION_CO2 ([FLAG] ) : Activation of CO2 downregulation (used for CMIP6 version 6.1.0-6.1.10) {OK_SECHIBA}
991	DOWNREGULATION_CO2 = y
992
993	# DOWNREGULATION_CO2_NEW ([FLAG] ) : Activation of CO2 downregulation (used for CMIP6 version 6.1.11 and later) {OK_SECHIBA}
994	DOWNREGULATION_CO2_NEW = n
995
996	# DOWNREGULATION_CO2_BASELEVEL ([ppm] ) : CO2 base level {DOWNREGULATION_CO2 or DOWNREGULATION_CO2_NEW}
997	DOWNREGULATION_CO2_BASELEVEL = 380.
998
999	# DOWNREGULATION_CO2_MINIMUM ([ppm] ) : CO2 value above which downregulation is taken into account {DOWNREGULATION_CO2_NEW}
1000	DOWNREGULATION_CO2_MINIMUM = 280.
1001
1002	# GB_REF ([s m-1] ) : Leaf bulk boundary layer resistance {}
1003	GB_REF = 1./25.
1004
1005	# MIN_VEGFRAC ([-] ) : Minimal fraction of mesh a vegetation type can occupy {OK_SECHIBA }
1006	MIN_VEGFRAC = 0.001
1007
1008	# STEMPDIAG_BID ([K]) : only needed for an initial LAI if there is no restart file {OK_SECHIBA }
1009	STEMPDIAG_BID = 280.
1010
1011	# LAI_LEVEL_DEPTH ([-] ) : {}
1012	LAI_LEVEL_DEPTH = 0.15
1013
1014	# Oi ([ubar] ) : Intercellular oxygen partial pressure {}
1015	Oi = 210000.
1016
1017	# TOO_LONG ([days] ) : longest sustainable time without regeneration (vernalization) {OK_STOMATE}
1018	TOO_LONG = 5.
1019
1020	# TAU_FIRE ([days] ) : Time scale for memory of the fire index (days). Validated for one year in the DGVM. {OK_STOMATE }
1021	TAU_FIRE = 30.
1022
1023	# LITTER_CRIT ([gC/m^2] ) : Critical litter quantity for fire {OK_STOMATE }
1024	LITTER_CRIT = 200.
1025
1026	# FIRE_RESIST_STRUCT ([-] ) : {OK_STOMATE }
1027	FIRE_RESIST_STRUCT = 0.5
1028
1029	# CO2FRAC ([-] ) : What fraction of a burned plant compartment goes into the atmosphere {OK_STOMATE }
1030	CO2FRAC = 0.95, 0.95, 0., 0.3, 0., 0., 0.95, 0.95
1031
1032	# BCFRAC_COEFF ([-] ) : {OK_STOMATE }
1033	BCFRAC_COEFF = 0.3, 1.3, 88.2
1034
1035	# FIREFRAC_COEFF ([-] ) : {OK_STOMATE }
1036	FIREFRAC_COEFF = 0.45, 0.8, 0.6, 0.13
1037
1038	# REF_GREFF ([1/year] ) : Asymptotic maximum mortality rate {OK_STOMATE }
1039	REF_GREFF = 0.035
1040
1041	# OK_MINRES ([FLAG]) : Do we try to reach a minimum reservoir even if we are severely stressed? {OK_STOMATE }
1042	OK_MINRES = y
1043
1044	# RESERVE_TIME_TREE ([days] ) : maximum time during which reserve is used (trees) {OK_STOMATE }
1045	RESERVE_TIME_TREE = 30.
1046
1047	# RESERVE_TIME_GRASS ([days] ) : maximum time during which reserve is used (grasses) {OK_STOMATE }
1048	RESERVE_TIME_GRASS = 20.
1049
1050	# F_FRUIT ([-] ) : Standard fruit allocation {OK_STOMATE }
1051	F_FRUIT = 0.1
1052
1053	# ALLOC_SAP_ABOVE_GRASS ([-] ) : fraction of sapwood allocation above ground {OK_STOMATE }
1054	ALLOC_SAP_ABOVE_GRASS = 1.0
1055
1056	# MIN_LTOLSR ([-] ) : extrema of leaf allocation fraction {OK_STOMATE }
1057	MIN_LTOLSR = 0.2
1058
1059	# MAX_LTOLSR ([-] ) : extrema of leaf allocation fraction {OK_STOMATE }
1060	MAX_LTOLSR = 0.5
1061
1062	# Z_NITROGEN ([m] ) : scaling depth for nitrogen limitation {OK_STOMATE}
1063	Z_NITROGEN = 0.2
1064
1065	# NLIM_TREF ([C] ) : {OK_STOMATE }
1066	NLIM_TREF = 25.
1067
1068	# PIPE_TUNE1 ([-] ) : crown area {OK_STOMATE }
1069	PIPE_TUNE1 = 100.0
1070
1071	# PIPE_TUNE2 ([-] ) : height {OK_STOMATE }
1072	PIPE_TUNE2 = 40.0
1073
1074	# PIPE_TUNE3 ([-] ) : height {OK_STOMATE }
1075	PIPE_TUNE3 = 0.5
1076
1077	# PIPE_TUNE4 ([-] ) : needed for stem diameter {OK_STOMATE }
1078	PIPE_TUNE4 = 0.3
1079
1080	# PIPE_DENSITY ([-] ) : Density {OK_STOMATE }
1081	PIPE_DENSITY = 2.e5
1082
1083	# PIPE_K1 ([-] ) : {OK_STOMATE }
1084	PIPE_K1 = 8.e3
1085
1086	# PIPE_TUNE_EXP_COEFF ([-] ) : pipe tune exponential coeff {OK_STOMATE }
1087	PIPE_TUNE_EXP_COEFF = 1.6
1088
1089	# PRECIP_CRIT ([mm/year] ) : minimum precip {OK_STOMATE }
1090	PRECIP_CRIT = 100.
1091
1092	# GDD_CRIT_ESTAB ([-] ) : minimum gdd for establishment of saplings {OK_STOMATE }
1093	GDD_CRIT_ESTAB = 150.
1094
1095	# FPC_CRIT ([-] ) : critical fpc, needed for light competition and establishment {OK_STOMATE }
1096	FPC_CRIT = 0.95
1097
1098	# ALPHA_GRASS ([-] ) : sapling characteristics : alpha's {OK_STOMATE }
1099	ALPHA_GRASS = 0.5
1100
1101	# ALPHA_TREE ([-] ) : sapling characteristics : alpha's {OK_STOMATE }
1102	ALPHA_TREE = 1.
1103
1104	# MASS_RATIO_HEART_SAP ([-] ) : mass ratio (heartwood+sapwood)/sapwood {OK_STOMATE }
1105	MASS_RATIO_HEART_SAP = 3.
1106
1107	# TAU_HUM_MONTH ([days] ) : time scales for phenology and other processes {OK_STOMATE }
1108	TAU_HUM_MONTH = 20.
1109
1110	# TAU_HUM_WEEK ([days] ) : time scales for phenology and other processes {OK_STOMATE }
1111	TAU_HUM_WEEK = 7.
1112
1113	# TAU_T2M_MONTH ([days] ) : time scales for phenology and other processes {OK_STOMATE }
1114	TAU_T2M_MONTH = 20.
1115
1116	# TAU_T2M_WEEK ([days] ) : time scales for phenology and other processes {OK_STOMATE }
1117	TAU_T2M_WEEK = 7.
1118
1119	# TAU_TSOIL_MONTH ([days] ) : time scales for phenology and other processes {OK_STOMATE }
1120	TAU_TSOIL_MONTH = 20.
1121
1122	# TAU_SOILHUM_MONTH ([days] ) : time scales for phenology and other processes {OK_STOMATE }
1123	TAU_SOILHUM_MONTH = 20.
1124
1125	# TAU_GPP_WEEK ([days] ) : time scales for phenology and other processes {OK_STOMATE }
1126	TAU_GPP_WEEK = 7.
1127
1128	# TAU_GDD ([days] ) : time scales for phenology and other processes {OK_STOMATE }
1129	TAU_GDD = 40.
1130
1131	# TAU_NGD ([days] ) : time scales for phenology and other processes {OK_STOMATE }
1132	TAU_NGD = 50.
1133
1134	# COEFF_TAU_LONGTERM ([days] ) : time scales for phenology and other processes {OK_STOMATE }
1135	COEFF_TAU_LONGTERM = 3.
1136
1137	# BM_SAPL_CARBRES ([-] ) : {OK_STOMATE }
1138	BM_SAPL_CARBRES = 5.
1139
1140	# BM_SAPL_SAPABOVE ([-] ) : {OK_STOMATE}
1141	BM_SAPL_SAPABOVE = 0.5
1142
1143	# BM_SAPL_HEARTABOVE ([-] ) : {OK_STOMATE }
1144	BM_SAPL_HEARTABOVE = 2.
1145
1146	# BM_SAPL_HEARTBELOW ([-] ) : {OK_STOMATE }
1147	BM_SAPL_HEARTBELOW = 2.
1148
1149	# INIT_SAPL_MASS_LEAF_NAT ([-] ) : {OK_STOMATE }
1150	INIT_SAPL_MASS_LEAF_NAT = 0.1
1151
1152	# INIT_SAPL_MASS_LEAF_AGRI ([-] ) : {OK_STOMATE }
1153	INIT_SAPL_MASS_LEAF_AGRI = 1.
1154
1155	# INIT_SAPL_MASS_CARBRES ([-] ) : {OK_STOMATE }
1156	INIT_SAPL_MASS_CARBRES = 5.
1157
1158	# INIT_SAPL_MASS_ROOT ([-] ) : {OK_STOMATE }
1159	INIT_SAPL_MASS_ROOT = 0.1
1160
1161	# INIT_SAPL_MASS_FRUIT ([-] ) : {OK_STOMATE }
1162	INIT_SAPL_MASS_FRUIT = 0.3
1163
1164	# CN_SAPL_INIT ([-] ) : {OK_STOMATE }
1165	CN_SAPL_INIT = 0.5
1166
1167	# MIGRATE_TREE ([m/year] ) : {OK_STOMATE }
1168	MIGRATE_TREE = 10000.
1169
1170	# MIGRATE_GRASS ([m/year] ) : {OK_STOMATE }
1171	MIGRATE_GRASS = 10000.
1172
1173	# LAI_INITMIN_TREE ([m^2/m^2] ) : {OK_STOMATE }
1174	LAI_INITMIN_TREE = 0.3
1175
1176	# LAI_INITMIN_GRASS ([m^2/m^2] ) : {OK_STOMATE }
1177	LAI_INITMIN_GRASS = 0.1
1178
1179	# DIA_COEFF ([-] ) : {OK_STOMATE }
1180	DIA_COEFF = 4., 0.5
1181
1182	# MAXDIA_COEFF ([-] ) : {OK_STOMATE }
1183	MAXDIA_COEFF = 100., 0.01
1184
1185	# BM_SAPL_LEAF ([-] ) : {OK_STOMATE }
1186	BM_SAPL_LEAF = 4., 4., 0.8, 5.
1187
1188	# METABOLIC_REF_FRAC ([-]) : {OK_STOMATE }
1189	METABOLIC_REF_FRAC = 0.85
1190
1191	# Z_DECOMP ([m] ) : scaling depth for soil activity {OK_STOMATE }
1192	Z_DECOMP = 0.2
1193
1194	# CN ([-] ) : C/N ratio {OK_STOMATE }
1195	CN = 40., 40., 40., 40., 40., 40., 40., 40.
1196
1197	# LC ([-] ) : Lignine/C ratio of the different plant parts {OK_STOMATE }
1198	LC = 0.22, 0.35, 0.35, 0.35, 0.35, 0.22, 0.22, 0.22
1199
1200	# FRAC_SOIL_STRUCT_AA ([-]) : frac_soil(istructural,iactive,iabove) {OK_STOMATE }
1201	FRAC_SOIL_STRUCT_AA = 0.55
1202
1203	# FRAC_SOIL_STRUCT_A ([-]) : frac_soil(istructural,iactive,ibelow) {OK_STOMATE }
1204	FRAC_SOIL_STRUCT_A = 0.45
1205
1206	# FRAC_SOIL_STRUCT_SA ([-] ) : frac_soil(istructural,islow,iabove) {OK_STOMATE}
1207	FRAC_SOIL_STRUCT_SA = 0.7
1208
1209	# FRAC_SOIL_STRUCT_SB ([-] ) : frac_soil(istructural,islow,ibelow) {OK_STOMATE }
1210	FRAC_SOIL_STRUCT_SB = 0.7
1211
1212	# FRAC_SOIL_METAB_AA ([-] ) : frac_soil(imetabolic,iactive,iabove) {OK_STOMATE }
1213	FRAC_SOIL_METAB_AA = 0.45
1214
1215	# FRAC_SOIL_METAB_AB ([-] ) : frac_soil(imetabolic,iactive,ibelow) {OK_STOMATE }
1216	FRAC_SOIL_METAB_AB = 0.45
1217
1218	# METABOLIC_LN_RATIO ([-] ) : {OK_STOMATE }
1219	METABOLIC_LN_RATIO = 0.018
1220
1221	# TAU_METABOLIC ([days] ) : {OK_STOMATE }
1222	TAU_METABOLIC = 0.066
1223
1224	# TAU_STRUCT ([days]) : {OK_STOMATE }
1225	TAU_STRUCT = 0.245
1226
1227	# SOIL_Q10 ([-]) : {OK_STOMATE }
1228	SOIL_Q10 = 0.69 (
1229
1230	# TSOIL_REF ([C] ) : {OK_STOMATE }
1231	TSOIL_REF = 30.
1232
1233	# LITTER_STRUCT_COEF ([-] ) : {OK_STOMATE }
1234	LITTER_STRUCT_COEF = 3.
1235
1236	# MOIST_COEFF ([-] ) : {OK_STOMATE }
1237	MOIST_COEFF = 1.1, 2.4, 0.29
1238
1239	# MOISTCONT_MIN ([-]) : minimum soil wetness to limit the heterotrophic respiration {OK_STOMATE }
1240	MOISTCONT_MIN = 0.25
1241
1242	# FRAC_TURNOVER_DAILY ([-]) : {OK_STOMATE }
1243	FRAC_TURNOVER_DAILY = 0.55
1244
1245	# TAX_MAX ([-] ) : maximum fraction of allocatable biomass used for maintenance respiration {OK_STOMATE }
1246	TAX_MAX = 0.8
1247
1248	# MIN_GROWTHINIT_TIME ([days] ) : minimum time since last beginning of a growing season {OK_STOMATE }
1249	MIN_GROWTHINIT_TIME = 300.
1250
1251	# MOIAVAIL_ALWAYS_TREE ([-] ) : moisture availability above which moisture tendency doesn't matter {OK_STOMATE }
1252	MOIAVAIL_ALWAYS_TREE = 1.0
1253
1254	# MOIAVAIL_ALWAYS_GRASS ([-] ) : moisture availability above which moisture tendency doesn't matter {OK_STOMATE }
1255	MOIAVAIL_ALWAYS_GRASS = 0.6
1256
1257	# T_ALWAYS_ADD ([C] ) : monthly temp. above which temp. tendency doesn't matter {OK_STOMATE }
1258	T_ALWAYS_ADD = 10.
1259
1260	# GDDNCD_REF ([-] ) : {OK_STOMATE }
1261	GDDNCD_REF = 603.
1262
1263	# GDDNCD_CURVE ([-] ) : {OK_STOMATE }
1264	GDDNCD_CURVE = 0.0091
1265
1266	# GDDNCD_OFFSET ([-] ) : {OK_STOMATE }
1267	GDDNCD_OFFSET = 64.
1268
1269	# BM_SAPL_RESCALE ([-] ) : {OK_STOMATE }
1270	BM_SAPL_RESCALE = 40.
1271
1272	# MAINT_RESP_MIN_VMAX ([-] ) : {OK_STOMATE }
1273	MAINT_RESP_MIN_VMAX = 0.3
1274
1275	# MAINT_RESP_COEFF ([-] ) : {OK_STOMATE }
1276	MAINT_RESP_COEFF = 1.4
1277
1278	# FRAC_CARB_AP ([-]) : frac carb coefficients from active pool: depends on clay content {OK_STOMATE }
1279	FRAC_CARB_AP = 0.004
1280
1281	# FRAC_CARB_SA ([-]) : frac_carb_coefficients from slow pool {OK_STOMATE }
1282	FRAC_CARB_SA = 0.42
1283
1284	# FRAC_CARB_SP ([-] ) : frac_carb_coefficients from slow pool {OK_STOMATE }
1285	FRAC_CARB_SP = 0.03
1286
1287	# FRAC_CARB_PA ([-]) : frac_carb_coefficients from passive pool {OK_STOMATE }
1288	FRAC_CARB_PA = 0.45
1289
1290	# FRAC_CARB_PS ([-]) : frac_carb_coefficients from passive pool {OK_STOMATE }
1291	FRAC_CARB_PS = 0.0
1292
1293	# ACTIVE_TO_PASS_CLAY_FRAC ([-] ) : {OK_STOMATE }
1294	ACTIVE_TO_PASS_CLAY_FRAC = 0.68
1295
1296	# CARBON_TAU_IACTIVE ( [days] ) : residence times in carbon pools {OK_STOMATE }
1297	CARBON_TAU_IACTIVE = 0.149
1298
1299	# CARBON_TAU_ISLOW ([days]) : residence times in carbon pools {OK_STOMATE }
1300	CARBON_TAU_ISLOW = 7.0
1301
1302	# CARBON_TAU_IPASSIVE ([days] ) : residence times in carbon pools {OK_STOMATE }
1303	CARBON_TAU_IPASSIVE = 300.
1304
1305	# FLUX_TOT_COEFF ([days] ) : {OK_STOMATE }
1306	FLUX_TOT_COEFF = 1.2, 1.4,.75
1307
1308	# NEW_TURNOVER_TIME_REF ([days] ) : {OK_STOMATE }
1309	NEW_TURNOVER_TIME_REF = 20.
1310
1311	# LEAF_AGE_CRIT_TREF ([days] ) : {OK_STOMATE }
1312	LEAF_AGE_CRIT_TREF = 20.
1313
1314	# LEAF_AGE_CRIT_COEFF ([-] ) : {OK_STOMATE }
1315	LEAF_AGE_CRIT_COEFF = 1.5, 0.75, 10.
1316
1317	# VMAX_OFFSET ([-] ) : offset (minimum relative vcmax) {OK_STOMATE }
1318	VMAX_OFFSET = 0.3
1319
1320	# LEAFAGE_FIRSTMAX ([-] ) : leaf age at which vmax attains vcmax_opt (in fraction of critical leaf age) {OK_STOMATE }
1321	LEAFAGE_FIRSTMAX = 0.03
1322
1323	# LEAFAGE_LASTMAX ([-] ) : leaf age at which vmax falls below vcmax_opt (in fraction of critical leaf age) {OK_STOMATE }
1324	LEAFAGE_LASTMAX = 0.5
1325
1326	# LEAFAGE_OLD ([-] ) : leaf age at which vmax attains its minimum (in fraction of critical leaf age) {OK_STOMATE }
1327	LEAFAGE_OLD = 1.
1328
1329	# GPPFRAC_DORMANCE ([-]) : rapport maximal GPP/GGP_max pour dormance {OK_STOMATE }
1330	GPPFRAC_DORMANCE = 0.2
1331
1332	# TAU_CLIMATOLOGY ([days]) : tau for "climatologic variables {OK_STOMATE }
1333	TAU_CLIMATOLOGY = 20
1334
1335	# HVC1 ([-] ) : parameters for herbivore activity {OK_STOMATE }
1336	HVC1 = 0.019
1337
1338	# HVC2 ([-] ) : parameters for herbivore activity {OK_STOMATE }
1339	HVC2 = 1.38
1340
1341	# LEAF_FRAC_HVC ([-] ) : parameters for herbivore activity {OK_STOMATE }
1342	LEAF_FRAC_HVC = 0.33
1343
1344	# TLONG_REF_MAX ([K] ) : maximum reference long term temperature {OK_STOMATE }
1345	TLONG_REF_MAX = 303.1
1346
1347	# TLONG_REF_MIN ([K] ) : minimum reference long term temperature {OK_STOMATE }
1348	TLONG_REF_MIN = 253.1
1349
1350	# NCD_MAX_YEAR ([days]) : {OK_STOMATE }
1351	NCD_MAX_YEAR = 3.
1352
1353	# GDD_THRESHOLD ([days] ) : {OK_STOMATE }
1354	GDD_THRESHOLD = 5.
1355
1356	# GREEN_AGE_EVER ([-] ) : {OK_STOMATE }
1357	GREEN_AGE_EVER = 2.
1358
1359	# GREEN_AGE_DEC ([-] ) : {OK_STOMATE }
1360	GREEN_AGE_DEC = 0.5
1361
1362	# ESTAB_MAX_TREE ([-] ) : Maximum tree establishment rate {OK_DGVM}
1363	ESTAB_MAX_TREE = 0.12
1364
1365	# ESTAB_MAX_GRASS ([-] ) : Maximum grass establishment rate {OK_DGVM}
1366	ESTAB_MAX_GRASS = 0.12
1367
1368	# ESTABLISH_SCAL_FACT ([-] ) : {OK_DGVM }
1369	ESTABLISH_SCAL_FACT = 5.
1370
1371	# MAX_TREE_COVERAGE ([-] ) : {OK_DGVM }
1372	MAX_TREE_COVERAGE = 0.98
1373
1374	# IND_0_ESTAB ([-] ) : {OK_DGVM }
1375	IND_0_ESTAB = 0.2
1376
1377	# ANNUAL_INCREASE ([FLAG]) : for diagnosis of fpc increase, compare today's fpc to last year's maximum (T) or to fpc of last time step (F)? {OK_DGVM}
1378	ANNUAL_INCREASE = y
1379
1380	# MIN_COVER ([-] ) : For trees, minimum fraction of crown area occupied {OK_DGVM}
1381	MIN_COVER = 0.05
1382
1383	# IND_0 ([-] ) : initial density of individuals {OK_DGVM}
1384	IND_0 = 0.02
1385
1386	# MIN_AVAIL ([-] ) : minimum availability {OK_DGVM}
1387	MIN_AVAIL = 0.01
1388
1389	# RIP_TIME_MIN ([year] ) : {OK_DGVM}
1390	RIP_TIME_MIN = 1.25
1391
1392	# NPP_LONGTERM_INIT ([gC/m^2/year]) : {OK_DGVM}
1393	NPP_LONGTERM_INIT = 10.
1394
1395	# EVERYWHERE_INIT ([-] ) : {OK_DGVM}
1396	EVERYWHERE_INIT = 0.05
1397
1398	# PRINTLEV ([0, 1, 2, 3, 4]) : Print level for text output {}
1399	PRINTLEV = 2
1400
1401	# PRINTLEV_modname ([0, 1, 2, 3, 4]) : Specific print level of text output for the module "modname". Default as PRINTLEV. {}
1402	PRINTLEV_modname = PRINTLEV
1403
1404	# DRY_SOIL_HEAT_CAPACITY ([J.m^{-3}.K^{-1}] ) : Dry soil Heat capacity of soils {OK_SECHIBA }
1405	DRY_SOIL_HEAT_CAPACITY = 1.80e+6
1406
1407	# DRY_SOIL_HEAT_COND ([W.m^{-2}.K^{-1}] ) : Dry soil Thermal Conductivity of soils {OK_SECHIBA}
1408	DRY_SOIL_HEAT_COND = 0.40
1409
1410	# SNOW_HEAT_COND ([W.m^{-2}.K^{-1}]) : Thermal Conductivity of snow {OK_SECHIBA }
1411	SNOW_HEAT_COND = 0.3
1412
1413	# SNOW_DENSITY ([-] ) : Snow density for the soil thermodynamics {OK_SECHIBA }
1414	SNOW_DENSITY = 330.0
1415
1416	# NOBIO_WATER_CAPAC_VOLUMETRI ([s/m^2]) : {}
1417	NOBIO_WATER_CAPAC_VOLUMETRI = 150.
1418
1419	# SECHIBA_QSINT ([m]) : Interception reservoir coefficient {OK_SECHIBA }
1420	SECHIBA_QSINT = 0.02
1421
1422	# OK_FREEZE ([FLAG]) : Activate the complet soil freezing scheme {OK_SECHIBA }
1423	OK_FREEZE = TRUE
1424
1425	# READ_REFTEMP ([FLAG]) : Initialize soil temperature using climatological temperature {}
1426	READ_REFTEMP = True/False depening on OK_FREEZE
1427
1428	# OK_FREEZE_THERMIX ([FLAG]) : Activate thermal part of the soil freezing scheme {}
1429	OK_FREEZE_THERMIX = True if OK_FREEZE else false
1430
1431	# OK_ECORR ([FLAG]) : Energy correction for freezing {OK_FREEZE_THERMIX}
1432	OK_ECORR = True if OK_FREEZE else false
1433
1434	# OK_FREEZE_THAW_LATENT_HEAT ([FLAG]) : Activate latent heat part of the soil freezing scheme {}
1435	OK_FREEZE_THAW_LATENT_HEAT = FALSE
1436
1437	# fr_dT ([K] ) : Freezing window {OK_SECHIBA}
1438	fr_dT = 2.0
1439
1440	# OK_FREEZE_CWRR ([FLAG]) : CWRR freezing scheme by I. Gouttevin {}
1441	OK_FREEZE_CWRR = True if OK_FREEZE else false
1442
1443	# OK_THERMODYNAMICAL_FREEZING ([FLAG]) : Calculate frozen fraction thermodynamically {OK_FREEZE_CWRR}
1444	OK_THERMODYNAMICAL_FREEZING = True
1445
1446	# CHECK_CWRR ([FLAG]) : Calculate diagnostics to check CWRR water balance {}
1447	CHECK_CWRR = n
1448
1449	# SOIL_FRACTIONS ([-]) : Areal fraction of the 13 soil USDA textures; the dominant one is selected, Loam by default {IMPOSE_SOILT}
1450	SOIL_FRACTIONS = 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0
1451
1452	# CLAY_FRACTION ([-]) : Fraction of the clay fraction (0-dim mode) {IMPOSE_SOIL}
1453	CLAY_FRACTION = 0.2 if Loam
1454
1455	# SAND_FRACTION ([-] ) : Fraction of the clay fraction (0-dim mode) {IMPOSE_SOIL}
1456	SAND_FRACTION = 0.4
1457
1458	# nvan ([-]) : nvan parameter from Van genutchen equations {IMPOSE_SOIL}
1459	nvan = 1.56 if Loam
1460
1461	# avan ([-]) : avan parameter from Van genutchen equations {IMPOSE_SOIL}
1462	avan = 0.0036 if Loam
1463
1464	# mcr ([-]) : residual soil moisture {IMPOSE_SOIL}
1465	mcr = 0.078 if Loam
1466
1467	# mcs ([-]) : saturation soil moisture {IMPOSE_SOIL}
1468	mcs = 0.43 if Loam
1469
1470	# ks ([mm/d]) : saturation conductivity {IMPOSE_SOIL}
1471	ks = 249.6 if Loam
1472
1473	# mcfc ([-]) : field capacity soil moisture {IMPOSE_SOIL}
1474	mcfc = 0.1654 if caclulated from default 5 parameters above
1475
1476	# mcw ([-]) : wilting point soil moisture {IMPOSE_SOIL}
1477	mcw = 0.0884 if caclulated from default 5 parameters above
1478
1479	# GET_SLOPE ([FLAG]) : Read slopes from file and do the interpolation {}
1480	GET_SLOPE = n
1481
1482	# REINF_SLOPE ([-]) : Fraction of reinfiltrated surface runoff {}
1483	REINF_SLOPE = 0.1
1484
1485	# SECHIBA_QSINT ([m]) : Interception reservoir coefficient {OK_SECHIBA }
1486	SECHIBA_QSINT = 0.02
1487
1488	# SECHIBA_ZCANOP ([m]) : Soil level used for canopy development (if STOMATE disactivated) {OK_SECHIBA and .NOT. OK_STOMATE }
1489	SECHIBA_ZCANOP = 0.5
1490
1491	# VEGET_UPDATE ([years]) : Update vegetation frequency: 0Y or 1Y {}
1492	VEGET_UPDATE = 0Y
1493
1494	# SECHIBA_VEGMAX ([-]) : Maximum vegetation distribution within the mesh (0-dim mode) {IMPOSE_VEG}
1495	SECHIBA_VEGMAX = 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8, 0.0, 0.0, 0.0
1496
1497	# SECHIBA_FRAC_NOBIO ([-]) : Fraction of other surface types within the mesh (0-dim mode) {IMPOSE_VEG}
1498	SECHIBA_FRAC_NOBIO = 0.0
1499
1500	# SECHIBA_LAI ([-]) : LAI for all vegetation types (0-dim mode) {IMPOSE_VEG}
1501	SECHIBA_LAI = 0., 8., 8., 4., 4.5, 4.5, 4., 4.5, 4., 2., 2., 2., 2.
1502
1503	# SLOWPROC_HEIGHT ([m]) : Height for all vegetation types {OK_SECHIBA}
1504	SLOWPROC_HEIGHT = 0., 30., 30., 20., 20., 20., 15., 15., 15., .5, .6, 1.0, 1.0
1505
1506	# LAI_FILE ([FILE]) : Name of file from which the vegetation map is to be read {LAI_MAP}
1507	LAI_FILE = lai2D.nc
1508
1509	# RENORM_LAI ([FLAG]) : flag to force LAI renormelization {LAI_MAP}
1510	RENORM_LAI = n
1511
1512	# VEGETATION_FILE ([FILE]) : Name of file from which the vegetation map is to be read {}
1513	VEGETATION_FILE = PFTmap.nc
1514
1515	# SPMIPEXP ([-]) : Types of alternative hydraulic parameters {}
1516	SPMIPEXP = 'texture'
1517
1518	# UNIF_CASE ([-]) : Types of uniform soil textures in SPMIP {SPMIPEXP='unif'}
1519	UNIF_CASE = 'b'
1520
1521	# SOILCLASS_FILE ([FILE]) : Name of file from which soil types are read {NOT(IMPOSE_VEG)}
1522	SOILCLASS_FILE = soils_param.nc
1523
1524	# PARAM_SPMIP_FILE ([FILE]) : Name of file from which soil parameter values are read {smipexp='maps'}
1525	PARAM_SPMIP_FILE = params_sp_mip.nc
1526
1527	# SLOPE_NOREINF ([%]) : Slope over which surface runoff does not reinfiltrate {}
1528	SLOPE_NOREINF = 0.5
1529
1530	# TOPOGRAPHY_FILE ([FILE]) : Name of file from which the topography map is to be read {}
1531	TOPOGRAPHY_FILE = cartepente2d_15min.nc
1532
1533	# WOODHARVEST_FILE ([FILE]) : Name of file from which the wood harvest will be read {DO_WOOD_HARVEST}
1534	WOODHARVEST_FILE = woodharvest.nc
1535
1536	# SOILALB_FILE ([FILE]) : Name of file from which the bare soil albedo {NOT(IMPOSE_AZE)}
1537	SOILALB_FILE = soils_param.nc
1538
1539	# ALB_BG_FILE ([FILE]) : Name of file from which the background albedo is read {ALB_BG_MODIS}
1540	ALB_BG_FILE = alb_bg.nc
1541
1542	# CDRAG_FROM_GCM ([FLAG]) : Keep cdrag coefficient from gcm. {OK_SECHIBA}
1543	CDRAG_FROM_GCM = y
1544
1545	# N_FERTIL_FILE (- ) : File name {CHEMISTRY_BVOC and NOx_FERTILIZERS_USE}
1546	N_FERTIL_FILE = orchidee_fertilizer_1995.nc
1547
1548	# N_FERTIL_FILE (-) : File name {CHEMISTRY_BVOC and NOx_FERTILIZERS_USE}
1549	N_FERTIL_FILE = orchidee_fertilizer_1995.nc
1550
1551	# ENERBIL_TSURF (Kelvin [K]) : Initial temperature if not found in restart {OK_SECHIBA}
1552	ENERBIL_TSURF = 280.
1553
1554	# ENERBIL_EVAPOT () : Initial Soil Potential Evaporation {OK_SECHIBA }
1555	ENERBIL_EVAPOT = 0.0
1556
1557	# BEDROCK_FLAG ([FLAG]) : Flag to consider bedrock at deeper layers. {}
1558	BEDROCK_FLAG = 0
1559
1560	# THERMOSOIL_TPRO (Kelvin [K]) : Initial soil temperature profile if not found in restart {OK_SECHIBA}
1561	THERMOSOIL_TPRO = 280.
1562
1563	# SOIL_REFTEMP_FILE ([FILE]) : File with climatological soil temperature {READ_REFTEMP}
1564	SOIL_REFTEMP_FILE = reftemp.nc
1565
1566	# DO_PONDS ([FLAG]) : Should we include ponds {}
1567	DO_PONDS = n
1568
1569	# FROZ_FRAC_CORR ([-]) : Coefficient for the frozen fraction correction {OK_FREEZE}
1570	FROZ_FRAC_CORR = 1.0
1571
1572	# MAX_FROZ_HYDRO ([-]) : Coefficient for the frozen fraction correction {OK_FREEZE}
1573	MAX_FROZ_HYDRO = 1.0
1574
1575	# SMTOT_CORR ([-]) : Coefficient for the frozen fraction correction {OK_FREEZE}
1576	SMTOT_CORR = 2.0
1577
1578	# DO_RSOIL ([FLAG]) : Should we reduce soil evaporation with a soil resistance {}
1579	DO_RSOIL = n
1580
1581	# OK_DYNROOT ([FLAG]) : Calculate dynamic root profile to optimize soil moisture usage {}
1582	OK_DYNROOT = n
1583
1584	# WETNESS_TRANSPIR_MAX ([-] ) : Soil moisture above which transpir is max, for each soil texture class {}
1585	WETNESS_TRANSPIR_MAX = 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8
1586
1587	# VWC_MIN_FOR_WET_ALB ([m3/m3] ) : Vol. wat. cont. above which albedo is cst {}
1588	VWC_MIN_FOR_WET_ALB = 0.25, 0.25, 0.25
1589
1590	# VWC_MAX_FOR_DRY_ALB ([m3/m3] ) : Vol. wat. cont. below which albedo is cst {}
1591	VWC_MAX_FOR_DRY_ALB = 0.1, 0.1, 0.1
1592
1593	# HYDROL_MOISTURE_CONTENT ([m3/m3]) : Soil moisture on each soil tile and levels {}
1594	HYDROL_MOISTURE_CONTENT = 0.3
1595
1596	# US_INIT ([-]) : US_NVM_NSTM_NSLM {}
1597	US_INIT = 0.0
1598
1599	# ZWT_FORCE ([m]) : Prescribed water depth, dimension nstm {}
1600	ZWT_FORCE = -9999. -9999. -9999.
1601
1602	# FREE_DRAIN_COEF ([-]) : Coefficient for free drainage at bottom, dimension nstm {}
1603	FREE_DRAIN_COEF = 1.0 1.0 1.0
1604
1605	# WATER_TO_INFILT ([mm]) : Water to be infiltrated on top of the soil {}
1606	WATER_TO_INFILT = 0.0
1607
1608	# EVAPNU_SOIL ([mm]) : Bare soil evap on each soil if not found in restart {}
1609	EVAPNU_SOIL = 0.0
1610
1611	# HYDROL_SNOW () : Initial snow mass if not found in restart {OK_SECHIBA}
1612	HYDROL_SNOW = 0.0
1613
1614	# HYDROL_SNOWAGE (***) : Initial snow age if not found in restart {OK_SECHIBA}
1615	HYDROL_SNOWAGE = 0.0
1616
1617	# HYDROL_SNOW_NOBIO ([mm]) : Initial snow amount on ice, lakes, etc. if not found in restart {OK_SECHIBA}
1618	HYDROL_SNOW_NOBIO = 0.0
1619
1620	# HYDROL_SNOW_NOBIO_AGE (***) : Initial snow age on ice, lakes, etc. if not found in restart {OK_SECHIBA}
1621	HYDROL_SNOW_NOBIO_AGE = 0.0
1622
1623	# HYDROL_QSV ([mm]) : Initial water on canopy if not found in restart {OK_SECHIBA}
1624	HYDROL_QSV = 0.0
1625
1626	# CWRR_NKS_N0 ([-]) : fitted value for relation log((n-n0)/(n_ref-n0)) {}
1627	CWRR_NKS_N0 = 0.0
1628
1629	# CWRR_NKS_POWER ([-]) : fitted value for relation log((n-n0)/(n_ref-n0)) {}
1630	CWRR_NKS_POWER = 0.0
1631
1632	# CWRR_AKS_A0 ([1/mm]) : fitted value for relation log((a-a0)/(a_ref-a0)) {}
1633	CWRR_AKS_A0 = 0.0
1634
1635	# CWRR_AKS_POWER ([-]) : fitted value for relation log((a-a0)/(a_ref-a0)) {}
1636	CWRR_AKS_POWER = 0.0
1637
1638	# KFACT_DECAY_RATE ([1/m]) : Factor for Ks decay with depth {}
1639	KFACT_DECAY_RATE = 2.0
1640
1641	# KFACT_STARTING_DEPTH ([m]) : Depth for compacted value of Ks {}
1642	KFACT_STARTING_DEPTH = 0.3
1643
1644	# KFACT_MAX ([-]) : Maximum Factor for Ks increase due to vegetation {}
1645	KFACT_MAX = 10.0
1646
1647	# KFACT_ROOT_CONST ([y/n]) : Set constant kfact_root in every soil layer. Otherwise kfact_root increase over soil depth in the rootzone. {}
1648	KFACT_ROOT_CONST = n
1649
1650	# DT_ROUTING ([seconds]) : Time step of the routing scheme {RIVER_ROUTING}
1651	DT_ROUTING = 86400.
1652
1653	# ROUTING_RIVERS ([-]) : Number of rivers {RIVER_ROUTING}
1654	ROUTING_RIVERS = 50
1655
1656	# DO_FLOODINFILT ([FLAG]) : Should floodplains reinfiltrate into the soil {RIVER_ROUTING}
1657	DO_FLOODINFILT = n
1658
1659	# DO_SWAMPS ([FLAG]) : Should we include swamp parameterization {RIVER_ROUTING}
1660	DO_SWAMPS = n
1661
1662	# DO_PONDS ([FLAG]) : Should we include ponds {RIVER_ROUTING}
1663	DO_PONDS = n
1664
1665	# SLOW_TCST ([days]) : Time constant for the slow reservoir {RIVER_ROUTING }
1666	SLOW_TCST = 25.0
1667
1668	# FAST_TCST ([days]) : Time constant for the fast reservoir {RIVER_ROUTING }
1669	FAST_TCST = 3.0
1670
1671	# STREAM_TCST ([days]) : Time constant for the stream reservoir {RIVER_ROUTING}
1672	STREAM_TCST = 0.24
1673
1674	# FLOOD_TCST ([days]) : Time constant for the flood reservoir {RIVER_ROUTING}
1675	FLOOD_TCST = 4.0
1676
1677	# SWAMP_CST ([-]) : Fraction of the river that flows back to swamps {RIVER_ROUTING}
1678	SWAMP_CST = 0.2
1679
1680	# FLOOD_BETA ([-] ) : Parameter to fix the shape of the floodplain {RIVER_ROUTING}
1681	FLOOD_BETA = 2.0
1682
1683	# POND_BETAP ([-] ) : Ratio of the basin surface intercepted by ponds and the maximum surface of ponds {RIVER_ROUTING}
1684	POND_BETAP = 0.5
1685
1686	# FLOOD_CRI ([mm] ) : Potential height for which all the basin is flooded {DO_FLOODPLAINS or DO_PONDS}
1687	FLOOD_CRI = 2000.
1688
1689	# POND_CRI ([mm] ) : Potential height for which all the basin is a pond {DO_FLOODPLAINS or DO_PONDS}
1690	POND_CRI = 2000.
1691
1692	# MAX_LAKE_RESERVOIR ([kg/m2(routing area)] ) : Maximum limit of water in lake_reservoir {RIVER_ROUTING}
1693	MAX_LAKE_RESERVOIR = 7000
1694
1695	# RIVER_DESC ([FLAG]) : Writes out a description of the rivers {RIVER_ROUTING}
1696	RIVER_DESC = n
1697
1698	# RIVER_DESC_FILE ([FILE]) : Filename in which we write the description of the rivers. If suffix is ".nc" a netCDF file is created {RIVER_DESC}
1699	RIVER_DESC_FILE = river_desc.nc
1700
1701	# ROUTING_FILE ([FILE]) : Name of file which contains the routing information {RIVER_ROUTING}
1702	ROUTING_FILE = routing.nc
1703
1704	# IRRIGATION_FILE ([FILE]) : Name of file which contains the map of irrigated areas {DO_IRRIGATION OR DO_FLOODPLAINS}
1705	IRRIGATION_FILE = floodplains.nc
1706
1707	# ROUTING_METHOD (character string) : Choice of routing module to be used {RIVER_ROUTING=T}
1708	ROUTING_METHOD = standard
1709
1710	# EPS_CARBON ([%] ) : Allowed error on carbon stock {SPINUP_ANALYTIC}
1711	EPS_CARBON = 0.01
1712
1713	# SPINUP_PERIOD ([years] ) : Period to calulcate equilibrium during spinup analytic {SPINUP_ANALYTIC}
1714	SPINUP_PERIOD = -1
1715
1716	# STOMATE_FORCING_NAME ([FILE]) : Name of STOMATE's forcing file {OK_STOMATE}
1717	STOMATE_FORCING_NAME = NONE
1718
1719	# STOMATE_FORCING_MEMSIZE ([MegaBytes]) : Size of STOMATE forcing data in memory {OK_STOMATE}
1720	STOMATE_FORCING_MEMSIZE = 50
1721
1722	# STOMATE_CFORCING_NAME ([FILE]) : Name of STOMATE's carbon forcing file {OK_STOMATE}
1723	STOMATE_CFORCING_NAME = NONE
1724
1725	# FORCESOIL_STEP_PER_YEAR ([days, months, year]) : Number of time steps per year for carbon spinup. {OK_STOMATE}
1726	FORCESOIL_STEP_PER_YEAR = 365
1727
1728	# FORCESOIL_NB_YEAR ([years]) : Number of years saved for carbon spinup. {OK_STOMATE}
1729	FORCESOIL_NB_YEAR = 1
1730
1731	# XIOS_ORCHIDEE_OK ([FLAG]) : Use XIOS for writing diagnostics file {}
1732	XIOS_ORCHIDEE_OK = y
1733
1734	# XIOS_INTERPOLATION ([FLAG]) : Actiave reading and intrepolation using XIOS {XIOS_ORCHIDEE_OK}
1735	XIOS_INTERPOLATION = n
1736
1737	# XIOS_REMAP_OUTPUT ([FLAG]) : Actiave remaping of diagnostic output files to regular grid {XIOS_ORCHIDEE_OK .AND. grid_type=unstructured}
1738	XIOS_REMAP_OUTPUT = True
1739
1740	# DT_STOMATE ([seconds]) : Time step of STOMATE and other slow processes {OK_STOMATE}
1741	DT_STOMATE = 86400.
1742

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

Download in other formats:

Original Format