[6455] | 1 | # $HeadURL: svn://forge.ipsl.jussieu.fr/orchidee/trunk/ORCHIDEE_OL/OOL_SEC_STO_FG2nd/PARAM/orchidee.def $ |
---|
| 2 | # $Date: 2023-05-12 13:04:57 +0200 (ven. 12 mai 2023) $ |
---|
| 3 | # $Revision: 8000 $ |
---|
[6084] | 4 | |
---|
[4979] | 5 | #************************************************************************** |
---|
[6084] | 6 | # Suggested parameter set up for ORCHIDEE trunk |
---|
[396] | 7 | # |
---|
[3528] | 8 | # See all optional parameters in modipsl/modeles/ORCHIDEE/orchidee.default |
---|
[4979] | 9 | # |
---|
| 10 | # Note when running with libIGCM, variables set here to _AUTOBLOCKER_ can |
---|
| 11 | # not be changed in run.def otherwise the simulation will stop. Variables |
---|
| 12 | # set to _AUTO_ can be changed from the comp.card but they can also be |
---|
| 13 | # changed directly in the run.def. A value set in run.def can never be |
---|
| 14 | # changed again by libIGCM. |
---|
[793] | 15 | #************************************************************************** |
---|
[2410] | 16 | |
---|
[6084] | 17 | # Parameters that should be removed once these issues are resolved |
---|
| 18 | #************************************************************************** |
---|
[4979] | 19 | |
---|
[6084] | 20 | # From a scientific point of view hydraulic architecture is much nicer |
---|
| 21 | # than a direct dependency on soil water. Julien is working on both |
---|
| 22 | # the nummerics and the science of the hydraulic architecture. To |
---|
| 23 | # avoid further delays with milestone 4.1 we decided to change the |
---|
| 24 | # settings. |
---|
| 25 | ENERGY_CONTROL=2 |
---|
| 26 | |
---|
[6455] | 27 | # Prescribe forest management to facilitate the interpretation of the |
---|
| 28 | # results when comparing versions. |
---|
[6084] | 29 | OK_READ_FM_MAP=n |
---|
| 30 | |
---|
[6455] | 31 | # Hack to use a higher threshold for small lai in the calculation of kB_m1 in condveg |
---|
| 32 | HACK_KB_M1 = y |
---|
[6084] | 33 | |
---|
[4979] | 34 | # Parameters related to the restart file and the start date |
---|
[793] | 35 | #************************************************************************** |
---|
[4979] | 36 | # Name of restart file for sechiba part of the model |
---|
| 37 | # (default NONE) |
---|
[2410] | 38 | SECHIBA_restart_in = _AUTOBLOCKER_ |
---|
[793] | 39 | |
---|
[4979] | 40 | # Name of restart file for stomate part of the model |
---|
| 41 | # (default NONE) |
---|
[2410] | 42 | STOMATE_RESTART_FILEIN = _AUTOBLOCKER_ |
---|
[793] | 43 | |
---|
[4979] | 44 | |
---|
| 45 | # Parameters related to the diagnostic output files |
---|
| 46 | #************************************************************************** |
---|
[2410] | 47 | # Use XIOS for writing diagnostics file |
---|
[6084] | 48 | # (defualt y) |
---|
| 49 | XIOS_ORCHIDEE_OK = y |
---|
[793] | 50 | |
---|
| 51 | |
---|
[4979] | 52 | # Parameters related to debugging and error checking |
---|
| 53 | #************************************************************************** |
---|
| 54 | # Mass balance and surface area checks |
---|
| 55 | # 1 - is recommended when running global long-term simulations. |
---|
| 56 | # Under this option, mass and area conservation are checked for all |
---|
| 57 | # biogeochemical processes but only at the highest level thus stomate.f90 |
---|
| 58 | # and stomate_lpj.f90. Although these checks are not very expensive in terms |
---|
| 59 | # of computer time, skipping the numerous lower level checks is expected to |
---|
| 60 | # save some time. Under this option the mass balance error is only written |
---|
| 61 | # to the history file. No information is provided in which subroutine the |
---|
| 62 | # problem occurred. |
---|
| 63 | # 2 - is recommended when developing and testing the model. Now the mass |
---|
| 64 | # balance is explicitly checked in stomate.f90, stomate_lpj.f90 and all |
---|
| 65 | # its subroutines. Under this option the mass balance error is written to |
---|
| 66 | # the history file and if the mass balance is not closed, the warning |
---|
| 67 | # message will indicate in which subroutine the problem likely |
---|
| 68 | # originated. |
---|
| 69 | # 3 - is recommended when having a problem with mass balance closure. The |
---|
| 70 | # mass balance is explicitly checked in stomate.f90, stomate_lpj.f90 and all |
---|
| 71 | # its subroutines. If a mass balance occurs, the model is stopped. |
---|
| 72 | # (default) |
---|
[6084] | 73 | ERR_ACT = 3 |
---|
[4979] | 74 | |
---|
| 75 | # If the model crashes and it is not clear in which subroutine the crash |
---|
| 76 | # occurs one could increase the print level for text output in out_orchidee_* |
---|
[3430] | 77 | # PRINTLEV = 1 : some output in initialization phase |
---|
| 78 | # PRINTLEV = 2 : more output in initialization phase and printing of new date |
---|
| 79 | # PRINTLEV = 3, 4,.. : more output also at each timestep |
---|
[4979] | 80 | # (default 2) |
---|
[6394] | 81 | PRINTLEV=1 |
---|
[3430] | 82 | |
---|
[4979] | 83 | # If the model crashes on an ipslerr there is a fair change that the pixel and |
---|
| 84 | # PFT where the crash occur will be printed in the out_orchidee_* file. Try to |
---|
| 85 | # re-run the model for a single pixel or a small region (for example, |
---|
| 86 | # 3 by 3). LIMIT_EAST, LIMIT_NORTH, LIMIT_SOUTH and LIMIT_WEST should |
---|
| 87 | # be set in the run.def |
---|
| 88 | |
---|
| 89 | # If it is known in which subroutine the model crashes, the local print level |
---|
| 90 | # of that subroutine could be increased. Additional write statements will then |
---|
| 91 | # be activated. These write statements print the values of specific variables. |
---|
| 92 | # (default 2) |
---|
| 93 | #PRINTLEV_sapiens_lcchange = 4 |
---|
| 94 | |
---|
| 95 | # When a local print level is set, the grid and pft number for which the |
---|
| 96 | # values needs to written should be set. This setting reduces reduces |
---|
| 97 | # the size and enhances the readability of the out_orchidee_* files. |
---|
| 98 | TEST_GRID = 1 |
---|
| 99 | TEST_PFT = 2 |
---|
| 100 | |
---|
| 101 | |
---|
[2410] | 102 | # Hydrology parameters |
---|
| 103 | #************************************************************************** |
---|
[6414] | 104 | # Activate river routing |
---|
| 105 | # default = y |
---|
| 106 | RIVER_ROUTING = _AUTO_ |
---|
[793] | 107 | |
---|
[6414] | 108 | # Choice of routing method: standard or simple |
---|
| 109 | ROUTING_METHOD = _AUTO_ |
---|
| 110 | |
---|
[3528] | 111 | # Activate creation of river_desc.nc file |
---|
[4979] | 112 | # RIVER_DESC will be activated only the first execution in the simulation |
---|
[3528] | 113 | RIVER_DESC = _AUTO_ |
---|
| 114 | |
---|
[4979] | 115 | # SOILTYPE_CLASSIF : Type of classification used for the map of soil types |
---|
| 116 | # (default zobler) |
---|
[6653] | 117 | SOILTYPE_CLASSIF = usda |
---|
[3015] | 118 | |
---|
[6653] | 119 | # SOILALB_FILE : File to read bare soil albedo |
---|
| 120 | SOILALB_FILE=soilcolor.nc |
---|
[3093] | 121 | |
---|
[4979] | 122 | # Parameters related to vegetation map |
---|
[2410] | 123 | #************************************************************************** |
---|
[4979] | 124 | # Update vegetation frequency |
---|
| 125 | # (default 0Y) |
---|
| 126 | VEGET_UPDATE = _AUTO_ |
---|
[793] | 127 | |
---|
[4979] | 128 | # Read lai map |
---|
| 129 | # (default n) |
---|
[6084] | 130 | LAI_MAP = n |
---|
[793] | 131 | |
---|
[4979] | 132 | # Prescribed vegetation |
---|
| 133 | # (default n) |
---|
| 134 | IMPOSE_VEG = n |
---|
[793] | 135 | |
---|
| 136 | |
---|
[4979] | 137 | # Parameters related to surface and thermal physical properties |
---|
| 138 | #************************************************************************ |
---|
| 139 | # ROUGH_DYN : Account for a dynamic roughness height (activation of Su et |
---|
| 140 | # al. parametrization) |
---|
| 141 | # (default y) |
---|
| 142 | ROUGH_DYN = y |
---|
[3093] | 143 | |
---|
[4979] | 144 | # OK_FREEZE : Activate the complet soil freezing scheme |
---|
| 145 | # (default y) |
---|
| 146 | OK_FREEZE = y |
---|
[3093] | 147 | |
---|
[4979] | 148 | # OK_EXPLICITSNOW : Activate explict snow scheme |
---|
| 149 | # (default y) |
---|
| 150 | OK_EXPLICITSNOW = y |
---|
[3418] | 151 | |
---|
[6653] | 152 | # NSNOW : Number of snow levels. Possible values are 3 or 12 |
---|
| 153 | NSNOW=12 |
---|
| 154 | |
---|
[4979] | 155 | # Controls several settings related to the energy budget and water stress. |
---|
| 156 | # 1 - DEFAULT uses the enerbil module in combination with the hydraulic |
---|
[6455] | 157 | # architecture (ok_hydrol_arch and ok_gs_feedback true, while |
---|
[4979] | 158 | # ok_mleb and ok_impose_canopy_structure are set to false). |
---|
| 159 | # 2 - option to use enerbil module and original water stress |
---|
| 160 | # (not hydraulic architecture) |
---|
| 161 | # 3 - The energy budget is calculated using the multi-layer energy scheme |
---|
| 162 | # with a single layer: ok_hydrol_arch, ok_gs_feedback, ok_impose_canopy_structure |
---|
| 163 | # and and ok_mleb all TRUE, but the energy budget is calculated for a single |
---|
| 164 | # layer (jnlvls=1,jnlvls_under=0,jnlvls_canopy=1,jnlvls_over=0). No mleb output, |
---|
| 165 | # ok_mleb_history_file is set to FALSE. |
---|
| 166 | # 4 - multi-layer energy budget: ok_hydrol_arch, ok_gs_feedback and ok_mleb all TRUE. |
---|
| 167 | # ok_impose_canopy_structure is False, and the energy budget is calculated for |
---|
[6455] | 168 | # multiple layers (jnlvls=29,jnlvls_under=10,jnlvls_canopy=10,jnlvls_over=9). |
---|
[4979] | 169 | # No mleb output, ok_mleb_history_file is set to FALSE. |
---|
| 170 | # 5 - user specific: user specific settings for these controls and layers will need |
---|
| 171 | # to be defined in the run.def by the user. |
---|
[6084] | 172 | #ENERGY_CONTROL=1 |
---|
[3418] | 173 | |
---|
[6455] | 174 | |
---|
[4979] | 175 | # Carbon related parameters |
---|
| 176 | #************************************************************************** |
---|
[6455] | 177 | # Analytic spinup |
---|
| 178 | # (default n) |
---|
| 179 | SPINUP_ANALYTIC = n |
---|
| 180 | SPINUP_PERIOD = 10 |
---|
| 181 | |
---|
[6084] | 182 | # Activate Stomate component |
---|
| 183 | # (default y) |
---|
| 184 | STOMATE_OK_STOMATE = _AUTOBLOCKER_ |
---|
| 185 | |
---|
[4979] | 186 | # As an alternative to DO_WOOD_HARVEST forest management can be accounted |
---|
[6455] | 187 | # for. DO_WOOD_HARVEST is not yet available in ORCHIDEE v4. |
---|
| 188 | # OK_READ_FM_MAP=y |
---|
[3418] | 189 | |
---|
[4173] | 190 | # Parameters related to nitrogen cycle |
---|
| 191 | #************************************************************************** |
---|
[4979] | 192 | # Update nitrogen input maps |
---|
| 193 | # (default 0Y) |
---|
[4173] | 194 | NINPUT_UPDATE = _AUTO_ |
---|
| 195 | |
---|
| 196 | # Impose the CN ratio of leaves |
---|
| 197 | STOMATE_IMPOSE_CN = _AUTO_ |
---|
| 198 | |
---|
[4979] | 199 | # Set STOMATE_READ_CN=y to read cnleaf_map.nc file if STOMATE_IMPOSE_CN=y |
---|
[4177] | 200 | STOMATE_READ_CN = n |
---|
[4173] | 201 | CNLEAF_VAR = CN_LEAF |
---|
| 202 | |
---|
[2769] | 203 | |
---|
[4173] | 204 | # File and variable name for nitrogen input files |
---|
| 205 | #************************************************************************** |
---|
| 206 | Nammonium_FILE = ndep_nhx.nc |
---|
| 207 | Nammonium_VAR = nhx |
---|
| 208 | |
---|
| 209 | Nnitrate_FILE = ndep_noy.nc |
---|
| 210 | Nnitrate_VAR = noy |
---|
| 211 | |
---|
| 212 | Nfert_FILE = NONE |
---|
| 213 | Nfert_VAR = nfer |
---|
| 214 | |
---|
| 215 | Nmanure_FILE = NONE |
---|
| 216 | Nmanure_VAR = Nmanure |
---|
| 217 | |
---|
| 218 | Nfert_cropland_FILE = nfert_cropland.nc |
---|
| 219 | Nfert_cropland_VAR = nfer |
---|
| 220 | |
---|
| 221 | Nmanure_cropland_FILE = nmanure_cropland.nc |
---|
| 222 | Nmanure_cropland_VAR = Nmanure |
---|
| 223 | |
---|
| 224 | Nfert_pasture_FILE = nfert_pasture.nc |
---|
| 225 | Nfert_pasture_VAR = Nfer |
---|
| 226 | |
---|
| 227 | Nmanure_pasture_FILE = nmanure_pasture.nc |
---|
| 228 | Nmanure_pasture_VAR = Nmanure |
---|
| 229 | |
---|
| 230 | Nbnf_FILE= bnf.nc |
---|
| 231 | Nbnf_VAR= BNF_MGN_PERM2_PERYR |
---|
| 232 | |
---|
| 233 | |
---|
[4979] | 234 | # Parameters describing the canopy structure |
---|
| 235 | #************************************************************************** |
---|
| 236 | # Number of diameter classes - enables to create complex canopies |
---|
| 237 | # (default = 1) |
---|
| 238 | NCIRC = 3 |
---|
| 239 | |
---|
| 240 | |
---|
| 241 | # Parameters related to chemistry bvoc |
---|
| 242 | #************************************************************************ |
---|
| 243 | # CHEMISTRY_OK_BVOC : Activate chemistry |
---|
| 244 | CHEMISTRY_BVOC = _AUTO_: DEFAULT=n |
---|
| 245 | |
---|
| 246 | # CO2 FOR BVOC - WILKINSON : CO2 inhibition effect for isoprene based on Wilkinson approach? |
---|
| 247 | CO2_FOR_BVOC_WILKINSON=n |
---|
| 248 | |
---|
| 249 | # CO2 FOR BVOC - POSSELL : CO2 inhibition effect for isoprene based on Possell approach? |
---|
| 250 | CO2_FOR_BVOC_POSSELL=n |
---|
[6455] | 251 | |
---|