Calculate thermal time using cardinal temperatures

thermalTime(weather, x_temp, y_temp, method = NULL)

Arguments

weather: WeaAna object
x_temp: The cardinal temperatures
y_temp: The effective thermal time
method: The method to calculate thermal time. The default method is ( maxt + mint ) / 2 - base. The three hour temperature methods will be used if method = '3hr'

Value

A data.frame with three columns: year, day and thermalTime.

Examples

met_file <- system.file("extdata/WeatherRecordsDemo1.met", package = "weaana")
records <- readWeatherRecords(met_file)
x_temp <- c(0, 26, 34)
y_temp <- c(0, 26, 0)
res <- thermalTime(records, x_temp, y_temp)
head(res)
#>   year day thermalTime
#> 1 1981   1     19.5500
#> 2 1981   2     20.8000
#> 3 1981   3     23.8875
#> 4 1981   4     22.6500
#> 5 1981   5     25.3000
#> 6 1981   6     24.0000
res <- thermalTime(records, x_temp, y_temp, method = "3hr")
head(res)
#>   year day thermalTime
#> 1 1981   1    19.55000
#> 2 1981   2    17.77214
#> 3 1981   3    14.70688
#> 4 1981   4    21.76661
#> 5 1981   5    21.36310
#> 6 1981   6    19.09683