R/wbt.Stull.R

In anacv/HeatStress: Calculate heat stress indices

#' Calculation of wet bulb temperature, following Stull's method.
#' 
#' Calculation of wet bulb temperature from temperature and relative humidity.
#' 
#' @param tas vector of air temperature in degC.
#' @param hurs vector of relative humidity in \%.
#' 
#' @return Wet bulb temperature in degC.
#' @author A.Casanueva (15.08.2016).
#' @details Formulation from Stull 2011, Journal of Applied Meteorology and Climatology.
#' 
#' @export
#' 
#' @examples \dontrun{ 
#' # load the meteorological variables for example data in Salamanca:
#' data("data_obs") 
#' wbt <- wbt.Stull(data_obs$tasmean, hurs=data_obs$hurs)
#' }
#' 

wbt.Stull <- function(tas,hurs){

  # assertion statements
  assertthat::assert_that(length(hurs)==length(tas), msg="Input vectors do not have the same length")
  
  # Constants
  c1 <- 0.151977
  c2 <- 8.313659
  c3 <- 1.676331
  c4 <- 0.00391838
  c5 <- 0.023101
  c6 <- 4.686035
  
  # calculation of the wet bulb temperature in degC (Stull 2011)
  wetbulb <- tas * atan(c1 * sqrt(hurs + c2)) + atan(tas + hurs) - atan(hurs - c3) + c4 * (hurs^(3/2)) * atan(c5 * hurs) - c6 
  
  return(wetbulb)

}