Changeset 4979 for CONFIG/UNIFORM/v6/LMDZOR_v6.4/GENERAL/PARAM/orchidee.def

Timestamp:

02/25/20 13:52:12 (5 years ago)

Author:

jgipsl

Message:

Added new LMDZOR configuration: This configuration is different from LMDZOR_v6.3 only for the version of ORCHIDEE. LMDZOR_v6.4 is adapted for ORCHIDEE-CN-CAN which soon will be the ORCHIDEE/trunk wheras LMDZOR_v6.3 is adapted for ORCHIDEE_3.

Location:

CONFIG/UNIFORM/v6/LMDZOR_v6.4

Files:

• . (copied) (copied from CONFIG/UNIFORM/v6/LMDZOR_v6.3)
• GENERAL/PARAM/orchidee.def (copied) (copied from CONFIG/UNIFORM/v6/LMDZOR_v6.3/GENERAL/PARAM/orchidee.def_CWRR) (3 diffs)

CONFIG/UNIFORM/v6/LMDZOR_v6.4/GENERAL/PARAM/orchidee.def

@@ +1 @@
+#**************************************************************************
+#           Suggested parameter set up for ORCHIDEE
 #
-#**************************************************************************
-#                    Parameter file list for ORCHIDEE
 # See all optional parameters in modipsl/modeles/ORCHIDEE/orchidee.default
-#**************************************************************************
-
-# Input and output
-#**************************************************************************
-# Restart the time from the GCM. 
-# default = n
-SECHIBA_reset_time = y
-
-# Name of restart to read for initial conditions
-# default = NONE
+#
+# Note when running with libIGCM, variables set here to _AUTOBLOCKER_ can 
+# not be changed in run.def otherwise the simulation will stop. Variables 
+# set to _AUTO_ can be changed from the comp.card but they can also be 
+# changed directly in the run.def. A value set in run.def can never be
+# changed again by libIGCM. 
+#**************************************************************************
+
+
+# Parameters related to the restart file and the start date
+#**************************************************************************
+# Name of restart file for sechiba part of the model
+# (default NONE)
 SECHIBA_restart_in = _AUTOBLOCKER_
 
-# Name of restart to read for initial conditions of STOMATE
-# default = NONE
+# Name of restart file for stomate part of the model
+# (default NONE)
 STOMATE_RESTART_FILEIN = _AUTOBLOCKER_
 
+
+# Parameters related to the diagnostic output files
+#**************************************************************************
 # Use XIOS for writing diagnostics file
-# default = y
-XIOS_ORCHIDEE_OK = y
-
-# Writefrequency in seconds in sechiba_history.nc
-# default = 0
-WRITE_STEP = 0
-
-# Print level for text output in out_orchidee_*
+# (default y)
+# All IOIPSL output are deactivated if XIOS_ORCHIDEE_OK=y
+XIOS_ORCHIDEE_OK = y 
+
+
+# Parameters related to debugging and error checking
+#**************************************************************************
+# Mass balance and surface area checks
+# 1 - is recommended when running global long-term simulations.
+# Under this option, mass and area conservation are checked for all
+# biogeochemical processes but only at the highest level thus stomate.f90
+# and stomate_lpj.f90. Although these checks are not very expensive in terms
+# of computer time, skipping the numerous lower level checks is expected to
+# save some time. Under this option the mass balance error is only written
+# to the history file. No information is provided in which subroutine the
+# problem occurred.
+# 2 - is recommended when developing and testing the model. Now the mass
+# balance is explicitly checked in stomate.f90, stomate_lpj.f90 and all
+# its subroutines. Under this option the mass balance error is written to
+# the history file and if the mass balance is not closed, the warning
+# message will indicate in which subroutine the problem likely
+# originated.
+# 3 - is recommended when having a problem with mass balance closure. The
+# mass balance is explicitly checked in stomate.f90, stomate_lpj.f90 and all
+# its subroutines. If a mass balance occurs, the model is stopped.
+# (default)
+ERR_ACT = 3
+
+# If the model crashes and it is not clear in which subroutine the crash
+# occurs one could increase the print level for text output in out_orchidee_*
 # PRINTLEV = 1 : some output in initialization phase
 # PRINTLEV = 2 : more output in initialization phase and printing of new date
 # PRINTLEV = 3, 4,.. : more output also at each timestep
+# (default 2)
 PRINTLEV=2
 
+# If the model crashes on an ipslerr there is a fair change that the pixel and
+# PFT where the crash occur will be printed in the out_orchidee_* file. Try to
+# re-run the model for a single pixel or a small region (for example,
+# 3 by 3). LIMIT_EAST, LIMIT_NORTH, LIMIT_SOUTH and LIMIT_WEST should
+# be set in the run.def 
+
+# If it is known in which subroutine the model crashes, the local print level
+# of that subroutine could be increased. Additional write statements will then
+# be activated. These write statements print the values of specific variables.
+# (default 2)
+#PRINTLEV_sapiens_lcchange = 4
+
+# When a local print level is set, the grid and pft number for which the
+# values needs to written should be set. This setting reduces reduces
+# the size and enhances the readability of the out_orchidee_* files.  
+TEST_GRID = 1
+TEST_PFT = 2
+
+
 # Hydrology parameters
 #**************************************************************************
-# Activate the multi-layer diffusion scheme adapted from CWRR
-# default = y
-HYDROL_CWRR = y
-
-# Activate river routing
-# default = y
-RIVER_ROUTING = y
+# Activate river routing. This version has not yet been tested with river routing.
+# (default y)
+RIVER_ROUTING = y 
 
 # Activate creation of river_desc.nc file
-# RIVER_DESC will be activated only the first execution in the simulation. 
+# RIVER_DESC will be activated only the first execution in the simulation
 RIVER_DESC = _AUTO_
 
-# SOILTYPE_CLASSIF :  Type of classification used for the map of soil types (default zobler)
+# SOILTYPE_CLASSIF :  Type of classification used for the map of soil types
+# (default zobler)
 SOILTYPE_CLASSIF = zobler 
-SOILALB_FILE     = soils_param.nc
-SOILCLASS_FILE   = soils_param.nc
-
-# DO_RSOIL : Reduce soil evaporation with a soil resistance (default n)
-DO_RSOIL = n
-
-# To make AVAN and NVAN constant with depth (default CWRR_AKS_A0=0, CWRR_AKS_POWER=0, CWRR_NKS_A0=0, CWRR_NKS_POWER=0 )
-CWRR_AKS_A0 = 0.
-CWRR_AKS_POWER = 0.
-CWRR_NKS_A0 = 0.
-CWRR_NKS_POWER = 0.
-
-#**************************************************************************
+
+
+# Parameters related to vegetation map
+#**************************************************************************
+# Update vegetation frequency
+# (default 0Y)
+VEGET_UPDATE = _AUTO_
+
+# Read lai map
+# (default n)
+LAI_MAP = _AUTO_ : DEFAULT = n
+
+# Prescribed vegetation
+# (default n)
+IMPOSE_VEG = n
+
+
+# Parameters related to surface and thermal physical properties
+#************************************************************************
+# ROUGH_DYN : Account for a dynamic roughness height (activation of Su et
+# al. parametrization)
+# (default y)
+ROUGH_DYN = y
+
+# OK_FREEZE :  Activate the complet soil freezing scheme
+# (default y)
+OK_FREEZE = y
+
+# OK_EXPLICITSNOW :  Activate explict snow scheme
+# (default y)
+OK_EXPLICITSNOW = y
+
+# Controls several settings related to the energy budget and water stress.
+# 1 - DEFAULT uses the enerbil module in combination with the hydraulic
+# architecture (ok_hydrol_arch and ok_gs_feedback true, while   
+# ok_mleb and ok_impose_canopy_structure are set to false).
+# 2 - option to use enerbil module and original water stress
+# (not hydraulic architecture)
+# 3 - The energy budget is calculated using the multi-layer energy scheme
+# with a single layer: ok_hydrol_arch, ok_gs_feedback, ok_impose_canopy_structure
+# and and ok_mleb all TRUE, but the energy budget is calculated for a single
+# layer (jnlvls=1,jnlvls_under=0,jnlvls_canopy=1,jnlvls_over=0). No mleb output,
+# ok_mleb_history_file is set to  FALSE.
+# 4 - multi-layer energy budget: ok_hydrol_arch, ok_gs_feedback and ok_mleb all TRUE.
+# ok_impose_canopy_structure is False, and the energy budget is calculated for
+# multiple layers (jnlvls=29,jnlvls_under=10,jnlvls_canopy=10,jnlvls_over=9).
+# No mleb output, ok_mleb_history_file is set to  FALSE.
+# 5 - user specific: user specific settings for these controls and layers will need
+# to be defined in the run.def by the user.
+ENERGY_CONTROL=1
+
+# Carbon related parameters
+#**************************************************************************
+# Analytic spinup
+# (default n)
+SPINUP_ANALYTIC = n 
+SPINUP_PERIOD = 10
+
 # Activate Stomate component
+# (default y)
 STOMATE_OK_STOMATE = _AUTOBLOCKER_ 
 
-# Activate calculations of CO2 according to Farqhuar and Ball
-# default = y
-STOMATE_OK_CO2 = y
-
-# CDRAG_from_GCM: Use cdrag coefficient from gcm (default=y if coupled to LMDZ)
-CDRAG_from_GCM = y
-
-# Interception reservoir coefficient.
-# default = 0.1
-SECHIBA_QSINT = 0.02
-
-# Use 15 PFTs. Note that PFTmap.nc must contains 15PFTs.
-# PFT=10, 14 and 15 share the same parametrization except for VCMAX and LAI_MAX specified below.  
-NVM=15
-PFT_TO_MTC=1,2,3,4,5,6,7,8,9,10,11,12,13,10,10
-PFT_NAME__10='temperate C3           grass      '
-PFT_NAME__14='tropical  C3           grass      '
-PFT_NAME__15='boreal    C3           grass      '
-
-# STRESS_GS, STRESS_GM, STRESS_VCMAX : Stress on GS, GM and VCMAX
-# Default STRESS_GS=1, STRESS_GM=1, STRESS_VCMAX=1 for all PFTs
-STRESS_GS=1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.
-STRESS_GM=1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.
-STRESS_VCMAX=1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.
-
-# Optimized parameters for albedo (default values)
-ALB_LEAF_NIR = 0.,0.23,0.18,0.18,0.20,0.24,0.15,0.26,0.20,0.24,0.27,0.28,0.26,0.24,0.24
-ALB_LEAF_VIS = 0.,0.04,0.04,0.04,0.04,0.03,0.03,0.03,0.03,0.06,0.06,0.06,0.06,0.06,0.06
-SNOWA_AGED_NIR = 0.50,0,0,0.10,0.37,0.08,0.16,0.17,0.27,0.44,0.44,0.44,0.44,0.44,0.44
-SNOWA_DEC_NIR = 0.13,0,0,0.10,0.10,0.16,0.04,0.07,0.08,0.12,0.12,0.12,0.12,0.12,0.12
-SNOWA_AGED_VIS = 0.74,0,0,0.08,0.24,0.07,0.18,0.18,0.33,0.57,0.57,0.57,0.57,0.57,0.57
-SNOWA_DEC_VIS = 0.21,0,0,0.14,0.08,0.17,0.05,0.06,0.09,0.15,0.15,0.15,0.15,0.15,0.15
-
-# Optimized parameters for photosynthetic capacity 
-# Default values except for PFT=15
-VCMAX25__02 = 45.0
-VCMAX25__03 = 45.0
-VCMAX25__04 = 35.0
-VCMAX25__05 = 40.0
-VCMAX25__06 = 50.0
-VCMAX25__07 = 45.0
-VCMAX25__08 = 35.0
-VCMAX25__09 = 35.0
-VCMAX25__10 = 50.0
-VCMAX25__11 = 50.0
-VCMAX25__12 = 60.0
-VCMAX25__13 = 60.0
-VCMAX25__14 = 50.0
-VCMAX25__15 = 40.0
-
-# Optimized parameters for LAI
-# Default values except for PFT=15
-LAI_MAX__02 = 7.0
-LAI_MAX__03 = 5.0
-LAI_MAX__04 = 5.0
-LAI_MAX__05 = 4.0
-LAI_MAX__06 = 5.0
-LAI_MAX__07 = 3.5
-LAI_MAX__08 = 4.0
-LAI_MAX__09 = 3.0
-LAI_MAX__10 = 2.5
-LAI_MAX__11 = 2.0
-LAI_MAX__12 = 5.0
-LAI_MAX__13 = 5.0
-LAI_MAX__14 = 2.5
-LAI_MAX__15 = 2.0
-
-# Activate downregulation for CO2
-DOWNREGULATION_CO2=y
-DOWNREGULATION_CO2_BASELEVEL=380.
-
-# Activate harvest of wood
-DO_WOOD_HARVEST=_AUTO_ : DEFAULT=y
-
-# Maintenance respiration slope C (default values except for PFT=14 and PFT=15)
-MAINT_RESP_SLOPE_C__02 = 0.12
-MAINT_RESP_SLOPE_C__03 = 0.12
-MAINT_RESP_SLOPE_C__04 = 0.16
-MAINT_RESP_SLOPE_C__05 = 0.16
-MAINT_RESP_SLOPE_C__06 = 0.16
-MAINT_RESP_SLOPE_C__07 = 0.25
-MAINT_RESP_SLOPE_C__08 = 0.25
-MAINT_RESP_SLOPE_C__09 = 0.25
-MAINT_RESP_SLOPE_C__10 = 0.16
-MAINT_RESP_SLOPE_C__11 = 0.12
-MAINT_RESP_SLOPE_C__12 = 0.16
-MAINT_RESP_SLOPE_C__13 = 0.12
-MAINT_RESP_SLOPE_C__14 = 0.12
-MAINT_RESP_SLOPE_C__15 = 0.25
-
-# Soil turnover (default CARBON_TAU=7, CARBON_TAU_IPASSOVE=300)
-CARBON_TAU_ISLOW = 7
-CARBON_TAU_IPASSIVE = 300
-
-# Take carbon from atmosphere if carbohydrate reserve too small, PFT dependent
-# Set to y for all pft except C4 grass (pft=11) and tropical C3 grass (pft=14)
-ALWAYS_INIT =  y, y, y, y, y, y, y, y, y, y, n, y, y, n, y
-
-# Do not use default fire module
-# Default FIRE_DISABLE=y
-FIRE_DISABLE = y
-
-# Parameters related to vegetation map
-#**************************************************************************
-# Update vegetation frequency
-# The maximum vegetation fractions will be update at this period
-# default = 1Y
-VEGET_UPDATE = _AUTO_
-
-# Read lai map (default n)
-LAI_MAP = _AUTO_ : DEFAULT = n
+# As an alternative to DO_WOOD_HARVEST forest management can be accounted
+# for. All forest will be managed by the same strategy unless a map is read
+# FM = 1 is unmanaged, FM = 2 is a thin and fell strategy, FM = 3 is a
+# coppice strategy and FM = 4 is a short rotation coppice strategy.
+# (default 2)
+FOREST_MANAGED_FORCED=2  
+
 
 # Parameters related to nitrogen cycle
 #**************************************************************************
-# Update nitrogen input maps (default 0Y)
+# Update nitrogen input maps
+# (default 0Y)
 NINPUT_UPDATE = _AUTO_
 
@@ -185 +173 @@
 STOMATE_IMPOSE_CN = _AUTO_
 
-# Set STOMATE_READ_CN=y to read leaf_cn.nc file if STOMATE_IMPOSE_CN=y
+# Set STOMATE_READ_CN=y to read cnleaf_map.nc file if STOMATE_IMPOSE_CN=y
 STOMATE_READ_CN = n
 CNLEAF_VAR = CN_LEAF
 
-# Parameters related to surface and thermal physical properties
-#************************************************************************
-# Read bare soil albedo from file with background MODIS data
-# Default ALB_BG_MODIS=y
-ALB_BG_MODIS = y
-ALB_BG_FILE  = alb_bg.nc
-
-# ROUGH_DYN : Account for a dynamic roughness height (activation of Su et al. parametrization) (default y)
-ROUGH_DYN=y
-
-# OK_FREEZE :  Activate the complet soil freezing scheme (default y)
-OK_FREEZE=y
-
-# DEPTH_MAX_T (m) :  Maximum depth of the soil thermodynamics (default 90)
-DEPTH_MAX_T=90
-
-# OK_EXPLICITSNOW :  Activate explict snow scheme (default y)
-OK_EXPLICITSNOW=y
- 
-# READ_REFTEMP :  Initialize soil temperature using climatological temperature (default y if OK_FREEZE=y)
-READ_REFTEMP=y
+
+# File and variable name for nitrogen input files
+#**************************************************************************
+Nammonium_FILE = ndep_nhx.nc
+Nammonium_VAR = nhx
+
+Nnitrate_FILE = ndep_noy.nc
+Nnitrate_VAR = noy
+
+Nfert_FILE = NONE
+Nfert_VAR = nfer
+
+Nmanure_FILE = NONE
+Nmanure_VAR = Nmanure
+
+Nfert_cropland_FILE = nfert_cropland.nc
+Nfert_cropland_VAR = nfer
+
+Nmanure_cropland_FILE = nmanure_cropland.nc
+Nmanure_cropland_VAR = Nmanure
+
+Nfert_pasture_FILE = nfert_pasture.nc
+Nfert_pasture_VAR = Nfer
+
+Nmanure_pasture_FILE = nmanure_pasture.nc
+Nmanure_pasture_VAR = Nmanure
+
+Nbnf_FILE= bnf.nc
+Nbnf_VAR= BNF_MGN_PERM2_PERYR
+
+
+# Parameters describing the canopy structure
+#**************************************************************************
+# Number of diameter classes - enables to create complex canopies
+# (default = 1)
+NCIRC = 3
+
+# Use recruiment. Note that this is an overall flag that makes use
+# of the variable RECRUITMENT_PFT which specifies whether a specific 
+# PFT grows recruits in canopy gaps. 
+OK_RECRUITMENT=y
+
+
+# Parameters that should be removed once these issues are resolved
+#**************************************************************************
+# This patch should be accepted or rejected and the flag should then
+# be removed
+LNVGRASSPATCH=y
 
 # Parameters related to chemistry bvoc 
@@ -222 +237 @@
 CO2_FOR_BVOC_POSSELL=n
 
-
-# Parameters below are different from default values of ORCHIDEE_2_0.
-# These values correspond to the previous default values and were used 
-# during developement of ORCHIDEE-CN
-#**************************************************************************
-# Growth respiration
-FRAC_GROWTHRESP__02 = 0.28
-FRAC_GROWTHRESP__03 = 0.28
-FRAC_GROWTHRESP__04 = 0.28
-FRAC_GROWTHRESP__05 = 0.28
-FRAC_GROWTHRESP__06 = 0.28
-FRAC_GROWTHRESP__07 = 0.28
-FRAC_GROWTHRESP__08 = 0.28
-FRAC_GROWTHRESP__09 = 0.28
-FRAC_GROWTHRESP__10 = 0.28
-FRAC_GROWTHRESP__11 = 0.28
-FRAC_GROWTHRESP__12 = 0.28
-FRAC_GROWTHRESP__13 = 0.28
-FRAC_GROWTHRESP__14 = 0.28
-FRAC_GROWTHRESP__15 = 0.28
-
-
-# Parameters below are different from default values of ORCHIDEE_2_0 
-# because they were tuned for the developements of ORCHIDEE-CN
-#**************************************************************************
-MOIAVAIL_ALWAYS_GRASS=0.4
-MOIAVAIL_ALWAYS_TREE=0.8
-TCST_SNOWA=3
-
-#PHENO_GDD_CRIT_C(K)
-PHENO_GDD_CRIT_C__00012=360
-PHENO_GDD_CRIT_C__00013=450
-
-#HUM_MIN_TIME(days)
-HUM_MIN_TIME__00012=50
-
-#GDD_SENESCENCE(K)
-GDD_SENESCENCE__00012=1600
-
-#ALLOC_MIN(nan)
-ALLOC_MIN__00007=0.35000
-ALLOC_MIN__00008=0.70
-ALLOC_MIN__00009=0.20
-
-
-# File and variable name for nitrogen input files
-# Following files are read if STOMATE_IMPOSE_CN=n
-#**************************************************************************
-Nammonium_FILE = ndep_nhx.nc
-Nammonium_VAR = nhx
-
-Nnitrate_FILE = ndep_noy.nc
-Nnitrate_VAR = noy
-
-Nfert_FILE = NONE
-Nfert_VAR = nfer
-
-Nmanure_FILE = NONE
-Nmanure_VAR = Nmanure
-
-Nfert_cropland_FILE = nfert_cropland.nc
-Nfert_cropland_VAR = nfer
-
-Nmanure_cropland_FILE = nmanure_cropland.nc
-Nmanure_cropland_VAR = Nmanure
-
-Nfert_pasture_FILE = nfert_pasture.nc
-Nfert_pasture_VAR = Nfer
-
-Nmanure_pasture_FILE = nmanure_pasture.nc
-Nmanure_pasture_VAR = Nmanure
-
-Nbnf_FILE= bnf.nc
-Nbnf_VAR= BNF_MGN_PERM2_PERYR
-
-