source: branches/publications/ORCHIDEE_gmd-2018-261/orchidee.default @ 6145

Last change on this file since 6145 was 3714, checked in by nicolas.vuichard, 8 years ago

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