source: branches/ORCHIDEE_2_2/ORCHIDEE/orchidee.default @ 7475

Last change on this file since 7475 was 7475, checked in by josefine.ghattas, 2 years ago

Update orchidee.default using the tool create_orchidee_defalut.sh tool. Note that manual modifications in orchidee.default are overwritten. Reported some of the changes done previously in orchidee.default to the fortran source code. Note also that all sections in !Config must be there to have the tool working. Therfore added some missing !Config If lines.

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