| 1 | ! ================================================================================================================================= |
|---|
| 2 | ! MODULE : grassland_constantes |
|---|
| 3 | ! |
|---|
| 4 | ! CONTACT : orchidee-help _at_ ipsl.jussieu.fr |
|---|
| 5 | ! |
|---|
| 6 | ! LICENCE : IPSL (2006) |
|---|
| 7 | ! This software is governed by the CeCILL licence see |
|---|
| 8 | ! ORCHIDEE/ORCHIDEE_CeCILL.LIC |
|---|
| 9 | ! |
|---|
| 10 | !>\BRIEF This module defined all constantes used in |
|---|
| 11 | !! grassland management module |
|---|
| 12 | !! |
|---|
| 13 | !!\n DESCRIPTION : None |
|---|
| 14 | !! |
|---|
| 15 | !! RECENT CHANGE(S) : None |
|---|
| 16 | !! |
|---|
| 17 | !! REFERENCE(S) : None |
|---|
| 18 | !! |
|---|
| 19 | !! \n |
|---|
| 20 | !_ |
|---|
| 21 | !================================================================================================================================ |
|---|
| 22 | MODULE grassland_constantes |
|---|
| 23 | |
|---|
| 24 | USE constantes |
|---|
| 25 | |
|---|
| 26 | LOGICAL, PARAMETER :: blabla_pasim = .FALSE. |
|---|
| 27 | LOGICAL, PARAMETER :: DEBUG_anne = .FALSE. |
|---|
| 28 | LOGICAL, PARAMETER :: BIG_DEBUG_anne = .FALSE. |
|---|
| 29 | INTEGER(i_std), PARAMETER :: taille = SELECTED_REAL_KIND(8) |
|---|
| 30 | INTEGER(i_std), PARAMETER :: taille_stdn = SELECTED_REAL_KIND(6,30) |
|---|
| 31 | INTEGER(i_std), PARAMETER :: nmaxprinc = 8760 |
|---|
| 32 | INTEGER(i_std), PARAMETER :: maxvar = 100 |
|---|
| 33 | INTEGER(i_std), PARAMETER :: NCUT = 10 |
|---|
| 34 | INTEGER(i_std), PARAMETER :: NSOILMAX = 6 |
|---|
| 35 | INTEGER(i_std), PARAMETER :: NFERT = 10 |
|---|
| 36 | INTEGER(i_std), PARAMETER :: NCANOPY = 10 |
|---|
| 37 | INTEGER(i_std), PARAMETER :: NGMEAN = 10 |
|---|
| 38 | INTEGER(i_std), PARAMETER :: NSTOCKING = 10 |
|---|
| 39 | INTEGER(i_std), PARAMETER :: NGAUSS = 5 |
|---|
| 40 | INTEGER(i_std), PARAMETER :: igrass = 10 |
|---|
| 41 | LOGICAL, PARAMETER :: Yearly_1snowzero = .FALSE. |
|---|
| 42 | REAL(r_std), PARAMETER :: minimum_vital = 1e-30 |
|---|
| 43 | REAL(r_std), PARAMETER :: DEFAULTMaxt = 0.02 |
|---|
| 44 | REAL(r_std), PARAMETER :: lambda = 2.49 |
|---|
| 45 | REAL(r_std), PARAMETER :: gamma = 0.065 |
|---|
| 46 | REAL(r_std), PARAMETER :: rho = 1.23 |
|---|
| 47 | REAL(R_std), PARAMETER :: alpha350 = 0.145 |
|---|
| 48 | REAL(r_std), PARAMETER :: devsecond = 0.77 |
|---|
| 49 | REAL(R_std), PARAMETER :: sigma = 5.67e-8 |
|---|
| 50 | LOGICAL , PARAMETER :: soilnh4pmobile = .TRUE. |
|---|
| 51 | REAL(R_std), PARAMETER :: nh3mol2nh3ug = 17.0e6 |
|---|
| 52 | REAL(R_std), PARAMETER :: qm2mol = 44.6 |
|---|
| 53 | REAL(R_std), PARAMETER :: nn2omol2kg = 0.028 |
|---|
| 54 | REAL(R_std), PARAMETER :: lf = 333750.0 |
|---|
| 55 | REAL(R_std), PARAMETER :: cp = 1.01e3 |
|---|
| 56 | REAL(R_std), PARAMETER :: rwt = 0 |
|---|
| 57 | REAL(R_std), PARAMETER :: dailysunshine = 0 |
|---|
| 58 | REAL(R_std), PARAMETER :: pa2atm = 101325.0 |
|---|
| 59 | ! konvertierung von umol co2 nach kg c (kg c/umol co2) |
|---|
| 60 | REAL(R_std), PARAMETER :: cumol2ckg = 12.0e-9 |
|---|
| 61 | ! convertion de mol n en kg n |
|---|
| 62 | REAL(R_std), PARAMETER :: nmol2kg = 0.014 |
|---|
| 63 | ! konvertierung von m nach mm (mm/m) |
|---|
| 64 | REAL(R_std), PARAMETER :: m2mm = 1000.0 |
|---|
| 65 | REAL(R_std), PARAMETER :: highresoutput = 0 |
|---|
| 66 | REAL(R_std), PARAMETER :: tfr = 273.15 |
|---|
| 67 | REAL(R_std), PARAMETER :: d0 = 2.15e6 |
|---|
| 68 | REAL(R_std), PARAMETER :: nstructfix = 0 |
|---|
| 69 | ! dichte von wasser @293.16k (kg h2o/m**3) |
|---|
| 70 | REAL(R_std), PARAMETER :: rhow = 998.206 |
|---|
| 71 | REAL(R_std), PARAMETER :: cw = 4.18e6 |
|---|
| 72 | REAL(R_std), PARAMETER :: kg2mg = 1.0e6 |
|---|
| 73 | ! anzahl sekunden pro tag (s/d) |
|---|
| 74 | REAL(R_std), PARAMETER :: d2s = 86400 |
|---|
| 75 | LOGICAL , PARAMETER :: unanaerob = .TRUE. |
|---|
| 76 | REAL(R_std), PARAMETER :: karman = 0.41 |
|---|
| 77 | REAL(R_std), PARAMETER :: dailyweatherdata = 0 |
|---|
| 78 | REAL(R_std), PARAMETER :: rgas = 8.3143 |
|---|
| 79 | REAL(R_std), PARAMETER :: solar = 1370 |
|---|
| 80 | REAL(R_std), PARAMETER :: infinit = 1.0e20 |
|---|
| 81 | REAL(R_std), PARAMETER :: grav = 9.81 |
|---|
| 82 | ! ! konvertierung von l nach m**3 (m**3/l) |
|---|
| 83 | REAL(R_std), PARAMETER :: l2m3 = 1.0e-3 |
|---|
| 84 | REAL(R_std), PARAMETER :: ug2kg = 1.0e-9 |
|---|
| 85 | REAL(R_std), PARAMETER :: cmol2kg = 0.012 |
|---|
| 86 | REAL(R_std), PARAMETER :: mw = 0.018 |
|---|
| 87 | REAL(R_std), PARAMETER :: mc = 28.5 |
|---|
| 88 | REAL(R_std), PARAMETER :: fcsh = 0.39 |
|---|
| 89 | REAL(R_std), PARAMETER :: ssta = 6.6 |
|---|
| 90 | REAL(R_std), PARAMETER :: fwnapo = 0.05 |
|---|
| 91 | REAL(R_std), PARAMETER :: mn = 62.0 |
|---|
| 92 | REAL(R_std), PARAMETER :: slam = 33.5 ! = 33.5 |
|---|
| 93 | REAL(R_std), PARAMETER :: fcr = 0.50 |
|---|
| 94 | REAL(R_std), PARAMETER :: tconstraintmax = 10000.0 |
|---|
| 95 | REAL(R_std), PARAMETER :: nage = 4.0 |
|---|
| 96 | ! parameter for the calculation of enteric methane emission (kg CH4 (kg life weight)-1 d-1) |
|---|
| 97 | REAL(r_std), PARAMETER :: aCH4 = 0.0002867 ! -0.0087 @equation constantes::aCH4 |
|---|
| 98 | ! parameter for the calculation of enteric methane emission (kg CH4 (kg life weight)-1 d-1) |
|---|
| 99 | REAL(r_std), PARAMETER :: bCH4 = 0.000045 ! 0.0541 @equation constantes::bCH4 |
|---|
| 100 | ! parameter for the calculation of enteric methane emission !!! |
|---|
| 101 | REAL(r_std), PARAMETER :: CH4toC = 0.75 |
|---|
| 102 | ! @equation constantes::CH4toC |
|---|
| 103 | ! parameter for transform the biomass(gC/m2) to wsh (kgDM/m2) |
|---|
| 104 | ! for grass C ratio in Dry Matter usually 0.45 for trees usually 0.5 |
|---|
| 105 | REAL(r_std), PARAMETER :: CtoDM = 0.45 |
|---|
| 106 | REAL(R_std), PARAMETER :: zeta = 10.0 |
|---|
| 107 | !facteur de conversion de MS a MO |
|---|
| 108 | REAL(R_std), PARAMETER :: dm2om = 0.9 |
|---|
| 109 | REAL(R_std), PARAMETER :: vcmaxadap = 1.0 |
|---|
| 110 | REAL(R_std), PARAMETER :: absorvl = 0.85 |
|---|
| 111 | ! fraction of net energy converted to metabolizable energy (-) |
|---|
| 112 | REAL(Taille), PARAMETER :: k_CH4 = 0.6 |
|---|
| 113 | ! flag d'activation de l'effet des temperatures elevees sur l'ingestion |
|---|
| 114 | LOGICAL , PARAMETER :: f_temperature_DMI = .TRUE. |
|---|
| 115 | ! TRUE : CH4 is calculated from N.Vuichard equation |
|---|
| 116 | LOGICAL , PARAMETER :: f_CH4_methode = .FALSE. |
|---|
| 117 | |
|---|
| 118 | REAL(R_std), PARAMETER :: devear = 0.52 |
|---|
| 119 | REAL(r_std), PARAMETER :: nanimaltotmax = 0.001 |
|---|
| 120 | REAL(r_std), PARAMETER :: tbase = 278.15 !278.0 |
|---|
| 121 | REAL(r_std), PARAMETER :: trep = 278.15 !278.0 |
|---|
| 122 | REAL(r_std), PARAMETER :: tasumrep = 225.0 |
|---|
| 123 | |
|---|
| 124 | ! variables for version 3.7 |
|---|
| 125 | REAL(r_std), PARAMETER :: l_min = 0.1 |
|---|
| 126 | REAL(r_std), PARAMETER :: age_init_1 = 0 |
|---|
| 127 | REAL(r_std), PARAMETER :: age_init_2 = 0 |
|---|
| 128 | REAL(r_std), PARAMETER :: age_init_3 = 0 |
|---|
| 129 | REAL(r_std), PARAMETER :: age_init_4 = 0 |
|---|
| 130 | |
|---|
| 131 | REAL(R_std), PARAMETER :: fligninstructinit = 0.54 |
|---|
| 132 | REAL(r_std), PARAMETER :: fligninresidue = 0.25 |
|---|
| 133 | |
|---|
| 134 | REAL(R_std), PARAMETER :: devstagemin = 0.6 |
|---|
| 135 | REAL(R_std), PARAMETER :: tgrowthmin = 45 !30.0 |
|---|
| 136 | REAL(R_std), PARAMETER :: tgrowthmax = 70.0 |
|---|
| 137 | REAL(R_std), PARAMETER :: misval = -99.999 |
|---|
| 138 | REAL(R_std), PARAMETER :: gmeansloperel = -0.002 |
|---|
| 139 | REAL(R_std), PARAMETER :: devstocking = 0.001 |
|---|
| 140 | ! original tseasonendmin = 250.0 is for France grassland |
|---|
| 141 | ! JC change it to 0 for global application |
|---|
| 142 | REAL(R_std), PARAMETER :: tseasonendmin = -1.0 |
|---|
| 143 | |
|---|
| 144 | ! constante for soil physics module |
|---|
| 145 | REAL(r_std), PARAMETER :: Temp_1csoil = 2.4E6 |
|---|
| 146 | REAL(r_std), PARAMETER :: Temp_1eta = 0.6 |
|---|
| 147 | REAL(r_std), PARAMETER :: sswmax = 10000.0 |
|---|
| 148 | REAL(r_std), PARAMETER :: tspring = 100.0 |
|---|
| 149 | REAL(r_std), PARAMETER :: tautomn = 300.0 |
|---|
| 150 | |
|---|
| 151 | ! constante for principal module |
|---|
| 152 | REAL(R_std), PARAMETER :: psifc = -612. |
|---|
| 153 | REAL(R_std), PARAMETER :: psipwp = -152905. |
|---|
| 154 | |
|---|
| 155 | ! constante for fertilisation module |
|---|
| 156 | REAL(r_std), PARAMETER :: tapplmist = 1.0 |
|---|
| 157 | REAL(r_std), PARAMETER :: fmistcmetabolic = 0.7 |
|---|
| 158 | REAL(r_std), PARAMETER :: tapplvg = 1.0 |
|---|
| 159 | REAL(r_std), PARAMETER :: asunshine = 0.19 |
|---|
| 160 | REAL(r_std), PARAMETER :: fvgcmetabolic = 0.8 |
|---|
| 161 | REAL(r_std), PARAMETER :: bsunshine = 0.56 |
|---|
| 162 | REAL(r_std), PARAMETER :: tapplka = 1.0 |
|---|
| 163 | REAL(r_std), PARAMETER :: c2nmist = 15.0 |
|---|
| 164 | REAL(r_std), PARAMETER :: fvgnurine = 0.61 |
|---|
| 165 | REAL(r_std), PARAMETER :: fkacmetabolic = 0.9 |
|---|
| 166 | REAL(r_std), PARAMETER :: c2nkot = 10.0 |
|---|
| 167 | REAL(r_std), PARAMETER :: fkanurine = 0.9 |
|---|
| 168 | ! si DMI/DMIpot < intake_tolerance alors on sort les animaux ou on les complemente |
|---|
| 169 | REAL(r_std), SAVE :: intake_tolerance = 0.80 |
|---|
| 170 | ! Maximum quantity of forage or concentrate to supplement animals when auto-supplementation (kg) |
|---|
| 171 | REAL(r_std), SAVE :: q_max_complement = 6 |
|---|
| 172 | |
|---|
| 173 | INTEGER, SAVE :: f_saturant = 0. |
|---|
| 174 | INTEGER, SAVE :: f_nonlimitant = 0. |
|---|
| 175 | INTEGER, SAVE :: f_autogestion = 0. !080110 AIG fin |
|---|
| 176 | INTEGER, SAVE :: f_complementation = 0. !11 janvier 2009 AIG |
|---|
| 177 | ! 0 : pas de fertilisation |
|---|
| 178 | ! 1 : fertilisation geree par le modele pour un niveau de satisfaction des besoins en N |
|---|
| 179 | INTEGER, SAVE :: f_fertilization = 0 |
|---|
| 180 | INTEGER, SAVE :: f_postauto = 0 |
|---|
| 181 | !071129 AIG |
|---|
| 182 | ! constante for auto management modifications by nicolas vuichard |
|---|
| 183 | ! INTEGER , PARAMETER :: n_somitermin = 2 |
|---|
| 184 | ! INTEGER , PARAMETER :: n_somitermax = 400 |
|---|
| 185 | ! REAL(Taille), PARAMETER :: nsatur_somerrormax = 0.01 |
|---|
| 186 | ! Minimum threshold of Wshtot, below which animals are kept indoors (kg.m-2), |
|---|
| 187 | ! either used when management is optimized by the model or for usual runs. |
|---|
| 188 | REAL(r_std), PARAMETER :: min_grazing = 0.11 |
|---|
| 189 | ! INTEGER , PARAMETER :: cte_nb_high_out = 1 |
|---|
| 190 | ! INTEGER , PARAMETER :: cte_nb_low_out = 1 |
|---|
| 191 | ! INTEGER , PARAMETER :: nb_last_year = 1 |
|---|
| 192 | !20/03/2009 AIG & MG |
|---|
| 193 | INTEGER, SAVE :: n_out = 3 |
|---|
| 194 | !Is the number of output for autogestion mangement run in import Yield file. |
|---|
| 195 | ! Dim1 : Fraction of grazed aera (F) Dim2: Ratio F/(1-F) Dim3: number of grazing days |
|---|
| 196 | |
|---|
| 197 | END MODULE grassland_constantes |
|---|
