Calculate thermal time using cardinal temperatures

Usage

thermal_time(
  mint,
  maxt,
  x_temp = fababean$get("phenology.thermal_time.x"),
  y_temp = fababean$get("phenology.thermal_time.y"),
  method = fababean$get("phenology.thermal_time.method"),
  ...
)

Arguments

mint: The minimum temperature
maxt: The maximum temperature
x_temp: The cardinal temperatures using to calculate thermal time. The default values are c(2, 30, 35) defined in fababean$get("phenology.thermal_time.x").
y_temp: The effective thermal time at the corresponding cardinal temperatures. The default values are c(0, 28, 0) defined in fababean$get("phenology.thermal_time.y").
method: The method to calculate thermal time. The default is "HourlySinPpAdjusted". Other options are "3hr" and NULL (mean temperature).
...: Additional arguments passed to the underlying tidyweather::thermal_time function.

Value

The thermal time.

Details

The daily thermal time is calculated by interpolating the effective thermal time at three hour intervals between the daily minimum and maximum temperatures, and then averaging these values to get the daily thermal time. For more information about the fababean phenology model, see the APSIM Next Generation fababean Documentation.

Examples

mint <- c(0, 10)
maxt <- c(30, 40)
date <- as.Date(c("2024-01-01", "2024-01-02"))
thermal_time(mint, maxt, date = date, lat = -27)
#> [1] 16.64386 24.88400