hollorol/RBBGCMuso
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Add files via upload

Browse Source
This commit is contained in:
Zoltán BARCZA 2019-01-23 12:36:36 +01:00 committed by GitHub
parent afccf21682
commit 897499d571
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
12 changed files with 22271 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

55
RBBGCMuso/inst/CO2_from1961.txt Normal file
View File

@ -0,0 +1,55 @@
1961 317.33
1962 318.08
1963 318.70
1964 319.36
1965 320.02
1966 321.09
1967 321.99
1968 322.93
1969 324.21
1970 325.24
1971 326.06
1972 327.18
1973 328.84
1974 329.73
1975 330.73
1976 331.83
1977 333.25
1978 334.60
1979 336.85
1980 338.69
1981 339.93
1982 341.13
1983 342.78
1984 344.42
1985 345.90
1986 347.15
1987 348.93
1988 351.48
1989 352.91
1990 354.19
1991 355.59
1992 356.37
1993 357.04
1994 358.88
1995 360.88
1996 362.64
1997 363.76
1998 366.63
1999 368.31
2000 369.48
2001 372.59
2002 374.37
2003 378.04
2004 380.88
2005 383.88
2006 385.64
2007 385.76
2008 386.13
2009 387.37
2010 389.85
2011 391.62
2012 393.82
2013 396.48
2014 398.61
2015 400.00

180
RBBGCMuso/inst/c3grass.epc Normal file
View File

@ -0,0 +1,180 @@
ECOPHYS FILE - C3 grass muso5
----------------------------------------------------------------------------------------
FLAGS
0 (flag) biome type flag (1 = WOODY 0 = NON-WOODY)
0 (flag) woody type flag (1 = EVERGREEN 0 = DECIDUOUS)
1 (flag) photosyn. type flag (1 = C3 PSN 0 = C4 PSN)
1 (flag) phenology flag (1 = MODEL PHENOLOGY 0 = USER-SPECIFIED PHENOLOGY)
0 (flag) transferGDD flag (1= transfer calc. from GDD 0 = transfer calc. from EPC)
1 (flag) q10 flag (1 = temperature dependent q10 value; 0= constans q10 value)
0 (flag) acclimation flag (1 = acclimation 0 = no acclimation)
1 (flag) CO2 conductance reduction flag (0: no effect, 1: multiplier)
1 (flag) use GSI index to calculate growing season (0: no GSI, 1: GSI)
0 (flag) soil temperature calculation method (0: Zheng, 1: DSSAT)
1 (flag) soil hydrological calculation method (0: Richards, 1: tipping DSSAT)
0 (int) discretization level of SWC calculation (0: low, 1: medium, 2: high)
0 (flag) photosynthesis calculation method (0: Farquhar, 1: DSSAT)
0 (flag) evapotranspiration calculation method (0: Penman-Montieth, 1: Priestly-Taylor)
0 (flag) radiation calculation method (0: SWabs, 1: Rn)
----------------------------------------------------------------------------------------
PLANT FUNCTIONING PARAMETERS
0 (yday) yearday to start new growth (when phenology flag = 0)
364 (yday) yearday to end litterfall (when phenology flag = 0)
0.5 (prop.) transfer growth period as fraction of growing season (when transferGDD_flag = 0)
0.5 (prop.) litterfall as fraction of growing season (when transferGDD_flag = 0)
5 (Celsius) base temperature
-9999 (Celsius) minimum temperature for growth displayed on current day (-9999: no T-dependence of allocation)
-9999 (Celsius) optimal1 temperature for growth displayed on current day (-9999: no T-dependence of allocation)
-9999 (Celsius) optimal2 temperature for growth displayed on current day (-9999: no T-dependence of allocation)
-9999 (Celsius) maxmimum temperature for growth displayed on current day (-9999: no T-dependence of allocation)
-9999 (Celsius) minimum temperature for carbon assimilation displayed on current day (-9999: no limitation)
-9999 (Celsius) optimal1 temperature for carbon assimilation displayed on current day (-9999: no limitation)
-9999 (Celsius) optimal2 temperature for carbon assimilation displayed on current day (-9999: no limitation)
-9999 (Celsius) maxmimum temperature for carbon assimilation displayed on current day (-9999: no limitation)
1.0 (1/yr) annual leaf and fine root turnover fraction
0.00 (1/yr) annual live wood turnover fraction
0.01 (1/yr) annual whole-plant mortality fraction
0.0 (1/yr) annual fire mortality fraction
36.6 (kgC/kgN) C:N of leaves
45.0 (kgC/kgN) C:N of leaf litter, after retranslocation
50.0 (kgC/kgN) C:N of fine roots
36.6 *(kgC/kgN) C:N of fruit
36.6 (kgC/kgN) C:N of soft stem
0.0 *(kgC/kgN) C:N of live wood
0.0 *(kgC/kgN) C:N of dead wood
0.68 (DIM) leaf litter labile proportion
0.23 (DIM) leaf litter cellulose proportion
0.34 (DIM) fine root labile proportion
0.44 (DIM) fine root cellulose proportion
0.68 *(DIM) fruit litter labile proportion
0.23 *(DIM) fruit litter cellulose proportion
0.68 (DIM) soft stem litter labile proportion
0.23 (DIM) soft stem litter cellulose proportion
0.00 *(DIM) dead wood cellulose proportion
0.01 (1/LAI/d) canopy water interception coefficient
0.63 (DIM) canopy light extinction coefficient
2.0 (g/MJ) potential radiation use efficiency
0.781 (DIM) radiation parameter1 (Jiang et al.2015)
-13.596 (DIM) radiation parameter2 (Jiang et al.2015)
2.0 (DIM) all-sided to projected leaf area ratio
2.0 (DIM) ratio of shaded SLA:sunlit SLA
0.2 (DIM) fraction of leaf N in Rubisco
0.03 (DIM) fraction of leaf N in PEP Carboxylase
0.004 (m/s) maximum stomatal conductance (projected area basis)
0.00006 (m/s) cuticular conductance (projected area basis)
0.039 (m/s) boundary layer conductance (projected area basis)
1.00 (prop) relative SWC (prop. to FC) to calc. soil moisture limit 1 (-9999: not used)
0.99 (prop) relative SWC (prop. to SAT) to calc. soil moisture limit 2 (-9999: not used)
-9999 (prop) relative PSI (prop. to FC) to calc. soil moisture limit 1 (-9999: not used)
-9999 (prop) relative PSI (prop. to SAT) to calc. soil moisture limit 2 (-9999: not used)
1000 (Pa) vapor pressure deficit: start of conductance reduction
5000 (Pa) vapor pressure deficit: complete conductance reduction
1.5 (m) maximum height of plant
250 *(kgC/m2) stem weight at which maximum height attended
1.0 (m) maximum depth of rooting zone
3.67 (DIM) root distribution parameter
0.4 (kgC/m2) rootlenght parameter 1 (estimated max root weight)
0.5 (prop.) rootlenght parameter 2 (slope)
0.3 (prop.) growth resp per unit of C grown
0.195 (kgC/kgN/d) maintenance respiration in kgC/day per kg of tissue N
0.1 (DIM) theoretical maximum prop. of non-structural and structural carbohydrates
0.24 (DIM) prop. of non-structural carbohydrates available for maintanance respiration
0.0048 (kgN/m2/yr) symbiotic+asymbiotic fixation of N
----------------------------------------------------------------------------------------
CROP SPECIFIC PARAMETERS
0 (DIM) number of phenophase of germination (from 1 to 7; 0: NO specific)
0 (DIM) number of phenophase of emergence (from 1 to 7; 0: NO specific)
0.5 (prop.) critical relative SWC (prop. to FC) in germination
-9999 (DIM) number of phenophase of photoperiodic slowing effect (from 1 to 7; -9999: NO effect)
20 (hour) critical photoslow daylength
0.005 (DIM) slope of relative photoslow development rate
-9999 (DIM) number of phenophase of vernalization (from 1 to 7; -9999: NO)
0 (Celsius) critical vernalization temperature 1
5 (Celsius) critical vernalization temperature 2
8 (Celsius) critical vernalization temperature 3
15 (Celsius) critical vernalization temperature 4
0.04 (DIM) slope of relative vernalization development rate
50 (n) required vernalization days (in vernalization development rate)
-9999 (DIM) number of flowering phenophase (from 1 to 7; -9999: NO effect)
35 (Celsius) critical flowering heat stress temperature 1
40 (Celsius) critical flowering heat stress temperature 2
0.2 (prop.) mortality parameter of flowering thermal stress
----------------------------------------------------------------------------------------
SENESCENCE AND SOIL PARAMETERS
0.1 (prop.) maximum senescence mortality coefficient of aboveground plant material
0.05 (prop.) maximum senescence mortality coefficient of belowground plant material
0.0 (prop.) maximum senescence mortality coefficient of non-structured plant material
2 (prop) effect of extreme high temperature on senesncene mortality
30 (Celsius) lower limit extreme high temperature effect on senescence mortality
40 (Celsius) upper limit extreme high temperature effect on senescence mortality
-9999 (Celsius) maximal lifetime of plant tissue (-9999: no genetically programmed senescence)
0.01 (prop.) turnover rate of wilted standing biomass to litter
0.047 (prop.) turnover rate of non-woody cut-down biomass to litter
0.01 (prop.) turnover rate of woody cut-down biomass to litter
17 (prop.) drought tolerance parameter (critical value of DSWS)
0.02 (prop.) denitrification rate per g of CO2 respiration of SOM
0.2 (prop.) nitrification coefficient 1
0.07 (prop.) nitrification coefficient 2
0.02 (prop.) coefficient of N2O emission of nitrification
0.5 (prop.) proprortion of NH4 flux of N-deposition
0.5 (prop.) NH4 mobilen proportion
1.0 (prop.) NO3 mobilen proportion
10 (m) e-folding depth of decomposition rate's depth scalar
0.001 (prop.) fraction of dissolved part of SOIL1 organic matter
0.001 (prop.) fraction of dissolved part of SOIL2 organic matter
0.001 (prop.) fraction of dissolved part of SOIL3 organic matter
0.001 (prop.) fraction of dissolved part of SOIL4 organic matter
0.6 (DIM) ratio of bare soil evaporation and pot.evaporation
----------------------------------------------------------------------------------------
RATE SCALARS
0.39 (DIM) respiration fractions for fluxes between compartments (l1s1)
0.55 (DIM) respiration fractions for fluxes between compartments (l2s2)
0.29 (DIM) respiration fractions for fluxes between compartments (l4s3)
0.28 (DIM) respiration fractions for fluxes between compartments (s1s2)
0.46 (DIM) respiration fractions for fluxes between compartments (s2s3)
0.55 (DIM) respiration fractions for fluxes between compartments (s3s4)
0.7 (DIM) rate constant scalar of labile litter pool
0.07 (DIM) rate constant scalar of cellulose litter pool
0.014 (DIM) rate constant scalar of lignin litter pool
0.07 (DIM) rate constant scalar of fast microbial recycling pool
0.014 (DIM) rate constant scalar of medium microbial recycling pool
0.0014 (DIM) rate constant scalar of slow microbial recycling pool
0.0001 (DIM) rate constant scalar of recalcitrant SOM (humus) pool
0.001 (DIM) rate constant scalar of physical fragmentation of coarse woody debris
----------------------------------------------------------------------------------------
GROWING SEASON PARAMETERS
5 (kg/m2) crit. amount of snow limiting photosyn.
20 (Celsius) limit1 (under:full constrained) of HEATSUM index
60 (Celsius) limit2 (above:unconstrained) of HEATSUM index
0 (Celsius) limit1 (under:full constrained) of TMIN index
5 (Celsius) limit2 (above:unconstrained) of TMIN index
4000 (Pa) limit1 (above:full constrained) of VPD index
1000 (Pa) limit2 (under:unconstrained) of VPD index
0 (s) limit1 (under:full constrained) of DAYLENGTH index
0 (s) limit2 (above:unconstrained) of DAYLENGTH index
10 (day) moving average (to avoid the effects of extreme events)
0.10 (dimless) GSI limit1 (greater that limit -> start of vegper)
0.01 (dimless) GSI limit2 (less that limit -> end of vegper)
----------------------------------------------------------------------------------------
CH4 PARAMETERS
212.5 (DIM) param1 for CH4 calculations (empirical function of BD)
1.81 (DIM) param2 for CH4 calculations (empirical function of BD)
-1.353 (DIM) param1 for CH4 calculations (empirical function of VWC)
0.2 (DIM) param2 for CH4 calculations (empirical function of VWC)
1.781 (DIM) param3 for CH4 calculations (empirical function of VWC)
6.786 (DIM) param4 for CH4 calculations (empirical function of VWC)
0.010 (DIM) param1 for CH4 calculations (empirical function of Tsoil)
----------------------------------------------------------------------------------------
PHENOLOGICAL (ALLOCATION) PARAMETERS (7 phenological phases)
phase1 phase2 phase3 phase4 phase5 phase6 phase7 (text) name of the phenophase
500 200 500 200 400 200 100 (Celsius) length of phenophase (GDD)
0.4 0.4 0.4 0.4 0.4 0.4 0.4 (ratio) leaf ALLOCATION
0.4 0.4 0.4 0.4 0.4 0.4 0.4 (ratio) fine root ALLOCATION
0.0 0.0 0.0 0.0 0.0 0.0 0.0 (ratio) fruit ALLOCATION
0.2 0.2 0.2 0.2 0.2 0.2 0.2 (ratio) soft stem ALLOCATION
0 0 0 0 0 0 0 (ratio) live woody stem ALLOCATION
0 0 0 0 0 0 0 (ratio) dead woody stem ALLOCATION
0 0 0 0 0 0 0 (ratio) live coarse root ALLOCATION
0 0 0 0 0 0 0 (ratio) dead coarse root ALLOCATION
49 49 49 49 49 49 49 (m2/kgC) canopy average specific leaf area (projected area basis)
0.37 0.37 0.37 0.37 0.37 0.37 0.37 (prop.) current growth proportion

BIN
RBBGCMuso/inst/cygwin1.dll Normal file
View File

Binary file not shown.

19714
RBBGCMuso/inst/hhs_1961-2014.mtc43 Normal file
View File

File diff suppressed because it is too large Load Diff

1464
RBBGCMuso/inst/hhs_2013-2016.mtc43 Normal file
View File

File diff suppressed because it is too large Load Diff

180
RBBGCMuso/inst/maize.epc Normal file
View File

@ -0,0 +1,180 @@
ECOPHYS FILE - MAIZE
----------------------------------------------------------------------------------------
FLAGS
0 (flag) biome type flag (1 = WOODY 0 = NON-WOODY)
0 (flag) woody type flag (1 = EVERGREEN 0 = DECIDUOUS)
0 (flag) photosyn. type flag (1 = C3 PSN 0 = C4 PSN)
0 (flag) phenology flag (1 = MODEL PHENOLOGY 0 = USER-SPECIFIED PHENOLOGY)
1 (flag) transferGDD flag (1= transfer calc. from GDD 0 = transfer calc. from EPC)
1 (flag) q10 flag (1 = temperature dependent q10 value; 0= constans q10 value)
0 (flag) acclimation flag (1 = acclimation 0 = no acclimation)
1 (flag) CO2 conductance reduction flag (0: no effect, 1: multiplier)
0 (flag) use GSI index to calculate growing season (0: no GSI, 1: GSI)
0 (flag) soil temperature calculation method (0: Zheng, 1: DSSAT)
1 (flag) soil hydrological calculation method (0: Richards, 1: tipping DSSAT)
0 (int) discretization level of SWC calculation (0: low, 1: medium, 2: high)
0 (flag) photosynthesis calculation method (0: Farquhar, 1: DSSAT)
0 (flag) evapotranspiration calculation method (0: Penman-Montieth, 1: Priestly-Taylor)
0 (flag) radiation calculation method (0: SWabs, 1: Rn)
----------------------------------------------------------------------------------------
PLANT FUNCTIONING PARAMETERS
0 (yday) yearday to start new growth (when phenology flag = 0)
364 (yday) yearday to end litterfall (when phenology flag = 0)
0.4 (prop.) transfer growth period as fraction of growing season
0. (prop.) litterfall as fraction of growing season
8 (Celsius) base temperature
-9999 (Celsius) minimum temperature for growth displayed on current day (-9999: no T-dependence of allocation)
-9999 (Celsius) optimal1 temperature for growth displayed on current day (-9999: no T-dependence of allocation)
-9999 (Celsius) optimal2 temperature for growth displayed on current day (-9999: no T-dependence of allocation)
-9999 (Celsius) maxmimum temperature for growth displayed on current day (-9999: no T-dependence of allocation)
-9999 (Celsius) minimum temperature for carbon assimilation displayed on current day (-9999: no limitation)
-9999 (Celsius) optimal1 temperature for carbon assimilation displayed on current day (-9999: no limitation)
-9999 (Celsius) optimal2 temperature for carbon assimilation displayed on current day (-9999: no limitation)
-9999 (Celsius) maxmimum temperature for carbon assimilation displayed on current day (-9999: no limitation)
1.0 (1/yr) annual leaf and fine root turnover fraction
0.00 (1/yr) annual live wood turnover fraction
0.0 (1/yr) annual whole-plant mortality fraction
0.0 (1/yr) annual fire mortality fraction
40.0 (kgC/kgN) C:N of leaves
65.0 (kgC/kgN) C:N of leaf litter, after retranslocation ion
42.0 (kgC/kgN) C:N of fine roots
50.0 *(kgC/kgN) C:N of fruit
50.0 (kgC/kgN) C:N of soft stem
0.0 (kgC/kgN) C:N of live wood
0.0 (kgC/kgN) C:N of dead wood
0.68 (DIM) leaf litter labile proportion
0.23 (DIM) leaf litter cellulose proportion
0.34 (DIM) fine root labile proportion
0.44 (DIM) fine root cellulose proportion
0.68 *(DIM) fruit litter labile proportion
0.23 *(DIM) fruit litter cellulose proportion
0.68 (DIM) soft stem litter labile proportion
0.23 (DIM) soft stem litter cellulose proportion
0.00 *(DIM) dead wood cellulose proportion
0.01 (1/LAI/d) canopy water interception coefficient
0.69 (DIM) canopy light extinction coefficient
2.0 (g/MJ) potential radiation use efficiency
0.781 (DIM) radiation parameter1 (Jiang et al.2015)
-13.596 (DIM) radiation parameter2 (Jiang et al.2015)
2.0 (DIM) all-sided to projected leaf area ratio
2.0 (DIM) ratio of shaded SLA:sunlit SLA
0.14 (DIM) fraction of leaf N in Rubisco CLM alapjan
0.0424 (DIM) fraction of leaf N in PEP Carboxylase
0.012 (m/s) maximum stomatal conductance (projected area basis)
0.00006 (m/s) cuticular conductance (projected area basis)
0.04 (m/s) boundary layer conductance (projected area basis)
0.99 (prop) relative SWC (prop. to FC) to calc. soil moisture limit 1 (-9999: not used)
1.0 (prop) relative SWC (prop. to SAT) to calc. soil moisture limit 2 (-9999: not used)
-9999 (prop) relative PSI (prop. to FC) to calc. soil moisture limit 1 (-9999: not used)
-9999 (prop) relative PSI (prop. to SAT) to calc. soil moisture limit 2 (-9999: not used)
1000 (Pa) vapor pressure deficit: start of conductance reduction
5000 (Pa) vapor pressure deficit: complete conductance reduction
2.0 (m) maximum height of plant
25 *(kgC/m2) stem weight at which maximum height attended
1.7 (m) (m) maximum depth of rooting zonemaximum depth of rooting zone
3.67 (DIM) root distribution parameter
0.4 (kgC/m2) rootlenght parameter 1 (estimated max root weight)
0.5 (prop.) rootlenght parameter 2 (slope)
0.3 (prop.) growth resp per unit of C grown
0.218 (kgC/kgN/d) maintenance respiration in kgC/day per kg of tissue N
0.1 (DIM) theoretical maximum prop. of non-structural and structural carbohydrates
0.3 (DIM) prop. of non-structural carbohydrates available for maintanance respiration
0.0005 (kgN/m2/yr) symbiotic+asymbiotic fixation of N
----------------------------------------------------------------------------------------
CROP SPECIFIC PARAMETERS
2 (DIM) number of the first phenophase when germination is complete (from 2 to 7; 0: NO specific)
3 (DIM) number of the first phenophase when emergence is complete (from 3 to 7; 0: NO specific)
0.5 (prop.) critical relative SWC (prop. to FC) in germination
-9999 (DIM) number of phenophase of photoperiodic slowing effect (from 1 to 7; -9999: NO)
16 (hour) critical photoslow daylength
0.005 (DIM) slope of relative photoslow development rate
-9999 (DIM) number of phenophase of vernalization (from 1 to 7; -9999: NO)
0 (Celsius) critical vernalization temperature 1
5 (Celsius) critical vernalization temperature 2
8 (Celsius) critical vernalization temperature 3
15 (Celsius) critical vernalization temperature 4
0.04 (DIM) slope of relative vernalization development rate
50 (n) required vernalization days (in vernalization development rate)
5 (DIM) number of flowering phenophase (from 1 to 7; -9999: NO)
35 (Celsius) critical flowering heat stress temperature 1
40 (Celsius) critical flowering heat stress temperature 2
0.2 (prop.) mortality parameter of flowering thermal stress
----------------------------------------------------------------------------------------
SENESCENCE AND SOIL PARAMETERS
0.01 (prop.) maximum senescence mortality coefficient of aboveground plant material
0.01 (prop.) maximum senescence mortality coefficient of belowground plant material
0.0 (prop.) maximum senescence mortality coefficient of non-structured plant material
2 (prop) effect of extreme high temperature on senesncene mortality
35 (Celsius) lower limit extreme high temperature effect on senescence mortality
40 (Celsius) upper limit extreme high temperature effect on senescence mortality
1000 (Celsius) maximal lifetime of plant tissue (-9999: no genetically programmed senescence)
0.01 (prop.) turnover rate of wilted standing biomass to litter
0.05 (prop.) turnover rate of non-woody cut-down biomass to litter
0.01 (prop.) turnover rate of woody cut-down biomass to litter
30 (prop.) drought tolerance parameter (critical value of DSWS)
0.02 (prop.) denitrification rate per g of CO2 respiration of SOM
0.5 (prop.) nitrification coefficient 1
0.4 (prop.) nitrification coefficient 2
0.02 (prop.) coefficient of N2O emission of nitrification
0.5 (prop.) proprortion of NH4 flux of N-deposition
0.5 (prop.) NH4 mobilen proportion
1.0 (prop.) NO3 mobilen proportion
10 (m) e-folding depth of decomposition rate's depth scalar
0.001 (prop.) fraction of dissolved part of SOIL1 organic matter
0.001 (prop.) fraction of dissolved part of SOIL2 organic matter
0.001 (prop.) fraction of dissolved part of SOIL3 organic matter
0.001 (prop.) fraction of dissolved part of SOIL4 organic matter
0.6 (DIM) ratio of bare soil evaporation and pot.evaporation
----------------------------------------------------------------------------------------
RATE SCALARS
0.39 (DIM) respiration fractions for fluxes between compartments (l1s1)
0.55 (DIM) respiration fractions for fluxes between compartments (l2s2)
0.29 (DIM) respiration fractions for fluxes between compartments (l4s3)
0.28 (DIM) respiration fractions for fluxes between compartments (s1s2)
0.46 (DIM) respiration fractions for fluxes between compartments (s2s3)
0.55 (DIM) respiration fractions for fluxes between compartments (s3s4)
0.7 (DIM) rate constant scalar of labile litter pool
0.07 (DIM) rate constant scalar of cellulose litter pool
0.014 (DIM) rate constant scalar of lignin litter pool
0.07 (DIM) rate constant scalar of fast microbial recycling pool
0.014 (DIM) rate constant scalar of medium microbial recycling pool
0.0014 (DIM) rate constant scalar of slow microbial recycling pool
0.0001 (DIM) rate constant scalar of recalcitrant SOM (humus) pool
0.001 (DIM) rate constant scalar of physical fragmentation of coarse woody debris
----------------------------------------------------------------------------------------
GROWING SEASON PARAMETERS
5 (kg/m2) crit. amount of snow limiting photosyn.
20 (Celsius) limit1 (under:full constrained) of HEATSUM index
60 (Celsius) limit2 (above:unconstrained) of HEATSUM index
0 (Celsius) limit1 (under:full constrained) of TMIN index
5 (Celsius) limit2 (above:unconstrained) of TMIN index
4000 (Pa) limit1 (above:full constrained) of VPD index
1000 (Pa) limit2 (under:unconstrained) of VPD index
0 (s) limit1 (under:full constrained) of DAYLENGTH index
0 (s) limit2 (above:unconstrained) of DAYLENGTH index
10 (day) moving average (to avoid the effects of extreme events)
0.10 (dimless) GSI limit1 (greater that limit -> start of vegper)
0.01 (dimless) GSI limit2 (less that limit -> end of vegper)
----------------------------------------------------------------------------------------
CH4 PARAMETERS
212.5 (DIM) param1 for CH4 calculations (empirical function of BD)
1.81 (DIM) param2 for CH4 calculations (empirical function of BD)
-1.353 (DIM) param1 for CH4 calculations (empirical function of VWC)
0.2 (DIM) param2 for CH4 calculations (empirical function of VWC)
1.781 (DIM) param3 for CH4 calculations (empirical function of VWC)
6.786 (DIM) param4 for CH4 calculations (empirical function of VWC)
0.010 (DIM) param1 for CH4 calculations (empirical function of Tsoil)
----------------------------------------------------------------------------------------
PHENOLOGICAL (ALLOCATION) PARAMETERS (7 phenological phases)
germination emergence leafunfold leafunfoldlin flowering grainfilling maturity (text) name of the phenophase
0 45 225 630 90 730 10000 (Celsius) length of phenophase (GDD)
0.0 0.5 0.60 0.30 0.05 0.0 0. (ratio) leaf ALLOCATION
0.0 0.5 0.35 0.25 0.20 0.15 0. (ratio) fine root ALLOCATION
0.0 0. 0.0 0.0 0.05 0.65 1. (ratio) fruit ALLOCATION
0.0 0. 0.05 0.45 0.70 0.20 0. (ratio) soft stem ALLOCATION
0 0 0 0 0 0 0 (ratio) live woody stem ALLOCATION
0 0 0 0 0 0 0 (ratio) dead woody stem ALLOCATION
0 0 0 0 0 0 0 (ratio) live coarse root ALLOCATION
0 0 0 0 0 0 0 (ratio) dead coarse root ALLOCATION
50 50 50 50 50 50 50 (m2/kgC) canopy average specific leaf area (projected area basis)
1.0 1.0 1.0 1.0 1.0 1.0 1.0 (prop.) current growth proportion

BIN
RBBGCMuso/inst/muso Normal file
View File

Binary file not shown.

BIN
RBBGCMuso/inst/muso.exe Normal file
View File

Binary file not shown.

229
RBBGCMuso/inst/n.ini Normal file
View File

@ -0,0 +1,229 @@
BBGC_MuSo simulation (missing data: -9999)
MET_INPUT
hhs_2013-2016.mtc43 (filename) met file name
4 (int) number of header lines in met file
365 (int) number of simdays in last simyear (truncated year: < 365)
RESTART
1 (flag) 1 = read restart; 0 = dont read restart
0 (flag) 1 = write restart; 0 = dont write restart
0 (flag) 1 = use restart metyear; 0 = reset metyear
hhs.endpoint (filename) name of the input restart file
hhs.endpoint (filename) name of the output restart file
TIME_DEFINE
4 (int) number of meteorological data years
4 (int) number of simulation years
2013 (int) first simulation year
0 (flag) 1 = spinup run; 0 = normal run
6000 (int) maximum number of spinup years
CLIM_CHANGE
0.0 (degC) - offset for Tmax
0.0 (degC) - offset for Tmin
1.0 (degC) - multiplier for PRCP
1.0 (degC) - multiplier for VPD
1.0 (degC) - multiplier for RAD
CO2_CONTROL
0 (flag) 0=constant; 1=vary with file
395.0 (ppm) constant atmospheric CO2 concentration
CO2_from2013.txt (filename) name of the CO2 file
NDEP_CONTROL
0 (flag) 0=constant; 1=vary with file
0.001400 (kgN/m2/yr) wet+dry atmospheric deposition of N
Ndep_from2013.txt (filename) name of the N-dep file
SITE
30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 (%) sand percentage by volume in rock-free soil
50.0 50.0 50.0 50.0 50.0 50.0 50.0 50.0 50.0 50.0 (%) silt percentage by volume in rock-free soil
7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 (dimless) soil pH
248.0 (m) site elevation
46.95 (degrees) site latitude (- for S.Hem.)
0.20 (DIM) site shortwave albedo
9.00 (Celsius) mean annual air temperature
10.15 (Celsius) mean annual air temperature range
20.00 (mm) maximum height of pond water
-9999 (dimless) measured runoff curve number
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (g/cm3) bulk density
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at saturation
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at field capacity
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at wilting point
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at hygroscopic water content
EPC_FILE
maize.epc (filename) EPC file name
W_STATE
0.0 (kg/m2) water stored in snowpack
1.0 (DIM) initial soil water as a proportion of field capacity
CN_STATE
0.001 (kgC/m2) first-year maximum leaf carbon
0.001 (kgC/m2) first-year maximum fine root carbon
0.001 (kgC/m2) first-year maximum fruit carbon
0.001 (kgC/m2) first-year maximum softstem carbon
0.001 (kgC/m2) first-year maximum live woody stem carbon
0.001 (kgC/m2) first-year maximum live coarse root carbon
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) coarse woody debris carbon
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, labile pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, unshielded cellulose pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, shielded cellulose pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, lignin pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, fast microbial recycling pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, medium microbial recycling pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, slow microbial recycling pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, recalcitrant SOM (slowest)
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgN/m2) litter nitrogen, labile pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgN/m2) soil mineralized nitrogen, NH4 pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgN/m2) soil mineralized nitrogen, NO3 pool
OUTPUT_CONTROL
hegyhatsal (filename) output prefix
1 (flag) writing daily output (0 = no; 1 = binary; 2 = ascii)
0 (flag) writing monthly average of daily output (0 = no; 1 = binary; 2 = ascii)
0 (flag) writing annual average of daily output (0 = no; 1 = binary; 2 = ascii)
0 (flag) writing annual output (0 = no; 1 = binary; 2 = ascii)
1 (flag) for on-screen progress indicator
DAILY_OUTPUT
39 number of daily output variables
49 1 tsoil[0]
171 2 evapotranspiration
159 3 soilw_evap
397 4 litr1c_total
398 5 litr2c_total
399 6 litr3c_total
400 7 litr4c_total
518 8 soil1c_total
519 9 soil2c_total
520 10 soil3c_total
521 11 soil4c_total
313 12 fruitc
552 13 HRV_transportC
2520 14 proj_lai
3017 15 total litterfall aboveground
3018 16 total litterfall belowground
2599 17 vwc[0] soilw layer1 first layer 0-2 cm
2600 18 vwc[1] the second 2-5 cm
2601 19 vwc[2] the third is at 5-10 cm
2602 20 vwc[3] the fourth is at 10-20 cm
2603 21 vwc[4] the fifth is at 20-50 cm
2604 22 vwc[5] the sixth is at 50-100
407 23 standing dead biomass
3006 24 daily_npp;
3005 25 daily_nep
3009 26 daily_gpp
3010 27 daily_mr
3011 28 daily_gr
3012 29 daily_hr
3051 30 litrc
3052 31 soilc
3037 32 cum_yieldC_HRV (kgC/m2)
74 33 GDD
1531 34 SUM of the soil mineral NH4 in the total soil [kgN/m2]
1532 35 SUM of the soil mineral NO3 in the total soil [kgN/m2]
3013 36 daily soil respiration
307 37 leafC
310 38 fine root C
316 39 soft stem C
ANNUAL_OUTPUT
12 number of annual output variables
0 onday
1 offday
3000 annprcp
3001 anntavg
3002 annrunoff
2765 ytd_maxplai
3024 cum_nee
3033 cum_Closs_THN_w
3035 cum_Closs_MOW
3037 cum_yieldC_HRV
3059 NH4_top10
3060 NO3_top10
-------------------
MANAGEMENT_SECTION
-------------------
PLANTING
1 (flag) do PLANTING? 0=no; 1=yes; filepath=reading from file
121 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) PLANTING day
0.05 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m) germination depth
7 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(n/m2) number of seedlings
380 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(g/1000n) weight of 1000-seed
40 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) C content of seed
100 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) emergence rate
THINNING
0 (flag) do THINNING? 0=no; 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) THINNING day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) thinning rate of woody plant material
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) thinning rate of non-woody plant material
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of stem
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of leaf
MOWING
0 (flag) do MOWING? 0=no; 1=yes; filepath=reading from file
0 (flag) mowing method? 0 - fixday method, 1 - fixvalue method
-9999 (m2/m2) fixed value of the LAI before MOWING (fixvalue method)
-9999 (m2/m2) fixed value of the LAI after MOWING (fixvalue method)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) MOWING day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m2/m2) value of the LAI after MOWING (fixday method)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of plant material
GRAZING
0 (flag) do GRAZING? 0=no, 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) first day of GRAZING
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) last day of GRAZING
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kg/LSU) weight equivalent of one unit
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(LSU/ha) animal stocking rate: Livestock Units per hectare
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kg dry matter/LSU) daily ingested dry matter
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) trampling effect (standing dead biome to litter)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) ratio of DM intake formed excrement
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) ratio of excrement returning to litter
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) carbon content of dry matter
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) N content of manure
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) C content of manure
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN2O-N:kgN) manure emission factor for direct N2O (T10.21)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/1000 kg animal mass/day) default N excretion rate (T10.19)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgCH4/LSU/yr*) manure emission factor for CH4 (T10.14)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgCH4/LSU/yr*) fermentation emission factor for CH4 (T10.11)
HARVESTING
1 (flag) do HARVESTING? 0=no; 1=yes; filepath=reading from file
275 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) HARVESTING day
0.03 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgC/m2) soft stem C content after HARVESTING
90 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of plant material
PLOUGHING
1 (flag) do PLOUGHING? 0=no; 1=yes; filepath=reading from file
110 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) PLOUGHING day
0.3 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m) PLOUGHING depth
0.1 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) dissolving coefficient of ploughed biome to litter
FERTILIZING
1 (flag) do FERTILIZING? 0=no; 1=yes; filepath=reading from file
90 140 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) FERTILIZING day
0.28 0.28 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m) fertilizing depth
120. 30. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN fertilizer/ha/day) amount of fertilizer
100 100 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) dry matter content of fertilizer
0.5 0.5 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/kg fertilizer) nitrate content of fertilizer
0.5 0.5 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/kg fertilizer) ammonium content of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/kg fertilizer) organic nitrogen content of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgC/kg fertilizer) organic carbon content of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) labile fraction of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) unshielded cellulose fraction of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) shielded cellulose fraction of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) lignin fraction of fertilizer
0.01 0.01 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN2O-N:kgN) emission factor for N-additions
IRRIGATION
0 (flag) do IRRIGATION? 0=no, 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) IRRIGATION day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgH2O/m2/day) amount of water
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) utilization coefficient of water
END_INIT

229
RBBGCMuso/inst/normal_grass.ini Normal file
View File

@ -0,0 +1,229 @@
BBGC_MuSo simulation (missing data: -9999)
MET_INPUT
hhs_2013-2016.mtc43 (filename) met file name
4 (int) number of header lines in met file
365 (int) number of simdays in last simyear (truncated year: < 365)
RESTART
1 (flag) 1 = read restart; 0 = dont read restart
0 (flag) 1 = write restart; 0 = dont write restart
0 (flag) 1 = use restart metyear; 0 = reset metyear
hhs.endpoint (filename) name of the input restart file
hhs.endpoint (filename) name of the output restart file
TIME_DEFINE
4 (int) number of meteorological data years
4 (int) number of simulation years
2013 (int) first simulation year
0 (flag) 1 = spinup run; 0 = normal run
6000 (int) maximum number of spinup years
CLIM_CHANGE
0.0 (degC) - offset for Tmax
0.0 (degC) - offset for Tmin
1.0 (degC) - multiplier for PRCP
1.0 (degC) - multiplier for VPD
1.0 (degC) - multiplier for RAD
CO2_CONTROL
0 (flag) 0=constant; 1=vary with file
395.0 (ppm) constant atmospheric CO2 concentration
CO2_from2013.txt (filename) name of the CO2 file
NDEP_CONTROL
0 (flag) 0=constant; 1=vary with file
0.001400 (kgN/m2/yr) wet+dry atmospheric deposition of N
Ndep_from2013.txt (filename) name of the N-dep file
SITE
30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 (%) sand percentage by volume in rock-free soil
50.0 50.0 50.0 50.0 50.0 50.0 50.0 50.0 50.0 50.0 (%) silt percentage by volume in rock-free soil
7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 (dimless) soil pH
248.0 (m) site elevation
46.95 (degrees) site latitude (- for S.Hem.)
0.20 (DIM) site shortwave albedo
9.00 (Celsius) mean annual air temperature
10.15 (Celsius) mean annual air temperature range
20.00 (mm) maximum height of pond water
-9999 (dimless) measured runoff curve number
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (g/cm3) bulk density
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at saturation
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at field capacity
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at wilting point
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at hygroscopic water content
EPC_FILE
maize.epc (filename) EPC file name
W_STATE
0.0 (kg/m2) water stored in snowpack
1.0 (DIM) initial soil water as a proportion of field capacity
CN_STATE
0.001 (kgC/m2) first-year maximum leaf carbon
0.001 (kgC/m2) first-year maximum fine root carbon
0.001 (kgC/m2) first-year maximum fruit carbon
0.001 (kgC/m2) first-year maximum softstem carbon
0.001 (kgC/m2) first-year maximum live woody stem carbon
0.001 (kgC/m2) first-year maximum live coarse root carbon
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) coarse woody debris carbon
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, labile pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, unshielded cellulose pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, shielded cellulose pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, lignin pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, fast microbial recycling pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, medium microbial recycling pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, slow microbial recycling pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, recalcitrant SOM (slowest)
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgN/m2) litter nitrogen, labile pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgN/m2) soil mineralized nitrogen, NH4 pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgN/m2) soil mineralized nitrogen, NO3 pool
OUTPUT_CONTROL
hegyhatsal (filename) output prefix
1 (flag) writing daily output (0 = no; 1 = binary; 2 = ascii)
0 (flag) writing monthly average of daily output (0 = no; 1 = binary; 2 = ascii)
0 (flag) writing annual average of daily output (0 = no; 1 = binary; 2 = ascii)
0 (flag) writing annual output (0 = no; 1 = binary; 2 = ascii)
1 (flag) for on-screen progress indicator
DAILY_OUTPUT
39 number of daily output variables
49 1 tsoil[0]
171 2 evapotranspiration
159 3 soilw_evap
397 4 litr1c_total
398 5 litr2c_total
399 6 litr3c_total
400 7 litr4c_total
518 8 soil1c_total
519 9 soil2c_total
520 10 soil3c_total
521 11 soil4c_total
313 12 fruitc
552 13 HRV_transportC
2520 14 proj_lai
3017 15 total litterfall aboveground
3018 16 total litterfall belowground
2599 17 vwc[0] soilw layer1 first layer 0-2 cm
2600 18 vwc[1] the second 2-5 cm
2601 19 vwc[2] the third is at 5-10 cm
2602 20 vwc[3] the fourth is at 10-20 cm
2603 21 vwc[4] the fifth is at 20-50 cm
2604 22 vwc[5] the sixth is at 50-100
407 23 standing dead biomass
3006 24 daily_npp
3005 25 daily_nep
3009 26 daily_gpp
3010 27 daily_mr
3011 28 daily_gr
3012 29 daily_hr
3051 30 litrc
3052 31 soilc
3037 32 cum_yieldC_HRV
74 33 GDD
1531 34 SUM of the soil mineral NH4 in the total soil [kgN/m2]
1532 35 SUM of the soil mineral NO3 in the total soil [kgN/m2]
3013 36 daily soil respiration
307 37 leafC
310 38 fine root C
316 39 soft stem C
ANNUAL_OUTPUT
12 number of annual output variables
0 onday
1 offday
3000 annprcp
3001 anntavg
3002 annrunoff
2765 ytd_maxplai
3024 cum_nee
3033 cum_Closs_THN_w
3035 cum_Closs_MOW
3037 cum_yieldC_HRV
3059 NH4_top10
3060 NO3_top10
-------------------
MANAGEMENT_SECTION
-------------------
PLANTING
1 (flag) do PLANTING? 0=no; 1=yes; filepath=reading from file
121 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) PLANTING day
0.05 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m) germination depth
7 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(n/m2) number of seedlings
380 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(g/1000n) weight of 1000-seed
40 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) C content of seed
100 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) emergence rate
THINNING
0 (flag) do THINNING? 0=no; 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) THINNING day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) thinning rate of woody plant material
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) thinning rate of non-woody plant material
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of stem
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of leaf
MOWING
0 (flag) do MOWING? 0=no; 1=yes; filepath=reading from file
0 (flag) mowing method? 0 - fixday method, 1 - fixvalue method
-9999 (m2/m2) fixed value of the LAI before MOWING (fixvalue method)
-9999 (m2/m2) fixed value of the LAI after MOWING (fixvalue method)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) MOWING day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m2/m2) value of the LAI after MOWING (fixday method)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of plant material
GRAZING
0 (flag) do GRAZING? 0=no, 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) first day of GRAZING
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) last day of GRAZING
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kg/LSU) weight equivalent of one unit
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(LSU/ha) animal stocking rate: Livestock Units per hectare
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kg dry matter/LSU) daily ingested dry matter
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) trampling effect (standing dead biome to litter)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) ratio of DM intake formed excrement
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) ratio of excrement returning to litter
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) carbon content of dry matter
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) N content of manure
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) C content of manure
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN2O-N:kgN) manure emission factor for direct N2O (T10.21)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/1000 kg animal mass/day) default N excretion rate (T10.19)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgCH4/LSU/yr*) manure emission factor for CH4 (T10.14)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgCH4/LSU/yr*) fermentation emission factor for CH4 (T10.11)
HARVESTING
0 (flag) do HARVESTING? 0=no; 1=yes; filepath=reading from file
275 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) HARVESTING day
0.03 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgC/m2) soft stem C content after HARVESTING
90 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of plant material
PLOUGHING
0 (flag) do PLOUGHING? 0=no; 1=yes; filepath=reading from file
110 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) PLOUGHING day
0.3 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m) PLOUGHING depth
0.1 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) dissolving coefficient of ploughed biome to litter
FERTILIZING
0 (flag) do FERTILIZING? 0=no; 1=yes; filepath=reading from file
90 140 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) FERTILIZING day
0.28 0.28 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m) fertilizing depth
120. 30. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN fertilizer/ha/day) amount of fertilizer
100 100 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) dry matter content of fertilizer
0.5 0.5 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/kg fertilizer) nitrate content of fertilizer
0.5 0.5 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/kg fertilizer) ammonium content of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/kg fertilizer) organic nitrogen content of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgC/kg fertilizer) organic carbon content of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) labile fraction of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) unshielded cellulose fraction of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) shielded cellulose fraction of fertilizer
0. 0. -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) lignin fraction of fertilizer
0.01 0.01 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN2O-N:kgN) emission factor for N-additions
IRRIGATION
0 (flag) do IRRIGATION? 0=no, 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) IRRIGATION day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgH2O/m2/day) amount of water
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) utilization coefficient of water
END_INIT

16
RBBGCMuso/inst/parameters.csv Normal file
View File

@ -0,0 +1,16 @@
NAME,INDEX,MIN,MAX
BASETEMP,25,3,9
WPM,36,0,1
EMERGENCE,170.61,0,1000
CN_lv,38,10,50
CN_li,39,32,70
CN_root,40,20,90
CN_fruit,41,0,200
CN_stem,42,0,200
CLEC,55,0.4,0.8
FLNR,61,0.05,0.8
STOMA,63,0.003,0.015
ROOTDEPTH,74,0.3,2.
SWCGERMIN,87,0.2,0.9
MAXLIFETIME,110,500,1500
NH4MOBILEPROP,120,0.05,0.7
1 NAME INDEX MIN MAX
2 BASETEMP 25 3 9
3 WPM 36 0 1
4 EMERGENCE 170.61 0 1000
5 CN_lv 38 10 50
6 CN_li 39 32 70
7 CN_root 40 20 90
8 CN_fruit 41 0 200
9 CN_stem 42 0 200
10 CLEC 55 0.4 0.8
11 FLNR 61 0.05 0.8
12 STOMA 63 0.003 0.015
13 ROOTDEPTH 74 0.3 2.
14 SWCGERMIN 87 0.2 0.9
15 MAXLIFETIME 110 500 1500
16 NH4MOBILEPROP 120 0.05 0.7

204
RBBGCMuso/inst/s.ini Normal file
View File

@ -0,0 +1,204 @@
BBGC_MuSo simulation (missing data: -9999)
MET_INPUT
hhs_1961-2014.mtc43 (filename) met file name
4 (int) number of header lines in met file
365 (int) number of simdays in last simyear (truncated year: < 365)
RESTART
0 (flag) 1 = read restart; 0 = dont read restart
1 (flag) 1 = write restart; 0 = dont write restart
0 (flag) 1 = use restart metyear; 0 = reset metyear
hhs.endpoint (filename) name of the input restart file
hhs.endpoint (filename) name of the output restart file
TIME_DEFINE
54 (int) number of meteorological data years
54 (int) number of simulation years
1961 (int) first simulation year
1 (flag) 1 = spinup run; 0 = normal run
6000 (int) maximum number of spinup years
CLIM_CHANGE
0.0 (degC) - offset for Tmax
0.0 (degC) - offset for Tmin
1.0 (degC) - multiplier for PRCP
1.0 (degC) - multiplier for VPD
1.0 (degC) - multiplier for RAD
CO2_CONTROL
0 (flag) 0=constant; 1=vary with file
290.0 (ppm) constant atmospheric CO2 concentration
CO2_from1961.txt (filename) name of the CO2 file
NDEP_CONTROL
0 (flag) 0=constant; 1=vary with file
0.000200 (kgN/m2/yr) wet+dry atmospheric deposition of N
Ndep_from1961.txt (filename) name of the N-dep file
SITE
30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 (%) sand percentage by volume in rock-free soil
50.0 50.0 50.0 50.0 50.0 50.0 50.0 50.0 50.0 50.0 (%) silt percentage by volume in rock-free soil
7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 (dimless) soil pH
248.0 (m) site elevation
46.95 (degrees) site latitude (- for S.Hem.)
0.20 (DIM) site shortwave albedo
9.00 (Celsius) mean annual air temperature
10.15 (Celsius) mean annual air temperature range
20.00 (mm) maximum height of pond water
-9999 (dimless) measured runoff curve number
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (g/cm3) bulk density
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at saturation
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at field capacity
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at wilting point
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 (m3/m3) SWC at hygroscopic water content
EPC_FILE
c3grass.epc (filename) EPC file name
W_STATE
0.0 (kg/m2) water stored in snowpack
1.0 (DIM) initial soil water as a proportion of field capacity
CN_STATE
0.001 (kgC/m2) first-year maximum leaf carbon
0.001 (kgC/m2) first-year maximum fine root carbon
0.001 (kgC/m2) first-year maximum fruit carbon
0.001 (kgC/m2) first-year maximum softstem carbon
0.001 (kgC/m2) first-year maximum live woody stem carbon
0.001 (kgC/m2) first-year maximum live coarse root carbon
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) coarse woody debris carbon
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, labile pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, unshielded cellulose pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, shielded cellulose pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) litter carbon, lignin pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, fast microbial recycling pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, medium microbial recycling pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, slow microbial recycling pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgC/m2) soil carbon, recalcitrant SOM (slowest)
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgN/m2) litter nitrogen, labile pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgN/m2) soil mineralized nitrogen, NH4 pool
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 (kgN/m2) soil mineralized nitrogen, NO3 pool
OUTPUT_CONTROL
spinup_ (filename) output prefix
0 (flag) writing daily output (0 = no; 1 = binary; 2 = ascii)
0 (flag) writing monthly average of daily output (0 = no; 1 = binary; 2 = ascii)
0 (flag) writing annual average of daily output (0 = no; 1 = binary; 2 = ascii)
0 (flag) writing annual output (0 = no; 1 = binary; 2 = ascii)
1 (flag) for on-screen progress indicator
DAILY_OUTPUT
14 number of daily output variables
49 tsoil_2cm
52 tsoil_20cm
2599 vwc_2cm
2602 vwc_20cm
2751 m_soilstress
671 vegc_to_SNSC
171 evapotransp
2520 proj_LAI
407 STDBc_above
3004 daily_n2o
3008 daily_nbp
3009 daily_gpp
3014 daily_tr
2502 n_actphen
ANNUAL_OUTPUT
12 number of annual output variables
0 onday
1 offday
3000 annprcp
3001 anntavg
3002 annrunoff
2765 ytd_maxplai
3024 cum_nee
3033 cum_Closs_THN_w
3035 cum_Closs_MOW
3037 cum_yieldC_HRV
3059 NH4_top10
3060 NO3_top10
-------------------
MANAGEMENT_SECTION
-------------------
PLANTING
0 (flag) do PLANTING? 0=no; 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) PLANTING day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m) germination depth
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(n/m2) number of seedlings
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(g/1000n) weight of 1000-seed
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) C content of seed
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) emergence rate
THINNING
0 (flag) do THINNING? 0=no; 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) THINNING day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) thinning rate of woody plant material
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) thinning rate of non-woody plant material
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of stem
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of leaf
MOWING
0 (flag) do MOWING? 0=no; 1=yes; filepath=reading from file
0 (flag) mowing method? 0 - fixday method, 1 - fixvalue method
-9999 (m2/m2) fixed value of the LAI before MOWING (fixvalue method)
-9999 (m2/m2) fixed value of the LAI after MOWING (fixvalue method)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) MOWING day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m2/m2) value of the LAI after MOWING (fixday method)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of plant material
GRAZING
0 (flag) do GRAZING? 0=no, 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) first day of GRAZING
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) last day of GRAZING
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kg/LSU) weight equivalent of one unit
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(LSU/ha) animal stocking rate: Livestock Units per hectare
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kg dry matter/LSU) daily ingested dry matter
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) trampling effect (standing dead biome to litter)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) ratio of DM intake formed excrement
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) ratio of excrement returning to litter
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) carbon content of dry matter
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) N content of manure
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) C content of manure
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN2O-N:kgN) manure emission factor for direct N2O (T10.21)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/1000 kg animal mass/day) default N excretion rate (T10.19)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgCH4/LSU/yr*) manure emission factor for CH4 (T10.14)
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgCH4/LSU/yr*) fermentation emission factor for CH4 (T10.11)
HARVESTING
0 (flag) do HARVESTING? 0=no, 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) HARVESTING day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgC/m2) soft stem C content after HARVESTING
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) transported part of plant material
PLOUGHING
0 (flag) do PLOUGHING? 0=no, 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) PLOUGHING day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m) PLOUGHING depth
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(prop) dissolving coefficient of ploughed biome to litter
FERTILIZING
0 (flag) do FERTILIZING? 0=no, 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) FERTILIZING day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(m) fertilizing depth
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN fertilizer/ha/day) amount of fertilizer
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) dry matter content of fertilizer
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/kg fertilizer) nitrate content of fertilizer
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/kg fertilizer) ammonium content of fertilizer
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN/kg fertilizer) organic nitrogen content of fertilizer
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgC/kg fertilizer) organic carbon content of fertilizer
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) labile fraction of fertilizer
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) unshielded cellulose fraction of fertilizer
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) shielded cellulose fraction of fertilizer
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) lignin fraction of fertilizer
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgN2O-N:kgN) emission factor for N-additions
IRRIGATION
0 (flag) do IRRIGATION? 0=no, 1=yes; filepath=reading from file
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(yday) IRRIGATION day
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(kgH2O/m2/day) amount of water
-9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 *(%) utilization coefficient of water
END_INIT