bring in merit data

workspace/alaska_prep.Rmd

+---
+title: "NHD Navigate"
+output: html_document
+---
+
+Project: GFv2.0  
+Script purpose: Get Alaska MERIT data ready for GF processing.
+Date: 12/24/2019  
+Author: [Dave Blodgett](dblodgett@usgs.gov)
+Notes: 
+
+```{r setup}
+library(nhdplusTools)
+library(dplyr)
+library(sf)
+library(hyRefactor)
+
+source("R/utils.R")
+
+data_paths <- jsonlite::read_json(file.path("cache", "data_paths.json"))
+
+ak_rivers <- read_sf(data_paths$merit_rivers)
+ak_cats <- read_sf(data_paths$merit_catchments)
+```
+

workspace/cache/data_paths.json

@@ -9,5 +9,7 @@
   },
   "fac": {
     "rpu_06a": "data/fdrfac/NHDPlusMS/NHDPlus06/NHDPlusFdrFac06a/fac"
-  }
+  },
+  "merit_catchments": "data/merged_AK_MERIT_Hydro/cat_pfaf_78_81_82_MERIT_Hydro_v07_Basins_v01.shp",
+  "merit_rivers": "data/merged_AK_MERIT_Hydro/riv_pfaf_78_81_82_MERIT_Hydro_v07_Basins_v01.shp"
 }

workspace/get_data.Rmd

@@ -111,6 +111,20 @@ out$fac[[paste0("rpu_", rpu)]] <- file.path(out_path, "fac")
 out_list<- c(out_list, out)
 ```
 
+```{r}
+merit_dir <- file.path(data_dir, "merged_AK_MERIT_Hydro")
+
+if(!dir.exists(merit_dir)) {
+  message("Talk to John Hammond to get the merged AK MERIT data.")
+} else {
+  out <- list(merit_catchments = file.path(merit_dir, 
+                                           "cat_pfaf_78_81_82_MERIT_Hydro_v07_Basins_v01.shp"),
+              merit_rivers = file.path(merit_dir, 
+                                       "riv_pfaf_78_81_82_MERIT_Hydro_v07_Basins_v01.shp"))
+  out_list <- c(out_list, out)
+}
+```
+
 Don't edit the following chunk.
 
 ```{r}