Context Navigation

orchidee.default @ 7432

Last change on this file since 7432 was 7432, checked in by agnes.ducharne, 2 years ago

Changes to make the IMPOSE_SOILT mode functional. This mode is not anymore dependent on IMPOSE_VEG, so we can impose soil
properties even when we read a vegetation map. A new output variable is added to export "ksref" before being vertically modified.

With IMPOSE_SOILT, we can either impose a certain USDA texture (via SOIL_FRACTIONS) and this propagates to all the soil parameters, but we can as well impose specific parameters in isolation (with a default txture as Loam). This committ was tested with various configs, and it also works if we use a restart file, and we can even change the soil texture map in this case.

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

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

Download in other formats:

Original Format