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

Last change on this file since 7444 was 7444, checked in by agnes.ducharne, 5 months ago

Finishes r7443 for WETNESS_TRANSPIR_MAX.

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