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

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

• 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 [-]

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.

• 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]

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

• 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]

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)

• 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]

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

• 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]

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.

• 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]

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 *

• 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) -

public static double KDIFF(double LAI, double BL, double SCP);

Calculates extinction coefficient for diffuse radiation

• 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]