source: branches/publications/ORCHIDEE-LEAK-r5919/orchidee.default @ 5925

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