In mikejohnson51/NWM: Provides access to the NOAA National Water Model output

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  fig.width = 6, 
  fig.height = 4
)

library(nwm)
library(rasterVis)

Define the AOI

Lets define the AOI:

#Define AOI
AOI = getAOI(list("Colorado Springs", 15, 15))

#Get some flowlines
nhd = getNHD(AOI)

# Visualize
AOI::check(AOI) %>% addPolylines(data = nhd, popup = paste0("COMID: ", nhd$comid), weight  = nhd$streamorde)

# Describe contexts
AOI %>% describe()

Channel Output Example

xx = system.time({

filelist = getFilelist(config = "medium_range",
                     date = "20180823",
                     type = "channel",
                     t =  c(0),
                     f = seq(3,33,3))

flows <- nwm::downloadNWM(AOI, filelist, param = "streamflow")
})

message(dim(flows$streamflow)[1], " streamflow records extracted in ", round(xx[3], 2), " seconds")

head(flows$streamflow, 5)

data = flows$streamflow
reach1 = unique(data$COMIDS)[1]
data_sub = data[data$COMIDS == reach1,]

plot(x = data_sub$DateTime, y = data_sub$streamflow, type = 'b',
      main = "National Water Model: Streamflow",
      xlab = "dateTime - UTC", 
      ylab = "Streamflow (cfs)")

Land output example

xx = system.time({

  filelist = getFilelist(config = "medium_range",
                         date = "20180823",
                         type = "land",
                         t =  c(0),
                         f = seq(3,33,3))

  et <- nwm::downloadNWM(AOI, filelist, param = "accet")
})

message(dim(et$accet)[3], " ",  dim(et$accet)[1] *dim(et$accet)[2],  " km2 land rasters extracted in ", round(xx[3], 2), " seconds")

data = et$accet[[1:4]]

rasterVis::levelplot(data, main = "National Water Model: Accumulated ET", 
                     names.attr = as.character(nwm::getGridTime(data)))

Forcing Examples

xx = system.time({

  filelist = getFilelist(config = "medium_range",
                         date = "20180823",
                         type = "forcing",
                         t =  c(0),
                         f = seq(3,33,3))

  prcp <- nwm::downloadNWM(AOI, filelist, param = "rainrate")
})

message(dim(prcp$rainrate)[3], " ",  dim(prcp$rainrate)[1] *dim(prcp$rainrate)[2],  " km2 land rasters extracted in ", round(xx[3], 2), " seconds")

data = prcp$rainrate[[1:4]]

rasterVis::levelplot(data, main = "National Water Model: Rainrate", 
                     names.attr = as.character(nwm::getGridTime(data)))