Calculate thermal time using the hourly temperature (non daily temperature)

thermalTimeHourly(timestamp, temperature, x_temp, y_temp)

Arguments

timestamp

The timestamp of weather records

temperature

The temperature

x_temp

The cardinal temperatures

y_temp

The effective thermal time

Value

A data frame with daily thermal time

Examples

met_file <- system.file("extdata/WeatherHourly.csv", package = "weaana")
hourly <- read.csv(met_file, as.is = TRUE) 

hourly$timestamp <- as.POSIXct(hourly$timestamp, format = "%Y-%m-%dT%H:%M:%SZ")
x_temp <- c(0, 20, 35)
y_temp <- c(0, 20, 0)
thermalTimeHourly(hourly$timestamp, hourly$temperature, x_temp, y_temp)
#> # A tibble: 43 × 2
#>    date       value
#>    <date>     <dbl>
#>  1 2017-12-20  16.3
#>  2 2017-12-21  14.0
#>  3 2017-12-22  12.9
#>  4 2017-12-23  12.7
#>  5 2017-12-24  15.3
#>  6 2017-12-25  14.4
#>  7 2017-12-26  12.4
#>  8 2017-12-27  12.1
#>  9 2017-12-28  11.6
#> 10 2017-12-29  13.9
#> # … with 33 more rows