The sim components in the package lintul5 implement the Lintul5 algorithms from (L5). They are designed to reproduce same results (up to rounding errors) as the original FST program.
The original model is subdivided into multiple sim components (which should be used in the following order to be as close as possible to the original algorithm).
Daylength calculation is performed by the transformer component AstronomicParametersTransformer
RadiationUseEfficiency component calculates the radiation use efficiency
The component calculates radiation use efficiency as well as Co2 and temperature correction factor. Calculations taken from original's CROPP() routine.
It takes inputs from the component Phenology.
It's outputs are used by Lintul5.
For documentation of the algorithms please consult the orginal documentation (L5)
Calculation of used temperature is made customizable by a variable cDayTempFactor. $$ \begin{eqnarray} DTemp & = & TMAX - cDayTempFactor \cdot (TMAX - TMIN) \\ \label{dtemp} & = & (1-cDayTempFactor) \cdot TMAX + (cDayTempFactor)\cdot TMIN \end{eqnarray} $$
As some parameters are given as interpolation tables, there have been added additional scale factors to make these values changeable for calibrations:
(L5) Joost Wolf, User guide for LINTUL5, Wageningen UR, Wageningen, 2012
Content Type | Name | Description | Data Type | Unit | Min Value | Max Value | Default Value |
---|---|---|---|---|---|---|---|
constant | cCO | atmosferic CO2 concentration | DOUBLE | ppm | - | - | 0.0 |
constant | cCOTB | table with correction factor of RUE as a function of atmospheric CO2 concentration | DOUBLEARRAY | 1 | - | - | - |
constant | cCOTableCo2 | CO2 concentration for correction of radiation use efficiency as a function of atmospheric (c.f. COTB) | DOUBLEARRAY | ppm | - | - | 40.0 360.0 720.0 1000.0 2000.0 |
constant | cCOTableFactor | Correction of radiation use efficiency as a function of atmospheric CO2 concentration (c.f. COTB) | DOUBLEARRAY | 1 | - | - | 0.0 1.0 1.35 1.5 1.5 |
constant | cDayTempFactor | Calculation factor for temperature means T = TMAX - f*(TMAX-TMIN). Default 0.25 gives day temp, 0.5 gives mean temp. See $\eqref{dtemp}$ | DOUBLE | 1 | - | - | 0.25 |
constant | cRUETB | radiation use efficiency as function of DVS | DOUBLEARRAY | g/MJ | - | - | - |
constant | cRUETableDVS | DVS for radiation use efficiency for biomass production (c.f. RUETB) | DOUBLEARRAY | 1 | - | - | 0.0 1.0 1.3 2.0 |
constant | cRUETableRUE | Radiation use efficiency for biomass production as function of DVS (c.f. RUETB) | DOUBLEARRAY | g/MJ | - | - | 3.0 3.0 3.0 0.4 |
constant | cScaleFactorRUE | Scales the y-values of RUETB (for sensitivity analysis / calibration) | DOUBLE | 1 | - | - | 1.0 |
constant | cTMNFTB | reduction factor of RUE as function of low minimum temperature | DOUBLEARRAY | 1 | - | - | - |
constant | cTMNFTableFactor | Reduction factor of RUE as function of low min. temp. (c.f. TMNFTB) | DOUBLEARRAY | 1 | - | - | 0.0 0.0 1.0 1.0 |
constant | cTMNFTableMinTemperature | Min Temp for reduction factor of RUE (c.f. TMNFTB) | DOUBLEARRAY | °C | - | - | -5.0 0.0 3.0 30.0 |
constant | cTMPFTB | reduction factor of RUE as function of mean daytime temperature | DOUBLEARRAY | 1 | - | - | - |
constant | cTMPFTableFactor | Reduction factor of RUE as function of mean daytime temp (c.f. TMPFTB) | DOUBLEARRAY | 1 | - | - | 0.0 0.0 0.6 1.0 1.0 0.0 0.0 |
constant | cTMPFTableMeanTemperature | Mean daytime temp for reduction factor of RUE (c.f. TMPFTB) | DOUBLEARRAY | °C | - | - | -1.0 0.0 10.0 15.0 30.0 35.0 40.0 |
input | iDVS | initial development stage of crop (from 0 to 2) | DOUBLE | 1 | - | - | 0.0 |
input | iTMAX | maximal air temperature during day (output of routine WEATHR) | DOUBLE | °C | - | - | 0.0 |
input | iTMIN | minimal air temperature during day (output of routine WEATHR) | DOUBLE | °C | - | - | 0.0 |
out | RTMCO | overall correction factor for RUE in dependence of both CO2 concentration and non-optimal daytime and minimal temperatures | DOUBLE | 1 | - | - | 0.0 |
out | RUE | radiation use efficiency | DOUBLE | g/MJ | - | - | 0.0 |