R/reformat_ERA.R

In AgronomicForecastingLab/GEFS_pSIMS:

#' Reformat ERA5 weather dataframe list
#'
#' Converts temp variables from K to C. Derives wind speed from wind u- and
#' v-components.
#'
#' @param data List of ERA5 weather data dataframes
#' @return A list of reformatted dataframes
#' @export
reformat_ERA_data <- function(data) {
  drops <- c("mx2t","mn2t", "d2m", "t2m", "tp", "ssrd", "u10", "v10")
  for (x in seq_len(length(data))) {
    data[[x]]$tmax <- units::drop_units(data[[x]]$mx2t) - 273.15
    data[[x]]$tmin <- units::drop_units(data[[x]]$mn2t) - 273.15
    data[[x]]$tdew <- units::drop_units(data[[x]]$d2m) - 273.15
    data[[x]]$tavg <- units::drop_units(data[[x]]$t2m)- 273.15
    data[[x]]$prcp <- units::drop_units(data[[x]]$tp) *1000 #m to mm
    data[[x]]$srad <- units::drop_units(data[[x]]$ssrd)/86400 #J/m2 to wm2
    data[[x]]$wind <- sqrt(units::drop_units(data[[x]]$u10)^2 + units::drop_units(data[[x]]$v10)^2)
  }
  data <- lapply(data, function(x) x[!(names(x) %in% drops)])
  return(data)
}