net.simplace.sim.components.models.gecros.GecrosFunctions

public class GecrosFunctions extends Object

Here you find functions that are part of the Gecros algorithm and are used by the related sim components. WIKI_START == Modifications == In the function TOTPT the DTR calculation was corrected. By setting the parameter legacymode to true one can switch to the original calculation method. == Reference: == (Y) Yin, X. and van Laar, H.H. (2005): Crop Systems Dynamics. An ecophysiological simulation model for genotype-by-environment interactions. Wageningen Academic Publishers, Wageningen, The Netherlands. WIKI_END

See Also:

Constructor Summary

Constructors

Constructor

Description

GecrosFunctions()
Method Summary

Modifier and Type

Method

Description

static double

BETAF(double DVR, double TE, double TX, double TI)

Calculates the dynamics of expected growth of sinks, based on the beta sigmoid growth equation

static double

KDIFF(double LAI, double BL, double SCP)

Calculates extinction coefficient for diffuse radiation

static double

PHENO(double DS, boolean SLP, double DDLP, double SPSP, double EPSP, double PSEN, double MTDV, double MTDR, double TDU)

Calculates phenological development rate.

static double

RLAIC(double DS, double SLA0, double RWLV, double LAI, double KN, double NLV, double RNLV, double SLNB, double RSLNB)

Calculates the daily increase of leaf area index (m2 leaf/m2 ground/day)

static double[]

RNACC(double FNSH, double NUPT, double RWST, double STEMNC, double LNCMIN, double RNCMIN, double LNC, double RNC, double NLV, double NRT, double WLV, double WRT, double DELT, double CB, double CX, double TM, double DS, double SEEDNC, double RWSO, double LNLV, double LNRT)

Calculates rate of N accumulation in organs

static double[]

SINKG(double DS, double SSG, double TOTC, double YG, double FD, double DCDR, double DCS, double DELT)

Calculates carbon demand for sink growth.

static double[]

TOTPT(double SC, double SINLD, double COSLD, double DAYL, double DSINBE, double DDTR, double TMAX, double TMIN, double DVP, double WNM, boolean C3C4, double LAI, double TLAI, double HT, double LWIDTH, double RD, double SD1, double RSS, double BLD, double KN, double KW, double SLN, double SLNT, double SLNN, double SLNMIN, double DWSUP, double CO2A, double LS, double EAJMAX, double XVN, double XJN, double THETA, double WCUL, double FVPD, boolean legacymode)

Calculates daily total gross photosynthesis and transpiration by performing a Gaussian integration over time.

static double

TUNIT(double DS, double TMAX, double TMIN, double DIF, double DAYL, double TBD, double TOD, double TCD, double TSEN)

Calculates the daily amount of thermal day

Methods inherited from class java.lang.Object
clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Constructor Details
- GecrosFunctions
  
  public GecrosFunctions()
Method Details
- TUNIT
  
  public static double TUNIT(double DS, double TMAX, double TMIN, double DIF, double DAYL, double TBD, double TOD, double TCD, double TSEN)
  
  Calculates the daily amount of thermal day
  
  Parameters:
  
  DS - Development stage
  
  TMAX - Daily maximum temperature oC
  
  TMIN - Daily minimum temperature oC
  
  DIF - Daytime plant-air temperature differential oC
  
  DAYL - Astronomic daylength (base = 0 degrees) h
  
  TBD - Base temperature for phenology oC
  
  TOD - Optimum temperature for phenology oC
  
  TCD - Ceiling temperature for phenology oC
  
  TSEN - Curvature for temperature response
  
  Returns:
  
  Daily thermal-day unit [-]
- PHENO
  
  public static double PHENO(double DS, boolean SLP, double DDLP, double SPSP, double EPSP, double PSEN, double MTDV, double MTDR, double TDU)
  
  Calculates phenological development rate.
  
  Parameters:
  
  DS - Development stage
  
  SLP - Crop type(true for short-day, false for long-day)
  
  DDLP - Daylength for photoperiodism h
  
  SPSP - DS for start of photoperiod-sensitive phase
  
  EPSP - DS for end of photoperiod-sensitive phase
  
  PSEN - Photoperiod sensitivity (+ for SD, - for LD) h-1
  
  MTDV - Minimum thermal days for vegetative phase d
  
  MTDR - Minimum thermal days for reproductive phase d
  
  TDU - Daily thermal-day unit
  
  Returns:
  
  Development rate [d-1]
- RNACC
  
  public static double[] RNACC(double FNSH, double NUPT, double RWST, double STEMNC, double LNCMIN, double RNCMIN, double LNC, double RNC, double NLV, double NRT, double WLV, double WRT, double DELT, double CB, double CX, double TM, double DS, double SEEDNC, double RWSO, double LNLV, double LNRT)
  
  Calculates rate of N accumulation in organs
  
  Parameters:
  
  FNSH - Fraction of new N partitioned to shoot
  
  NUPT - Nitrogen uptake at a time step gN/m2/d
  
  RWST - Rate of stem weight g/m2/d
  
  STEMNC - Nitrogen concentration in stem gN/g
  
  LNCMIN - Minimum N concentration in leaf gN/g
  
  RNCMIN - Minimum N concentration in root gN/g
  
  LNC - Nitrogen concentration in leaf gN/g
  
  RNC - Nitrogen concentration in root gN/g
  
  NLV - Canopy (green)leaf N content gN/m2
  
  NRT - (living)root N content gN/m2
  
  WLV - Canopy (green)leaf weight g/m2
  
  WRT - (living)Root weight g/m2
  
  DELT - Time step of simulation d
  
  CB - Factor for initial N concent. of seed-fill
  
  CX - Factor for final N concent. of seed-fill
  
  TM - DS when transition from CB to CX is fastest
  
  DS - Development stage
  
  SEEDNC - Standard seed N concentration gN/g
  
  RWSO - growth rate of seed g/m2/d
  
  LNLV - Loss rate of NLV due to senescence gN/m2/d
  
  LNRT - Loss rate of NRT due to senescence gN/m2/d
  
  Returns:
  
  rate of N accumulation in root[gN/m2/d] / rate of N accumulation in stem [gN/m2/d] / rate of N accumulation in leaf [gN/m2/d] / Positive value of rate in leaf [gN/m2/d] /rate of N accumulation in seed(storage organ) [gN/m2/d]
- RLAIC
  
  public static double RLAIC(double DS, double SLA0, double RWLV, double LAI, double KN, double NLV, double RNLV, double SLNB, double RSLNB)
  
  Calculates the daily increase of leaf area index (m2 leaf/m2 ground/day)
  
  Parameters:
  
  DS - Development stage
  
  SLA0 - Specific leaf area constant m2 g-1
  
  RWLV - Rate of increment in leaf weight g m-2 d-1
  
  LAI - Leaf area index m2 m-2
  
  KN - Leaf nitrogen extinction coefficient m2 m-2
  
  NLV - Total leaf nitrogen content in a canopy g m-2
  
  RNLV - Rate of increment in NLV g m-2 d-1
  
  SLNB - Nitrogen content of bottom leaves g m-2
  
  RSLNB - Rate of increment in SLNB g m-2 d-1
  
  Returns:
  
  Rate of increment in leaf area index [m2 m-2d-1]
- BETAF
  
  public static double BETAF(double DVR, double TE, double TX, double TI)
  
  Calculates the dynamics of expected growth of sinks, based on the beta sigmoid growth equation
  
  Parameters:
  
  DVR - Development rate d-1
  
  TE - Stage at which sink growth stops
  
  TX - Stage at which sink growth rate is maximal
  
  TI - Stage of a day considered
  
  Returns:
  
  Relative expected growth of a sink at a day [d-1]
- SINKG
  
  public static double[] SINKG(double DS, double SSG, double TOTC, double YG, double FD, double DCDR, double DCS, double DELT)
  
  Calculates carbon demand for sink growth.
  
  Parameters:
  
  DS - Development stage
  
  SSG - Stage at which sink growth starts
  
  TOTC - Total carbon in a sink at end of its growth g C/m2
  
  YG - Growth efficiency g C/g C
  
  FD - Relative expected growth of a sink at a day d-1
  
  DCDR - Shortfall of C demand in previous days g C/m2
  
  DCS - Daily C supply for sink growth g C/m2/d
  
  DELT - Time step of integration d
  
  Returns:
  
  C demand of the current day [g C/m2/d] / Daily C demand for sink growth [g C/m2/d] / Flow of current assimilated C to sink [g C/m2/d]
- TOTPT
  
  public static double[] TOTPT(double SC, double SINLD, double COSLD, double DAYL, double DSINBE, double DDTR, double TMAX, double TMIN, double DVP, double WNM, boolean C3C4, double LAI, double TLAI, double HT, double LWIDTH, double RD, double SD1, double RSS, double BLD, double KN, double KW, double SLN, double SLNT, double SLNN, double SLNMIN, double DWSUP, double CO2A, double LS, double EAJMAX, double XVN, double XJN, double THETA, double WCUL, double FVPD, boolean legacymode)
  
  Calculates daily total gross photosynthesis and transpiration by performing a Gaussian integration over time. At five different times of the day, temperature and radiation are computed to determine assimilation and transpiration whereafter integration takes place. Resulting values are stored in the SimVariables PPCAN Potential canopy CO2 assimilation g m-2 d-1 O * APCANS Actual standing-canopy CO2 assimilation g m-2 d-1 O * APCANN APCANS with small plant-N increment g m-2 d-1 O * APCAN Actual canopy CO2 assimilation g m-2 d-1 O * PTCAN Potential canopy transpiration mm d-1 O * ATCAN Actual canopy transpiration mm d-1 O * PESOIL Potential soil evaporation mm d-1 O * AESOIL Actual soil evaporation mm d-1 O * DIFS Daytime average soil-air temp. difference oC O * DIFSU Daytime aver. sunlit leaf-air temp. diff. oC O * DIFSH Daytime aver. shaded leaf-air temp. diff. oC O * DAPAR Daily PAR absorbed by crop canopy J m-2 d-1 O *
  
  Parameters:
  
  SC - Solar constant J m-2 s-1
  
  SINLD - Seasonal offset of sine of solar height
  
  COSLD - Amplitude of sine of solar height
  
  DAYL - Astronomic daylength (base = 0 degrees) h
  
  DSINBE - Daily total of effective solar height s d-1
  
  DDTR - Daily global radiation J m-2 d-1
  
  TMAX - Daily maximum temperature oC
  
  TMIN - Daily minimum temperature oC
  
  DVP - Vapour pressure kPa
  
  WNM - daily average wind speed (>=0.1 m/s) m s-1
  
  C3C4 - Crop type (=true for C3, false for C4 crops)
  
  LAI - (green)Leaf area index m2 m-2
  
  TLAI - Total Leaf area index m2 m-2
  
  HT - Plant height m
  
  LWIDTH - Leaf width m
  
  RD - Rooting depth cm
  
  SD1 - Depth of evaporative upper soil layer cm
  
  RSS - Soil resistance,equivalent to leaf stomata s m-1
  
  BLD - Leaf angle from horizontal degree
  
  KN - Leaf nitrogen extinction coefficient m2 m-2
  
  KW - Windspeed extinction coefficient in canopy m2 m-2
  
  SLN - Average leaf nitrogen content in canopy g m-2
  
  SLNT - Top-leaf nitrogen content g m-2
  
  SLNN - Value of SLNT with small plant-N increment g m-2
  
  SLNMIN - Minimum or base SLNT for photosynthesis g m-2
  
  DWSUP - Daily water supply for evapotranspiration
  
  CO2A - Ambient CO2 concentration ml m-3
  
  LS - Lodging severity
  
  EAJMAX - Energy of activation for Jmax J mol-1
  
  XVN - Slope of linearity between Vcmax & leaf N umol/g/s
  
  XJN - Slope of linearity between Jmax & leaf N umol/g/s
  
  THETA - Convexity for light response of e-transport
  
  WCUL - Water content of the upper soil layer m3 m-3
  
  FVPD - Slope for linear effect of VPD on Ci/Ca (kPa)-1
  
  legacymode - If true, DTR is calculated with the old faulty formula from (Y) -
- KDIFF
  
  public static double KDIFF(double LAI, double BL, double SCP)
  
  Calculates extinction coefficient for diffuse radiation
  
  Parameters:
  
  LAI - Total leaf area index [ m2 m-2]
  
  BL - Leaf angle (from horizontal) [radians]
  
  SCP - Leaf scattering coefficient
  
  Returns:
  
  Diffuse radiation extinction coefficient [m2 m-2]

Class GecrosFunctions

Constructor Summary

Method Summary

Methods inherited from class java.lang.Object

Constructor Details

GecrosFunctions

Method Details

TUNIT

PHENO

RNACC

RLAIC

BETAF

SINKG

TOTPT

KDIFF