R/utci.r

In alfcrisci/rBiometeo: Biometeorological Functions in R

#' utci
#'
#' Calculate Universal Thermal Climate Index ( UTCI) in degC.
#'
#' @param t  numeric Air temperature in degC.
#' @param rh numeric Air Relative humidity.
#' @param wind numeric Windspeed in meters per second.
#' @param tr  numeric  Mean radiant temperature in degC.
#' @return UTCI index n degC
#' @references Brode P, Jendritzky G, Fiala D and Havenith G, 2011, The Universal Thermal Climate Index UTCI in Operational Use, International Journal of Biometeorology.
#' @author  Istituto per la Bioeconomia CNR Firenze Italy  Alfonso Crisci \email{alfonso.crisci@@ibe.cnr.it}
#' @keywords  UTCI 
#' 
#' @export
#'
#'
#'
#'
utci <-function(t,rh,wind,tr=NULL) {
                         if(is.null(tr)) {tr=t} 
                         ct$assign("ta", as.array(t))
                         ct$assign("rh", as.array(rh))
                         ct$assign("wind", as.array(wind))
                         ct$assign("tmrt", as.array(tr))
                         ct$eval("var res=[]; for(var i=0, len=ta.length; i < len; i++){ res[i]=UTCI(ta[i],rh[i],wind[i],tmrt[i])};")
                         res=ct$get("res")
                         res=ifelse(is.na(t),NA,res)
                         return(ifelse(res==9999,NA,res))
}