#' @export
subset_nomads <- function(comids = NULL) {
comids = as.numeric(comids)
mapping = read.csv("data/current_nc/map.csv", stringsAsFactors = FALSE)
files = unique(sapply(comids, function(id) subset(mapping, max >= id & min <= id)$filename))
res = list()
for (file in files) {
res_subset <- fst::read.fst(paste0('data/current_nc/', file))
res_new = dplyr::filter(res_subset, COMID %in% comids)
res = dplyr::bind_rows(res, res_new)
}
res$agency_code = "NOAA NWS"
return(res)
}
subset_nomads_sql <- function(comids = NULL) {
comids = as.numeric(comids)
mapping = read.csv("data/current_nc/map.csv", stringsAsFactors = FALSE)
files = unique(sapply(comids, function(id) subset(mapping, max >= id & min <= id)$num))
res = list()
for (file in files) {
res_new<- dbGetQuery(con, paste("SELECT * FROM", paste0("flow_", file)))
res = dplyr::bind_rows(res, res_new)
}
res <- res %>%
reshape2::melt(id.vars=c("comid")) %>%
dplyr::mutate(variable = lubridate::ymd_h(variable)) %>%
dplyr::rename(dateTime = variable, Q_cfs = value) %>%
dplyr::arrange(comid)
res$agency_code = "NOAA NWS"
return(res)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.