station-gap.R

library(tidyhydat)
library(dplyr)
library(purrr)

rl_stations <- realtime_stations(prov_terr_state_loc = "BC")
hy_stations <- hy_stations(prov_terr_state_loc = "BC")

## all stations that are both in hydat and have a realtime record
bc_stations <- rl_stations %>%
  inner_join(hy_stations, by = "STATION_NUMBER")

find_station_gap <- function(station_number) {

  data.frame(
    station_number = station_number,
    station_name = hy_stations(station_number)$STATION_NAME,
    earliest_rl_date = as.Date(min(realtime_dd(station_number)$Date)),
    latest_hydat_date = max(hy_daily(station_number)$Date)
  )
}

station_gaps <- map_df(bc_stations$STATION_NUMBER, find_station_gap)

write.csv(station_gaps,
          file = paste0("station-gap-", as.Date(hy_version()$Date), ".csv"),
          row.names = FALSE)