Added fullcats for CONUS and HI here

workspace/00_get_data.Rmd

@@ -19,6 +19,7 @@ library(tidyr)
 library(dplyr)
 library(sf)
 
+source("R/utils.R")
 
 if(!dir.exists("data")) {dir.create("data")}
 if(!dir.exists("bin")) {dir.create("bin")}
@@ -274,6 +275,25 @@ if(!file.exists(waterbodies_path)) {
 out_list <- c(out_list, list(waterbodies_path = waterbodies_path))
 ```
 
+```{r NHDPlusV2 Full cats}
+fullcat_path <- file.path(nhdplus_dir, "nhdcat_full.rds")
+islandcat_path <- file.path(data_paths$islands_dir, "nhdcat_full.rds")
+
+# Create full cat dataset
+if(!file.exists(fullcat_path)){
+  
+  cat_tab <- cat_rpu(staged_nhd$catchment, nhdplus_gdb)
+  saveRDS(cat_tab, fullcat_path)
+  
+  island_cats <- file.path(data_paths$islands_dir, "NHDPlusNationalData/nhdplus_catchment.rds")
+  island_tab <- cat_rpu(island_cats, data_paths$islands_gdb)
+  saveRDS(island_tab, islandcat_path)
+}
+  
+out_list <- c(out_list, list(fullcats_table = fullcat_path)) 
+out_list <- c(out_list, list(islandcats_table = islandcat_path))
+```
+
 ```{r NHDPlusV2 FDR_FAC}
 # NHDPlus FDR/FAC grids available by raster processing unit