R/update_wx_stns.R

In adamdsmith/nrsmisc: Miscellaneous functions for work

#' Update list of US weather stations
#'
#' Meant for developer use only

update_wx_stns <- function() {

  # # Don't use GHCN stations for now (except those in US and major territories)
  # ### GHCN STATIONS
  # ghcn_dat <- "C:/Users/adsmith/Downloads/ghcnd-stations.txt"
  # ghcn <- readr::read_fwf(ghcn_dat,
  #                         readr::fwf_empty(ghcn_dat,
  #                                          col_names = c("GHCN", "lat", "lon", "elev_m", "name")),
  #                         col_types = "cddnc---") %>%
  #   mutate(elev_ft = round(elev_m * 3.28084)) %>%
  #   select(GHCN, name, lat, lon, elev_ft)

  ### CO-OP STATIONS
  coop_dat <- "https://www.ncdc.noaa.gov/homr/file/coop-stations.txt"
  coop <- readr::read_fwf(coop_dat,
                          readr::fwf_empty(coop_dat)) %>%
    mutate(WBAN = NA_character_, FAA = NA_character_) %>%
    select(.data$WBAN, COOP = .data$X2, .data$FAA, GHCN = .data$X4, name = .data$X5,
           lat = .data$X12, lon = .data$X13, elev_ft = .data$X15) %>%
    # Needn't save GHCN for now
    select(-.data$GHCN)

  ### ASOS STATIONS
  asos_dat <- "https://www.ncdc.noaa.gov/homr/file/asos-stations.txt"
  asos <- readr::read_fwf(asos_dat,
                          readr::fwf_empty(asos_dat, skip = 2),
                          skip = 2) %>%
    select(WBAN = .data$X2, COOP = .data$X3, FAA = .data$X4, name = .data$X5,
           lat = .data$X11, lon = .data$X12, elev_ft = .data$X13) %>%
    filter(!is.na(.data$WBAN)) %>%
    mutate(ASOS = TRUE)
  # Manually fix one missing station FAA ID (19 May 2018)
  asos[grep("CHERRY POINT MCAS ASOS", asos$name), "FAA"] <- "NKT"

  ### AWOS STATIONS
  awos_dat <- "https://www.ncdc.noaa.gov/homr/file/awos-stations.txt"
  awos <- readr::read_fwf(awos_dat,
                          readr::fwf_empty(awos_dat, skip = 2),
                          skip = 2) %>%
    select(WBAN = .data$X2, COOP = .data$X3, FAA = .data$X4, name = .data$X5,
           lat = .data$X9, lon = .data$X10, elev_ft = .data$X11) %>%
    filter(!is.na(.data$WBAN),
           !.data$WBAN %in% asos$WBAN) # If ASOS too, keep ASOS

  # Reduce AWOS stations already in COOP
  awos <- awos %>%
    filter(!.data$COOP %in% coop$COOP)

  # Remove COOPs that are also ASOS (prefer to keep ASOS)
  coop <- coop %>%
    filter(!.data$COOP %in% asos$COOP)

  # Consolidate
  all_wx <- bind_rows(coop, awos, asos) %>%
    mutate(ASOS = ifelse(is.na(.data$ASOS), FALSE, TRUE),
           id = ifelse(!is.na(.data$FAA), .data$FAA,
                       ifelse(!is.na(.data$WBAN), .data$WBAN, .data$COOP))) %>%
    select(.data$id, .data$name, .data$lon, .data$lat, .data$ASOS)

  saveRDS(all_wx, "./inst/extdata/wx_stations.rds")

}