R/processWTH.R

In allogamous/EnvRtype: Envirotyping Tools for Crop Research and Genetic Improvment

#'@title  Basic Processing Tools for Enriching get_weather outputs
#'
#'
#' @description A wraper of \code{param_padiation()}, \code{param_atmospheric()}, and \code{param_temperature()} for evaluating get_weather() datasets. Calculates a series of parameters based on the \code{get_weather()} object.
#' @author Germano Costa Neto
#'
#' @param env.data data.frame. A \code{get_weather()} output.
#' @param PREC character. Indicates the column of precipitation.
#' @param Tdew character. Indicates the column of dew/frost.
#' @param Tmin character. Indicates the column of minimum temperature.
#' @param Tmax character. Indicates the column of maximum temperature.]
#' @param Tbase1 numeric. Minimum cardinal value for temperature base for phenological development (Celsius).
#' @param Tbase2 numeric. Maximum cardinal value for temperature base for phenological development (Celsius).
#' @param Topt1 numeric. Lower temperature bound for phenological development (Celsius).
#' @param Topt2 numeric. Upper temperature bound for phenological development (Celsius).
#' @param RH character. Indicates the column of reative humidity.
#' @param Rad character. Indicates the column of solar radiation. This parameter can be calculated from \code{Param_Radiation()}.
#' @param G numeric. Flux of heat conducted into the ground. Default is 0.
#' @param alpha Alpha of Priestley & Taylor's (1972) equation. Default is 1.26, which fit data from most sources.
#' @param Alt numeric. Elevation above sea level (meters)
#'
#' @return
#' Returns a \code{get_wheather()} object with an additional set of parameters calculated from the nasapower data.
#'
#' @details
#' This function is old. P requires a dataframe with all parameters listed above. If any is missing, an error will be returned.
#' The estimated variables are:
#' \itemize{
#'  \item n: Actual duration of sunshine (hour)
#'  \item N: Daylight hours (hour)
#'  \item RTA: Extraterrestrial radiation (MJ/m^2/day)
#'  \item SRAD: Solar radiation (MJ/m^2/day)
#'  \item SPV: Slope of saturation vapour pressure curve (kPa.Celsius)
#'  \item VPD: Vapour pressure deficit (kPa)
#'  \item ETP: Potential Evapotranspiration (mm.day)
#'  \item PEPT: Deficit by Precipitation (mm.day)
#'  \item GDD: Growing Degree Day (oC/day)
#'  \item FRUE: Effect of temperature on radiation use efficiency (from 0 to 1)
#'  \item T2M_RANGE: Daily Temperature Range (oC day)
#' }
#'
#' @examples
#'\dontrun{
#' ### Fetching weather information from NASA-POWER
#' env.data = get_weather(lat = -13.05, lon = -56.05)
#'
#' ### Returning calculated parameters merged to the \code{get_weather()} dataframe
#' processWTH(env.data)
#'}
#' @export

processWTH <- function(env.data,PREC=NULL, Tdew=NULL,
                       Tmin=NULL, Tmax=NULL, RH=NULL,
                       Rad=NULL, G=NULL, Alt=600,
                       Tbase1=9,Tbase2=45, Topt1=26,Topt2=32,alpha=1.26){

  # computing radiation paramters
  env.data <-param_radiation(env.data=env.data, merge = TRUE)
  # computing atmospheric paramters
  env.data <- param_atmospheric(env.data=env.data,
                                PREC=PREC,Tdew=Tdew,Tmin = Tmin, Tmax=Tmax,
                                Rad = Rad, G = G, Alt = Alt, alpha=alpha,
                                merge = TRUE)
  # computing temperature parameters
  env.data <- param_temperature(env.data=env.data,Tbase1=Tbase1,Tmax=Tmax,Tmin=Tmin,
                                Tbase2=Tbase2,Topt1=Topt1, Topt2=Topt2,merge = TRUE)

  return(env.data)
}