Address sites with multiple flow measurements

[lindsayplatt]

There are some sites that have multiple timeseries. How do we switch from one to the next? And determine which stat to use?

sites_with_mult_flow_ts <- c("01574500", "03292555", "06903900", "08188590")

dv_stat <- 
  dataRetrieval::readNWISstat(
    siteNumbers = sites_with_mult_flow_ts,
    parameterCd = "00060",
    statType = paste0("P", c("05","10","20","25","50","75","80","90","95")),
    statReportType="daily"
  ) %>% 
  filter(month_nu == 9, day_nu >= 12 & day_nu <= 19)

dv_stat %>% select(site_no, day_nu, ts_id, loc_web_ds) %>% arrange(site_no, day_nu)

dv_data <- 
  dataRetrieval::readNWISdata(
    service = "dv",
    siteNumbers = problem_sites,
    parameterCd = "00060",
    startDate = "2018-09-12",
    endDate = "2018-09-19"
  ) %>% dataRetrieval::renameNWISColumns()