Calculate temperatures and thermal stress from weather data (daily minimum and maximum temperatures).

1 2 3 | ```
thermalStressSeasonal(criticalTemp, dailyWeather, trialData, trialLocs, startEnd = c("START","END"), zenith=96)
thermalStressDaily(criticalTemp, Tmax, Tmin, sunr, TminNext)
temperature(Time, Tmax, Tmin, sunr, TminNext)
``` |

`criticalTemp` |
temperature above which thermal stress is experienced |

`dailyWeather` |
data frame containing daily weather data |

`trialData` |
data frame containing trial data, with start (START) and end (END) dates of the growing season |

`trialLocs` |
data frame containing trial locations |

`startEnd` |
character vector of two, indicating which column to select for start and end dates from trialData |

`zenith` |
zenith angle, see ?sunrise |

`Tmax` |
single number or vector with daily maximum temperature values |

`Tmin` |
single number or vector with daily minimum temperature values |

`sunr` |
single number or vector with daily time of sunrise |

`TminNext` |
single number or vector with minimum temperature values for the next day |

`Time` |
single number or vector with (local mean) time of day in hours |

single number or vector with thermal stress duration (in hours) or temperature values

Jacob van Etten

1 2 3 4 5 6 7 8 9 10 11 12 | ```
# Illustration of the thermal stress concept
s <- sunrise(25, 6, 1990, -74.3, 40.9, zenith=(90+50/60))
temp <- temperature(0:24, 31, 15, s, 15)
plot(ts(temp,0,24)) # Modelled course of temperature during the day
abline(30,0) # Line with critical temperature of 30 degrees Celsius
# Now let's calculate for how long the temperature exceeded 30 degrees this day
thermalStressDaily(criticalTemp=30, Tmax=31, Tmin=15, sunr=s, TminNext=15)
# Example with thermalStressSeasonal()
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.