Context Navigation

orchidee.default @ 6412

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

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

Download in other formats:

Original Format