Context Navigation

orchidee.default @ 8612

Last change on this file since 8612 was 4719, checked in by albert.jornet, 7 years ago

Merge: from revisions [4491:4695/trunk/ORCHIDEE]

Merge done in [4671:4718/perso/albert.jornet/MICT_MERGE]

File size: 74.7 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	# SPRED_PREC ([-]) : Spread the precipitation. {[-]}
54	SPRED_PREC = Half of the forcing time step or uniform, depending on dt_force and dt_sechiba
55
56	# ATM_CO2 ([ppm]) : Value for atm CO2 {[-]}
57	ATM_CO2 = 350.
58
59	# ALLOW_WEATHERGEN ([FLAG]) : Allow weather generator to create data {[-]}
60	ALLOW_WEATHERGEN = n
61
62	# DT_WEATHGEN ([seconds]) : Calling frequency of weather generator {ALLOW_WEATHERGEN}
63	DT_WEATHGEN = 1800.
64
65	# LIMIT_WEST ([Degrees] ) : Western limit of region {[-]}
66	LIMIT_WEST = -180.
67
68	# LIMIT_EAST ([Degrees] ) : Eastern limit of region {[-]}
69	LIMIT_EAST = 180.
70
71	# LIMIT_NORTH ([Degrees]) : Northern limit of region {[-]}
72	LIMIT_NORTH = 90.
73
74	# LIMIT_SOUTH ([Degrees]) : Southern limit of region {[-]}
75	LIMIT_SOUTH = -90.
76
77	# MERID_RES ([Degrees]) : North-South Resolution {ALLOW_WEATHERGEN}
78	MERID_RES = 2.
79
80	# ZONAL_RES ([Degrees] ) : East-West Resolution {ALLOW_WEATHERGEN}
81	ZONAL_RES = 2.
82
83	# HEIGHT_LEV1 ([m]) : Height at which T and Q are given {offline mode}
84	HEIGHT_LEV1 = 2.0
85
86	# HEIGHT_LEVW ([m]) : Height at which the wind is given {offline mode}
87	HEIGHT_LEVW = 10.0
88
89	# NBUFF (-) : Number of time steps of data to buffer between each reading of the forcing file {OFF_LINE}
90	NBUFF = 15
91
92	# IPPREC ([-] ) : Use prescribed values {ALLOW_WEATHERGEN}
93	IPPREC = 0
94
95	# WEATHGEN_PRECIP_EXACT ([FLAG]) : Exact monthly precipitation {ALLOW_WEATHERGEN}
96	WEATHGEN_PRECIP_EXACT = n
97
98	# DUMP_WEATHER ([FLAG]) : Write weather from generator into a forcing file {ALLOW_WEATHERGEN }
99	DUMP_WEATHER = n
100
101	# DUMP_WEATHER_FILE ([FILE]) : Name of the file that contains the weather from generator {DUMP_WEATHER}
102	DUMP_WEATHER_FILE = weather_dump.nc
103
104	# DUMP_WEATHER_GATHERED ([FLAG]) : Dump weather data on gathered grid {DUMP_WEATHER}
105	DUMP_WEATHER_GATHERED = y
106
107	# HEIGHT_LEV1_DUMP ([m]) : {DUMP_WEATHER}
108	HEIGHT_LEV1_DUMP = 10.
109
110	#*******************************************************************************************
111	# ORCHIDEE parameters
112	#*******************************************************************************************
113
114	# SOILTYPE_CLASSIF ([-]) : Type of classification used for the map of soil types {!IMPOSE_VEG}
115	SOILTYPE_CLASSIF = zobler
116
117	# RIVER_ROUTING ([FLAG]) : Decides if we route the water or not {OK_SECHIBA}
118	RIVER_ROUTING = n
119
120	# HYDROL_CWRR ([FLAG]) : Allows to switch on the multilayer hydrology of CWRR {OK_SECHIBA}
121	HYDROL_CWRR = n
122
123	# DO_IRRIGATION ([FLAG]) : Should we compute an irrigation flux {RIVER_ROUTING }
124	DO_IRRIGATION = n
125
126	# DO_FLOODPLAINS ([FLAG] ) : Should we include floodplains {RIVER_ROUTING }
127	DO_FLOODPLAINS = n
128
129	# CHECK_WATERBAL ([FLAG] ) : Should we check the global water balance {OK_SECHIBA}
130	CHECK_WATERBAL = n
131
132	# OK_EXPLICITSNOW ([FLAG]) : Activate explict snow scheme {OK_SECHIBA}
133	OK_EXPLICITSNOW = FALSE
134
135	# STOMATE_OK_STOMATE ([FLAG]) : Activate STOMATE? {OK_SECHIBA}
136	STOMATE_OK_STOMATE = n
137
138	# STOMATE_OK_CO2 ([FLAG]) : Activate CO2? {OK_SECHIBA }
139	STOMATE_OK_CO2 = y if OK_STOMATE else n
140
141	# DO_WOOD_HARVEST ([FLAG]) : Activate Wood Harvest ? {OK_STOMATE}
142	DO_WOOD_HARVEST = n
143
144	# STOMATE_OK_DGVM ([FLAG]) : Activate DGVM? {OK_STOMATE}
145	STOMATE_OK_DGVM = n
146
147	# CHEMISTRY_BVOC ([FLAG]) : Activate calculations for BVOC {OK_SECHIBA}
148	CHEMISTRY_BVOC = n
149
150	# CHEMISTRY_LEAFAGE ([FLAG]) : Activate LEAFAGE? {CHEMISTRY_BVOC}
151	CHEMISTRY_LEAFAGE = n
152
153	# CANOPY_EXTINCTION ([FLAG]) : Use canopy radiative transfer model? {CHEMISTRY_BVOC }
154	CANOPY_EXTINCTION = n
155
156	# CANOPY_MULTILAYER ([FLAG]) : Use canopy radiative transfer model with multi-layers {CANOPY_EXTINCTION }
157	CANOPY_MULTILAYER = n
158
159	# NOx_RAIN_PULSE ([FLAG]) : Calculate NOx emissions with pulse? {CHEMISTRY_BVOC }
160	NOx_RAIN_PULSE = n
161
162	# NOx_BBG_FERTIL ([FLAG]) : Calculate NOx emissions with bbg fertilizing effect? {CHEMISTRY_BVOC }
163	NOx_BBG_FERTIL = n
164
165	# NOx_FERTILIZERS_USE ([FLAG] ) : Calculate NOx emissions with fertilizers use? {CHEMISTRY_BVOC }
166	NOx_FERTILIZERS_USE = n
167
168	# NVM ([-]) : number of PFTs {OK_SECHIBA or OK_STOMATE}
169	NVM = 13
170
171	# IMPOSE_PARAM ([FLAG]) : Do you impose the values of the parameters? {OK_SECHIBA or OK_STOMATE}
172	IMPOSE_PARAM = y
173
174	# THERMOSOIL_NBLEV (m) : Number of soil level {HDYROL_CWRR=FALSE}
175	THERMOSOIL_NBLEV = Number of soil level
176
177	# DEPTH_MAX_T (m) : Maximum depth of the soil thermodynamics {}
178	DEPTH_MAX_T = 10.0
179
180	# DEPTH_MAX_H (m) : Maximum depth of soil moisture {}
181	DEPTH_MAX_H = 2.0 or 4.0 depending on hydrol_cwrr
182
183	# DEPTH_TOPTHICK (m) : Thickness of upper most Layer {}
184	DEPTH_TOPTHICK = 9.77517107e-04
185
186	# DEPTH_CSTTHICK (m) : Depth at which constant layer thickness start {}
187	DEPTH_CSTTHICK = DEPTH_MAX_H
188
189	# REFINEBOTTOM (-) : Depth at which the hydrology layers will be refined towards the bottom. {}
190	REFINEBOTTOM = .FALSE.
191
192	# DEPTH_GEOM (m) : Depth at which we resume geometrical increases for temperature {}
193	DEPTH_GEOM = DEPTH_MAX_H
194
195	# RATIO_GEOM_BELOW (-) : Ratio of the geometrical series defining the thickness below DEPTH_GEOM {}
196	RATIO_GEOM_BELOW = 2
197
198	# ALMA_OUTPUT ([FLAG]) : Should the output follow the ALMA convention {OK_SECHIBA}
199	ALMA_OUTPUT = n
200
201	# OUTPUT_FILE ([FILE]) : Name of file in which the output is going to be written {OK_SECHIBA}
202	OUTPUT_FILE = sechiba_history.nc
203
204	# WRITE_STEP ([seconds]) : Frequency in seconds for sechiba_history.nc file with IOIPSL {OK_SECHIBA, NOT XIOS_ORCHIDEE_OK}
205	WRITE_STEP = 86400.
206
207	# SECHIBA_HISTLEVEL ([-]) : SECHIBA history output level (0..10) {OK_SECHIBA and HF}
208	SECHIBA_HISTLEVEL = 5
209
210	# SECHIBA_HISTFILE2 ([FLAG]) : Flag to switch on histfile 2 for SECHIBA (hi-frequency ?) {OK_SECHIBA}
211	SECHIBA_HISTFILE2 = n
212
213	# WRITE_STEP2 ([seconds]) : Frequency in seconds at which to WRITE output {SECHIBA_HISTFILE2}
214	WRITE_STEP2 = 1800.0
215
216	# SECHIBA_OUTPUT_FILE2 ([FILE]) : Name of file in which the output number 2 is going to be written {SECHIBA_HISTFILE2}
217	SECHIBA_OUTPUT_FILE2 = sechiba_out_2.nc
218
219	# SECHIBA_HISTLEVEL2 ([-] ) : SECHIBA history 2 output level (0..10) {SECHIBA_HISTFILE2}
220	SECHIBA_HISTLEVEL2 = 1
221
222	# STOMATE_OUTPUT_FILE ([FILE]) : Name of file in which STOMATE's output is going to be written {OK_STOMATE}
223	STOMATE_OUTPUT_FILE = stomate_history.nc
224
225	# STOMATE_HIST_DT ([days]) : STOMATE history time step {OK_STOMATE}
226	STOMATE_HIST_DT = 10.
227
228	# STOMATE_IPCC_OUTPUT_FILE ([FILE]) : Name of file in which STOMATE's output is going to be written {OK_STOMATE}
229	STOMATE_IPCC_OUTPUT_FILE = stomate_ipcc_history.nc
230
231	# STOMATE_IPCC_HIST_DT ([days]) : STOMATE IPCC history time step {OK_STOMATE}
232	STOMATE_IPCC_HIST_DT = 0.
233
234	# OK_HISTSYNC ([FLAG]) : Syncronize and write IOIPSL output files at each time step {}
235	OK_HISTSYNC = FALSE
236
237	# STOMATE_HISTLEVEL ([-]) : STOMATE history output level (0..10) {OK_STOMATE}
238	STOMATE_HISTLEVEL = 10
239
240	# SECHIBA_restart_in ([FILE]) : Name of restart to READ for initial conditions {OK_SECHIBA }
241	SECHIBA_restart_in = NONE
242
243	# SECHIBA_rest_out ([FILE]) : Name of restart files to be created by SECHIBA {OK_SECHIBA}
244	SECHIBA_rest_out = sechiba_rest_out.nc
245
246	# STOMATE_RESTART_FILEIN ([FILE]) : Name of restart to READ for initial conditions of STOMATE {STOMATE_OK_STOMATE}
247	STOMATE_RESTART_FILEIN = NONE
248
249	# STOMATE_RESTART_FILEOUT ([FILE]) : Name of restart files to be created by STOMATE {STOMATE_OK_STOMATE}
250	STOMATE_RESTART_FILEOUT = stomate_rest_out.nc
251
252	# FORCE_CO2_VEG ([FLAG]) : Flag to force the value of atmospheric CO2 for vegetation. {Only in coupled mode}
253	FORCE_CO2_VEG = FALSE
254
255	# ATM_CO2 ([ppm]) : Value for atm CO2 {FORCE_CO2_VEG (only in coupled mode)}
256	ATM_CO2 = 350.
257
258	# TAU_OUTFLOW ([days]) : Number of days over which the coastal- and riverflow will be distributed {Only in coupled mode}
259	TAU_OUTFLOW = 0
260
261	# ECCENTRICITY ([-]) : Use prescribed values {ALLOW_WEATHERGEN}
262	ECCENTRICITY = 0.016724
263
264	# PERIHELIE ([-]) : Use prescribed values {ALLOW_WEATHERGEN}
265	PERIHELIE = 102.04
266
267	# OBLIQUITY ([Degrees]) : Use prescribed values {ALLOW_WEATHERGEN}
268	OBLIQUITY = 23.446
269
270	# PFT_TO_MTC ([-]) : correspondance array linking a PFT to MTC {OK_SECHIBA or OK_STOMATE}
271	PFT_TO_MTC = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13
272
273	# PFT_NAME ([-]) : Name of a PFT {OK_SECHIBA or OK_STOMATE}
274	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
275
276	# LEAF_TAB ([-] ) : leaf type : 1 {OK_STOMATE}
277	LEAF_TAB = 4, 1, 1, 2, 1, 1, 2, 1, 2, 3, 3, 3, 3
278
279	# PHENO_MODEL ([-] ) : which phenology model is used? (tabulated) {OK_STOMATE}
280	PHENO_MODEL = none, none, moi, none, none, ncdgdd, none, ncdgdd, ngd, moigdd, moigdd, moigdd, moigdd
281
282	# SECHIBA_LAI ([m^2/m^2]) : laimax for maximum lai(see also type of lai interpolation) {OK_SECHIBA or IMPOSE_VEG}
283	SECHIBA_LAI = 0., 8., 8., 4., 4.5, 4.5, 4., 4.5, 4., 2., 2., 2., 2.
284
285	# LLAIMIN ([m^2/m^2]) : laimin for minimum lai(see also type of lai interpolation) {OK_SECHIBA or IMPOSE_VEG}
286	LLAIMIN = 0., 8., 0., 4., 4.5, 0., 4., 0., 0., 0., 0., 0., 0.
287
288	# SLOWPROC_HEIGHT ([m] ) : prescribed height of vegetation {OK_SECHIBA}
289	SLOWPROC_HEIGHT = 0., 30., 30., 20., 20., 20., 15., 15., 15., .5, .6, 1., 1.
290
291	# Z0_OVER_HEIGHT ([-] ) : factor to calculate roughness height from height of canopy {OK_SECHIBA}
292	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
293
294	# RATIO_Z0M_Z0H ([-]) : Ratio between z0m and z0h {OK_SECHIBA}
295	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
296
297	# TYPE_OF_LAI ([-]) : Type of behaviour of the LAI evolution algorithm {OK_SECHIBA}
298	TYPE_OF_LAI = inter, inter, inter, inter, inter, inter, inter, inter, inter, inter, inter, inter, inter
299
300	# NATURAL ([BOOLEAN]) : natural? {OK_SECHIBA, OK_STOMATE}
301	NATURAL = y, y, y, y, y, y, y, y, y, y, y, n, n
302
303	# IS_C4 ([BOOLEAN]) : flag for C4 vegetation types {OK_SECHIBA or OK_STOMATE}
304	IS_C4 = n, n, n, n, n, n, n, n, n, n, n, y, n, y
305
306	# VCMAX_FIX ([micromol/m^2/s] ) : values used for vcmax when STOMATE is not activated {OK_SECHIBA and NOT(OK_STOMATE)}
307	VCMAX_FIX = 0., 40., 50., 30., 35., 40.,30., 40., 35., 60., 60., 70., 70.
308
309	# DOWNREG_CO2 ([-]) : coefficient for CO2 downregulation (unitless) {OK_CO2}
310	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
311
312	# E_KmC ([J mol-1]) : Energy of activation for KmC {OK_CO2}
313	E_KmC = -9999., 79430., 79430., 79430., 79430., 79430., 79430., 79430., 79430., 79430., 79430., 79430., 79430.
314
315	# E_KmO ([J mol-1]) : Energy of activation for KmO {OK_CO2}
316	E_KmO = -9999., 36380., 36380., 36380., 36380., 36380., 36380., 36380., 36380., 36380., 36380., 36380., 36380.
317
318	# E_Sco ([J mol-1]) : Energy of activation for Sco {OK_CO2}
319	E_Sco = -9999., -24460., -24460., -24460., -24460., -24460., -24460., -24460., -24460., -24460., -24460., -24460., -24460.
320
321	# E_gamma_star ([J mol-1]) : Energy of activation for gamma_star {OK_CO2}
322	E_gamma_star = -9999., 37830., 37830., 37830., 37830., 37830., 37830., 37830., 37830., 37830., 37830., 37830., 37830.
323
324	# E_Vcmax ([J mol-1]) : Energy of activation for Vcmax {OK_CO2}
325	E_Vcmax = -9999., 71513., 71513., 71513., 71513., 71513., 71513., 71513., 71513., 71513., 67300., 71513., 67300.
326
327	# E_Jmax ([J mol-1]) : Energy of activation for Jmax {OK_CO2}
328	E_Jmax = -9999., 49884., 49884., 49884., 49884., 49884., 49884., 49884., 49884., 49884., 77900., 49884., 77900.
329
330	# aSV ([J K-1 mol-1]) : a coefficient of the linear regression (a+bT) defining the Entropy term for Vcmax {OK_CO2}
331	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
332
333	# bSV ([J K-1 mol-1 Â°C-1]) : b coefficient of the linear regression (a+bT) defining the Entropy term for Vcmax {OK_CO2}
334	bSV = -9999., -1.07, -1.07, -1.07, -1.07, -1.07, -1.07, -1.07, -1.07, -1.07, 0., -1.07, 0.
335
336	# TPHOTO_MIN ([-]) : minimum photosynthesis temperature (deg C) {OK_STOMATE}
337	TPHOTO_MIN = -9999., -4., -4., -4., -4.,-4.,-4., -4., -4., -4., -4., -4., -4.
338
339	# TPHOTO_MAX ([-]) : maximum photosynthesis temperature (deg C) {OK_STOMATE}
340	TPHOTO_MAX = -9999., 55., 55., 55., 55., 55., 55., 55., 55., 55., 55., 55., 55.
341
342	# aSJ ([J K-1 mol-1]) : a coefficient of the linear regression (a+bT) defining the Entropy term for Jmax {OK_CO2}
343	aSJ = -9999., 659.70, 659.70, 659.70, 659.70, 659.70, 659.70, 659.70, 659.70, 659.70, 630., 659.70, 630.
344
345	# bSJ ([J K-1 mol-1 Â°C-1]) : b coefficient of the linear regression (a+bT) defining the Entropy term for Jmax {OK_CO2}
346	bSJ = -9999., -0.75, -0.75, -0.75, -0.75, -0.75, -0.75, -0.75, -0.75, -0.75, 0., -0.75, 0.
347
348	# D_Vcmax ([J mol-1]) : Energy of deactivation for Vcmax {OK_CO2}
349	D_Vcmax = -9999., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 192000., 200000., 192000.
350
351	# D_Jmax ([J mol-1]) : Energy of deactivation for Jmax {OK_CO2}
352	D_Jmax = -9999., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 200000., 192000., 200000., 192000.
353
354	# E_gm ([J mol-1] ) : Energy of activation for gm {OK_CO2 }
355	E_gm = -9999., 49600., 49600., 49600., 49600., 49600., 49600., 49600., 49600., 49600., -9999., 49600., -9999.
356
357	# S_gm ([J K-1 mol-1] ) : Entropy term for gm {OK_CO2 }
358	S_gm = -9999., 1400., 1400., 1400., 1400., 1400., 1400., 1400., 1400., 1400., -9999., 1400., -9999.
359
360	# D_gm ([J mol-1] ) : Energy of deactivation for gm {OK_CO2 }
361	D_gm = -9999., 437400., 437400., 437400., 437400., 437400., 437400., 437400., 437400., 437400., -9999., 437400., -9999.
362
363	# E_Rd ([J mol-1]) : Energy of activation for Rd {OK_CO2}
364	E_Rd = -9999., 46390., 46390., 46390., 46390., 46390., 46390., 46390., 46390., 46390., 46390., 46390., 46390.
365
366	# VCMAX25 ([micromol/m^2/s]) : Maximum rate of Rubisco activity-limited carboxylation at 25Â°C {OK_STOMATE}
367	VCMAX25 = -9999., 65., 65., 35., 45., 55., 35., 45., 35., 70., 70., 70., 70.
368
369	# ARJV ([mu mol e- (mu mol CO2)-1]) : a coefficient of the linear regression (a+bT) defining the Jmax25/Vcmax25 ratio {OK_STOMATE}
370	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
371
372	# 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}
373	BRJV = -9999., -0.035, -0.035, -0.035, -0.035, -0.035, -0.035, -0.035, -0.035, -0.035, 0., -0.035, 0.
374
375	# KmC25 ([ubar]) : MichaelisâMenten constant of Rubisco for CO2 at 25Â°C {OK_CO2}
376	KmC25 = -9999., 404.9, 404.9, 404.9, 404.9, 404.9, 404.9, 404.9, 404.9, 404.9, 650., 404.9, 650.
377
378	# KmO25 ([ubar]) : MichaelisâMenten constant of Rubisco for O2 at 25Â°C {OK_CO2}
379	KmO25 = -9999., 278400., 278400., 278400., 278400., 278400., 278400., 278400., 278400., 278400., 450000., 278400., 450000.
380
381	# Sco25 ([bar bar-1]) : Relative CO2 /O2 specificity factor for Rubisco at 25ÃÂ°C {OK_CO2}
382	Sco25 = -9999., 2800., 2800., 2800., 2800., 2800., 2800., 2800., 2800., 2800., 2590., 2800., 2590.
383
384	# gm25 ([mol m-2 s-1 bar-1] ) : Mesophyll diffusion conductance at 25ÃÂÃÂ°C {OK_CO2 }
385	gm25 = -9999., 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, -9999., 0.4, -9999.
386
387	# gamma_star25 ([ubar]) : Ci-based CO2 compensation point in the absence of Rd at 25Â°C (ubar) {OK_CO2}
388	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
389
390	# a1 ([-]) : Empirical factor involved in the calculation of fvpd {OK_CO2}
391	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
392
393	# b1 ([-]) : Empirical factor involved in the calculation of fvpd {OK_CO2}
394	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
395
396	# g0 ([mol mâ2 sâ1 barâ1]) : Residual stomatal conductance when irradiance approaches zero {OK_CO2}
397	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
398
399	# h_protons ([mol mol-1]) : Number of protons required to produce one ATP {OK_CO2}
400	h_protons = -9999., 4., 4., 4., 4., 4., 4., 4., 4., 4., 4., 4., 4.
401
402	# fpsir ([-]) : Fraction of PSII eâ transport rate partitioned to the C4 cycle {OK_CO2}
403	fpsir = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.4, -9999., 0.4
404
405	# fQ ([-]) : Fraction of electrons at reduced plastoquinone that follow the Q-cycle {OK_CO2}
406	fQ = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 1., -9999., 1.
407
408	# fpseudo ([-]) : Fraction of electrons at PSI that follow pseudocyclic transport {OK_CO2}
409	fpseudo = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.1, -9999., 0.1
410
411	# kp ([mol mâ2 sâ1 barâ1]) : Initial carboxylation efficiency of the PEP carboxylase {OK_CO2}
412	kp = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.7, -9999., 0.7
413
414	# alpha ([-]) : Fraction of PSII activity in the bundle sheath {OK_CO2}
415	alpha = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.1, -9999., 0.1
416
417	# gbs ([mol mâ2 sâ1 barâ1]) : Bundle-sheath conductance {OK_CO2}
418	gbs = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.003, -9999., 0.003
419
420	# theta ([â]) : Convexity factor for response of J to irradiance {OK_CO2}
421	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
422
423	# alpha_LL ([mol eâ (mol photon)â1]) : Conversion efficiency of absorbed light into J at strictly limiting light {OK_CO2}
424	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
425
426	# STRESS_VCMAX ([-]) : Stress on vcmax {OK_SECHIBA or OK_STOMATE}
427	STRESS_VCMAX = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.
428
429	# STRESS_GS ([-]) : Stress on gs {OK_SECHIBA or OK_STOMATE}
430	STRESS_GS = 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.
431
432	# STRESS_GM ([-]) : Stress on gm {OK_SECHIBA or OK_STOMATE}
433	STRESS_GM = 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.
434
435	# EXT_COEFF ([-]) : extinction coefficient of the Monsi&Seaki relationship (1953) {OK_SECHIBA or OK_STOMATE}
436	EXT_COEFF = .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5
437
438	# EXT_COEFF_VEGETFRAC ([-]) : extinction coefficient used for the calculation of the bare soil fraction {OK_SECHIBA or OK_STOMATE}
439	EXT_COEFF_VEGETFRAC = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.
440
441	# HYDROL_HUMCSTE ([m]) : Root profile {OK_SECHIBA}
442	HYDROL_HUMCSTE = humcste_ref2m or humcste_ref4m depending on zmaxh
443
444	# PREF_SOIL_VEG ([-] ) : The soil tile number for each vegetation {OK_SECHIBA or OK_STOMATE}
445	PREF_SOIL_VEG = 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3
446
447	# RSTRUCT_CONST ([s/m]) : Structural resistance {OK_SECHIBA}
448	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
449
450	# KZERO ([kg/m^2/s]) : A vegetation dependent constant used in the calculation of the surface resistance. {OK_SECHIBA}
451	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
452
453	# RVEG_PFT ([-]) : Artificial parameter to increase or decrease canopy resistance. {OK_SECHIBA}
454	RVEG_PFT = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.
455
456	# WMAX_VEG ([kg/m^3]) : Maximum field capacity for each of the vegetations (Temporary): max quantity of water {OK_SECHIBA}
457	WMAX_VEG = 150., 150., 150., 150., 150., 150., 150.,150., 150., 150., 150., 150., 150.
458
459	# PERCENT_THROUGHFALL_PFT ([%]) : Percent by PFT of precip that is not intercepted by the canopy. Default value depend on run mode. {OK_SECHIBA}
460	PERCENT_THROUGHFALL_PFT = Case offline+CWRR [0. 0. 0....] else [30. 30. 30.....]
461
462	# SNOWA_AGED_VIS ([-]) : Minimum snow albedo value for each vegetation type after aging (dirty old snow), visible albedo {OK_SECHIBA}
463	SNOWA_AGED_VIS = 0.5, 0., 0., 0.15, 0.14, 0.14, 0.15, 0.14, 0.22, 0.35, 0.35, 0.35, 0.35
464
465	# SNOWA_AGED_NIR ([-]) : Minimum snow albedo value for each vegetation type after aging (dirty old snow), near infrared albedo {OK_SECHIBA}
466	SNOWA_AGED_NIR = 0.35, 0., 0., 0.14, 0.14, 0.14, 0.14, 0.14, 0.14, 0.18, 0.18, 0.18, 0.18
467
468	# 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}
469	SNOWA_DEC_VIS = 0.45, 0., 0., 0.1, 0.06, 0.11, 0.10, 0.11, 0.18, 0.60, 0.60, 0.60, 0.60
470
471	# 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}
472	SNOWA_DEC_NIR = 0.45, 0., 0., 0.06, 0.06, 0.11, 0.06, 0.11, 0.11, 0.52 ,0.52, 0.52, 0.52
473
474	# ALB_LEAF_VIS ([-]) : leaf albedo of vegetation type, visible albedo {OK_SECHIBA}
475	ALB_LEAF_VIS = .0, .0397, .0474, .0386, .0484, .0411, .041, .0541, .0435, .0524, .0508, .0509, .0606
476
477	# ALB_LEAF_NIR ([-]) : leaf albedo of vegetation type, near infrared albedo {OK_SECHIBA}
478	ALB_LEAF_NIR = .0, .227, .214, .193, .208, .244, .177, .218, .213, .252, .265, .272, .244
479
480	# ISO_ACTIVITY ([-]) : Biogenic activity for each age class : isoprene {CHEMISTRY_BVOC}
481	ISO_ACTIVITY = 0.5, 1.5, 1.5, 0.5
482
483	# METHANOL_ACTIVITY ([-]) : Isoprene emission factor for each age class : methanol {CHEMISTRY_BVOC}
484	METHANOL_ACTIVITY = 1., 1., 0.5, 0.5
485
486	# EM_FACTOR_ISOPRENE ([ugC/g/h] ) : Isoprene emission factor {CHEMISTRY_BVOC}
487	EM_FACTOR_ISOPRENE = 0., 24., 24., 8., 16., 45., 8., 18., 0.5, 12., 18., 5., 5.
488
489	# EM_FACTOR_MONOTERPENE ([ugC/g/h] ) : Monoterpene emission factor {CHEMISTRY_BVOC }
490	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
491
492	# C_LDF_MONO ([]) : Monoterpenes fraction dependancy to light {CHEMISTRY_BVOC}
493	C_LDF_MONO = 0.6
494
495	# C_LDF_SESQ ([]) : Sesquiterpenes fraction dependancy to light {CHEMISTRY_BVOC}
496	C_LDF_SESQ = 0.5
497
498	# C_LDF_METH ([]) : Methanol fraction dependancy to light {CHEMISTRY_BVOC}
499	C_LDF_METH = 0.8
500
501	# C_LDF_ACET ([]) : Acetone fraction dependancy to light {CHEMISTRY_BVOC}
502	C_LDF_ACET = 0.2
503
504	# EM_FACTOR_APINENE ([ugC/g/h] ) : Alfa pinene emission factor {CHEMISTRY_BVOC }
505	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
506
507	# EM_FACTOR_BPINENE ([ugC/g/h] ) : Beta pinene emission factor {CHEMISTRY_BVOC }
508	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
509
510	# EM_FACTOR_LIMONENE ([ugC/g/h] ) : Limonene emission factor {CHEMISTRY_BVOC}
511	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
512
513	# EM_FACTOR_MYRCENE ([ugC/g/h] ) : Myrcene emission factor {CHEMISTRY_BVOC}
514	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
515
516	# EM_FACTOR_SABINENE ([ugC/g/h] ) : Sabinene emission factor {CHEMISTRY_BVOC}
517	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
518
519	# EM_FACTOR_CAMPHENE ([ugC/g/h] ) : Camphene emission factor {CHEMISTRY_BVOC}
520	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
521
522	# EM_FACTOR_3CARENE ([ugC/g/h] ) : 3-Carene emission factor {CHEMISTRY_BVOC}
523	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
524
525	# EM_FACTOR_TBOCIMENE ([ugC/g/h] ) : T-beta-ocimene emission factor {CHEMISTRY_BVOC}
526	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
527
528	# EM_FACTOR_OTHERMONOT ([ugC/g/h] ) : Other monoterpenes emission factor {CHEMISTRY_BVOC}
529	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
530
531	# EM_FACTOR_SESQUITERP ([ugC/g/h] ) : Sesquiterpenes emission factor {CHEMISTRY_BVOC}
532	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
533
534	# C_BETA_MONO ([]) : Monoterpenes temperature dependency coefficient {CHEMISTRY_BVOC}
535	C_BETA_MONO = 0.1
536
537	# C_BETA_SESQ ([]) : Sesquiterpenes temperature dependency coefficient {CHEMISTRY_BVOC}
538	C_BETA_SESQ = 0.17
539
540	# C_BETA_METH ([]) : Methanol temperature dependency coefficient {CHEMISTRY_BVOC}
541	C_BETA_METH = 0.08
542
543	# C_BETA_ACET ([]) : Acetone temperature dependency coefficient {CHEMISTRY_BVOC}
544	C_BETA_ACET = 0.1
545
546	# C_BETA_OXYVOC ([]) : Other oxygenated BVOC temperature dependency coefficient {CHEMISTRY_BVOC}
547	C_BETA_OXYVOC = 0.13
548
549	# EM_FACTOR_ORVOC ([ugC/g/h] ) : ORVOC emissions factor {CHEMISTRY_BVOC }
550	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
551
552	# EM_FACTOR_OVOC ([ugC/g/h] ) : OVOC emissions factor {CHEMISTRY_BVOC}
553	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
554
555	# EM_FACTOR_MBO ([ugC/g/h] ) : MBO emissions factor {CHEMISTRY_BVOC }
556	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
557
558	# EM_FACTOR_METHANOL ([ugC/g/h] ) : Methanol emissions factor {CHEMISTRY_BVOC }
559	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.
560
561	# EM_FACTOR_ACETONE ([ugC/g/h] ) : Acetone emissions factor {CHEMISTRY_BVOC }
562	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
563
564	# EM_FACTOR_ACETAL ([ugC/g/h] ) : Acetaldehyde emissions factor {CHEMISTRY_BVOC}
565	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
566
567	# EM_FACTOR_FORMAL ([ugC/g/h] ) : Formaldehyde emissions factor {CHEMISTRY_BVOC }
568	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
569
570	# EM_FACTOR_ACETIC ([ugC/g/h] ) : Acetic Acid emissions factor {CHEMISTRY_BVOC }
571	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
572
573	# EM_FACTOR_FORMIC ([ugC/g/h] ) : Formic Acid emissions factor {CHEMISTRY_BVOC}
574	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
575
576	# EM_FACTOR_NO_WET ([ngN/m^2/s]) : NOx emissions factor wet soil emissions and exponential dependancy factor {CHEMISTRY_BVOC}
577	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
578
579	# EM_FACTOR_NO_DRY ([ngN/m^2/s] ) : NOx emissions factor dry soil emissions and exponential dependancy factor {CHEMISTRY_BVOC}
580	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
581
582	# LARCH ([-] ) : Larcher 1991 SAI/LAI ratio {CHEMISTRY_BVOC }
583	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
584
585	# SLA ([m^2/gC]) : specif leaf area {OK_STOMATE}
586	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
587
588	# AVAILABILITY_FACT ([-] ) : Calculate dynamic mortality in lpj_gap, pft dependent parameter {OK_STOMATE }
589	AVAILABILITY_FACT = -9999., 0.14, 0.14, 0.10, 0.10, 0.10, 0.05, 0.05, 0.05, -9999., -9999., -9999., -9999.
590
591	# R0 ([-] ) : Standard root allocation {OK_STOMATE }
592	R0 = -9999., .30, .30, .30, .30, .30, .30, .30, .30, .30, .30, .30, .30
593
594	# S0 ([-] ) : Standard sapwood allocation {OK_STOMATE }
595	S0 = -9999., .25, .25, .30, .30, .30, .30, .30, .30, .30, .30, .30, .30
596
597	# FRAC_GROWTHRESP ([-]) : fraction of GPP which is lost as growth respiration {OK_STOMATE}
598	FRAC_GROWTHRESP = -9999., .28, .28, .28, .28, .28, .28, .28, .28, .28, .28, .28, .28
599
600	# MAINT_RESP_SLOPE_C ([-]) : slope of maintenance respiration coefficient (1/K), constant c of aT^2+bT+c , tabulated {OK_STOMATE}
601	MAINT_RESP_SLOPE_C = -9999., .20, .20, .16, .16, .16, .16, .16, .16, .16, .12, .16, .12
602
603	# MAINT_RESP_SLOPE_B ([-]) : slope of maintenance respiration coefficient (1/K), constant b of aT^2+bT+c , tabulated {OK_STOMATE}
604	MAINT_RESP_SLOPE_B = -9999., .0, .0, .0, .0, .0, .0, .0, .0, -.00133, .0, -.00133, .0
605
606	# MAINT_RESP_SLOPE_A ([-]) : slope of maintenance respiration coefficient (1/K), constant a of aT^2+bT+c , tabulated {OK_STOMATE}
607	MAINT_RESP_SLOPE_A = -9999., .0, .0, .0, .0, .0, .0, .0, .0, .0, .0, .0, .0
608
609	# CM_ZERO_LEAF ([g/g/day]) : maintenance respiration coefficient at 0 deg C, for leaves, tabulated {OK_STOMATE}
610	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
611
612	# CM_ZERO_SAPABOVE ([g/g/day]) : maintenance respiration coefficient at 0 deg C,for sapwood above, tabulated {OK_STOMATE}
613	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
614
615	# CM_ZERO_SAPBELOW ([g/g/day]) : maintenance respiration coefficient at 0 deg C, for sapwood below, tabulated {OK_STOMATE}
616	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
617
618	# CM_ZERO_HEARTABOVE ([g/g/day]) : maintenance respiration coefficient at 0 deg C, for heartwood above, tabulated {OK_STOMATE }
619	CM_ZERO_HEARTABOVE = -9999., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.
620
621	# CM_ZERO_HEARTBELOW ([g/g/day] ) : maintenance respiration coefficient at 0 deg C,for heartwood below, tabulated {OK_STOMATE }
622	CM_ZERO_HEARTBELOW = -9999., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.
623
624	# CM_ZERO_ROOT ([g/g/day] ) : maintenance respiration coefficient at 0 deg C, for roots, tabulated {OK_STOMATE}
625	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
626
627	# CM_ZERO_FRUIT ([g/g/day] ) : maintenance respiration coefficient at 0 deg C, for fruits, tabulated {OK_STOMATE}
628	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
629
630	# CM_ZERO_CARBRES ([g/g/day] ) : maintenance respiration coefficient at 0 deg C, for carbohydrate reserve, tabulated {OK_STOMATE}
631	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
632
633	# FLAM ([-]) : flamability: critical fraction of water holding capacity {OK_STOMATE}
634	FLAM = -9999., .15, .25, .25, .25, .25, .25, .25, .25, .25, .25, .35, .35
635
636	# RESIST ([-]) : fire resistance {OK_STOMATE}
637	RESIST = -9999., .95, .90, .12, .50, .12, .12, .12, .12, .0, .0, .0, .0
638
639	# COEFF_LCCHANGE_1 ([-]) : Coeff of biomass export for the year {OK_STOMATE}
640	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
641
642	# COEFF_LCCHANGE_10 ([-]) : Coeff of biomass export for the decade {OK_STOMATE}
643	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
644
645	# COEFF_LCCHANGE_100 ([-]) : Coeff of biomass export for the century {OK_STOMATE}
646	COEFF_LCCHANGE_100 = -9999., 0., 0., 0.104, 0.104, 0.104, 0.104, 0.104, 0.104, 0.104, 0., 0.104, 0.
647
648	# LAI_MAX_TO_HAPPY ([-]) : threshold of LAI below which plant uses carbohydrate reserves {OK_STOMATE}
649	LAI_MAX_TO_HAPPY = -9999., .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5
650
651	# LAI_MAX ([m^2/m^2]) : maximum LAI, PFT-specific {OK_STOMATE}
652	LAI_MAX = -9999., 7., 7., 5., 5., 5., 4.5, 4.5, 3.0, 2.5, 2.5, 5.,5.
653
654	# PHENO_TYPE ([-]) : type of phenology, 0 {OK_STOMATE}
655	PHENO_TYPE = 0, 1, 3, 1, 1, 2, 1, 2, 2, 4, 4, 2, 3
656
657	# PHENO_GDD_CRIT_C ([-]) : critical gdd, tabulated (C), constant c of aT^2+bT+c {OK_STOMATE}
658	PHENO_GDD_CRIT_C = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 270., 400., 125., 400.
659
660	# PHENO_GDD_CRIT_B ([-]) : critical gdd, tabulated (C), constant b of aT^2+bT+c {OK_STOMATE}
661	PHENO_GDD_CRIT_B = -9999., -9999., -9999., -9999., -9999., -9999., -9999.,-9999., -9999., 6.25, 0., 0., 0.
662
663	# PHENO_GDD_CRIT_A ([-]) : critical gdd, tabulated (C), constant a of aT^2+bT+c {OK_STOMATE}
664	PHENO_GDD_CRIT_A = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0.03125, 0., 0., 0.
665
666	# PHENO_MOIGDD_T_CRIT ([C]) : Average temperature threashold for C4 grass used in pheno_moigdd {OK_STOMATE}
667	PHENO_MOIGDD_T_CRIT = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 22.0, -9999., -9999.
668
669	# NGD_CRIT ([days]) : critical ngd, tabulated. Threshold -5 degrees {OK_STOMATE}
670	NGD_CRIT = -9999., -9999., -9999., -9999., -9999., -9999., -9999., 0., -9999., -9999., -9999., -9999., -9999.
671
672	# NCDGDD_TEMP ([C] ) : critical temperature for the ncd vs. gdd function in phenology {OK_STOMATE}
673	NCDGDD_TEMP = -9999., -9999., -9999., -9999., -9999., 5., -9999., 0., -9999., -9999., -9999., -9999., -9999.
674
675	# HUM_FRAC ([%]) : critical humidity (relative to min/max) for phenology {OK_STOMATE}
676	HUM_FRAC = -9999., -9999., .5, -9999., -9999., -9999., -9999., -9999., -9999., .5, .5, .5,.5
677
678	# HUM_MIN_TIME ([days]) : minimum time elapsed since moisture minimum {OK_STOMATE}
679	HUM_MIN_TIME = -9999., -9999., 50., -9999., -9999., -9999., -9999., -9999., -9999., 35., 35., 75., 75.
680
681	# TAU_SAP ([days]) : sapwood -> heartwood conversion time {OK_STOMATE}
682	TAU_SAP = -9999., 730., 730., 730., 730., 730., 730., 730., 730., -9999., -9999., -9999., -9999.
683
684	# TAU_LEAFINIT ([days]) : time to attain the initial foliage using the carbohydrate reserve {OK_STOMATE}
685	TAU_LEAFINIT = -9999., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10.
686
687	# TAU_FRUIT ([days]) : fruit lifetime {OK_STOMATE}
688	TAU_FRUIT = -9999., 90., 90., 90., 90., 90., 90., 90., 90., -9999., -9999., -9999., -9999.
689
690	# ECUREUIL ([-]) : fraction of primary leaf and root allocation put into reserve {OK_STOMATE}
691	ECUREUIL = -9999., .0, 1., .0, .0, 1., .0, 1., 1., 1., 1., 1., 1.
692
693	# ALLOC_MIN ([-]) : minimum allocation above/below {OK_STOMATE}
694	ALLOC_MIN = -9999., 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, -9999., -9999., -9999., -9999.
695
696	# ALLOC_MAX ([-]) : maximum allocation above/below {OK_STOMATE}
697	ALLOC_MAX = -9999., 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, -9999., -9999., -9999., -9999.
698
699	# DEMI_ALLOC ([-]) : mean allocation above/below {OK_STOMATE}
700	DEMI_ALLOC = -9999., 5., 5., 5., 5., 5., 5., 5., 5., -9999., -9999., -9999., -9999.
701
702	# LEAFLIFE_TAB ([years]) : leaf longevity {OK_STOMATE}
703	LEAFLIFE_TAB = -9999., .5, 2., .33, 1., 2., .33, 2., 2., 2., 2., 2., 2.
704
705	# LEAFFALL ([days]) : length of death of leaves, tabulated {OK_STOMATE}
706	LEAFFALL = -9999., -9999., 10., -9999., -9999., 10., -9999., 10., 10., 10., 10., 10., 10.
707
708	# LEAFAGECRIT ([days]) : critical leaf age, tabulated {OK_STOMATE}
709	LEAFAGECRIT = -9999., 730., 180., 910., 730., 180., 910., 180., 180., 120., 120., 90., 90.
710
711	# SENESCENCE_TYPE ([-]) : type of senescence, tabulated {OK_STOMATE}
712	SENESCENCE_TYPE = none, none, dry, none, none, cold, none, cold, cold, mixed, mixed, mixed, mixed
713
714	# SENESCENCE_HUM ([-] ) : critical relative moisture availability for senescence {OK_STOMATE}
715	SENESCENCE_HUM = -9999., -9999., .3, -9999., -9999., -9999., -9999., -9999., -9999., .2, .2, .3, .2
716
717	# NOSENESCENCE_HUM ([-]) : relative moisture availability above which there is no humidity-related senescence {OK_STOMATE}
718	NOSENESCENCE_HUM = -9999., -9999., .8, -9999., -9999., -9999., -9999., -9999., -9999., .3, .3, .3, .3
719
720	# MAX_TURNOVER_TIME ([days]) : maximum turnover time for grasse {OK_STOMATE}
721	MAX_TURNOVER_TIME = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 80., 80., 80., 80.
722
723	# MIN_TURNOVER_TIME ([days]) : minimum turnover time for grasse {OK_STOMATE}
724	MIN_TURNOVER_TIME = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 10., 10., 10., 10.
725
726	# MIN_LEAF_AGE_FOR_SENESCENCE ([days] ) : minimum leaf age to allow senescence g {OK_STOMATE}
727	MIN_LEAF_AGE_FOR_SENESCENCE = -9999., -9999., 90., -9999., -9999., 90., -9999., 60., 60., 30., 30., 30., 30.
728
729	# SENESCENCE_TEMP_C ([-]) : critical temperature for senescence (C), constant c of aT^2+bT+c, tabulated {OK_STOMATE}
730	SENESCENCE_TEMP_C = -9999., -9999., -9999., -9999., -9999., 12., -9999., 7., 2., -1.375, 5., 5., 10.
731
732	# SENESCENCE_TEMP_B ([-]) : critical temperature for senescence (C), constant b of aT^2+bT+c ,tabulated {OK_STOMATE }
733	SENESCENCE_TEMP_B = -9999., -9999., -9999., -9999., -9999., 0., -9999., 0., 0., .1, 0., 0., 0.
734
735	# SENESCENCE_TEMP_A ([-] ) : critical temperature for senescence (C), constant a of aT^2+bT+c , tabulated {OK_STOMATE}
736	SENESCENCE_TEMP_A = -9999., -9999., -9999., -9999., -9999., 0., -9999., 0., 0.,.00375, 0., 0., 0.
737
738	# GDD_SENESCENCE ([days] ) : minimum gdd to allow senescence of crops {OK_STOMATE}
739	GDD_SENESCENCE = -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., -9999., 950., 4000.
740
741	# RESIDENCE_TIME ([years]) : residence time of trees {OK_DGVM and NOT(LPJ_GAP_CONST_MORT)}
742	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
743
744	# TMIN_CRIT ([C]) : critical tmin, tabulated {OK_STOMATE}
745	TMIN_CRIT = -9999., 0.0, 0.0, -30.0, -14.0, -30.0, -45.0, -45.0, -9999., -9999., -9999., -9999., -9999.
746
747	# TCM_CRIT ([C]) : critical tcm, tabulated {OK_STOMATE}
748	TCM_CRIT = -9999., -9999., -9999., 5.0, 15.5, 15.5, -8.0, -8.0, -8.0, -9999., -9999., -9999., -9999.
749
750	# HERBIVORES ([FLAG]) : herbivores allowed? {OK_STOMATE }
751	HERBIVORES = n
752
753	# TREAT_EXPANSION ([FLAG]) : treat expansion of PFTs across a grid cell? {OK_STOMATE }
754	TREAT_EXPANSION = n
755
756	# LPJ_GAP_CONST_MORT ([FLAG]) : Constant mortality {OK_STOMATE AND NOT OK_DGVM}
757	LPJ_GAP_CONST_MORT = y/n depending on OK_DGVM
758
759	# HARVEST_AGRI ([FLAG]) : Harvest model for agricultural PFTs. {OK_STOMATE }
760	HARVEST_AGRI = y
761
762	# FIRE_DISABLE ([FLAG]) : no fire allowed {OK_STOMATE }
763	FIRE_DISABLE = n
764
765	# SPINUP_ANALYTIC (BOOLEAN ) : Activation of the analytic resolution of the spinup. {OK_STOMATE}
766	SPINUP_ANALYTIC = n
767
768	# AGRICULTURE ([FLAG]) : agriculture allowed? {OK_SECHIBA or OK_STOMATE}
769	AGRICULTURE = y
770
771	# IMPOSE_VEG ([FLAG]) : Should the vegetation be prescribed ? {OK_SECHIBA or OK_STOMATE}
772	IMPOSE_VEG = n
773
774	# IMPOSE_SOILT ([FLAG]) : Should the soil type be prescribed ? {IMPOSE_VEG}
775	IMPOSE_SOILT = n
776
777	# LAI_MAP ([FLAG]) : Read the LAI map {OK_SECHIBA or OK_STOMATE}
778	LAI_MAP = n
779
780	# MAP_PFT_FORMAT ([FLAG]) : Read a land use vegetation map on PFT format {OK_SECHIBA or OK_STOMATE}
781	MAP_PFT_FORMAT = y
782
783	# VEGET_REINIT ([FLAG] ) : booleen to indicate that a new LAND USE file will be used. {MAP_PFT_FORMAT}
784	VEGET_REINIT = y
785
786	# VEGET_YEAR ([FLAG] ) : Year of the vegetation map to be read {MAP_PFT_FORMAT}
787	VEGET_YEAR = 1
788
789	# MAXMASS_SNOW ([kg/m^2] ) : The maximum mass of a snow {OK_SECHIBA or HYDROL_CWRR}
790	MAXMASS_SNOW = 3000.
791
792	# SNOWCRI ([kg/m^2] ) : Sets the amount above which only sublimation occures {OK_SECHIBA or HYDROL_CWRR}
793	SNOWCRI = 1.5
794
795	# MIN_WIND ([m/s]) : Minimum wind speed {OK_SECHIBA}
796	MIN_WIND = 0.1
797
798	# MAX_SNOW_AGE ([days?]) : Maximum period of snow aging {OK_SECHIBA}
799	MAX_SNOW_AGE = 50.
800
801	# SNOW_TRANS ([m] ) : Transformation time constant for snow {OK_SECHIBA}
802	SNOW_TRANS = 0.2
803
804	# OK_NUDGE_MC ([FLAG]) : Activate nudging of soil moisture {HYDROL_CWRR}
805	OK_NUDGE_MC = n
806
807	# NUDGE_TAU_MC ([-]) : Relaxation time for nudging of soil moisture expressed in fraction of the day {OK_NUDGE_MC}
808	NUDGE_TAU_MC = 1
809
810	# OK_NUDGE_SNOW ([FLAG]) : Activate nudging of snow variables {HYDROL_CWRR}
811	OK_NUDGE_SNOW = n
812
813	# NUDGE_TAU_SNOW ([-]) : Relaxation time for nudging of snow variables {OK_NUDGE_SNOW}
814	NUDGE_TAU_SNOW = 1
815
816	# NUDGE_INTERPOL_WITH_XIOS ([FLAG]) : Activate reading and interpolation with XIOS for nudging fields {OK_NUDGE_MC or OK_NUDGE_SNOW}
817	NUDGE_INTERPOL_WITH_XIOS = n
818
819	# HEIGHT_DISPLACEMENT ([m] ) : Magic number which relates the height to the displacement height. {OK_SECHIBA }
820	HEIGHT_DISPLACEMENT = 0.75
821
822	# Z0_BARE ([m] ) : bare soil roughness length {OK_SECHIBA }
823	Z0_BARE = 0.01
824
825	# Z0_ICE ([m] ) : ice roughness length {OK_SECHIBA }
826	Z0_ICE = 0.001
827
828	# TCST_SNOWA ([days]) : Time constant of the albedo decay of snow {OK_SECHIBA }
829	TCST_SNOWA = 10.0
830
831	# SNOWCRI_ALB ([cm] ) : Critical value for computation of snow albedo {OK_SECHIBA}
832	SNOWCRI_ALB = 10.
833
834	# VIS_DRY ([-] ) : The correspondance table for the soil color numbers and their albedo {OK_SECHIBA }
835	VIS_DRY = 0.24, 0.22, 0.20, 0.18, 0.16, 0.14, 0.12, 0.10, 0.27
836
837	# NIR_DRY ([-] ) : The correspondance table for the soil color numbers and their albedo {OK_SECHIBA }
838	NIR_DRY = 0.48, 0.44, 0.40, 0.36, 0.32, 0.28, 0.24, 0.20, 0.55
839
840	# VIS_WET ([-] ) : The correspondance table for the soil color numbers and their albedo {OK_SECHIBA }
841	VIS_WET = 0.12, 0.11, 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.15
842
843	# NIR_WET ([-] ) : The correspondance table for the soil color numbers and their albedo {OK_SECHIBA }
844	NIR_WET = 0.24, 0.22, 0.20, 0.18, 0.16, 0.14, 0.12, 0.10, 0.31
845
846	# ALBSOIL_VIS ([-] ) : {OK_SECHIBA }
847	ALBSOIL_VIS = 0.18, 0.16, 0.16, 0.15, 0.12, 0.105, 0.09, 0.075, 0.25
848
849	# ALBSOIL_NIR ([-] ) : {OK_SECHIBA }
850	ALBSOIL_NIR = 0.36, 0.34, 0.34, 0.33, 0.30, 0.25, 0.20, 0.15, 0.45
851
852	# ALB_DEADLEAF ([-] ) : albedo of dead leaves, VIS+NIR {OK_SECHIBA }
853	ALB_DEADLEAF = 0.12, 0.35
854
855	# ALB_ICE ([-] ) : albedo of ice, VIS+NIR {OK_SECHIBA}
856	ALB_ICE = 0.60, 0.20
857
858	# CONDVEG_SNOWA ([-]) : The snow albedo used by SECHIBA {OK_SECHIBA}
859	CONDVEG_SNOWA = 1.E+20
860
861	# ALB_BARE_MODEL ([FLAG]) : Switch bare soil albedo dependent (if TRUE) on soil wetness {OK_SECHIBA}
862	ALB_BARE_MODEL = n
863
864	# ALB_BG_MODIS ([FLAG]) : Read bare soil albedo from file with background MODIS data {OK_SECHIBA}
865	ALB_BG_MODIS = n
866
867	# IMPOSE_AZE ([FLAG]) : Should the surface parameters be prescribed {OK_SECHIBA}
868	IMPOSE_AZE = n
869
870	# CONDVEG_Z0 ([m]) : Surface roughness {IMPOSE_AZE}
871	CONDVEG_Z0 = 0.15
872
873	# ROUGHHEIGHT ([m] ) : Height to be added to the height of the first level {IMPOSE_AZE}
874	ROUGHHEIGHT = 0.0
875
876	# CONDVEG_ALBVIS ([-]) : SW visible albedo for the surface {IMPOSE_AZE}
877	CONDVEG_ALBVIS = 0.25
878
879	# CONDVEG_ALBNIR ([-] ) : SW near infrared albedo for the surface {IMPOSE_AZE}
880	CONDVEG_ALBNIR = 0.25
881
882	# CONDVEG_EMIS ([-] ) : Emissivity of the surface for LW radiation {IMPOSE_AZE}
883	CONDVEG_EMIS = 1.0
884
885	# ROUGH_DYN ([FLAG]) : Account for a dynamic roughness height {OK_SECHIBA}
886	ROUGH_DYN = y
887
888	# C1 ([-] ) : Constant used in the formulation of the ratio of {ROUGH_DYN}
889	C1 = 0.32
890
891	# C2 ([-] ) : Constant used in the formulation of the ratio of {ROUGH_DYN}
892	C2 = 0.264
893
894	# C3 ([-] ) : Constant used in the formulation of the ratio of {ROUGH_DYN}
895	C3 = 15.1
896
897	# Cdrag_foliage ([-] ) : Drag coefficient of the foliage {ROUGH_DYN}
898	Cdrag_foliage = 0.2
899
900	# Ct ([-] ) : Heat transfer coefficient of the leaf {ROUGH_DYN}
901	Ct = 0.01
902
903	# Prandtl ([-] ) : Prandtl number used in the calculation of Ct* {ROUGH_DYN}
904	Prandtl = 0.71
905
906	# xansmax ([-] ) : maximum snow albedo {OK_SECHIBA}
907	xansmax = 0.85
908
909	# xansmin ([-] ) : minimum snow albedo {OK_SECHIBA}
910	xansmin = 0.50
911
912	# xans_todry ([S-1] ) : albedo decay rate for the dry snow {OK_SECHIBA}
913	xans_todry = 0.008
914
915	# xans_t ([S-1] ) : albedo decay rate for the wet snow {OK_SECHIBA}
916	xans_t = 0.24
917
918	# xrhosmax ([-] ) : maximum snow density {OK_SECHIBA}
919	xrhosmax = 750
920
921	# xwsnowholdmax1 ([-] ) : snow holding capacity 1 {OK_SECHIBA}
922	xwsnowholdmax1 = 0.03
923
924	# xwsnowholdmax2 ([-] ) : snow holding capacity 2 {OK_SECHIBA}
925	xwsnowholdmax2 = 0.10
926
927	# xsnowrhohold ([kg/m3] ) : snow density {OK_SECHIBA}
928	xsnowrhohold = 200.0
929
930	# ZSNOWTHRMCOND1 ([W/m/K] ) : Thermal conductivity Coef 1 {OK_SECHIBA}
931	ZSNOWTHRMCOND1 = 0.02
932
933	# ZSNOWTHRMCOND2 ([W m5/(kg2 K)] ) : Thermal conductivity Coef 2 {OK_SECHIBA}
934	ZSNOWTHRMCOND2 = 2.5E-6
935
936	# ZSNOWTHRMCOND_AVAP ([W/m/K] ) : Thermal conductivity Coef 1 water vapor {OK_SECHIBA}
937	ZSNOWTHRMCOND_AVAP = -0.06023
938
939	# ZSNOWTHRMCOND_BVAP ([W/m] ) : Thermal conductivity Coef 2 water vapor {OK_SECHIBA}
940	ZSNOWTHRMCOND_BVAP = -2.5425
941
942	# ZSNOWTHRMCOND_CVAP ([K] ) : Thermal conductivity Coef 3 water vapor {OK_SECHIBA}
943	ZSNOWTHRMCOND_CVAP = -289.99
944
945	# ZSNOWCMPCT_RHOD ([kg/m3]) : Snow compaction coefficent {OK_SECHIBA}
946	ZSNOWCMPCT_RHOD = 150.0
947
948	# ZSNOWCMPCT_ACM ([1/s]) : Coefficent for the thermal conductivity {OK_SECHIBA}
949	ZSNOWCMPCT_ACM = 2.8e-6
950
951	# ZSNOWCMPCT_BCM ([1/K]) : Coefficent for the thermal conductivity {OK_SECHIBA}
952	ZSNOWCMPCT_BCM = 0.04
953
954	# ZSNOWCMPCT_CCM ([m3/kg] ) : Coefficent for the thermal conductivity {OK_SECHIBA}
955	ZSNOWCMPCT_CCM = 460.
956
957	# ZSNOWCMPCT_V0 ([Pa/s]) : Vapor coefficent for the thermal conductivity {OK_SECHIBA}
958	ZSNOWCMPCT_V0 = 3.7e7
959
960	# ZSNOWCMPCT_VT ([1/K]) : Vapor coefficent for the thermal conductivity {OK_SECHIBA}
961	ZSNOWCMPCT_VT = 0.081
962
963	# ZSNOWCMPCT_VR ([m3/kg]) : Vapor coefficent for the thermal conductivity {OK_SECHIBA}
964	ZSNOWCMPCT_VR = 0.018
965
966	# CB ([-] ) : Constant of the Louis scheme {OK_SECHIBA}
967	CB = 5.0
968
969	# CC ([-] ) : Constant of the Louis scheme {OK_SECHIBA}
970	CC = 5.0
971
972	# CD ([-] ) : Constant of the Louis scheme {OK_SECHIBA}
973	CD = 5.0
974
975	# RAYT_CSTE ([W.m^{-2}] ) : Constant in the computation of surface resistance {OK_SECHIBA}
976	RAYT_CSTE = 125
977
978	# DEFC_PLUS ([K.W^{-1}] ) : Constant in the computation of surface resistance {OK_SECHIBA}
979	DEFC_PLUS = 23.E-3
980
981	# DEFC_MULT ([K.W^{-1}] ) : Constant in the computation of surface resistance {OK_SECHIBA}
982	DEFC_MULT = 1.5
983
984	# NLAI ([-] ) : Number of LAI levels {OK_SECHIBA}
985	NLAI = 20
986
987	# LAIMAX ([m^2/m^2] ) : Maximum LAI {OK_SECHIBA}
988	LAIMAX =
989
990	# DEW_VEG_POLY_COEFF ([-] ) : coefficients of the polynome of degree 5 for the dew {OK_SECHIBA}
991	DEW_VEG_POLY_COEFF = 0.887773, 0.205673, 0.110112, 0.014843, 0.000824, 0.000017
992
993	# DOWNREGULATION_CO2 ([FLAG] ) : Activation of CO2 downregulation {OK_SECHIBA}
994	DOWNREGULATION_CO2 = n
995
996	# DOWNREGULATION_CO2_BASELEVEL ([ppm] ) : CO2 base level {OK_SECHIBA }
997	DOWNREGULATION_CO2_BASELEVEL = 280.
998
999	# GB_REF ([s m-1] ) : Leaf bulk boundary layer resistance {OK_CO2}
1000	GB_REF = 1./25.
1001
1002	# CLAYFRACTION_DEFAULT ([-] ) : default fraction of clay {OK_SECHIBA }
1003	CLAYFRACTION_DEFAULT = 0.2
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 ([-] ) : {OK_CO2}
1012	LAI_LEVEL_DEPTH = 0.15
1013
1014	# Oi ([ubar] ) : Intercellular oxygen partial pressure {OK_CO2}
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	# ALWAYS_INIT ([-] ) : take carbon from atmosphere if carbohydrate reserve too small? {OK_STOMATE }
1249	ALWAYS_INIT = n
1250
1251	# MIN_GROWTHINIT_TIME ([days] ) : minimum time since last beginning of a growing season {OK_STOMATE }
1252	MIN_GROWTHINIT_TIME = 300.
1253
1254	# MOIAVAIL_ALWAYS_TREE ([-] ) : moisture availability above which moisture tendency doesn't matter {OK_STOMATE }
1255	MOIAVAIL_ALWAYS_TREE = 1.0
1256
1257	# MOIAVAIL_ALWAYS_GRASS ([-] ) : moisture availability above which moisture tendency doesn't matter {OK_STOMATE }
1258	MOIAVAIL_ALWAYS_GRASS = 0.6
1259
1260	# T_ALWAYS_ADD ([C] ) : monthly temp. above which temp. tendency doesn't matter {OK_STOMATE }
1261	T_ALWAYS_ADD = 10.
1262
1263	# GDDNCD_REF ([-] ) : {OK_STOMATE }
1264	GDDNCD_REF = 603.
1265
1266	# GDDNCD_CURVE ([-] ) : {OK_STOMATE }
1267	GDDNCD_CURVE = 0.0091
1268
1269	# GDDNCD_OFFSET ([-] ) : {OK_STOMATE }
1270	GDDNCD_OFFSET = 64.
1271
1272	# BM_SAPL_RESCALE ([-] ) : {OK_STOMATE }
1273	BM_SAPL_RESCALE = 40.
1274
1275	# MAINT_RESP_MIN_VMAX ([-] ) : {OK_STOMATE }
1276	MAINT_RESP_MIN_VMAX = 0.3
1277
1278	# MAINT_RESP_COEFF ([-] ) : {OK_STOMATE }
1279	MAINT_RESP_COEFF = 1.4
1280
1281	# FRAC_CARB_AP ([-]) : frac carb coefficients from active pool: depends on clay content {OK_STOMATE }
1282	FRAC_CARB_AP = 0.004
1283
1284	# FRAC_CARB_SA ([-]) : frac_carb_coefficients from slow pool {OK_STOMATE }
1285	FRAC_CARB_SA = 0.42
1286
1287	# FRAC_CARB_SP ([-] ) : frac_carb_coefficients from slow pool {OK_STOMATE }
1288	FRAC_CARB_SP = 0.03
1289
1290	# FRAC_CARB_PA ([-]) : frac_carb_coefficients from passive pool {OK_STOMATE }
1291	FRAC_CARB_PA = 0.45
1292
1293	# FRAC_CARB_PS ([-]) : frac_carb_coefficients from passive pool {OK_STOMATE }
1294	FRAC_CARB_PS = 0.0
1295
1296	# ACTIVE_TO_PASS_CLAY_FRAC ([-] ) : {OK_STOMATE }
1297	ACTIVE_TO_PASS_CLAY_FRAC = 0.68
1298
1299	# CARBON_TAU_IACTIVE ( [days] ) : residence times in carbon pools {OK_STOMATE }
1300	CARBON_TAU_IACTIVE = 0.149
1301
1302	# CARBON_TAU_ISLOW ([days]) : residence times in carbon pools {OK_STOMATE }
1303	CARBON_TAU_ISLOW = 5.48
1304
1305	# CARBON_TAU_IPASSIVE ([days] ) : residence times in carbon pools {OK_STOMATE }
1306	CARBON_TAU_IPASSIVE = 241.
1307
1308	# FLUX_TOT_COEFF ([days] ) : {OK_STOMATE }
1309	FLUX_TOT_COEFF = 1.2, 1.4,.75
1310
1311	# NEW_TURNOVER_TIME_REF ([days] ) : {OK_STOMATE }
1312	NEW_TURNOVER_TIME_REF = 20.
1313
1314	# LEAF_AGE_CRIT_TREF ([days] ) : {OK_STOMATE }
1315	LEAF_AGE_CRIT_TREF = 20.
1316
1317	# LEAF_AGE_CRIT_COEFF ([-] ) : {OK_STOMATE }
1318	LEAF_AGE_CRIT_COEFF = 1.5, 0.75, 10.
1319
1320	# VMAX_OFFSET ([-] ) : offset (minimum relative vcmax) {OK_STOMATE }
1321	VMAX_OFFSET = 0.3
1322
1323	# LEAFAGE_FIRSTMAX ([-] ) : leaf age at which vmax attains vcmax_opt (in fraction of critical leaf age) {OK_STOMATE }
1324	LEAFAGE_FIRSTMAX = 0.03
1325
1326	# LEAFAGE_LASTMAX ([-] ) : leaf age at which vmax falls below vcmax_opt (in fraction of critical leaf age) {OK_STOMATE }
1327	LEAFAGE_LASTMAX = 0.5
1328
1329	# LEAFAGE_OLD ([-] ) : leaf age at which vmax attains its minimum (in fraction of critical leaf age) {OK_STOMATE }
1330	LEAFAGE_OLD = 1.
1331
1332	# GPPFRAC_DORMANCE ([-]) : rapport maximal GPP/GGP_max pour dormance {OK_STOMATE }
1333	GPPFRAC_DORMANCE = 0.2
1334
1335	# TAU_CLIMATOLOGY ([days]) : tau for "climatologic variables {OK_STOMATE }
1336	TAU_CLIMATOLOGY = 20
1337
1338	# HVC1 ([-] ) : parameters for herbivore activity {OK_STOMATE }
1339	HVC1 = 0.019
1340
1341	# HVC2 ([-] ) : parameters for herbivore activity {OK_STOMATE }
1342	HVC2 = 1.38
1343
1344	# LEAF_FRAC_HVC ([-] ) : parameters for herbivore activity {OK_STOMATE }
1345	LEAF_FRAC_HVC = 0.33
1346
1347	# TLONG_REF_MAX ([K] ) : maximum reference long term temperature {OK_STOMATE }
1348	TLONG_REF_MAX = 303.1
1349
1350	# TLONG_REF_MIN ([K] ) : minimum reference long term temperature {OK_STOMATE }
1351	TLONG_REF_MIN = 253.1
1352
1353	# NCD_MAX_YEAR ([days]) : {OK_STOMATE }
1354	NCD_MAX_YEAR = 3.
1355
1356	# GDD_THRESHOLD ([days] ) : {OK_STOMATE }
1357	GDD_THRESHOLD = 5.
1358
1359	# GREEN_AGE_EVER ([-] ) : {OK_STOMATE }
1360	GREEN_AGE_EVER = 2.
1361
1362	# GREEN_AGE_DEC ([-] ) : {OK_STOMATE }
1363	GREEN_AGE_DEC = 0.5
1364
1365	# ESTAB_MAX_TREE ([-] ) : Maximum tree establishment rate {OK_DGVM}
1366	ESTAB_MAX_TREE = 0.12
1367
1368	# ESTAB_MAX_GRASS ([-] ) : Maximum grass establishment rate {OK_DGVM}
1369	ESTAB_MAX_GRASS = 0.12
1370
1371	# ESTABLISH_SCAL_FACT ([-] ) : {OK_DGVM }
1372	ESTABLISH_SCAL_FACT = 5.
1373
1374	# MAX_TREE_COVERAGE ([-] ) : {OK_DGVM }
1375	MAX_TREE_COVERAGE = 0.98
1376
1377	# IND_0_ESTAB ([-] ) : {OK_DGVM }
1378	IND_0_ESTAB = 0.2
1379
1380	# 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}
1381	ANNUAL_INCREASE = y
1382
1383	# MIN_COVER ([-] ) : For trees, minimum fraction of crown area occupied {OK_DGVM}
1384	MIN_COVER = 0.05
1385
1386	# IND_0 ([-] ) : initial density of individuals {OK_DGVM}
1387	IND_0 = 0.02
1388
1389	# MIN_AVAIL ([-] ) : minimum availability {OK_DGVM}
1390	MIN_AVAIL = 0.01
1391
1392	# RIP_TIME_MIN ([year] ) : {OK_DGVM}
1393	RIP_TIME_MIN = 1.25
1394
1395	# NPP_LONGTERM_INIT ([gC/m^2/year]) : {OK_DGVM}
1396	NPP_LONGTERM_INIT = 10.
1397
1398	# EVERYWHERE_INIT ([-] ) : {OK_DGVM}
1399	EVERYWHERE_INIT = 0.05
1400
1401	# PRINTLEV ([0, 1, 2, 3, 4]) : Print level for text output {}
1402	PRINTLEV = 2
1403
1404	# PRINTLEV_modname ([0, 1, 2, 3, 4]) : Specific print level of text output for the module "modname". Default as PRINTLEV. {}
1405	PRINTLEV_modname = PRINTLEV
1406
1407	# DRY_SOIL_HEAT_CAPACITY ([J.m^{-3}.K^{-1}] ) : Dry soil Heat capacity of soils {OK_SECHIBA }
1408	DRY_SOIL_HEAT_CAPACITY = 1.80e+6
1409
1410	# DRY_SOIL_HEAT_COND ([W.m^{-2}.K^{-1}] ) : Dry soil Thermal Conductivity of soils {OK_SECHIBA}
1411	DRY_SOIL_HEAT_COND = 0.40
1412
1413	# WET_SOIL_HEAT_CAPACITY ([J.m^{-3}.K^{-1}]) : Wet soil Heat capacity of soils {OK_SECHIBA}
1414	WET_SOIL_HEAT_CAPACITY = 3.03e+6
1415
1416	# WET_SOIL_HEAT_COND ([W.m^{-2}.K^{-1}]) : Wet soil Thermal Conductivity of soils {OK_SECHIBA }
1417	WET_SOIL_HEAT_COND = 1.89
1418
1419	# SNOW_HEAT_COND ([W.m^{-2}.K^{-1}]) : Thermal Conductivity of snow {OK_SECHIBA }
1420	SNOW_HEAT_COND = 0.3
1421
1422	# SNOW_DENSITY ([-] ) : Snow density for the soil thermodynamics {OK_SECHIBA }
1423	SNOW_DENSITY = 330.0
1424
1425	# NOBIO_WATER_CAPAC_VOLUMETRI ([s/m^2]) : {OK_SECHIBA and .NOT.(HYDROL_CWRR)}
1426	NOBIO_WATER_CAPAC_VOLUMETRI = 150.
1427
1428	# SECHIBA_QSINT ([m]) : Interception reservoir coefficient {OK_SECHIBA }
1429	SECHIBA_QSINT = 0.1
1430
1431	# CHOISNEL_DIFF_MIN ([kg/m^2/dt]) : Diffusion constant for the slow regime {OK_SECHIBA and .NOT.(HYDROL_CWRR)}
1432	CHOISNEL_DIFF_MIN = 0.001
1433
1434	# CHOISNEL_DIFF_MAX ([kg/m^2/dt]) : Diffusion constant for the fast regime {OK_SECHIBA and .NOT.(HYDROL_CWRR)}
1435	CHOISNEL_DIFF_MAX = 0.1
1436
1437	# CHOISNEL_DIFF_EXP ([-]) : The exponential in the diffusion law {OK_SECHIBA and .NOT.(HYDROL_CWRR)}
1438	CHOISNEL_DIFF_EXP = 1.5
1439
1440	# CHOISNEL_RSOL_CSTE ([s/m^2]) : Constant in the computation of resistance for bare soil evaporation {OK_SECHIBA and .NOT.(HYDROL_CWRR)}
1441	CHOISNEL_RSOL_CSTE = 33.E3
1442
1443	# HCRIT_LITTER ([m]) : Scaling depth for litter humidity {OK_SECHIBA and .NOT.(HYDROL_CWRR) }
1444	HCRIT_LITTER = 0.08
1445
1446	# OK_FREEZE ([FLAG]) : Activate the complet soil freezing scheme {OK_SECHIBA }
1447	OK_FREEZE = FALSE
1448
1449	# READ_REFTEMP ([FLAG]) : Initialize soil temperature using climatological temperature {}
1450	READ_REFTEMP = True/False depening on OK_FREEZE
1451
1452	# OK_FREEZE_THERMIX ([FLAG]) : Activate thermal part of the soil freezing scheme {}
1453	OK_FREEZE_THERMIX = True if OK_FREEZE else false
1454
1455	# OK_ECORR ([FLAG]) : Energy correction for freezing {OK_FREEZE_THERMIX}
1456	OK_ECORR = True if OK_FREEZE else false
1457
1458	# POROS ([-] ) : Soil porosity {OK_SECHIBA}
1459	POROS = 0.41
1460
1461	# fr_dT ([K] ) : Freezing window {OK_SECHIBA}
1462	fr_dT = 2.0
1463
1464	# OK_FREEZE_CWRR ([FLAG]) : CWRR freezing scheme by I. Gouttevin {}
1465	OK_FREEZE_CWRR = True if OK_FREEZE else false
1466
1467	# OK_THERMODYNAMICAL_FREEZING ([FLAG]) : Calculate frozen fraction thermodynamically {HYDROL_CWRR .AND. OK_FREEZE_CWRR}
1468	OK_THERMODYNAMICAL_FREEZING = True
1469
1470	# CHECK_CWRR ([FLAG]) : Check detailed CWRR water balance {HYDROL_CWRR}
1471	CHECK_CWRR = n
1472
1473	# CHECK_CWRR2 ([FLAG]) : Caluculate diagnostics to check CWRR water balance {HYDROL_CWRR2}
1474	CHECK_CWRR2 = n
1475
1476	# VEGET_UPDATE ([years]) : Update vegetation frequency {MAP_PFT_FORMAT}
1477	VEGET_UPDATE = 0Y
1478
1479	# SECHIBA_ZCANOP ([m]) : Soil level used for canopy development (if STOMATE disactivated) {OK_SECHIBA and .NOT. OK_STOMATE }
1480	SECHIBA_ZCANOP = 0.5
1481
1482	# SECHIBA_QSINT ([m]) : Interception reservoir coefficient {OK_SECHIBA }
1483	SECHIBA_QSINT = 0.1
1484
1485	# SECHIBA_VEGMAX ([-]) : Maximum vegetation distribution within the mesh (0-dim mode) {IMPOSE_VEG}
1486	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
1487
1488	# SECHIBA_FRAC_NOBIO ([-]) : Fraction of other surface types within the mesh (0-dim mode) {IMPOSE_VEG}
1489	SECHIBA_FRAC_NOBIO = 0.0
1490
1491	# SECHIBA_LAI ([-]) : LAI for all vegetation types (0-dim mode) {IMPOSE_VEG}
1492	SECHIBA_LAI = 0., 8., 8., 4., 4.5, 4.5, 4., 4.5, 4., 2., 2., 2., 2.
1493
1494	# SOIL_FRACTIONS ([-]) : Fraction of the 3 soil types (0-dim mode) {IMPOSE_VEG and IMPOSE_SOILT}
1495	SOIL_FRACTIONS = -9999._sechiba
1496
1497	# CLAY_FRACTION ([-] ) : Fraction of the clay fraction (0-dim mode) {IMPOSE_VEG and IMPOSE_SOIL}
1498	CLAY_FRACTION = 0.2
1499
1500	# REINF_SLOPE ([-]) : Slope coef for reinfiltration {IMPOSE_VEG}
1501	REINF_SLOPE = 0.1
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	# GET_SLOPE ([FLAG]) : Read slopes from file and do the interpolation {}
1507	GET_SLOPE = n
1508
1509	# LAI_FILE ([FILE]) : Name of file from which the vegetation map is to be read {LAI_MAP}
1510	LAI_FILE = lai2D.nc
1511
1512	# RENORM_LAI ([FLAG]) : flag to force LAI renormelization {LAI_MAP}
1513	RENORM_LAI = n
1514
1515	# VEGETATION_FILE ([FILE]) : Name of file from which the vegetation map is to be read {MAP_PFT_FORMAT}
1516	VEGETATION_FILE = PFTmap.nc
1517
1518	# VEGETATION_FILE ([FILE]) : Name of file from which the vegetation map is to be read {NOT(IMPOSE_VEG) and NOT(MAP_PFT_FORMAT)}
1519	VEGETATION_FILE = carteveg5km.nc
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	# SLOPE_NOREINF ([-]) : See slope_noreinf above {}
1525	SLOPE_NOREINF = 0.5
1526
1527	# TOPOGRAPHY_FILE ([FILE]) : Name of file from which the topography map is to be read {}
1528	TOPOGRAPHY_FILE = cartepente2d_15min.nc
1529
1530	# WOODHARVEST_FILE ([FILE]) : Name of file from which the wood harvest will be read {DO_WOOD_HARVEST}
1531	WOODHARVEST_FILE = woodharvest.nc
1532
1533	# SOILALB_FILE ([FILE]) : Name of file from which the bare soil albedo {NOT(IMPOSE_AZE)}
1534	SOILALB_FILE = soils_param.nc
1535
1536	# ALB_BG_FILE ([FILE]) : Name of file from which the background albedo is read {}
1537	ALB_BG_FILE = alb_bg.nc
1538
1539	# CDRAG_FROM_GCM ([FLAG]) : Keep cdrag coefficient from gcm. {OK_SECHIBA}
1540	CDRAG_FROM_GCM = y
1541
1542	# N_FERTIL_FILE (- ) : File name {CHEMISTRY_BVOC and NOx_FERTILIZERS_USE}
1543	N_FERTIL_FILE = orchidee_fertilizer_1995.nc
1544
1545	# N_FERTIL_FILE (-) : File name {CHEMISTRY_BVOC and NOx_FERTILIZERS_USE}
1546	N_FERTIL_FILE = orchidee_fertilizer_1995.nc
1547
1548	# ENERBIL_TSURF (Kelvin [K]) : Initial temperature if not found in restart {OK_SECHIBA}
1549	ENERBIL_TSURF = 280.
1550
1551	# ENERBIL_EVAPOT () : Initial Soil Potential Evaporation {OK_SECHIBA }
1552	ENERBIL_EVAPOT = 0.0
1553
1554	# BEDROCK_FLAG ([FLAG]) : Flag to consider bedrock at deeper layers. {}
1555	BEDROCK_FLAG = 0
1556
1557	# THERMOSOIL_TPRO (Kelvin [K]) : Initial soil temperature profile if not found in restart {OK_SECHIBA}
1558	THERMOSOIL_TPRO = 280.
1559
1560	# SOIL_REFTEMP_FILE ([FILE]) : File with climatological soil temperature {READ_REFTEMP}
1561	SOIL_REFTEMP_FILE = reftemp.nc
1562
1563	# HYDROL_OK_HDIFF ([FLAG]) : do horizontal diffusion? {OK_SECHIBA and .NOT.(HYDROL_CWRR) }
1564	HYDROL_OK_HDIFF = n
1565
1566	# HYDROL_SNOW ([kg/m^2]) : Initial snow mass if not found in restart {OK_SECHIBA}
1567	HYDROL_SNOW = 0.0
1568
1569	# HYDROL_SNOWAGE ([days]) : Initial snow age if not found in restart {OK_SECHIBA }
1570	HYDROL_SNOWAGE = 0.0
1571
1572	# HYDROL_SNOW_NOBIO ([m]) : Initial snow amount on ice, lakes, etc. if not found in restart {OK_SECHIBA }
1573	HYDROL_SNOW_NOBIO = 0.0
1574
1575	# HYDROL_SNOW_NOBIO_AGE ([days]) : Initial snow age on ice, lakes, etc. if not found in restart {OK_SECHIBA }
1576	HYDROL_SNOW_NOBIO_AGE = 0.0
1577
1578	# HYDROL_HUMR ([-]) : Initial soil moisture stress if not found in restart {OK_SECHIBA }
1579	HYDROL_HUMR = 1.0
1580
1581	# HYDROL_BQSB ([kg/m^2]) : Initial restart deep soil moisture if not found in restart {OK_SECHIBA }
1582	HYDROL_BQSB = 999999.
1583
1584	# HYDROL_GQSB ([kg/m^2]) : Initial upper soil moisture if not found in restart {OK_SECHIBA }
1585	HYDROL_GQSB = 0.0
1586
1587	# HYDROL_DSG ([m]) : Initial upper reservoir depth if not found in restart {OK_SECHIBA }
1588	HYDROL_DSG = 0.0
1589
1590	# HYDROL_QSV ([kg/m^2]) : Initial water on canopy if not found in restart {OK_SECHIBA }
1591	HYDROL_QSV = 0.0
1592
1593	# HYDROL_DSP ([m]) : Initial dry soil above upper reservoir if not found in restart {OK_SECHIBA }
1594	HYDROL_DSP = 999999.
1595
1596	# HYDROL_TAU_HDIFF ([seconds]) : time scale (s) for horizontal diffusion of water {HYDROL_OK_HDIFF}
1597	HYDROL_TAU_HDIFF = 86400.
1598
1599	# DO_PONDS ([FLAG]) : Should we include ponds {HYDROL_CWRR}
1600	DO_PONDS = n
1601
1602	# FROZ_FRAC_CORR ([-]) : Coefficient for the frozen fraction correction {HYDROL_CWRR and OK_FREEZE}
1603	FROZ_FRAC_CORR = 1.0
1604
1605	# MAX_FROZ_HYDRO ([-]) : Coefficient for the frozen fraction correction {HYDROL_CWRR and OK_FREEZE}
1606	MAX_FROZ_HYDRO = 1.0
1607
1608	# SMTOT_CORR ([-]) : Coefficient for the frozen fraction correction {HYDROL_CWRR and OK_FREEZE}
1609	SMTOT_CORR = 2.0
1610
1611	# DO_RSOIL ([FLAG]) : Should we reduce soil evaporation with a soil resistance {HYDROL_CWRR}
1612	DO_RSOIL = n
1613
1614	# OK_DYNROOT ([FLAG]) : Calculate dynamic root profile to optimize soil moisture usage {HYDROL_CWRR}
1615	OK_DYNROOT = n
1616
1617	# CWRR_N_VANGENUCHTEN ([-]) : Van genuchten coefficient n {HYDROL_CWRR}
1618	CWRR_N_VANGENUCHTEN = 1.89, 1.56, 1.31
1619
1620	# CWRR_A_VANGENUCHTEN ([1/mm] ) : Van genuchten coefficient a {HYDROL_CWRR}
1621	CWRR_A_VANGENUCHTEN = 0.0075, 0.0036, 0.0019
1622
1623	# VWC_RESIDUAL ([m3/m3] ) : Residual soil water content {HYDROL_CWRR}
1624	VWC_RESIDUAL = 0.065, 0.078, 0.095
1625
1626	# VWC_SAT ([m3/m3] ) : Saturated soil water content {HYDROL_CWRR}
1627	VWC_SAT = 0.41, 0.43, 0.41
1628
1629	# CWRR_KS ([mm/d] ) : Hydraulic conductivity Saturation {HYDROL_CWRR }
1630	CWRR_KS = 1060.8, 249.6, 62.4
1631
1632	# WETNESS_TRANSPIR_MAX ([-] ) : Soil moisture above which transpir is max {HYDROL_CWRR}
1633	WETNESS_TRANSPIR_MAX = 0.5, 0.5, 0.5
1634
1635	# VWC_FC ([m3/m3] ) : Volumetric water content field capacity {HYDROL_CWRR}
1636	VWC_FC = 0.32, 0.32, 0.32
1637
1638	# VWC_WP ([m3/m3] ) : Volumetric water content Wilting pt {HYDROL_CWRR}
1639	VWC_WP = 0.10, 0.10, 0.10
1640
1641	# VWC_MIN_FOR_WET_ALB ([m3/m3] ) : Vol. wat. cont. above which albedo is cst {HYDROL_CWRR}
1642	VWC_MIN_FOR_WET_ALB = 0.25, 0.25, 0.25
1643
1644	# VWC_MAX_FOR_DRY_ALB ([m3/m3] ) : Vol. wat. cont. below which albedo is cst {HYDROL_CWRR}
1645	VWC_MAX_FOR_DRY_ALB = 0.1, 0.1, 0.1
1646
1647	# HYDROL_MOISTURE_CONTENT ([m3/m3]) : Soil moisture on each soil tile and levels {HYDROL_CWRR }
1648	HYDROL_MOISTURE_CONTENT = 0.3
1649
1650	# US_INIT ([-]) : US_NVM_NSTM_NSLM {HYDROL_CWRR }
1651	US_INIT = 0.0
1652
1653	# ZWT_FORCE ([m]) : Prescribed water depth, dimension nstm {HYDROL_CWRR }
1654	ZWT_FORCE = -9999. -9999. -9999.
1655
1656	# FREE_DRAIN_COEF ([-]) : Coefficient for free drainage at bottom, dimension nstm {HYDROL_CWRR }
1657	FREE_DRAIN_COEF = 1.0 1.0 1.0
1658
1659	# WATER_TO_INFILT ([mm]) : Water to be infiltrated on top of the soil {HYDROL_CWRR }
1660	WATER_TO_INFILT = 0.0
1661
1662	# EVAPNU_SOIL ([mm]) : Bare soil evap on each soil if not found in restart {HYDROL_CWRR }
1663	EVAPNU_SOIL = 0.0
1664
1665	# HYDROL_SNOW () : Initial snow mass if not found in restart {OK_SECHIBA}
1666	HYDROL_SNOW = 0.0
1667
1668	# HYDROL_SNOWAGE (***) : Initial snow age if not found in restart {OK_SECHIBA}
1669	HYDROL_SNOWAGE = 0.0
1670
1671	# HYDROL_SNOW_NOBIO ([mm]) : Initial snow amount on ice, lakes, etc. if not found in restart {OK_SECHIBA}
1672	HYDROL_SNOW_NOBIO = 0.0
1673
1674	# HYDROL_SNOW_NOBIO_AGE (***) : Initial snow age on ice, lakes, etc. if not found in restart {OK_SECHIBA}
1675	HYDROL_SNOW_NOBIO_AGE = 0.0
1676
1677	# HYDROL_QSV ([mm]) : Initial water on canopy if not found in restart {OK_SECHIBA}
1678	HYDROL_QSV = 0.0
1679
1680	# CWRR_NKS_N0 ([-]) : fitted value for relation log((n-n0)/(n_ref-n0)) {HYDROL_CWRR }
1681	CWRR_NKS_N0 = 0.95
1682
1683	# CWRR_NKS_POWER ([-]) : fitted value for relation log((n-n0)/(n_ref-n0)) {HYDROL_CWRR }
1684	CWRR_NKS_POWER = 0.34
1685
1686	# CWRR_AKS_A0 ([1/mm]) : fitted value for relation log((a-a0)/(a_ref-a0)) {HYDROL_CWRR }
1687	CWRR_AKS_A0 = 0.00012
1688
1689	# CWRR_AKS_POWER ([-]) : fitted value for relation log((a-a0)/(a_ref-a0)) {HYDROL_CWRR }
1690	CWRR_AKS_POWER = 0.53
1691
1692	# KFACT_DECAY_RATE ([1/m]) : Factor for Ks decay with depth {HYDROL_CWRR }
1693	KFACT_DECAY_RATE = 2.0
1694
1695	# KFACT_STARTING_DEPTH ([m]) : Depth for compacted value of Ks {HYDROL_CWRR }
1696	KFACT_STARTING_DEPTH = 0.3
1697
1698	# KFACT_MAX ([-]) : Maximum Factor for Ks increase due to vegetation {HYDROL_CWRR }
1699	KFACT_MAX = 10.0
1700
1701	# DT_ROUTING ([seconds]) : Time step of the routing scheme {RIVER_ROUTING}
1702	DT_ROUTING = 86400.
1703
1704	# ROUTING_RIVERS ([-]) : Number of rivers {RIVER_ROUTING}
1705	ROUTING_RIVERS = 50
1706
1707	# DO_FLOODINFILT ([FLAG]) : Should floodplains reinfiltrate into the soil {RIVER_ROUTING}
1708	DO_FLOODINFILT = n
1709
1710	# DO_SWAMPS ([FLAG]) : Should we include swamp parameterization {RIVER_ROUTING}
1711	DO_SWAMPS = n
1712
1713	# DO_PONDS ([FLAG]) : Should we include ponds {RIVER_ROUTING}
1714	DO_PONDS = n
1715
1716	# SLOW_TCST ([days]) : Time constant for the slow reservoir {RIVER_ROUTING }
1717	SLOW_TCST = n
1718
1719	# FAST_TCST ([days]) : Time constant for the fast reservoir {RIVER_ROUTING }
1720	FAST_TCST = fast_tcst_cwrr or fast_tcst_chois depending on flag HYDROL_CWRR
1721
1722	# STREAM_TCST ([days]) : Time constant for the stream reservoir {RIVER_ROUTING}
1723	STREAM_TCST = stream_tcst_cwrr or stream_tcst_chois depending on flag HYDROL_CWRR
1724
1725	# FLOOD_TCST ([days]) : Time constant for the flood reservoir {RIVER_ROUTING}
1726	FLOOD_TCST = 4.0
1727
1728	# SWAMP_CST ([-]) : Fraction of the river that flows back to swamps {RIVER_ROUTING}
1729	SWAMP_CST = 0.2
1730
1731	# FLOOD_BETA ([-] ) : Parameter to fix the shape of the floodplain {RIVER_ROUTING}
1732	FLOOD_BETA = 2.0
1733
1734	# POND_BETAP ([-] ) : Ratio of the basin surface intercepted by ponds and the maximum surface of ponds {RIVER_ROUTING}
1735	POND_BETAP = 0.5
1736
1737	# FLOOD_CRI ([mm] ) : Potential height for which all the basin is flooded {DO_FLOODPLAINS or DO_PONDS}
1738	FLOOD_CRI = 2000.
1739
1740	# POND_CRI ([mm] ) : Potential height for which all the basin is a pond {DO_FLOODPLAINS or DO_PONDS}
1741	POND_CRI = 2000.
1742
1743	# MAX_LAKE_RESERVOIR ([kg/m2(routing area)] ) : Maximum limit of water in lake_reservoir {RIVER_ROUTING}
1744	MAX_LAKE_RESERVOIR = 7000
1745
1746	# RIVER_DESC ([FLAG]) : Writes out a description of the rivers {RIVER_ROUTING}
1747	RIVER_DESC = n
1748
1749	# 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}
1750	RIVER_DESC_FILE = river_desc.nc
1751
1752	# ROUTING_FILE ([FILE]) : Name of file which contains the routing information {RIVER_ROUTING}
1753	ROUTING_FILE = routing.nc
1754
1755	# IRRIGATION_FILE ([FILE]) : Name of file which contains the map of irrigated areas {DO_IRRIGATION OR DO_FLOODPLAINS}
1756	IRRIGATION_FILE = floodplains.nc
1757
1758	# EPS_CARBON ([%] ) : Allowed error on carbon stock {SPINUP_ANALYTIC}
1759	EPS_CARBON = 0.01
1760
1761	# SPINUP_PERIOD ([years] ) : Period to calulcate equilibrium during spinup analytic {SPINUP_ANALYTIC}
1762	SPINUP_PERIOD = -1
1763
1764	# STOMATE_FORCING_NAME ([FILE]) : Name of STOMATE's forcing file {OK_STOMATE}
1765	STOMATE_FORCING_NAME = NONE
1766
1767	# STOMATE_FORCING_MEMSIZE ([MegaBytes]) : Size of STOMATE forcing data in memory {OK_STOMATE}
1768	STOMATE_FORCING_MEMSIZE = 50
1769
1770	# STOMATE_CFORCING_NAME ([FILE]) : Name of STOMATE's carbon forcing file {OK_STOMATE}
1771	STOMATE_CFORCING_NAME = NONE
1772
1773	# FORCESOIL_STEP_PER_YEAR ([days, months, year]) : Number of time steps per year for carbon spinup. {OK_STOMATE}
1774	FORCESOIL_STEP_PER_YEAR = 365
1775
1776	# FORCESOIL_NB_YEAR ([years]) : Number of years saved for carbon spinup. {OK_STOMATE}
1777	FORCESOIL_NB_YEAR = 1
1778
1779	# XIOS_ORCHIDEE_OK ([FLAG]) : Use XIOS for writing diagnostics file {}
1780	XIOS_ORCHIDEE_OK = y
1781

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

Download in other formats:

Original Format