source: branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes.f90 @ 7432

Last change on this file since 7432 was 7432, checked in by agnes.ducharne, 2 years ago

Changes to make the IMPOSE_SOILT mode functional. This mode is not anymore dependent on IMPOSE_VEG, so we can impose soil
properties even when we read a vegetation map. A new output variable is added to export "ksref" before being vertically modified.

With IMPOSE_SOILT, we can either impose a certain USDA texture (via SOIL_FRACTIONS) and this propagates to all the soil parameters, but we can as well impose specific parameters in isolation (with a default txture as Loam). This committ was tested with various configs, and it also works if we use a restart file, and we can even change the soil texture map in this case.

  • Property svn:keywords set to Date Revision
File size: 72.8 KB
Line 
1! =================================================================================================================================
2! MODULE       : constantes
3!
4! CONTACT      : orchidee-help _at_ listes.ipsl.fr
5!
6! LICENCE      : IPSL (2006)
7! This software is governed by the CeCILL licence see ORCHIDEE/ORCHIDEE_CeCILL.LIC
8!
9!>\BRIEF        "constantes" module contains subroutines to initialize most of the exernalized parameters. This module
10!!              also make a use to the module constantes_var where the parameters are declared.
11!!
12!!\n DESCRIPTION: This module contains subroutines to initialize most of the exernalized parameters. This module
13!!                also make a use to the module constantes_var where the parameters are declared.\n
14!!                This module can be used to acces the subroutines and the constantes. The constantes declarations
15!!                can also be used seperatly with "USE constantes_var".
16!!
17!! RECENT CHANGE(S): Didier Solyga : This module contains now all the externalized parameters of ORCHIDEE
18!!                   listed by modules which are not pft-dependent 
19!!                   Josefine Ghattas 2013 : The declaration part has been extracted and moved to module constates_var
20!!
21!! REFERENCE(S) :
22!! - Louis, Jean-Francois (1979), A parametric model of vertical eddy fluxes in the atmosphere.
23!! Boundary Layer Meteorology, 187-202.
24!!
25!! SVN          :
26!! $HeadURL: $
27!! $Date$
28!! $Revision$
29!! \n
30!_ ================================================================================================================================
31
32MODULE constantes
33
34  USE constantes_var
35  USE defprec
36  USE ioipsl_para, ONLY : getin_p, ipslerr_p
37  USE mod_orchidee_para_var, ONLY : numout
38  USE time, ONLY : one_day, dt_sechiba
39
40  IMPLICIT NONE
41 
42CONTAINS
43
44
45!! ================================================================================================================================
46!! SUBROUTINE   : activate_sub_models
47!!
48!>\BRIEF         This subroutine reads the flags in the configuration file to
49!! activate some sub-models like routing, irrigation, fire, herbivory, ... 
50!!
51!! DESCRIPTION  : None
52!!
53!! RECENT CHANGE(S): None
54!!
55!! MAIN OUTPUT VARIABLE(S): None
56!!
57!! REFERENCE(S) : None
58!!
59!! FLOWCHART    : None
60!! \n
61!_ ================================================================================================================================
62
63  SUBROUTINE activate_sub_models()
64
65    IMPLICIT NONE
66
67    !! 0. Variables and parameters declaration
68
69    !! 0.4 Local variables
70
71    !_ ================================================================================================================================
72
73    IF (ok_stomate) THEN
74
75       !Config Key   = HERBIVORES
76       !Config Desc  = herbivores allowed?
77       !Config If    = OK_STOMATE
78       !Config Def   = n
79       !Config Help  = With this variable, you can determine
80       !Config         if herbivores are activated
81       !Config Units = [FLAG]
82       CALL getin_p('HERBIVORES', ok_herbivores)
83       !
84       !Config Key   = TREAT_EXPANSION
85       !Config Desc  = treat expansion of PFTs across a grid cell?
86       !Config If    = OK_STOMATE
87       !Config Def   = n
88       !Config Help  = With this variable, you can determine
89       !Config         whether we treat expansion of PFTs across a
90       !Config         grid cell.
91       !Config Units = [FLAG]
92       CALL getin_p('TREAT_EXPANSION', treat_expansion)
93
94       !Config Key   = LPJ_GAP_CONST_MORT
95       !Config Desc  = Constant mortality
96       !Config If    = OK_STOMATE AND NOT OK_DGVM
97       !Config Def   = y/n depending on OK_DGVM
98       !Config Help  = set to TRUE if constant mortality is to be activated
99       !Config         
100       !Config Units = [FLAG]
101
102       ! Set Default value different if DGVM is activated.
103       IF ( ok_dgvm ) THEN
104          lpj_gap_const_mort=.FALSE.
105       ELSE
106          lpj_gap_const_mort=.TRUE.
107       END IF
108       CALL getin_p('LPJ_GAP_CONST_MORT', lpj_gap_const_mort)
109
110       IF (ok_dgvm .AND. lpj_gap_const_mort) THEN
111          CALL ipslerr_p(1,"activate_sub_models","Both OK_DGVM and LPJ_GAP_CONST_MORT are activated.",&
112               "This combination is possible but unusual","The simulation will continue with these flags activated." )
113       ELSEIF (.NOT. ok_dgvm  .AND. .NOT. lpj_gap_const_mort) THEN
114           CALL ipslerr_p(3,"activate_sub_models", &
115                "The combination of OK_DGVM=false and LPJ_GAP_CONST_MORT=false is not operational in this version", &
116                "Some parts of the code should first be revised.","" )
117       END IF
118
119       !Config Key   = HARVEST_AGRI
120       !Config Desc  = Harvest model for agricultural PFTs.
121       !Config If    = OK_STOMATE
122       !Config Def   = y
123       !Config Help  = Compute harvest above ground biomass for agriculture.
124       !Config         Change daily turnover.
125       !Config Units = [FLAG]
126       CALL getin_p('HARVEST_AGRI', harvest_agri)
127       !
128       !Config Key   = FIRE_DISABLE
129       !Config Desc  = no fire allowed
130       !Config If    = OK_STOMATE
131       !Config Def   = y
132       !Config Help  = With this variable, you can allow or not
133       !Config         the estimation of CO2 lost by fire
134       !Config Units = [FLAG]
135       CALL getin_p('FIRE_DISABLE', disable_fire)
136       !
137       !Config Key   = SPINUP_ANALYTIC
138       !Config Desc  = Activation of the analytic resolution of the spinup.
139       !Config If    = OK_STOMATE
140       !Config Def   = n
141       !Config Help  = Activate this option if you want to solve the spinup by the Gauss-Jordan method.
142       !Config Units = BOOLEAN   
143       CALL getin_p('SPINUP_ANALYTIC',spinup_analytic)
144
145    ENDIF
146
147    !
148    ! Check consistency (see later)
149    !
150!!$        IF(.NOT.(ok_routing) .AND. (doirrigation .OR. dofloodplains)) THEN
151!!$           CALL ipslerr_p(2,'activate_sub_models', &
152!!$               &     'Problem :you tried to activate the irrigation and floodplains without activating the routing',&
153!!$               &     'Are you sure ?', &
154!!$               &     '(check your parameters).')
155!!$        ENDIF
156
157!!$        IF(.NOT.(ok_stomate) .AND. (ok_herbivores .OR. treat_expansion .OR. lpj_gap_const_mort &
158!!$            & .OR. harvest_agri .OR. disable_fire)) THEN
159!!$          CALL ipslerr_p(2,'activate_sub_models', &
160!!$               &     'Problem : try to activate the following options : herbivory, treat_expansion, fire,',&
161!!$               &     'harvest_agri and constant mortality without stomate activated.',&
162!!$               &     '(check your parameters).')
163!!$        ENDIF
164
165
166  END SUBROUTINE activate_sub_models
167
168!! ================================================================================================================================
169!! SUBROUTINE   : veget_config
170!!
171!>\BRIEF         This subroutine reads the flags controlling the configuration for
172!! the vegetation : impose_veg, veget_mpa, lai_map, etc...       
173!!
174!! DESCRIPTION  : None
175!!
176!! RECENT CHANGE(S): None
177!!
178!! MAIN OUTPUT VARIABLE(S):
179!!
180!! REFERENCE(S) :
181!!
182!! FLOWCHART    :
183!! \n
184!_ ================================================================================================================================
185
186  SUBROUTINE veget_config
187
188    IMPLICIT NONE
189
190    !! 0. Variables and parameters declaration
191
192    !! 0.4 Local variables 
193
194    !_ ================================================================================================================================
195
196    !Config Key   = AGRICULTURE
197    !Config Desc  = agriculture allowed?
198    !Config If    = OK_SECHIBA or OK_STOMATE
199    !Config Def   = y
200    !Config Help  = With this variable, you can determine
201    !Config         whether agriculture is allowed
202    !Config Units = [FLAG]
203    CALL getin_p('AGRICULTURE', agriculture)
204    !
205    !Config Key   = IMPOSE_VEG
206    !Config Desc  = Should the vegetation be prescribed ?
207    !Config If    = OK_SECHIBA or OK_STOMATE
208    !Config Def   = n
209    !Config Help  = This flag allows the user to impose a vegetation distribution
210    !Config         and its characteristics. It is espacially interesting for 0D
211    !Config         simulations. On the globe it does not make too much sense as
212    !Config         it imposes the same vegetation everywhere
213    !Config Units = [FLAG]
214    CALL getin_p('IMPOSE_VEG', impveg)
215
216    !IF (impveg) THEN
217       !Config Key   = IMPOSE_SOILT
218       !Config Desc  = Should the soil type be prescribed ?
219       !Config Def   = n
220       !Config If    = IMPOSE_VEG
221       !Config Help  = This flag allows the user to impose a soil type distribution.
222       !Config         It is espacially interesting for 0D
223       !Config         simulations. On the globe it does not make too much sense as
224       !Config         it imposes the same soil everywhere
225       !Config Units = [FLAG]
226       CALL getin_p('IMPOSE_SOILT', impsoilt)     
227    !ENDIF
228
229    !Config Key   = LAI_MAP
230    !Config Desc  = Read the LAI map
231    !Config If    = OK_SECHIBA or OK_STOMATE
232    !Config Def   = n
233    !Config Help  = It is possible to read a 12 month LAI map which will
234    !Config         then be interpolated to daily values as needed.
235    !Config Units = [FLAG]
236    CALL getin_p('LAI_MAP',read_lai)
237
238    !Config Key   = VEGET_REINIT
239    !Config Desc  = Reset veget_year counter (obsolet)
240    !Config If    = VEGET_UPDATE > 0Y
241    !Config Def   = y
242    !Config Help  = The parameter is used to bypass veget_year count
243    !Config         and reinitialize it with VEGET_YEAR parameter.
244    !Config         Then it is possible to change LAND USE file.
245    !Config Units = [FLAG]
246    CALL getin_p('VEGET_REINIT', veget_reinit)
247   
248    !Config Key   = VEGETMAP_RESET
249    !Config Desc  = Flag to change vegetation map without activating LAND USE change for carbon fluxes and reset carbon related variables to zero
250    !Config If    =
251    !Config Def   = n
252    !Config Help  = Use this option to change vegetation map while keeping VEGET_UPDATE=0Y
253    !Config Units = [FLAG]
254    CALL getin_p('VEGETMAP_RESET', vegetmap_reset)
255
256
257    !Config Key   = VEGET_YEAR
258    !Config Desc  = Year of the vegetation map to be read
259    !Config If    =
260    !Config Def   = 1
261    !Config Help  = First year for land use vegetation (2D map by pft).
262    !Config         If VEGET_YEAR is set to 0, this means there is no time axis.
263    !Config Units = [FLAG]
264    CALL getin_p('VEGET_YEAR', veget_year_orig)
265   
266
267!!$        ! DS : Add warning in case of a wrong configuration (need to be discussed)
268!!$        ! 4.
269!!$        IF ( .NOT.(impveg) .AND. impsoilt) THEN
270!!$           CALL ipslerr_p(2,'veget_config', &
271!!$               &     'Problem : try to activate impose_soilt without activating impose_veg.',&
272!!$               &     'Are you sure ?', &
273!!$               &     '(check your parameters).')
274!!$        ENDIF
275!!$
276
277  END SUBROUTINE veget_config
278
279
280!! ================================================================================================================================
281!! SUBROUTINE   : veget_config
282!!
283!>\BRIEF         This subroutine reads in the configuration file the imposed values of the parameters for all SECHIBA modules. 
284!!
285!! DESCRIPTION  : None
286!!
287!! RECENT CHANGE(S): None
288!!
289!! MAIN OUTPUT VARIABLE(S):
290!!
291!! REFERENCE(S) :
292!!
293!! FLOWCHART    :
294!! \n
295!_ ================================================================================================================================
296
297  SUBROUTINE config_sechiba_parameters
298
299    IMPLICIT NONE
300
301    !! 0. Variables and parameters declaration
302
303    !! 0.4 Local variables
304    REAL(r_std) :: nudge_tau_mc     !! Temporary variable read from run.def
305    REAL(r_std) :: nudge_tau_snow   !! Temporary variable read from run.def
306
307    !_ ================================================================================================================================
308
309    ! Global : parameters used by many modules
310    CALL getin_p('TESTPFT',testpft)
311
312    !
313    !Config Key   = MAXMASS_SNOW
314    !Config Desc  = The maximum mass of a snow
315    !Config If    = OK_SECHIBA
316    !Config Def   = 3000.
317    !Config Help  =
318    !Config Units = [kg/m^2] 
319    CALL getin_p('MAXMASS_SNOW',maxmass_snow)
320    !
321    !Config Key   = SNOWCRI
322    !Config Desc  = Sets the amount above which only sublimation occures
323    !Config If    = OK_SECHIBA
324    !Config Def   = 1.5
325    !Config Help  =
326    !Config Units = [kg/m^2] 
327    CALL getin_p('SNOWCRI',snowcri)
328    !
329    !! Initialization of sneige
330    sneige = snowcri/mille
331    !
332    !Config Key   = MIN_WIND
333    !Config Desc  = Minimum wind speed
334    !Config If    = OK_SECHIBA
335    !Config Def   = 0.1
336    !Config Help  =
337    !Config Units = [m/s]
338    CALL getin_p('MIN_WIND',min_wind)
339    !
340    !Config Key   = MAX_SNOW_AGE
341    !Config Desc  = Maximum period of snow aging
342    !Config If    = OK_SECHIBA
343    !Config Def   = 50.
344    !Config Help  =
345    !Config Units = [days?]
346    CALL getin_p('MAX_SNOW_AGE',max_snow_age)
347    !
348    !Config Key   = SNOW_TRANS
349    !Config Desc  = Transformation time constant for snow
350    !Config If    = OK_SECHIBA
351    !Config Def   = 0.2
352    !Config Help  = optimized on 04/07/2016
353    !Config Units = [m]   
354    CALL getin_p('SNOW_TRANS',snow_trans)
355
356   
357    !Config Key   = OK_NUDGE_MC
358    !Config Desc  = Activate nudging of soil moisture
359    !Config Def   = n
360    !Config If    =
361    !Config Help  =
362    !Config Units = [FLAG]
363    ok_nudge_mc = .FALSE.
364    CALL getin_p('OK_NUDGE_MC', ok_nudge_mc)
365
366    !Config Key   = NUDGE_TAU_MC
367    !Config Desc  = Relaxation time for nudging of soil moisture expressed in fraction of the day
368    !Config Def   = 1
369    !Config If    = OK_NUDGE_MC
370    !Config Help  =
371    !Config Units = [-]
372    nudge_tau_mc = 1.0
373    CALL getin_p('NUDGE_TAU_MC', nudge_tau_mc)
374    IF (nudge_tau_mc < dt_sechiba/one_day) CALL ipslerr_p(3, 'hydrol_initialize', &
375         'NUDGE_TAU_MC is smaller than the time step in sechiba which is not allowed.', &
376         'Set NUDGE_TAU_MC higher or equal to dt_sechiba/one_day','')
377    ! Calculate alpha to be used in hydrol
378    alpha_nudge_mc = dt_sechiba/(one_day*nudge_tau_mc)
379    IF (printlev>=2) WRITE(numout, *) 'ok_nudge_mc, nudge_tau_mc, alpha_nudge_mc =', &
380         ok_nudge_mc, nudge_tau_mc, alpha_nudge_mc
381
382
383    !Config Key   = OK_NUDGE_SNOW
384    !Config Desc  = Activate nudging of snow variables
385    !Config Def   = n
386    !Config If    =
387    !Config Help  =
388    !Config Units = [FLAG]
389    ok_nudge_snow = .FALSE.
390    CALL getin_p('OK_NUDGE_SNOW', ok_nudge_snow)
391
392    !Config Key   = NUDGE_TAU_SNOW
393    !Config Desc  = Relaxation time for nudging of snow variables
394    !Config Def   = 1
395    !Config If    = OK_NUDGE_SNOW
396    !Config Help  =
397    !Config Units = [-]
398    nudge_tau_snow = 1.0
399    CALL getin_p('NUDGE_TAU_SNOW', nudge_tau_snow)
400    IF (nudge_tau_snow < dt_sechiba/one_day) CALL ipslerr_p(3, 'hydrol_initialize', &
401         'NUDGE_TAU_SNOW is smaller than the time step in sechiba which is not allowed.', &
402         'Set NUDGE_TAU_SNOW higher or equal to dt_sechiba/one_day','')
403    ! Calculate alpha to be used in hydrol
404    alpha_nudge_snow = dt_sechiba/(one_day*nudge_tau_snow)
405    IF (printlev>=2) WRITE(numout, *) 'ok_nudge_snow, nudge_tau_snow, alpha_nudge_snow =', &
406         ok_nudge_snow, nudge_tau_snow, alpha_nudge_snow
407
408
409    !Config Key   = NUDGE_INTERPOL_WITH_XIOS
410    !Config Desc  = Activate reading and interpolation with XIOS for nudging fields
411    !Config Def   = n
412    !Config If    = OK_NUDGE_MC or OK_NUDGE_SNOW
413    !Config Help  =
414    !Config Units = [FLAG]
415    nudge_interpol_with_xios = .FALSE.
416    CALL getin_p('NUDGE_INTERPOL_WITH_XIOS', nudge_interpol_with_xios)
417
418    !-
419    ! condveg
420    !-
421    !
422    !Config Key   = HEIGHT_DISPLACEMENT
423    !Config Desc  = Magic number which relates the height to the displacement height.
424    !Config If    = OK_SECHIBA
425    !Config Def   = 0.75
426    !Config Help  =
427    !Config Units = [m] 
428    CALL getin_p('HEIGHT_DISPLACEMENT',height_displacement)
429    !
430    !Config Key   = Z0_BARE
431    !Config Desc  = bare soil roughness length
432    !Config If    = OK_SECHIBA
433    !Config Def   = 0.01
434    !Config Help  =
435    !Config Units = [m]   
436    CALL getin_p('Z0_BARE',z0_bare)
437    !
438    !Config Key   = Z0_ICE
439    !Config Desc  = ice roughness length
440    !Config If    = OK_SECHIBA
441    !Config Def   = 0.001
442    !Config Help  =
443    !Config Units = [m]   
444    CALL getin_p('Z0_ICE',z0_ice)
445    !
446    !Config Key   = TCST_SNOWA
447    !Config Desc  = Time constant of the albedo decay of snow
448    !Config If    = OK_SECHIBA
449    !Config Def   = 10.0
450    !Config Help  = optimized on 04/07/2016
451    !Config Units = [days]
452    CALL getin_p('TCST_SNOWA',tcst_snowa)
453    !
454    !Config Key   = SNOWCRI_ALB
455    !Config Desc  = Critical value for computation of snow albedo
456    !Config If    = OK_SECHIBA
457    !Config Def   = 10.
458    !Config Help  =
459    !Config Units = [cm] 
460    CALL getin_p('SNOWCRI_ALB',snowcri_alb)
461    !
462    !
463    !Config Key   = VIS_DRY
464    !Config Desc  = The correspondance table for the soil color numbers and their albedo
465    !Config If    = OK_SECHIBA
466    !Config Def   = 0.24, 0.22, 0.20, 0.18, 0.16, 0.14, 0.12, 0.10, 0.27
467    !Config Help  =
468    !Config Units = [-] 
469    CALL getin_p('VIS_DRY',vis_dry)
470    !
471    !Config Key   = NIR_DRY
472    !Config Desc  = The correspondance table for the soil color numbers and their albedo
473    !Config If    = OK_SECHIBA
474    !Config Def   = 0.48, 0.44, 0.40, 0.36, 0.32, 0.28, 0.24, 0.20, 0.55
475    !Config Help  =
476    !Config Units = [-]   
477    CALL getin_p('NIR_DRY',nir_dry)
478    !
479    !Config Key   = VIS_WET
480    !Config Desc  = The correspondance table for the soil color numbers and their albedo
481    !Config If    = OK_SECHIBA 
482    !Config Def   = 0.12, 0.11, 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.15
483    !Config Help  =
484    !Config Units = [-]   
485    CALL getin_p('VIS_WET',vis_wet)
486    !
487    !Config Key   = NIR_WET
488    !Config Desc  = The correspondance table for the soil color numbers and their albedo
489    !Config If    = OK_SECHIBA
490    !Config Def   = 0.24, 0.22, 0.20, 0.18, 0.16, 0.14, 0.12, 0.10, 0.31
491    !Config Help  =
492    !Config Units = [-]   
493    CALL getin_p('NIR_WET',nir_wet)
494    !
495    !Config Key   = ALBSOIL_VIS
496    !Config Desc  =
497    !Config If    = OK_SECHIBA
498    !Config Def   = 0.18, 0.16, 0.16, 0.15, 0.12, 0.105, 0.09, 0.075, 0.25
499    !Config Help  =
500    !Config Units = [-] 
501    CALL getin_p('ALBSOIL_VIS',albsoil_vis)
502    !
503    !Config Key   = ALBSOIL_NIR
504    !Config Desc  =
505    !Config If    = OK_SECHIBA
506    !Config Def   = 0.36, 0.34, 0.34, 0.33, 0.30, 0.25, 0.20, 0.15, 0.45
507    !Config Help  =
508    !Config Units = [-] 
509    CALL getin_p('ALBSOIL_NIR',albsoil_nir)
510    !-
511    !
512    !Config Key   = ALB_DEADLEAF
513    !Config Desc  = albedo of dead leaves, VIS+NIR
514    !Config If    = OK_SECHIBA
515    !Config Def   = 0.12, 0.35
516    !Config Help  =
517    !Config Units = [-]     
518    CALL getin_p('ALB_DEADLEAF',alb_deadleaf)
519    !
520    !Config Key   = ALB_ICE
521    !Config Desc  = albedo of ice, VIS+NIR
522    !Config If    = OK_SECHIBA
523    !Config Def   = 0.60, 0.20
524    !Config Help  =
525    !Config Units = [-] 
526    CALL getin_p('ALB_ICE',alb_ice)
527    !
528    ! Get the fixed snow albedo if needed
529    !
530    !Config Key   = CONDVEG_SNOWA
531    !Config Desc  = The snow albedo used by SECHIBA
532    !Config Def   = 1.E+20
533    !Config if    = OK_SECHIBA
534    !Config Help  = This option allows the user to impose a snow albedo.
535    !Config         Default behaviour is to use the model of snow albedo
536    !Config         developed by Chalita (1993).
537    !Config Units = [-]
538    CALL getin_p('CONDVEG_SNOWA',fixed_snow_albedo)
539    !
540    !Config Key   = ALB_BARE_MODEL
541    !Config Desc  = Switch bare soil albedo dependent (if TRUE) on soil wetness
542    !Config Def   = n
543    !Config if    = OK_SECHIBA
544    !Config Help  = If TRUE, the model for bare soil albedo is the old formulation.
545    !Config         Then it depend on the soil dry or wetness. If FALSE, it is the
546    !Config         new computation that is taken, it is the mean of soil albedo.
547    !Config Units = [FLAG]
548    CALL getin_p('ALB_BARE_MODEL',alb_bare_model)
549    !
550    !Config Key   = ALB_BG_MODIS
551    !Config Desc  = Read bare soil albedo from file with background MODIS data
552    !Config Def   = y
553    !Config if    = OK_SECHIBA
554    !Config Help  = If TRUE, the bare soil albedo is read from file
555    !Config         based on background MODIS data. 
556    !Config         If FALSE, computaion depends on ALB_BARE_MODEL
557    !Config Units = [FLAG]
558    CALL getin_p('ALB_BG_MODIS',alb_bg_modis)
559    !
560    !Config Key   = IMPOSE_AZE
561    !Config Desc  = Should the surface parameters be prescribed
562    !Config Def   = n
563    !Config if    = OK_SECHIBA
564    !Config Help  = This flag allows the user to impose the surface parameters
565    !Config         (Albedo Roughness and Emissivity). It is espacially interesting for 0D
566    !Config         simulations. On the globe it does not make too much sense as
567    !Config         it imposes the same vegetation everywhere
568    !Config Units = [FLAG]
569    CALL getin_p('IMPOSE_AZE',impaze)
570    !
571    IF(impaze) THEN
572       !
573       !Config Key   = CONDVEG_Z0
574       !Config Desc  = Surface roughness
575       !Config Def   = 0.15
576       !Config If    = IMPOSE_AZE
577       !Config Help  = Surface rougness to be used on the point if a 0-dim version
578       !Config         of SECHIBA is used. Look at the description of the forcing 
579       !Config         data for the correct value.
580       !Config Units = [m]
581       CALL getin_p('CONDVEG_Z0', z0_scal) 
582       !
583       !Config Key   = ROUGHHEIGHT
584       !Config Desc  = Height to be added to the height of the first level
585       !Config Def   = 0.0
586       !Config If    = IMPOSE_AZE
587       !Config Help  = ORCHIDEE assumes that the atmospheric level height is counted
588       !Config         from the zero wind level. Thus to take into account the roughness
589       !Config         of tall vegetation we need to correct this by a certain fraction
590       !Config         of the vegetation height. This is called the roughness height in
591       !Config         ORCHIDEE talk.
592       !Config Units = [m]
593       CALL getin_p('ROUGHHEIGHT', roughheight_scal)
594       !
595       !Config Key   = CONDVEG_ALBVIS
596       !Config Desc  = SW visible albedo for the surface
597       !Config Def   = 0.25
598       !Config If    = IMPOSE_AZE
599       !Config Help  = Surface albedo in visible wavelengths to be used
600       !Config         on the point if a 0-dim version of SECHIBA is used.
601       !Config         Look at the description of the forcing data for
602       !Config         the correct value.
603       !Config Units = [-]
604       CALL getin_p('CONDVEG_ALBVIS', albedo_scal(ivis))
605       !
606       !Config Key   = CONDVEG_ALBNIR
607       !Config Desc  = SW near infrared albedo for the surface
608       !Config Def   = 0.25
609       !Config If    = IMPOSE_AZE
610       !Config Help  = Surface albedo in near infrared wavelengths to be used
611       !Config         on the point if a 0-dim version of SECHIBA is used.
612       !Config         Look at the description of the forcing data for
613       !Config         the correct value.
614       !Config Units = [-] 
615       CALL getin_p('CONDVEG_ALBNIR', albedo_scal(inir))
616       !
617       !Config Key   = CONDVEG_EMIS
618       !Config Desc  = Emissivity of the surface for LW radiation
619       !Config Def   = 1.0
620       !Config If    = IMPOSE_AZE
621       !Config Help  = The surface emissivity used for compution the LE emission
622       !Config         of the surface in a 0-dim version. Values range between
623       !Config         0.97 and 1.. The GCM uses 0.98.
624       !Config Units = [-]
625       CALL getin_p('CONDVEG_EMIS', emis_scal)
626    ENDIF
627
628    CALL getin_p('NEW_WATSTRESS',new_watstress)
629    IF(new_watstress) THEN
630       CALL getin_p('ALPHA_WATSTRESS',alpha_watstress)
631    ENDIF
632
633    !
634    !Config Key   = ROUGH_DYN
635    !Config Desc  = Account for a dynamic roughness height
636    !Config Def   = y
637    !Config if    = OK_SECHIBA
638    !Config Help  = If this flag is set to true (y) then the roughness
639    !Config         height is computed dynamically, varying with LAI
640    !Config Units = [FLAG]
641    CALL getin_p('ROUGH_DYN',rough_dyn)
642
643    IF ( rough_dyn ) THEN
644       !
645       !Config Key   = C1
646       !Config Desc  = Constant used in the formulation of the ratio of
647       !Config         the ratio of friction velocity to the wind speed
648       !Config         at the canopy top
649       !Config         See Ershadi et al. (2015) for more info
650       !Config Def   = 0.32
651       !Config If    = ROUGH_DYN
652       !Config Help  =
653       !Config Units = [-]
654       CALL getin_p('C1', c1)
655       !
656       !Config Key   = C2
657       !Config Desc  = Constant used in the formulation of the ratio of
658       !Config         the ratio of friction velocity to the wind speed
659       !Config         at the canopy top
660       !Config         See Ershadi et al. (2015) for more info
661       !Config Def   = 0.264
662       !Config If    = ROUGH_DYN
663       !Config Help  =
664       !Config Units = [-]
665       CALL getin_p('C2', c2)
666       !
667       !Config Key   = C3
668       !Config Desc  = Constant used in the formulation of the ratio of
669       !Config         the ratio of friction velocity to the wind speed
670       !Config         at the canopy top
671       !Config         See Ershadi et al. (2015) for more info
672       !Config Def   = 15.1
673       !Config If    = ROUGH_DYN
674       !Config Help  =
675       !Config Units = [-]
676       CALL getin_p('C3', c3)
677       !
678       !Config Key   = Cdrag_foliage
679       !Config Desc  = Drag coefficient of the foliage
680       !Config         See Ershadi et al. (2015) and Su et al. (2001)
681       !Config         for more info
682       !Config Def   = 0.2
683       !Config If    = ROUGH_DYN
684       !Config Help  =
685       !Config Units = [-]
686       CALL getin_p('CDRAG_FOLIAGE', Cdrag_foliage)
687       !
688       !Config Key   = Ct
689       !Config Desc  = Heat transfer coefficient of the leaf
690       !Config         See Ershadi et al. (2015) and Su et al. (2001)
691       !Config         for more info
692       !Config Def   = 0.01
693       !Config If    = ROUGH_DYN
694       !Config Help  =
695       !Config Units = [-]
696       CALL getin_p('CT', Ct)
697       !
698       !Config Key   = Prandtl
699       !Config Desc  = Prandtl number used in the calculation of Ct*
700       !Config         See Su et al. (2001) for more info
701       !Config Def   = 0.71
702       !Config If    = ROUGH_DYN
703       !Config Help  =
704       !Config Units = [-]
705       CALL getin_p('PRANDTL', Prandtl)
706    ENDIF
707    !-
708    ! Variables related to the explicitsnow module
709    !-
710    !Config Key = xansmax
711    !Config Desc = maximum snow albedo
712    !Config If = OK_SECHIBA
713    !Config Def = 0.85
714    !Config Help =
715    !Config Units = [-]
716    CALL getin_p('XANSMAX',xansmax)
717    !
718    !Config Key = xansmin
719    !Config Desc = minimum snow albedo
720    !Config If = OK_SECHIBA
721    !Config Def = 0.50
722    !Config Help =
723    !Config Units = [-]
724    CALL getin_p('XANSMIN',xansmin)
725    !
726    !Config Key = xans_todry
727    !Config Desc = albedo decay rate for the dry snow
728    !Config If = OK_SECHIBA
729    !Config Def = 0.008
730    !Config Help =
731    !Config Units = [S-1]
732    CALL getin_p('XANSDRY',xans_todry)
733    !
734    !Config Key = xans_t
735    !Config Desc = albedo decay rate for the wet snow
736    !Config If = OK_SECHIBA
737    !Config Def = 0.24
738    !Config Help =
739    !Config Units = [S-1]
740    CALL getin_p('XANS_T',xans_t)
741
742    !Config Key = xrhosmax
743    !Config Desc = maximum snow density
744    !Config If = OK_SECHIBA
745    !Config Def = 750
746    !Config Help =
747    !Config Units = [-]
748    CALL getin_p('XRHOSMAX',xrhosmax)
749    !
750    !Config Key = xwsnowholdmax1
751    !Config Desc = snow holding capacity 1
752    !Config If = OK_SECHIBA
753    !Config Def = 0.03
754    !Config Help =
755    !Config Units = [-]
756    CALL getin_p('XWSNOWHOLDMAX1',xwsnowholdmax1)
757    !
758    !Config Key = xwsnowholdmax2
759    !Config Desc = snow holding capacity 2
760    !Config If = OK_SECHIBA
761    !Config Def = 0.10
762    !Config Help =
763    !Config Units = [-]
764    CALL getin_p('XWSNOWHOLDMAX2',xwsnowholdmax2)
765    !
766    !Config Key = xsnowrhohold
767    !Config Desc = snow density
768    !Config If = OK_SECHIBA
769    !Config Def = 200.0
770    !Config Help =
771    !Config Units = [kg/m3]
772    CALL getin_p('XSNOWRHOHOLD',xsnowrhohold)
773    !
774    !Config Key = ZSNOWTHRMCOND1
775    !Config Desc = Thermal conductivity Coef 1
776    !Config If = OK_SECHIBA
777    !Config Def = 0.02
778    !Config Help =
779    !Config Units = [W/m/K]
780    CALL getin_p('ZSNOWTHRMCOND1',ZSNOWTHRMCOND1)
781    !
782    !Config Key = ZSNOWTHRMCOND2
783    !Config Desc = Thermal conductivity Coef 2
784    !Config If = OK_SECHIBA
785    !Config Def = 2.5E-6
786    !Config Help =
787    !Config Units = [W m5/(kg2 K)]
788    CALL getin_p('ZSNOWTHRMCOND2',ZSNOWTHRMCOND2)
789    !
790    !Config Key = ZSNOWTHRMCOND_AVAP
791    !Config Desc = Thermal conductivity Coef 1 water vapor
792    !Config If = OK_SECHIBA
793    !Config Def = -0.06023
794    !Config Help =
795    !Config Units = [W/m/K]
796    CALL getin_p('ZSNOWTHRMCOND_AVAP',ZSNOWTHRMCOND_AVAP)
797    !
798    !Config Key = ZSNOWTHRMCOND_BVAP
799    !Config Desc = Thermal conductivity Coef 2 water vapor
800    !Config If = OK_SECHIBA
801    !Config Def = -2.5425
802    !Config Help =
803    !Config Units = [W/m]
804    CALL getin_p('ZSNOWTHRMCOND_BVAP',ZSNOWTHRMCOND_BVAP)
805    !
806    !Config Key = ZSNOWTHRMCOND_CVAP
807    !Config Desc = Thermal conductivity Coef 3 water vapor
808    !Config If = OK_SECHIBA
809    !Config Def = -289.99
810    !Config Help =
811    !Config Units = [K]
812    CALL getin_p('ZSNOWTHRMCOND_CVAP',ZSNOWTHRMCOND_CVAP)
813
814    !Snow compaction factors
815    !Config Key = ZSNOWCMPCT_RHOD
816    !Config Desc = Snow compaction coefficent
817    !Config If = OK_SECHIBA
818    !Config Def = 150.0
819    !Config Help =
820    !Config Units = [kg/m3]
821    CALL getin_p('ZSNOWCMPCT_RHOD',ZSNOWCMPCT_RHOD)
822
823    !Config Key = ZSNOWCMPCT_ACM
824    !Config Desc = Coefficent for the thermal conductivity
825    !Config If = OK_SECHIBA
826    !Config Def = 2.8e-6
827    !Config Help =
828    !Config Units = [1/s]
829    CALL getin_p('ZSNOWCMPCT_ACM',ZSNOWCMPCT_ACM)
830
831    !Config Key = ZSNOWCMPCT_BCM
832    !Config Desc = Coefficent for the thermal conductivity
833    !Config If = OK_SECHIBA
834    !Config Def = 0.04
835    !Config Help =
836    !Config Units = [1/K]
837    CALL getin_p('ZSNOWCMPCT_BCM',ZSNOWCMPCT_BCM)
838
839    !Config Key = ZSNOWCMPCT_CCM
840    !Config Desc = Coefficent for the thermal conductivity
841    !Config If = OK_SECHIBA
842    !Config Def = 460.
843    !Config Help =
844    !Config Units = [m3/kg]
845    CALL getin_p('ZSNOWCMPCT_CCM',ZSNOWCMPCT_CCM)
846
847    !Config Key = ZSNOWCMPCT_V0
848    !Config Desc = Vapor coefficent for the thermal conductivity
849    !Config If = OK_SECHIBA
850    !Config Def = 3.7e7
851    !Config Help =
852    !Config Units = [Pa/s]
853    CALL getin_p('ZSNOWCMPCT_V0',ZSNOWCMPCT_V0)
854
855    !Config Key = ZSNOWCMPCT_VT
856    !Config Desc = Vapor coefficent for the thermal conductivity
857    !Config If = OK_SECHIBA
858    !Config Def = 0.081
859    !Config Help =
860    !Config Units = [1/K]
861    CALL getin_p('ZSNOWCMPCT_VT',ZSNOWCMPCT_VT)
862
863    !Config Key = ZSNOWCMPCT_VR
864    !Config Desc = Vapor coefficent for the thermal conductivity
865    !Config If = OK_SECHIBA
866    !Config Def = 0.018
867    !Config Help =
868    !Config Units = [m3/kg]
869    CALL getin_p('ZSNOWCMPCT_VR',ZSNOWCMPCT_VR)
870
871
872    !Surface resistance
873    !
874    !Config Key = CB
875    !Config Desc = Constant of the Louis scheme
876    !Config If = OK_SECHIBA
877    !Config Def = 5.0
878    !Config Help =
879    !Config Units = [-]
880    CALL getin_p('CB',cb)
881    !
882    !Config Key = CC
883    !Config Desc = Constant of the Louis scheme
884    !Config If = OK_SECHIBA
885    !Config Def = 5.0
886    !Config Help =
887    !Config Units = [-]
888    CALL getin_p('CC',cc)
889    !
890    !Config Key = CD
891    !Config Desc = Constant of the Louis scheme
892    !Config If = OK_SECHIBA
893    !Config Def = 5.0
894    !Config Help =
895    !Config Units = [-]
896    CALL getin_p('CD',cd)
897    !
898    !Config Key = RAYT_CSTE
899    !Config Desc = Constant in the computation of surface resistance 
900    !Config If = OK_SECHIBA
901    !Config Def = 125
902    !Config Help =
903    !Config Units = [W.m^{-2}]
904    CALL getin_p('RAYT_CSTE',rayt_cste)
905    !
906    !Config Key = DEFC_PLUS
907    !Config Desc = Constant in the computation of surface resistance 
908    !Config If = OK_SECHIBA
909    !Config Def = 23.E-3
910    !Config Help =
911    !Config Units = [K.W^{-1}]
912    CALL getin_p('DEFC_PLUS',defc_plus)
913    !
914    !Config Key = DEFC_MULT
915    !Config Desc = Constant in the computation of surface resistance 
916    !Config If = OK_SECHIBA
917    !Config Def = 1.5
918    !Config Help =
919    !Config Units = [K.W^{-1}]
920    CALL getin_p('DEFC_MULT',defc_mult)
921    !
922
923    !
924    !-
925    ! diffuco
926    !-
927    !
928    !Config Key   = NLAI
929    !Config Desc  = Number of LAI levels
930    !Config If    = OK_SECHIBA
931    !Config Def   = 20
932    !Config Help  =
933    !Config Units = [-] 
934    CALL getin_p('NLAI',nlai)
935    !
936    !Config Key   = LAIMAX
937    !Config Desc  = Maximum LAI
938    !Config If    = OK_SECHIBA
939    !Config Def   =
940    !Config Help  =
941    !Config Units = [m^2/m^2]   
942    CALL getin_p('LAIMAX',laimax)
943    !
944    !Config Key   = DEW_VEG_POLY_COEFF
945    !Config Desc  = coefficients of the polynome of degree 5 for the dew
946    !Config If    = OK_SECHIBA
947    !Config Def   = 0.887773, 0.205673, 0.110112, 0.014843, 0.000824, 0.000017
948    !Config Help  =
949    !Config Units = [-]   
950    CALL getin_p('DEW_VEG_POLY_COEFF',dew_veg_poly_coeff)
951    !
952    !Config Key   = DOWNREGULATION_CO2
953    !Config Desc  = Activation of CO2 downregulation (used for CMIP6 version 6.1.0-6.1.10)
954    !Config If    = OK_SECHIBA
955    !Config Def   = y
956    !Config Help  =
957    !Config Units = [FLAG]   
958    CALL getin_p('DOWNREGULATION_CO2',downregulation_co2)
959
960    !Config Key   = DOWNREGULATION_CO2_NEW
961    !Config Desc  = Activation of CO2 downregulation (used for CMIP6 version 6.1.11 and later)
962    !Config If    = OK_SECHIBA
963    !Config Def   = n
964    !Config Help  = See also information in the ticket
965    !Config         http://forge.ipsl.jussieu.fr/orchidee/ticket/641
966    !Config Units = [FLAG]   
967    CALL getin_p('DOWNREGULATION_CO2_NEW',downregulation_co2_new)
968   
969    IF (downregulation_co2_new .AND. downregulation_co2) THEN
970       downregulation_co2=.FALSE.
971       CALL ipslerr_p(2,"config_sechiba_parameters",&
972            "Both DOWNREGULATION_CO2 and DOWNREGULATION_CO2_NEW were set to TRUE.",&
973            "DOWNREGULATION_CO2_NEW will be kept to TRUE.", &
974            "DOWNREGULATION_CO2 will be set to FALSE.")
975    END IF
976
977    !Config Key   = DOWNREGULATION_CO2_BASELEVEL
978    !Config Desc  = CO2 base level
979    !Config If    = DOWNREGULATION_CO2 or DOWNREGULATION_CO2_NEW
980    !Config Def   = 380.
981    !Config Help  =
982    !Config Units = [ppm]   
983    CALL getin_p('DOWNREGULATION_CO2_BASELEVEL',downregulation_co2_baselevel)
984   
985    !Config Key   = DOWNREGULATION_CO2_MINIMUM
986    !Config Desc  = CO2 value above which downregulation is taken into account
987    !Config If    = DOWNREGULATION_CO2_NEW
988    !Config Def   = 280.
989    !Config Help  =
990    !Config Units = [ppm]   
991    CALL getin_p('DOWNREGULATION_CO2_MINIMUM',downregulation_co2_minimum)
992 
993    !Config Key   = GB_REF
994    !Config Desc  = Leaf bulk boundary layer resistance
995    !Config If    =
996    !Config Def   = 1./25.
997    !Config Help  =
998    !Config Units = [s m-1]   
999    CALL getin_p('GB_REF',gb_ref)
1000
1001
1002    !-
1003    ! slowproc
1004    !-
1005    !
1006    !Config Key   = MIN_VEGFRAC
1007    !Config Desc  = Minimal fraction of mesh a vegetation type can occupy
1008    !Config If    = OK_SECHIBA
1009    !Config Def   = 0.001
1010    !Config Help  =
1011    !Config Units = [-] 
1012    CALL getin_p('MIN_VEGFRAC',min_vegfrac)
1013    !
1014    !Config Key   = STEMPDIAG_BID
1015    !Config Desc  = only needed for an initial LAI if there is no restart file
1016    !Config If    = OK_SECHIBA
1017    !Config Def   = 280.
1018    !Config Help  =
1019    !Config Units = [K]
1020    CALL getin_p('STEMPDIAG_BID',stempdiag_bid)
1021    !
1022
1023  END SUBROUTINE config_sechiba_parameters
1024
1025
1026!! ================================================================================================================================
1027!! SUBROUTINE   : config_co2_parameters
1028!!
1029!>\BRIEF        This subroutine reads in the configuration file all the parameters when impose_param=TRUE
1030!!
1031!! DESCRIPTION  : None
1032!!
1033!! RECENT CHANGE(S): None
1034!!
1035!! MAIN OUTPUT VARIABLE(S): None
1036!!
1037!! REFERENCE(S) :
1038!!
1039!! FLOWCHART    :
1040!! \n
1041!_ ================================================================================================================================
1042
1043  SUBROUTINE config_co2_parameters
1044
1045    IMPLICIT NONE
1046
1047    !! 0. Variables and parameters declaration
1048
1049    !! 0.4 Local variables
1050
1051    !_ ================================================================================================================================
1052
1053    !
1054    !Config Key   = LAI_LEVEL_DEPTH
1055    !Config Desc  =
1056    !Config If    =
1057    !Config Def   = 0.15
1058    !Config Help  =
1059    !Config Units = [-] 
1060    CALL getin_p('LAI_LEVEL_DEPTH',lai_level_depth)
1061    !
1062    !Config Key   = Oi
1063    !Config Desc  = Intercellular oxygen partial pressure
1064    !Config If    =
1065    !Config Def   = 210000.
1066    !Config Help  = See Legend of Figure 6 of Yin et al. (2009)
1067    !Config Units = [ubar] 
1068    CALL getin_p('Oi',Oi)
1069
1070
1071  END SUBROUTINE config_co2_parameters
1072
1073
1074!! ================================================================================================================================
1075!! SUBROUTINE   : config_stomate_parameters
1076!!
1077!>\BRIEF        This subroutine reads in the configuration file all the parameters
1078!! needed when stomate is activated (ie : when OK_STOMATE is set to true).
1079!!
1080!! DESCRIPTION  : None
1081!!
1082!! RECENT CHANGE(S): None
1083!!
1084!! MAIN OUTPUT VARIABLE(S):
1085!!
1086!! REFERENCE(S) :
1087!!
1088!! FLOWCHART    :
1089!! \n
1090!_ ================================================================================================================================
1091
1092  SUBROUTINE config_stomate_parameters
1093
1094    IMPLICIT NONE
1095
1096    !! 0. Variables and parameters declaration
1097
1098    !! 0.4 Local variables   
1099
1100
1101    !_ ================================================================================================================================
1102
1103    !-
1104    ! constraints_parameters
1105    !-
1106    !
1107    !Config Key   = TOO_LONG
1108    !Config Desc  = longest sustainable time without regeneration (vernalization)
1109    !Config If    = OK_STOMATE
1110    !Config Def   = 5.
1111    !Config Help  =
1112    !Config Units = [days]   
1113    CALL getin_p('TOO_LONG',too_long)
1114
1115    !-
1116    ! fire parameters
1117    !-
1118    !
1119    !Config Key   = TAU_FIRE
1120    !Config Desc  = Time scale for memory of the fire index (days). Validated for one year in the DGVM.
1121    !Config If    = OK_STOMATE
1122    !Config Def   = 30.
1123    !Config Help  =
1124    !Config Units = [days]   
1125    CALL getin_p('TAU_FIRE',tau_fire)
1126    !
1127    !Config Key   = LITTER_CRIT
1128    !Config Desc  = Critical litter quantity for fire
1129    !Config If    = OK_STOMATE
1130    !Config Def   = 200.
1131    !Config Help  =
1132    !Config Units = [gC/m^2] 
1133    CALL getin_p('LITTER_CRIT',litter_crit)
1134    !
1135    !Config Key   = FIRE_RESIST_STRUCT
1136    !Config Desc  =
1137    !Config If    = OK_STOMATE
1138    !Config Def   = 0.5
1139    !Config Help  =
1140    !Config Units = [-] 
1141    CALL getin_p('FIRE_RESIST_STRUCT',fire_resist_struct)
1142    !
1143    !
1144    !Config Key   = CO2FRAC
1145    !Config Desc  = What fraction of a burned plant compartment goes into the atmosphere
1146    !Config If    = OK_STOMATE
1147    !Config Def   = 0.95, 0.95, 0., 0.3, 0., 0., 0.95, 0.95
1148    !Config Help  =
1149    !Config Units = [-] 
1150    CALL getin_p('CO2FRAC',co2frac)
1151    !
1152    !Config Key   = BCFRAC_COEFF
1153    !Config Desc  =
1154    !Config If    = OK_STOMATE
1155    !Config Def   = 0.3, 1.3, 88.2
1156    !Config Help  =
1157    !Config Units = [-] 
1158    CALL getin_p('BCFRAC_COEFF',bcfrac_coeff)
1159    !
1160    !Config Key   = FIREFRAC_COEFF
1161    !Config Desc  =
1162    !Config If    = OK_STOMATE
1163    !Config Def   = 0.45, 0.8, 0.6, 0.13
1164    !Config Help  =
1165    !Config Units = [-]   
1166    CALL getin_p('FIREFRAC_COEFF',firefrac_coeff)
1167
1168    !Config Key   = REF_GREFF
1169    !Config Desc  = Asymptotic maximum mortality rate
1170    !Config If    = OK_STOMATE
1171    !Config Def   = 0.035
1172    !Config Help  = Set asymptotic maximum mortality rate from Sitch 2003
1173    !Config         (they use 0.01) (year^{-1})
1174    !Config Units = [1/year] 
1175    CALL getin_p('REF_GREFF',ref_greff)
1176    !-
1177    ! allocation parameters
1178    !-
1179    !
1180    !Config Key   = OK_MINRES
1181    !Config Desc  = Do we try to reach a minimum reservoir even if we are severely stressed?
1182    !Config If    = OK_STOMATE
1183    !Config Def   = y
1184    !Config Help  =
1185    !Config Units = [FLAG]
1186    CALL getin_p('OK_MINRES',ok_minres)
1187    !
1188    !Config Key   = RESERVE_TIME_TREE
1189    !Config Desc  = maximum time during which reserve is used (trees)
1190    !Config If    = OK_STOMATE
1191    !Config Def   = 30.
1192    !Config Help  =
1193    !Config Units = [days]   
1194    CALL getin_p('RESERVE_TIME_TREE',reserve_time_tree)
1195    !
1196    !Config Key   = RESERVE_TIME_GRASS
1197    !Config Desc  = maximum time during which reserve is used (grasses)
1198    !Config If    = OK_STOMATE
1199    !Config Def   = 20.
1200    !Config Help  =
1201    !Config Units = [days]   
1202    CALL getin_p('RESERVE_TIME_GRASS',reserve_time_grass)
1203    !
1204    !Config Key   = F_FRUIT
1205    !Config Desc  = Standard fruit allocation
1206    !Config If    = OK_STOMATE
1207    !Config Def   = 0.1
1208    !Config Help  =
1209    !Config Units = [-]   
1210    CALL getin_p('F_FRUIT',f_fruit)
1211    !
1212    !Config Key   = ALLOC_SAP_ABOVE_GRASS
1213    !Config Desc  = fraction of sapwood allocation above ground
1214    !Config If    = OK_STOMATE
1215    !Config Def   = 1.0
1216    !Config Help  =
1217    !Config Units = [-]   
1218    CALL getin_p('ALLOC_SAP_ABOVE_GRASS',alloc_sap_above_grass)
1219    !
1220    !Config Key   = MIN_LTOLSR
1221    !Config Desc  = extrema of leaf allocation fraction
1222    !Config If    = OK_STOMATE
1223    !Config Def   = 0.2
1224    !Config Help  =
1225    !Config Units = [-]   
1226    CALL getin_p('MIN_LTOLSR',min_LtoLSR)
1227    !
1228    !Config Key   = MAX_LTOLSR
1229    !Config Desc  = extrema of leaf allocation fraction
1230    !Config If    = OK_STOMATE
1231    !Config Def   = 0.5
1232    !Config Help  =
1233    !Config Units = [-]   
1234    CALL getin_p('MAX_LTOLSR',max_LtoLSR)
1235    !
1236    !Config Key   = Z_NITROGEN
1237    !Config Desc  = scaling depth for nitrogen limitation
1238    !Config If    = OK_STOMATE
1239    !Config Def   = 0.2
1240    !Config Help  =
1241    !Config Units = [m] 
1242    CALL getin_p('Z_NITROGEN',z_nitrogen)
1243    !
1244    !Config Key   = NLIM_TREF
1245    !Config Desc  =
1246    !Config If    = OK_STOMATE
1247    !Config Def   = 25.
1248    !Config Help  =
1249    !Config Units = [C] 
1250    CALL getin_p('NLIM_TREF',Nlim_tref) 
1251
1252    !-
1253    ! data parameters
1254    !-
1255    !
1256    !Config Key   = PIPE_TUNE1
1257    !Config Desc  = crown area = pipe_tune1. stem diameter**(1.6) (Reinicke's theory)
1258    !Config If    = OK_STOMATE
1259    !Config Def   = 100.0
1260    !Config Help  =
1261    !Config Units = [-]   
1262    CALL getin_p('PIPE_TUNE1',pipe_tune1)
1263    !
1264    !Config Key   = PIPE_TUNE2
1265    !Config Desc  = height=pipe_tune2 * diameter**pipe_tune3
1266    !Config If    = OK_STOMATE
1267    !Config Def   = 40.0
1268    !Config Help  =
1269    !Config Units = [-]     
1270    CALL getin_p('PIPE_TUNE2',pipe_tune2) 
1271    !
1272    !Config Key   = PIPE_TUNE3
1273    !Config Desc  = height=pipe_tune2 * diameter**pipe_tune3
1274    !Config If    = OK_STOMATE
1275    !Config Def   = 0.5
1276    !Config Help  =
1277    !Config Units = [-]   
1278    CALL getin_p('PIPE_TUNE3',pipe_tune3)
1279    !
1280    !Config Key   = PIPE_TUNE4
1281    !Config Desc  = needed for stem diameter
1282    !Config If    = OK_STOMATE
1283    !Config Def   = 0.3
1284    !Config Help  =
1285    !Config Units = [-] 
1286    CALL getin_p('PIPE_TUNE4',pipe_tune4)
1287    !
1288    !Config Key   = PIPE_DENSITY
1289    !Config Desc  = Density
1290    !Config If    = OK_STOMATE
1291    !Config Def   = 2.e5
1292    !Config Help  =
1293    !Config Units = [-] 
1294    CALL getin_p('PIPE_DENSITY',pipe_density)
1295    !
1296    !Config Key   = PIPE_K1
1297    !Config Desc  =
1298    !Config If    = OK_STOMATE
1299    !Config Def   = 8.e3
1300    !Config Help  =
1301    !Config Units = [-]   
1302    CALL getin_p('PIPE_K1',pipe_k1)
1303    !
1304    !Config Key   = PIPE_TUNE_EXP_COEFF
1305    !Config Desc  = pipe tune exponential coeff
1306    !Config If    = OK_STOMATE
1307    !Config Def   = 1.6
1308    !Config Help  =
1309    !Config Units = [-]   
1310    CALL getin_p('PIPE_TUNE_EXP_COEFF',pipe_tune_exp_coeff)
1311    !
1312    !
1313    !Config Key   = PRECIP_CRIT
1314    !Config Desc  = minimum precip
1315    !Config If    = OK_STOMATE
1316    !Config Def   = 100.
1317    !Config Help  =
1318    !Config Units = [mm/year] 
1319    CALL getin_p('PRECIP_CRIT',precip_crit)
1320    !
1321    !Config Key   = GDD_CRIT_ESTAB
1322    !Config Desc  = minimum gdd for establishment of saplings
1323    !Config If    = OK_STOMATE
1324    !Config Def   = 150.
1325    !Config Help  =
1326    !Config Units = [-] 
1327    CALL getin_p('GDD_CRIT_ESTAB',gdd_crit_estab)
1328    !
1329    !Config Key   = FPC_CRIT
1330    !Config Desc  = critical fpc, needed for light competition and establishment
1331    !Config If    = OK_STOMATE
1332    !Config Def   = 0.95
1333    !Config Help  =
1334    !Config Units = [-] 
1335    CALL getin_p('FPC_CRIT',fpc_crit)
1336    !
1337    !Config Key   = ALPHA_GRASS
1338    !Config Desc  = sapling characteristics : alpha's
1339    !Config If    = OK_STOMATE
1340    !Config Def   = 0.5
1341    !Config Help  =
1342    !Config Units = [-]   
1343    CALL getin_p('ALPHA_GRASS',alpha_grass)
1344    !
1345    !Config Key   = ALPHA_TREE
1346    !Config Desc  = sapling characteristics : alpha's
1347    !Config If    = OK_STOMATE
1348    !Config Def   = 1.
1349    !Config Help  =
1350    !Config Units = [-]   
1351    CALL getin_p('ALPHA_TREE',alpha_tree)
1352    !-
1353    !
1354    !Config Key   = MASS_RATIO_HEART_SAP
1355    !Config Desc  = mass ratio (heartwood+sapwood)/sapwood
1356    !Config If    = OK_STOMATE
1357    !Config Def   = 3.
1358    !Config Help  =
1359    !Config Units = [-]   
1360    CALL getin_p('MASS_RATIO_HEART_SAP',mass_ratio_heart_sap)
1361    !
1362    !Config Key   = TAU_HUM_MONTH
1363    !Config Desc  = time scales for phenology and other processes
1364    !Config If    = OK_STOMATE
1365    !Config Def   = 20.
1366    !Config Help  =
1367    !Config Units = [days] 
1368    CALL getin_p('TAU_HUM_MONTH',tau_hum_month)
1369    !
1370    !Config Key   = TAU_HUM_WEEK
1371    !Config Desc  = time scales for phenology and other processes
1372    !Config If    = OK_STOMATE
1373    !Config Def   = 7.
1374    !Config Help  =
1375    !Config Units = [days]   
1376    CALL getin_p('TAU_HUM_WEEK',tau_hum_week)
1377    !
1378    !Config Key   = TAU_T2M_MONTH
1379    !Config Desc  = time scales for phenology and other processes
1380    !Config If    = OK_STOMATE
1381    !Config Def   = 20.
1382    !Config Help  =
1383    !Config Units = [days]     
1384    CALL getin_p('TAU_T2M_MONTH',tau_t2m_month)
1385    !
1386    !Config Key   = TAU_T2M_WEEK
1387    !Config Desc  = time scales for phenology and other processes
1388    !Config If    = OK_STOMATE
1389    !Config Def   = 7.
1390    !Config Help  =
1391    !Config Units = [days]   
1392    CALL getin_p('TAU_T2M_WEEK',tau_t2m_week)
1393    !
1394    !Config Key   = TAU_TSOIL_MONTH
1395    !Config Desc  = time scales for phenology and other processes
1396    !Config If    = OK_STOMATE
1397    !Config Def   = 20.
1398    !Config Help  =
1399    !Config Units = [days]     
1400    CALL getin_p('TAU_TSOIL_MONTH',tau_tsoil_month)
1401    !
1402    !Config Key   = TAU_SOILHUM_MONTH
1403    !Config Desc  = time scales for phenology and other processes
1404    !Config If    = OK_STOMATE
1405    !Config Def   = 20.
1406    !Config Help  =
1407    !Config Units = [days]   
1408    CALL getin_p('TAU_SOILHUM_MONTH',tau_soilhum_month)
1409    !
1410    !Config Key   = TAU_GPP_WEEK
1411    !Config Desc  = time scales for phenology and other processes
1412    !Config If    = OK_STOMATE
1413    !Config Def   = 7.
1414    !Config Help  =
1415    !Config Units = [days]   
1416    CALL getin_p('TAU_GPP_WEEK',tau_gpp_week)
1417    !
1418    !Config Key   = TAU_GDD
1419    !Config Desc  = time scales for phenology and other processes
1420    !Config If    = OK_STOMATE
1421    !Config Def   = 40.
1422    !Config Help  =
1423    !Config Units = [days]   
1424    CALL getin_p('TAU_GDD',tau_gdd)
1425    !
1426    !Config Key   = TAU_NGD
1427    !Config Desc  = time scales for phenology and other processes
1428    !Config If    = OK_STOMATE
1429    !Config Def   = 50.
1430    !Config Help  =
1431    !Config Units = [days]   
1432    CALL getin_p('TAU_NGD',tau_ngd)
1433    !
1434    !Config Key   = COEFF_TAU_LONGTERM
1435    !Config Desc  = time scales for phenology and other processes
1436    !Config If    = OK_STOMATE
1437    !Config Def   = 3.
1438    !Config Help  =
1439    !Config Units = [days]   
1440    CALL getin_p('COEFF_TAU_LONGTERM',coeff_tau_longterm)
1441    !-
1442    !
1443    !Config Key   = BM_SAPL_CARBRES
1444    !Config Desc  =
1445    !Config If    = OK_STOMATE
1446    !Config Def   = 5.
1447    !Config Help  =
1448    !Config Units = [-]   
1449    CALL getin_p('BM_SAPL_CARBRES',bm_sapl_carbres)
1450    !
1451    !Config Key   = BM_SAPL_SAPABOVE
1452    !Config Desc  =
1453    !Config If    = OK_STOMATE
1454    !Config Def   = 0.5
1455    !Config Help  =
1456    !Config Units = [-]   
1457    CALL getin_p('BM_SAPL_SAPABOVE',bm_sapl_sapabove)
1458    !
1459    !Config Key   = BM_SAPL_HEARTABOVE
1460    !Config Desc  =
1461    !Config If    = OK_STOMATE
1462    !Config Def   = 2.
1463    !Config Help  =
1464    !Config Units = [-]   
1465    CALL getin_p('BM_SAPL_HEARTABOVE',bm_sapl_heartabove)
1466    !
1467    !Config Key   = BM_SAPL_HEARTBELOW
1468    !Config Desc  =
1469    !Config If    = OK_STOMATE
1470    !Config Def   = 2.
1471    !Config Help  =
1472    !Config Units = [-]   
1473    CALL getin_p('BM_SAPL_HEARTBELOW',bm_sapl_heartbelow)
1474    !
1475    !Config Key   = INIT_SAPL_MASS_LEAF_NAT
1476    !Config Desc  =
1477    !Config If    = OK_STOMATE
1478    !Config Def   = 0.1
1479    !Config Help  =
1480    !Config Units = [-]   
1481    CALL getin_p('INIT_SAPL_MASS_LEAF_NAT',init_sapl_mass_leaf_nat)
1482    !
1483    !Config Key   = INIT_SAPL_MASS_LEAF_AGRI
1484    !Config Desc  =
1485    !Config If    = OK_STOMATE
1486    !Config Def   = 1.
1487    !Config Help  =
1488    !Config Units = [-]   
1489    CALL getin_p('INIT_SAPL_MASS_LEAF_AGRI',init_sapl_mass_leaf_agri)
1490    !
1491    !Config Key   = INIT_SAPL_MASS_CARBRES
1492    !Config Desc  =
1493    !Config If    = OK_STOMATE
1494    !Config Def   = 5.
1495    !Config Help  =
1496    !Config Units = [-]   
1497    CALL getin_p('INIT_SAPL_MASS_CARBRES',init_sapl_mass_carbres)
1498    !
1499    !Config Key   = INIT_SAPL_MASS_ROOT
1500    !Config Desc  =
1501    !Config If    = OK_STOMATE
1502    !Config Def   = 0.1
1503    !Config Help  =
1504    !Config Units = [-]   
1505    CALL getin_p('INIT_SAPL_MASS_ROOT',init_sapl_mass_root)
1506    !
1507    !Config Key   = INIT_SAPL_MASS_FRUIT
1508    !Config Desc  =
1509    !Config If    = OK_STOMATE
1510    !Config Def   = 0.3
1511    !Config Help  =
1512    !Config Units = [-]   
1513    CALL getin_p('INIT_SAPL_MASS_FRUIT',init_sapl_mass_fruit)
1514    !
1515    !Config Key   = CN_SAPL_INIT
1516    !Config Desc  =
1517    !Config If    = OK_STOMATE
1518    !Config Def   = 0.5
1519    !Config Help  =
1520    !Config Units = [-]   
1521    CALL getin_p('CN_SAPL_INIT',cn_sapl_init)
1522    !
1523    !Config Key   = MIGRATE_TREE
1524    !Config Desc  =
1525    !Config If    = OK_STOMATE
1526    !Config Def   = 10000.
1527    !Config Help  =
1528    !Config Units = [m/year]   
1529    CALL getin_p('MIGRATE_TREE',migrate_tree)
1530    !
1531    !Config Key   = MIGRATE_GRASS
1532    !Config Desc  =
1533    !Config If    = OK_STOMATE
1534    !Config Def   = 10000.
1535    !Config Help  =
1536    !Config Units = [m/year]   
1537    CALL getin_p('MIGRATE_GRASS',migrate_grass)
1538    !
1539    !Config Key   = LAI_INITMIN_TREE
1540    !Config Desc  =
1541    !Config If    = OK_STOMATE
1542    !Config Def   = 0.3
1543    !Config Help  =
1544    !Config Units = [m^2/m^2] 
1545    CALL getin_p('LAI_INITMIN_TREE',lai_initmin_tree)
1546    !
1547    !Config Key   = LAI_INITMIN_GRASS
1548    !Config Desc  =
1549    !Config If    = OK_STOMATE
1550    !Config Def   = 0.1
1551    !Config Help  =
1552    !Config Units = [m^2/m^2]   
1553    CALL getin_p('LAI_INITMIN_GRASS',lai_initmin_grass)
1554    !
1555    !Config Key   = DIA_COEFF
1556    !Config Desc  =
1557    !Config If    = OK_STOMATE
1558    !Config Def   = 4., 0.5
1559    !Config Help  =
1560    !Config Units = [-]   
1561    CALL getin_p('DIA_COEFF',dia_coeff)
1562    !
1563    !Config Key   = MAXDIA_COEFF
1564    !Config Desc  =
1565    !Config If    = OK_STOMATE
1566    !Config Def   = 100., 0.01
1567    !Config Help  =
1568    !Config Units = [-]   
1569    CALL getin_p('MAXDIA_COEFF',maxdia_coeff)
1570    !
1571    !Config Key   = BM_SAPL_LEAF
1572    !Config Desc  =
1573    !Config If    = OK_STOMATE
1574    !Config Def   = 4., 4., 0.8, 5.
1575    !Config Help  =
1576    !Config Units = [-] 
1577    CALL getin_p('BM_SAPL_LEAF',bm_sapl_leaf)
1578
1579    !-
1580    ! litter parameters
1581    !-
1582    !
1583    !Config Key   = METABOLIC_REF_FRAC
1584    !Config Desc  =
1585    !Config If    = OK_STOMATE
1586    !Config Def   = 0.85 
1587    !Config Help  =
1588    !Config Units = [-]
1589    CALL getin_p('METABOLIC_REF_FRAC',metabolic_ref_frac)
1590    !
1591    !Config Key   = Z_DECOMP
1592    !Config Desc  = scaling depth for soil activity
1593    !Config If    = OK_STOMATE
1594    !Config Def   = 0.2
1595    !Config Help  =
1596    !Config Units = [m]   
1597    CALL getin_p('Z_DECOMP',z_decomp)
1598    !
1599    !Config Key   = CN
1600    !Config Desc  = C/N ratio
1601    !Config If    = OK_STOMATE
1602    !Config Def   = 40., 40., 40., 40., 40., 40., 40., 40.
1603    !Config Help  =
1604    !Config Units = [-] 
1605    CALL getin_p('CN',CN)
1606    !
1607    !Config Key   = LC
1608    !Config Desc  = Lignine/C ratio of the different plant parts
1609    !Config If    = OK_STOMATE
1610    !Config Def   = 0.22, 0.35, 0.35, 0.35, 0.35, 0.22, 0.22, 0.22
1611    !Config Help  =
1612    !Config Units = [-]   
1613    CALL getin_p('LC',LC)
1614    !
1615    !Config Key   = FRAC_SOIL_STRUCT_AA
1616    !Config Desc  = frac_soil(istructural,iactive,iabove)
1617    !Config If    = OK_STOMATE
1618    !Config Def   = 0.55
1619    !Config Help  =
1620    !Config Units = [-]
1621    CALL getin_p('FRAC_SOIL_STRUCT_AA',frac_soil_struct_aa)
1622    !
1623    !Config Key   = FRAC_SOIL_STRUCT_A
1624    !Config Desc  = frac_soil(istructural,iactive,ibelow)
1625    !Config If    = OK_STOMATE
1626    !Config Def   = 0.45
1627    !Config Help  =
1628    !Config Units = [-]
1629    CALL getin_p('FRAC_SOIL_STRUCT_AB',frac_soil_struct_ab)
1630    !
1631    !Config Key   = FRAC_SOIL_STRUCT_SA
1632    !Config Desc  = frac_soil(istructural,islow,iabove)
1633    !Config If    = OK_STOMATE
1634    !Config Def   = 0.7 
1635    !Config Help  =
1636    !Config Units = [-]   
1637    CALL getin_p('FRAC_SOIL_STRUCT_SA',frac_soil_struct_sa)
1638    !
1639    !Config Key   = FRAC_SOIL_STRUCT_SB
1640    !Config Desc  = frac_soil(istructural,islow,ibelow)
1641    !Config If    = OK_STOMATE
1642    !Config Def   = 0.7 
1643    !Config Help  =
1644    !Config Units = [-]   
1645    CALL getin_p('FRAC_SOIL_STRUCT_SB',frac_soil_struct_sb)
1646    !
1647    !Config Key   = FRAC_SOIL_METAB_AA
1648    !Config Desc  = frac_soil(imetabolic,iactive,iabove)
1649    !Config If    = OK_STOMATE
1650    !Config Def   = 0.45
1651    !Config Help  =
1652    !Config Units = [-]   
1653    CALL getin_p('FRAC_SOIL_METAB_AA',frac_soil_metab_aa)
1654    !
1655    !Config Key   = FRAC_SOIL_METAB_AB
1656    !Config Desc  = frac_soil(imetabolic,iactive,ibelow)
1657    !Config If    = OK_STOMATE
1658    !Config Def   = 0.45 
1659    !Config Help  =
1660    !Config Units = [-]   
1661    CALL getin_p('FRAC_SOIL_METAB_AB',frac_soil_metab_ab)
1662    !
1663    !
1664    !Config Key   = METABOLIC_LN_RATIO
1665    !Config Desc  =
1666    !Config If    = OK_STOMATE
1667    !Config Def   = 0.018 
1668    !Config Help  =
1669    !Config Units = [-]   
1670    CALL getin_p('METABOLIC_LN_RATIO',metabolic_LN_ratio) 
1671    !
1672    !Config Key   = TAU_METABOLIC
1673    !Config Desc  =
1674    !Config If    = OK_STOMATE
1675    !Config Def   = 0.066
1676    !Config Help  =
1677    !Config Units = [days]
1678    CALL getin_p('TAU_METABOLIC',tau_metabolic)
1679    !
1680    !Config Key   = TAU_STRUCT
1681    !Config Desc  =
1682    !Config If    = OK_STOMATE
1683    !Config Def   = 0.245
1684    !Config Help  =
1685    !Config Units = [days]
1686    CALL getin_p('TAU_STRUCT',tau_struct)
1687    !
1688    !Config Key   = SOIL_Q10
1689    !Config Desc  =
1690    !Config If    = OK_STOMATE
1691    !Config Def   = 0.69 (=ln2)
1692    !Config Help  =
1693    !Config Units = [-]
1694    CALL getin_p('SOIL_Q10',soil_Q10)
1695    !
1696    !Config Key   = TSOIL_REF
1697    !Config Desc  =
1698    !Config If    = OK_STOMATE
1699    !Config Def   = 30.
1700    !Config Help  =
1701    !Config Units = [C]   
1702    CALL getin_p('TSOIL_REF',tsoil_ref)
1703    !
1704    !Config Key   = LITTER_STRUCT_COEF
1705    !Config Desc  =
1706    !Config If    = OK_STOMATE
1707    !Config Def   = 3.
1708    !Config Help  =
1709    !Config Units = [-]   
1710    CALL getin_p('LITTER_STRUCT_COEF',litter_struct_coef)
1711    !
1712    !Config Key   = MOIST_COEFF
1713    !Config Desc  =
1714    !Config If    = OK_STOMATE
1715    !Config Def   = 1.1, 2.4, 0.29
1716    !Config Help  =
1717    !Config Units = [-]   
1718    CALL getin_p('MOIST_COEFF',moist_coeff)
1719    !
1720    !Config Key   = MOISTCONT_MIN
1721    !Config Desc  = minimum soil wetness to limit the heterotrophic respiration
1722    !Config If    = OK_STOMATE
1723    !Config Def   = 0.25
1724    !Config Help  =
1725    !Config Units = [-]
1726    CALL getin_p('MOISTCONT_MIN',moistcont_min)
1727
1728    !-
1729    ! lpj parameters
1730    !-
1731    !
1732    !Config Key   = FRAC_TURNOVER_DAILY
1733    !Config Desc  =
1734    !Config If    = OK_STOMATE
1735    !Config Def   = 0.55
1736    !Config Help  =
1737    !Config Units = [-]
1738    CALL getin_p('FRAC_TURNOVER_DAILY',frac_turnover_daily)   
1739
1740    !-
1741    ! npp parameters
1742    !-
1743    !
1744    !Config Key   = TAX_MAX
1745    !Config Desc  = maximum fraction of allocatable biomass used for maintenance respiration
1746    !Config If    = OK_STOMATE
1747    !Config Def   = 0.8
1748    !Config Help  =
1749    !Config Units = [-]   
1750    CALL getin_p('TAX_MAX',tax_max) 
1751
1752    !-
1753    ! phenology parameters
1754    !-
1755    !Config Key   = MIN_GROWTHINIT_TIME
1756    !Config Desc  = minimum time since last beginning of a growing season
1757    !Config If    = OK_STOMATE
1758    !Config Def   = 300.
1759    !Config Help  =
1760    !Config Units = [days] 
1761    CALL getin_p('MIN_GROWTHINIT_TIME',min_growthinit_time)
1762    !
1763    !Config Key   = MOIAVAIL_ALWAYS_TREE
1764    !Config Desc  = moisture availability above which moisture tendency doesn't matter
1765    !Config If    = OK_STOMATE
1766    !Config Def   = 1.0
1767    !Config Help  =
1768    !Config Units = [-]   
1769    CALL getin_p('MOIAVAIL_ALWAYS_TREE',moiavail_always_tree)
1770    !
1771    !Config Key   = MOIAVAIL_ALWAYS_GRASS
1772    !Config Desc  = moisture availability above which moisture tendency doesn't matter
1773    !Config If    = OK_STOMATE
1774    !Config Def   = 0.6
1775    !Config Help  =
1776    !Config Units = [-]   
1777    CALL getin_p('MOIAVAIL_ALWAYS_GRASS',moiavail_always_grass)
1778    !
1779    !Config Key   = T_ALWAYS_ADD
1780    !Config Desc  = monthly temp. above which temp. tendency doesn't matter
1781    !Config If    = OK_STOMATE
1782    !Config Def   = 10.
1783    !Config Help  =
1784    !Config Units = [C]   
1785    CALL getin_p('T_ALWAYS_ADD',t_always_add)
1786    !
1787    !
1788    !Config Key   = GDDNCD_REF
1789    !Config Desc  =
1790    !Config If    = OK_STOMATE
1791    !Config Def   = 603.
1792    !Config Help  =
1793    !Config Units = [-]   
1794    CALL getin_p('GDDNCD_REF',gddncd_ref)
1795    !
1796    !Config Key   = GDDNCD_CURVE
1797    !Config Desc  =
1798    !Config If    = OK_STOMATE
1799    !Config Def   = 0.0091
1800    !Config Help  =
1801    !Config Units = [-] 
1802    CALL getin_p('GDDNCD_CURVE',gddncd_curve)
1803    !
1804    !Config Key   = GDDNCD_OFFSET
1805    !Config Desc  =
1806    !Config If    = OK_STOMATE
1807    !Config Def   = 64.
1808    !Config Help  =
1809    !Config Units = [-] 
1810    CALL getin_p('GDDNCD_OFFSET',gddncd_offset)
1811    !-
1812    ! prescribe parameters
1813    !-
1814    !
1815    !Config Key   = BM_SAPL_RESCALE
1816    !Config Desc  =
1817    !Config If    = OK_STOMATE
1818    !Config Def   = 40.
1819    !Config Help  =
1820    !Config Units = [-] 
1821    CALL getin_p('BM_SAPL_RESCALE',bm_sapl_rescale)
1822
1823    !-
1824    ! respiration parameters
1825    !-
1826    !
1827    !Config Key   = MAINT_RESP_MIN_VMAX
1828    !Config Desc  =
1829    !Config If    = OK_STOMATE
1830    !Config Def   = 0.3
1831    !Config Help  =
1832    !Config Units = [-] 
1833    CALL getin_p('MAINT_RESP_MIN_VMAX',maint_resp_min_vmax) 
1834    !
1835    !Config Key   = MAINT_RESP_COEFF
1836    !Config Desc  =
1837    !Config If    = OK_STOMATE
1838    !Config Def   = 1.4
1839    !Config Help  =
1840    !Config Units = [-]
1841    CALL getin_p('MAINT_RESP_COEFF',maint_resp_coeff)
1842
1843    !-
1844    ! soilcarbon parameters
1845    !-
1846    !
1847    !Config Key   = FRAC_CARB_AP
1848    !Config Desc  = frac carb coefficients from active pool: depends on clay content
1849    !Config if    = OK_STOMATE
1850    !Config Def   = 0.004
1851    !Config Help  = fraction of the active pool going into the passive pool
1852    !Config Units = [-]
1853    CALL getin_p('FRAC_CARB_AP',frac_carb_ap) 
1854    !
1855    !Config Key   = FRAC_CARB_SA
1856    !Config Desc  = frac_carb_coefficients from slow pool
1857    !Config if    = OK_STOMATE
1858    !Config Def   = 0.42
1859    !Config Help  = fraction of the slow pool going into the active pool
1860    !Config Units = [-]
1861    CALL getin_p('FRAC_CARB_SA',frac_carb_sa)
1862    !
1863    !Config Key   = FRAC_CARB_SP
1864    !Config Desc  = frac_carb_coefficients from slow pool
1865    !Config if    = OK_STOMATE
1866    !Config Def   = 0.03
1867    !Config Help  = fraction of the slow pool going into the passive pool
1868    !Config Units = [-]
1869    CALL getin_p('FRAC_CARB_SP',frac_carb_sp)
1870    !
1871    !Config Key   = FRAC_CARB_PA
1872    !Config Desc  = frac_carb_coefficients from passive pool
1873    !Config if    = OK_STOMATE
1874    !Config Def   = 0.45
1875    !Config Help  = fraction of the passive pool going into the active pool
1876    !Config Units = [-]
1877    CALL getin_p('FRAC_CARB_PA',frac_carb_pa)
1878    !
1879    !Config Key   = FRAC_CARB_PS
1880    !Config Desc  = frac_carb_coefficients from passive pool
1881    !Config if    = OK_STOMATE
1882    !Config Def   = 0.0
1883    !Config Help  = fraction of the passive pool going into the slow pool
1884    !Config Units = [-]
1885    CALL getin_p('FRAC_CARB_PS',frac_carb_ps)
1886    !
1887    !Config Key   = ACTIVE_TO_PASS_CLAY_FRAC
1888    !Config Desc  =
1889    !Config if    = OK_STOMATE
1890    !Config Def   = 0.68 
1891    !Config Help  =
1892    !Config Units = [-]
1893    CALL getin_p('ACTIVE_TO_PASS_CLAY_FRAC',active_to_pass_clay_frac)
1894    !
1895    !Config Key   = CARBON_TAU_IACTIVE
1896    !Config Desc  = residence times in carbon pools
1897    !Config if    = OK_STOMATE
1898    !Config Def   = 0.149
1899    !Config Help  =
1900    !Config Units =  [days]
1901    CALL getin_p('CARBON_TAU_IACTIVE',carbon_tau_iactive)
1902    !
1903    !Config Key   = CARBON_TAU_ISLOW
1904    !Config Desc  = residence times in carbon pools
1905    !Config if    = OK_STOMATE
1906    !Config Def   = 7.0
1907    !Config Help  =
1908    !Config Units = [days]
1909    CALL getin_p('CARBON_TAU_ISLOW',carbon_tau_islow)
1910    !
1911    !Config Key   = CARBON_TAU_IPASSIVE
1912    !Config Desc  = residence times in carbon pools
1913    !Config if    = OK_STOMATE
1914    !Config Def   = 300.
1915    !Config Help  = residence time in the passive pool
1916    !Config Units = [days]
1917    CALL getin_p('CARBON_TAU_IPASSIVE',carbon_tau_ipassive)
1918    !
1919    !Config Key   = FLUX_TOT_COEFF
1920    !Config Desc  =
1921    !Config if    = OK_STOMATE
1922    !Config Def   = 1.2, 1.4,.75
1923    !Config Help  =
1924    !Config Units = [days]
1925    CALL getin_p('FLUX_TOT_COEFF',flux_tot_coeff)
1926
1927    !-
1928    ! turnover parameters
1929    !-
1930    !
1931    !Config Key   = NEW_TURNOVER_TIME_REF
1932    !Config Desc  =
1933    !Config If    = OK_STOMATE
1934    !Config Def   = 20.
1935    !Config Help  =
1936    !Config Units = [days] 
1937    CALL getin_p('NEW_TURNOVER_TIME_REF',new_turnover_time_ref)
1938
1939    !Config Key   = LEAF_AGE_CRIT_TREF
1940    !Config Desc  =
1941    !Config If    = OK_STOMATE
1942    !Config Def   = 20.
1943    !Config Help  =
1944    !Config Units = [days] 
1945    CALL getin_p('LEAF_AGE_CRIT_TREF',leaf_age_crit_tref)
1946    !
1947    !Config Key   = LEAF_AGE_CRIT_COEFF
1948    !Config Desc  =
1949    !Config If    = OK_STOMATE
1950    !Config Def   = 1.5, 0.75, 10.
1951    !Config Help  =
1952    !Config Units = [-]
1953    CALL getin_p('LEAF_AGE_CRIT_COEFF',leaf_age_crit_coeff)
1954
1955    !-
1956    ! vmax parameters
1957    !-
1958    !
1959    !Config Key   = VMAX_OFFSET
1960    !Config Desc  = offset (minimum relative vcmax)
1961    !Config If    = OK_STOMATE
1962    !Config Def   = 0.3
1963    !Config Help  = offset (minimum vcmax/vmax_opt)
1964    !Config Units = [-] 
1965    CALL getin_p('VMAX_OFFSET',vmax_offset)
1966    !
1967    !Config Key   = LEAFAGE_FIRSTMAX
1968    !Config Desc  = leaf age at which vmax attains vcmax_opt (in fraction of critical leaf age)
1969    !Config If    = OK_STOMATE
1970    !Config Def   = 0.03
1971    !Config Help  = relative leaf age at which vmax attains vcmax_opt
1972    !Config Units = [-]
1973    CALL getin_p('LEAFAGE_FIRSTMAX',leafage_firstmax)
1974    !
1975    !Config Key   = LEAFAGE_LASTMAX
1976    !Config Desc  = leaf age at which vmax falls below vcmax_opt (in fraction of critical leaf age)
1977    !Config If    = OK_STOMATE
1978    !Config Def   = 0.5
1979    !Config Help  = relative leaf age at which vmax falls below vcmax_opt
1980    !Config Units = [-] 
1981    CALL getin_p('LEAFAGE_LASTMAX',leafage_lastmax)
1982    !
1983    !Config Key   = LEAFAGE_OLD
1984    !Config Desc  = leaf age at which vmax attains its minimum (in fraction of critical leaf age)
1985    !Config If    = OK_STOMATE
1986    !Config Def   = 1.
1987    !Config Help  = relative leaf age at which vmax attains its minimum
1988    !Config Units = [-] 
1989    CALL getin_p('LEAFAGE_OLD',leafage_old)
1990
1991    !-
1992    ! season parameters
1993    !-
1994    !
1995    !Config Key   = GPPFRAC_DORMANCE
1996    !Config Desc  = rapport maximal GPP/GGP_max pour dormance
1997    !Config If    = OK_STOMATE
1998    !Config Def   = 0.2
1999    !Config Help  =
2000    !Config Units = [-]
2001    CALL getin_p('GPPFRAC_DORMANCE',gppfrac_dormance)
2002    !
2003    !Config Key   = TAU_CLIMATOLOGY
2004    !Config Desc  = tau for "climatologic variables
2005    !Config If    = OK_STOMATE
2006    !Config Def   = 20
2007    !Config Help  =
2008    !Config Units = [days]
2009    CALL getin_p('TAU_CLIMATOLOGY',tau_climatology)
2010    !
2011    !Config Key   = HVC1
2012    !Config Desc  = parameters for herbivore activity
2013    !Config If    = OK_STOMATE
2014    !Config Def   = 0.019
2015    !Config Help  =
2016    !Config Units = [-] 
2017    CALL getin_p('HVC1',hvc1)
2018    !
2019    !Config Key   = HVC2
2020    !Config Desc  = parameters for herbivore activity
2021    !Config If    = OK_STOMATE
2022    !Config Def   = 1.38
2023    !Config Help  =
2024    !Config Units = [-] 
2025    CALL getin_p('HVC2',hvc2)
2026    !
2027    !Config Key   = LEAF_FRAC_HVC
2028    !Config Desc  = parameters for herbivore activity
2029    !Config If    = OK_STOMATE
2030    !Config Def   = 0.33
2031    !Config Help  =
2032    !Config Units = [-]
2033    CALL getin_p('LEAF_FRAC_HVC',leaf_frac_hvc)
2034    !
2035    !Config Key   = TLONG_REF_MAX
2036    !Config Desc  = maximum reference long term temperature
2037    !Config If    = OK_STOMATE
2038    !Config Def   = 303.1
2039    !Config Help  =
2040    !Config Units = [K] 
2041    CALL getin_p('TLONG_REF_MAX',tlong_ref_max)
2042    !
2043    !Config Key   = TLONG_REF_MIN
2044    !Config Desc  = minimum reference long term temperature
2045    !Config If    = OK_STOMATE
2046    !Config Def   = 253.1
2047    !Config Help  =
2048    !Config Units = [K] 
2049    CALL getin_p('TLONG_REF_MIN',tlong_ref_min)
2050    !
2051    !Config Key   = NCD_MAX_YEAR
2052    !Config Desc  =
2053    !Config If    = OK_STOMATE
2054    !Config Def   = 3.
2055    !Config Help  = NCD : Number of Chilling Days
2056    !Config Units = [days]
2057    CALL getin_p('NCD_MAX_YEAR',ncd_max_year)
2058    !
2059    !Config Key   = GDD_THRESHOLD
2060    !Config Desc  =
2061    !Config If    = OK_STOMATE
2062    !Config Def   = 5.
2063    !Config Help  = GDD : Growing-Degree-Day
2064    !Config Units = [days]
2065    CALL getin_p('GDD_THRESHOLD',gdd_threshold)
2066    !
2067    !Config Key   = GREEN_AGE_EVER
2068    !Config Desc  =
2069    !Config If    = OK_STOMATE
2070    !Config Def   = 2.
2071    !Config Help  =
2072    !Config Units = [-] 
2073    CALL getin_p('GREEN_AGE_EVER',green_age_ever)
2074    !
2075    !Config Key   = GREEN_AGE_DEC
2076    !Config Desc  =
2077    !Config If    = OK_STOMATE
2078    !Config Def   = 0.5
2079    !Config Help  =
2080    !Config Units = [-]
2081    CALL getin_p('GREEN_AGE_DEC',green_age_dec)
2082
2083  END SUBROUTINE config_stomate_parameters
2084
2085!! ================================================================================================================================
2086!! SUBROUTINE   : config_dgvm_parameters
2087!!
2088!>\BRIEF        This subroutine reads in the configuration file all the parameters
2089!! needed when the DGVM model is activated (ie : when ok_dgvm is set to true).
2090!!
2091!! DESCRIPTION  : None
2092!!
2093!! RECENT CHANGE(S): None
2094!!
2095!! MAIN OUTPUT VARIABLE(S):
2096!!
2097!! REFERENCE(S) :
2098!!
2099!! FLOWCHART    :
2100!! \n
2101!_ ================================================================================================================================
2102
2103  SUBROUTINE config_dgvm_parameters   
2104
2105    IMPLICIT NONE
2106
2107    !! 0. Variables and parameters declaration
2108
2109    !! 0.4 Local variables
2110
2111    !_ ================================================================================================================================   
2112
2113    !-
2114    ! establish parameters
2115    !-
2116    !
2117    !Config Key   = ESTAB_MAX_TREE
2118    !Config Desc  = Maximum tree establishment rate
2119    !Config If    = OK_DGVM
2120    !Config Def   = 0.12
2121    !Config Help  =
2122    !Config Units = [-]   
2123    CALL getin_p('ESTAB_MAX_TREE',estab_max_tree)
2124    !
2125    !Config Key   = ESTAB_MAX_GRASS
2126    !Config Desc  = Maximum grass establishment rate
2127    !Config If    = OK_DGVM
2128    !Config Def   = 0.12
2129    !Config Help  =
2130    !Config Units = [-] 
2131    CALL getin_p('ESTAB_MAX_GRASS',estab_max_grass)
2132    !
2133    !Config Key   = ESTABLISH_SCAL_FACT
2134    !Config Desc  =
2135    !Config If    = OK_DGVM
2136    !Config Def   = 5.
2137    !Config Help  =
2138    !Config Units = [-]
2139    CALL getin_p('ESTABLISH_SCAL_FACT',establish_scal_fact)
2140    !
2141    !Config Key   = MAX_TREE_COVERAGE
2142    !Config Desc  =
2143    !Config If    = OK_DGVM
2144    !Config Def   = 0.98
2145    !Config Help  =
2146    !Config Units = [-]
2147    CALL getin_p('MAX_TREE_COVERAGE',max_tree_coverage)
2148    !
2149    !Config Key   = IND_0_ESTAB
2150    !Config Desc  =
2151    !Config If    = OK_DGVM
2152    !Config Def   = 0.2
2153    !Config Help  =
2154    !Config Units = [-] 
2155    CALL getin_p('IND_0_ESTAB',ind_0_estab)
2156
2157    !-
2158    ! light parameters
2159    !-
2160    !
2161    !Config Key   = ANNUAL_INCREASE
2162    !Config Desc  = for diagnosis of fpc increase, compare today's fpc to last year's maximum (T) or to fpc of last time step (F)?
2163    !Config If    = OK_DGVM
2164    !Config Def   = y
2165    !Config Help  =
2166    !Config Units = [FLAG]
2167    CALL getin_p('ANNUAL_INCREASE',annual_increase)
2168    !
2169    !Config Key   = MIN_COVER
2170    !Config Desc  = For trees, minimum fraction of crown area occupied
2171    !Config If    = OK_DGVM
2172    !Config Def   = 0.05
2173    !Config Help  =
2174    !Config Units = [-] 
2175    CALL getin_p('MIN_COVER',min_cover)
2176
2177    !-
2178    ! pftinout parameters
2179    !
2180    !Config Key   = IND_0
2181    !Config Desc  = initial density of individuals
2182    !Config If    = OK_DGVM
2183    !Config Def   = 0.02
2184    !Config Help  =
2185    !Config Units = [-] 
2186    CALL getin_p('IND_0',ind_0)
2187    !
2188    !Config Key   = MIN_AVAIL
2189    !Config Desc  = minimum availability
2190    !Config If    = OK_DGVM
2191    !Config Def   = 0.01
2192    !Config Help  =
2193    !Config Units = [-] 
2194    CALL getin_p('MIN_AVAIL',min_avail)
2195    !
2196    !Config Key   = RIP_TIME_MIN
2197    !Config Desc  =
2198    !Config If    = OK_DGVM
2199    !Config Def   = 1.25
2200    !Config Help  =
2201    !Config Units = [year] 
2202    CALL getin_p('RIP_TIME_MIN',RIP_time_min)
2203    !
2204    !Config Key   = NPP_LONGTERM_INIT
2205    !Config Desc  =
2206    !Config If    = OK_DGVM
2207    !Config Def   = 10.
2208    !Config Help  =
2209    !Config Units = [gC/m^2/year]
2210    CALL getin_p('NPP_LONGTERM_INIT',npp_longterm_init)
2211    !
2212    !Config Key   = EVERYWHERE_INIT
2213    !Config Desc  =
2214    !Config If    = OK_DGVM
2215    !Config Def   = 0.05
2216    !Config Help  =
2217    !Config Units = [-]
2218    CALL getin_p('EVERYWHERE_INIT',everywhere_init)
2219
2220
2221  END SUBROUTINE config_dgvm_parameters
2222
2223
2224!! ================================================================================================================================
2225!! FUNCTION   : get_printlev
2226!!
2227!>\BRIEF        Read global PRINTLEV parmeter and local PRINTLEV_modname
2228!!
2229!! DESCRIPTION  : The first time this function is called the parameter PRINTLEV is read from run.def file.
2230!!                It is stored in the variable named printlev which is declared in constantes_var.f90. printlev
2231!!                can be accesed each module in ORCHIDEE which makes use of constantes_var module.
2232!!
2233!!                This function also reads the parameter PRINTLEV_modname for run.def file. modname is the
2234!!                intent(in) character string to this function. If the variable is set in run.def file, the corresponding
2235!!                value is returned. Otherwise the value of printlev is returnd as default.
2236!!
2237!! RECENT CHANGE(S): None
2238!!
2239!! MAIN OUTPUT VARIABLE(S): The local output level for the module set as intent(in) argument.
2240!!
2241!! REFERENCE(S) :
2242!!
2243!! FLOWCHART    :
2244!! \n
2245!_ ================================================================================================================================
2246
2247  FUNCTION get_printlev ( modname )
2248
2249    !! 0.1 Input arguments
2250    CHARACTER(LEN=*), INTENT(IN) :: modname
2251
2252    !! 0.2 Returned variable
2253    INTEGER(i_std)               :: get_printlev
2254
2255    !! 0.3 Local variables
2256    LOGICAL, SAVE :: first=.TRUE.
2257
2258    !_ ================================================================================================================================
2259
2260    !! 1.0  Read the global PRINTLEV from run.def. This is only done at first call to this function.
2261    IF (first) THEN
2262       !Config Key   = PRINTLEV
2263       !Config Desc  = Print level for text output
2264       !Config If    =
2265       !Config Help  = Possible values are:
2266       !Config         0    No output,
2267       !Config         1    Minimum writing for long simulations,
2268       !Config         2    More basic information for long simulations,
2269       !Config         3    First debug level,
2270       !Config         4    Higher debug level
2271       !Config Def   = 2
2272       !Config Units = [0, 1, 2, 3, 4]
2273       ! Default value is set in constantes_var
2274       CALL getin_p('PRINTLEV',printlev)
2275       first=.FALSE.
2276
2277       !Config Key   = PRINTLEV_modname
2278       !Config Desc  = Specific print level of text output for the module "modname". Default as PRINTLEV.
2279       !Config Def   = PRINTLEV
2280       !Config If    =
2281       !Config Help  = Use this option to activate a different level of text output
2282       !Config         for a specific module. This can be activated for several modules
2283       !Config         at the same time. Use for example PRINTLEV_sechiba.
2284       !Config Units = [0, 1, 2, 3, 4]
2285    END IF
2286
2287    ! Set default value as the standard printlev
2288    get_printlev=printlev
2289    ! Read optional value from run.def file
2290    CALL getin_p('PRINTLEV_'//modname, get_printlev)
2291
2292  END FUNCTION get_printlev
2293
2294
2295END MODULE constantes
Note: See TracBrowser for help on using the repository browser.