local_eli/look_at_precipitation.R

In MazamaScience/RAWSmet: Download and Process data from Remote Automatic Weather Stations

setSpatialDataDir("~/Data/Spatial")
tsObj <- cefa_createRawsObject(nwsID = 451702)
start <- MazamaCoreUtils::parseDatetime(20051101, timezone = "UTC")
end <- MazamaCoreUtils::parseDatetime(20051110, timezone = "UTC")
tsObj$data$date <- tsObj$data$datetime
filtered <- dplyr::filter(a$data, date >= start & date < end)
filtered$precipPerHour <- filtered$precipitation / filtered$precipDuration
filtered <- filtered %>% dplyr::select(c("date", "precipitation", "precipDuration", "precipPerHour"))

openair::timePlot(filtered, pollutant = "precipPerHour")