diff --git a/workspace/01_HI_prep.Rmd b/workspace/01_HI_prep.Rmd
index d46f07446fd4b145490e7d619b7c485fb2394a22..f7a9564fa6a10b35d7d817d5d828820306b0ef7c 100644
--- a/workspace/01_HI_prep.Rmd
+++ b/workspace/01_HI_prep.Rmd
@@ -129,7 +129,7 @@ if(needs_layer(src_gpkg, TE)) {
 
 ```{r NID POIs}
 # Derive NID locations for HI domain
-if(needs_layer(src_gpkg, NID)) {
+if(needs_layer(src_gpkg, NID_layer)) {
   message("formatting NID locations...")
   
   # Download NID file from website for 2019
@@ -148,8 +148,8 @@ if(needs_layer(src_gpkg, NID)) {
   NID_data$COMID <- nhd$ComID[nhd_index]
   
   # write to geopackage
-  write_sf(NID_data, src_gpkg, NID)
+  write_sf(NID_data, src_gpkg, NID_layer)
 } else {
-  NID_data <- read_sf(src_gpkg, NID)
+  NID_data <- read_sf(src_gpkg, NID_layer)
 }
 ```
diff --git a/workspace/02_HI_navigate.Rmd b/workspace/02_HI_navigate.Rmd
index fdb8f2a854d665c7acb15012afbc019bb361e00b..f91b909a93cb9a7d56a36c90957de4b90e2658a3 100644
--- a/workspace/02_HI_navigate.Rmd
+++ b/workspace/02_HI_navigate.Rmd
@@ -103,7 +103,7 @@ if(needs_layer(nav_gpkg, nhd_flowline)) {
 ```
 
 ```{r gageLoc_pois}
-if(needs_layer(nav_gpkg, gage_events)) {
+if(needs_layer(nav_gpkg, gage_events_layer)) {
   
   # Reads Streamstats gages/POIs
   ss_basins <- read_sf(file.path(data_paths$islands_dir, src_gpkg), "streamstats_basins") %>%
@@ -174,7 +174,7 @@ if(needs_layer(nav_gpkg, gage_events)) {
     filter(if(numCom == 1) TRUE else nrecords == max(nrecords, na.rm = T)) %>%
     select(COMID, Station)
 
-  write_sf(filter(gages_Data, Station %in% gages_Data$Station), nav_gpkg, gage_events)
+  write_sf(filter(gages_Data, Station %in% gages_Data$Station), nav_gpkg, gage_events_layer)
   
   # Create POIs 
   gage_POIs <- POI_creation(gage_COMIDs, filter(nhd_prep, dend == 1), "Gages")
@@ -244,7 +244,7 @@ if(all(is.na(nav_POIs$Type_Conf))) {
 if(all(is.na(nav_POIs$Type_NID))) {
   
   # Read in NID shapefile
-  NID_COMIDs <- read_sf(src_gpkg, NID) %>%
+  NID_COMIDs <- read_sf(src_gpkg, NID_layer) %>%
     select(COMID, NIDID) %>% 
     filter(COMID %in% filter(nhd_prep, dend == 1)$COMID) %>%
     st_drop_geometry() %>%
@@ -314,7 +314,7 @@ if(all(is.na(nav_POIs$Type_WBOut))) {
 
 ```{r}
 # Derive final POI set 
-if(needs_layer(nav_gpkg, inc_pois)) {
+if(needs_layer(nav_gpkg, inc_pois_layer)) {
   
   unCon_POIs <- filter(st_drop_geometry(filter(nhd_prep, minNet == 1)), COMID %in% nav_POIs$COMID, AreaSqKM == 0)
   
@@ -331,21 +331,21 @@ if(needs_layer(nav_gpkg, inc_pois)) {
     inc_POIs <- bind_rows(tmp_POIs_fixed, select(poi_fix$new_POIs, -c(oldPOI, to_com)))
     
     # Write out fixes
-    write_sf(new_POIs, nav_gpkg, poi_moved)
+    write_sf(new_POIs, nav_gpkg, poi_moved_layer)
     
     # write out final POIs
-    write_sf(inc_POIs, nav_gpkg, inc_pois)
+    write_sf(inc_POIs, nav_gpkg, inc_pois_layer)
   } else {
     # If no fixes designate as NA
     poi_fix <- NA
     # write out final POIs
-    write_sf(nav_POIs, nav_gpkg, inc_pois)
+    write_sf(nav_POIs, nav_gpkg, inc_pois_layer)
   }
 
 } else {
   # Need all three sets for attribution below
-  inc_POIs <- read_sf(nav_gpkg, inc_pois)
-  new_POIs <- if(layer_exists(nav_gpkg, poi_moved)) read_sf(nav_gpkg, poi_moved) else (NA)
+  inc_POIs <- read_sf(nav_gpkg, inc_pois_layer)
+  new_POIs <- if(layer_exists(nav_gpkg, poi_moved_layer)) read_sf(nav_gpkg, poi_moved_layer) else (NA)
   xWalk <- select(st_drop_geometry(new_POIs), COMID, oldPOI)
   unCon_POIs <- filter(st_drop_geometry(filter(nhd_prep, minNet == 1)), COMID %in% nav_POIs$COMID, AreaSqKM == 0)
 }
@@ -380,7 +380,7 @@ if(needs_layer(nav_gpkg, nsegments_layer)) {
   write_sf(nsegments, nav_gpkg, nsegments_layer)
 } else {
   # Read in NHDPlusV2 flowline simple features and filter by vector processing unit (VPU)
-  final_POIs <- read_sf(nav_gpkg, pois_all)
+  final_POIs <- read_sf(nav_gpkg, pois_all_layer)
   nhd_Final <- read_sf(nav_gpkg, nhd_flowline)
   nsegments <- read_sf(nav_gpkg, nsegments_layer)
 }
diff --git a/workspace/02_NHD_navigate.Rmd b/workspace/02_NHD_navigate.Rmd
index ae3cbe11bd6fb3d4b8c783c64899b1258a4675e0..028f88447e6c4a12207df7cb877c05d22e80cf9e 100644
--- a/workspace/02_NHD_navigate.Rmd
+++ b/workspace/02_NHD_navigate.Rmd
@@ -281,7 +281,7 @@ if(all(is.na(tmp_POIs$Type_WBOut))) {
 
 ```{r Final POIs}
 # Derive final POI set ----------------------
-if(needs_layer(nav_gpkg, pois_all)) {
+if(needs_layer(nav_gpkg, pois_all_layer)) {
 
   unCon_POIs <- filter(st_drop_geometry(filter(nhd, minNet == 1)), COMID %in% tmp_POIs$COMID, AreaSqKM == 0)
   
@@ -297,22 +297,22 @@ if(needs_layer(nav_gpkg, pois_all)) {
     # bind together
     final_POIs <- bind_rows(tmp_POIs_fixed, select(poi_fix$new_POIs, -c(oldPOI, to_com)))
     # Write out fixes
-    write_sf(new_POIs, nav_gpkg, poi_moved)
-    write_sf(xWalk, nav_gpkg, poi_xWalk)
+    write_sf(new_POIs, nav_gpkg, poi_moved_layer)
+    write_sf(xWalk, nav_gpkg, poi_xwalk_layer)
     # write out final POIs
-    write_sf(final_POIs, nav_gpkg, pois_all)
+    write_sf(final_POIs, nav_gpkg, pois_all_layer)
   } else {
     # If no fixes designate as NA
     poi_fix <- NA
     # write out final POIs
-    write_sf(tmp_POIs, nav_gpkg, pois_all)
+    write_sf(tmp_POIs, nav_gpkg, pois_all_layer)
   }
   
 } else {
   # Need all three sets for attribution below
-  final_POIs <- read_sf(nav_gpkg, pois_all)
-  new_POIs <- if(layer_exists(nav_gpkg, poi_moved)) read_sf(nav_gpkg, poi_moved) else (NA)
-  xWalk <- if(layer_exists(nav_gpkg, poi_xWalk)) read_sf(nav_gpkg, poi_xWalk) else (NA)
+  final_POIs <- read_sf(nav_gpkg, pois_all_layer)
+  new_POIs <- if(layer_exists(nav_gpkg, poi_moved_layer)) read_sf(nav_gpkg, poi_moved_layer) else (NA)
+  xWalk <- if(layer_exists(nav_gpkg, poi_xwalk_layer)) read_sf(nav_gpkg, poi_xwalk_layer) else (NA)
   unCon_POIs <- filter(st_drop_geometry(filter(nhd, minNet == 1)), COMID %in% tmp_POIs$COMID, AreaSqKM == 0)
 }
 ```
@@ -345,7 +345,7 @@ if(needs_layer(nav_gpkg, nsegments_layer)) {
   write_sf(nsegments, nav_gpkg, nsegments_layer)
 } else {
   # Read in NHDPlusV2 flowline simple features and filter by vector processing unit (VPU)
-  final_POIs <- read_sf(nav_gpkg, pois_all)
+  final_POIs <- read_sf(nav_gpkg, pois_all_layer)
   nhd_Final <- read_sf(nav_gpkg, nhd_flowline)
   nsegments <- read_sf(nav_gpkg, nsegments_layer)
 }
diff --git a/workspace/05_hyRefactor_flines.Rmd b/workspace/05_hyRefactor_flines.Rmd
index af817b30e9b451c283e8a3a4224d7a9a99204e1c..ab4223db7d17093a7ddfe6e6acf30b72d1aa7459 100644
--- a/workspace/05_hyRefactor_flines.Rmd
+++ b/workspace/05_hyRefactor_flines.Rmd
@@ -93,7 +93,7 @@ POIs <- read_sf(nav_gpkg, final_poi_layer) %>%
 
 # Gages that were collapsed during navigate; might not correspond to their Gage_info COMIDs;
 #   We know DAR for those are acceptable, so will keep those out of the events generation.
-POIs_mv <- read_sf(nav_gpkg, poi_moved) %>% 
+POIs_mv <- read_sf(nav_gpkg, poi_moved_layer) %>% 
   filter(!is.na(Type_Gages))
 
 # Also need to avoid modification to flowlines immediately downstream of POIs
@@ -217,7 +217,7 @@ lookup_table <- tibble::tibble(NHDPlusV2_COMID = unique(as.character(refactor_lo
 
 if(!file.exists(lookup_table_file)) {
   readr::write_csv(lookup_table, lookup_table_file)
-  write_sf(lookup_table, out_gpkg, lookup_VPU)
+  write_sf(lookup_table, out_gpkg, lookup_VPU_layer)
 }
 ```
 
diff --git a/workspace/07-1_hyRefactor_cats.Rmd b/workspace/07-1_hyRefactor_cats.Rmd
index b9709d059dabf21c9945d452545a3e7dd969cc41..86b574f5f463374549397b32fa33ffbf0ba5f480 100644
--- a/workspace/07-1_hyRefactor_cats.Rmd
+++ b/workspace/07-1_hyRefactor_cats.Rmd
@@ -37,8 +37,9 @@ if(needs_layer(out_gpkg, divide_layer)) {
   # refactored has the original COMIDs and other attributes
   refactored <- st_transform(read_sf(out_gpkg, refactored_layer), crs) 
   
+  rpu_string <- paste0("rpu_", rpu_code)
   
-  fdr <- raster::raster(data_paths$fdr[[rpu]])
+  fdr <- raster::raster(data_paths$fdr[[rpu_string]])
   crs <- st_crs(fdr)
 
   cats <- st_transform(read_sf(out_gpkg, nhd_catchment), crs) 
@@ -48,8 +49,8 @@ if(needs_layer(out_gpkg, divide_layer)) {
   # Write out catchments to use in 06_non_dend
   write_sf(cats, out_gpkg, nhd_catchment)
   
-  fdr <- raster::raster(data_paths$fdr[[rpu]])
-  fac <- raster::raster(data_paths$fac[[rpu]])
+  fdr <- raster::raster(data_paths$fdr[[rpu_string]])
+  fac <- raster::raster(data_paths$fac[[rpu_string]])
   crs <- raster::crs(fac)
 
   divides <- reconcile_catchment_divides(catchment = cats,
diff --git a/workspace/07-2_aggregate_cats.Rmd b/workspace/07-2_aggregate_cats.Rmd
index 6e0a46bd2217f5bf3ae8e60c79a829a536b5214c..33b1f25ca2d725b6b8df8cb9efae68ac12ac17fe 100644
--- a/workspace/07-2_aggregate_cats.Rmd
+++ b/workspace/07-2_aggregate_cats.Rmd
@@ -274,5 +274,5 @@ lookup_table <- tibble::tibble(NHDPlusV2_COMID = unique(as.character(refactor_lo
 
 
 readr::write_csv(lookup_table, lookup_table_file)
-write_sf(lookup_table, out_gpkg, lookup_VPU)
+write_sf(lookup_table, out_gpkg, lookup_VPU_layer)
 ```
diff --git a/workspace/08_hyRefactor_HI.Rmd b/workspace/08_hyRefactor_HI.Rmd
index ce8f4b01b63be3d906804926c78c4107590181dd..6626bed74d6804e2b43e56fb37af3dea1a5dce5f 100644
--- a/workspace/08_hyRefactor_HI.Rmd
+++ b/workspace/08_hyRefactor_HI.Rmd
@@ -109,7 +109,7 @@ if(needs_layer(out_gpkg, outlets_layer)) {
   write_sf(outlets, out_gpkg, outlets_layer)
 }
 
-events <- read_sf(nav_gpkg, gage_events) %>%
+events <- read_sf(nav_gpkg, gage_events_layer) %>%
   inner_join(select(st_drop_geometry(nhd), AreaSqKM, LENGTHKM, COMID, FromMeas, ToMeas), by = "COMID") %>%
   filter(REACH_meas - FromMeas > 5 & AreaSqKM > 1 & ToMeas - REACH_meas > 5) %>% # & LENGTHKM > (combine_meters / 500)) %>%
   select(Station, COMID, REACHCODE, REACH_meas) %>%
@@ -174,12 +174,15 @@ lookup_table <- tibble::tibble(NHDPlusV2_COMID = unique(as.character(refactor_lo
 
 if(!file.exists(lookup_table_file)) {
   readr::write_csv(lookup_table, lookup_table_file)
-  write_sf(lookup_table, out_gpkg, lookup_VPU)
+  write_sf(lookup_table, out_gpkg, lookup_VPU_layer)
 }
 ```
 
 ```{r reconcile_catchments}
-fdr <- raster::raster(data_paths$fdr[[rpu]])
+
+rpu_string <- paste0("rpu_", rpu_code)
+
+fdr <- raster::raster(data_paths$fdr[[rpu_string]])
 crs <- st_crs(fdr)
 # reconciled have no ID applied in refactoring flowlines
 reconciled <- st_transform(read_sf(out_gpkg, reconciled_layer), crs) 
@@ -194,7 +197,7 @@ cats <- sf::st_make_valid(cats)
 
 if(needs_layer(out_gpkg, divide_layer)) {
   
-  fac <- raster::raster(data_paths$fac[[rpu]])
+  fac <- raster::raster(data_paths$fac[[rpu_string]])
   crs <- raster::crs(fac)
 
   divides <- reconcile_catchment_divides(catchment = cats,
@@ -320,6 +323,6 @@ lookup_table <- tibble::tibble(NHDPlusV2_COMID = unique(as.character(refactor_lo
 
 
 readr::write_csv(lookup_table, lookup_table_file)
-write_sf(lookup_table, out_gpkg, lookup_VPU)
+write_sf(lookup_table, out_gpkg, lookup_VPU_layer)
 
 ```
\ No newline at end of file
diff --git a/workspace/R/NHD_navigate.R b/workspace/R/NHD_navigate.R
index 5217376f4ef7160b363b9a69d04a88af5da6cb50..c257a7a78745a496b4dc9729d852bdfb5aae9b5b 100644
--- a/workspace/R/NHD_navigate.R
+++ b/workspace/R/NHD_navigate.R
@@ -527,8 +527,8 @@ POI_move_down<-function(out_gpkg, POIs, Seg, Seg2, exp, DA_diff){
     inner_join(select(st_drop_geometry(merged_POIs), COMID, oldPOI)) %>%
     select(COMID, oldPOI, Type_HUC12, Type_Gages, Type_TE, Type_NID, Type_WBIn, Type_WBOut, Type_Conf) %>%
     mutate(to_com = COMID) %>%
-    st_compatibalize(., read_sf(out_gpkg, poi_moved))
-  st_write(changed_POIs, out_gpkg, poi_moved, append = TRUE) # write_sf not appending?
+    st_compatibalize(., read_sf(out_gpkg, poi_moved_layer))
+  st_write(changed_POIs, out_gpkg, poi_moved_layer, append = TRUE) # write_sf not appending?
   
   # Format POI moves to table
   xWalk <- select(st_drop_geometry(segsub_dir), Move, COMID) %>%
diff --git a/workspace/R/config.R b/workspace/R/config.R
index bd670b10907adcadda9bee118ca5538d9643f896..7d3090ab28a641ac12b9107b13b2dcbdbc835560 100644
--- a/workspace/R/config.R
+++ b/workspace/R/config.R
@@ -46,131 +46,74 @@ options(scipen = 9999)
 # Source Data
 crs <- 5070
 data_paths <- jsonlite::read_json(file.path("cache", "data_paths.json"))
-huc12_xWalk <- file.path(data_paths$nhdplus_dir, "CrosswalkTable_NHDplus_HU12.csv")
-HUC12 <- file.path(data_paths$nhdplus_dir, "HUC12.rds")
 
-# Defined at various points
+# not used?
+# huc12_xWalk <- file.path(data_paths$nhdplus_dir, "CrosswalkTable_NHDplus_HU12.csv")
+# HUC12 <- file.path(data_paths$nhdplus_dir, "HUC12.rds")
+
+nav_gpkg <- file.path("cache", paste0("reference_", VPU,".gpkg"))
+
+# Defined at various points for HI workflow
 raw_wbd <- paste0("HU12_", VPU) # Raw HUC12
 gage_loc <- paste0("GageLoc_", VPU) # HI GageLoc
+
+### Layer Names used broadly ###
 TE <- paste0("TE_", VPU) # Thermo-electric 
-NID <- paste0("NID_", VPU) # Raw NID points
-gage_events <- paste0("gage_events_", VPU) # gages with flowline index information for possible splitting
+NID_layer <- paste0("NID_", VPU) # Raw NID points
+
+gage_events_layer <- paste0("gage_events_", VPU) # gages with flowline index information for possible splitting
 
-# Defined during NHD_Navigate.Rmd
-nav_gpkg <- file.path("cache", paste0("reference_", VPU,".gpkg"))
 nhd_flowline <- "nhd_flowline"
 nhd_catchment <- "nhd_catchment"
+
 xwalk_layer <- paste0("HUC12_nhd_", VPU) # HUC12 - nhdcat crosswalk, built in Nav for VPU 20
 nav_poi_layer <- paste0("pois_", VPU) # Rolling Nav POI layer added to/modified througout nav workflow
 WBs_layer <-  paste0("WB_", VPU) # Waterbodies within VPU
-poi_moved <- paste0("POIs_mv_", VPU) # POIs moved from original COMID assignment
-inc_pois <- paste0("inc_POIs_", VPU) # layer of POIs used to derive first cut of segments
+poi_moved_layer <- paste0("POIs_mv_", VPU) # POIs moved from original COMID assignment
 nsegments_layer <- paste0("nsegment_", VPU) # Minimally-sufficient network dissolved by POI_ID
-pois_all <- paste0("POIs_", VPU) # All POIs binded together
-poi_xWalk <- paste0("poi_xWalk_", VPU) # POIs that changed COMIDS during the navigate part of the workflow
+pois_all_layer <- paste0("POIs_", VPU) # All POIs binded together
+poi_xwalk_layer <- paste0("poi_xwalk_layer_", VPU) # POIs that changed COMIDS during the navigate part of the workflow
 final_poi_layer <- paste0("final_POIS_", VPU)
 
-# Defined during Refactor
-rpu <- paste0("rpu_", rpu_code)
+# only used in HI?
+inc_pois_layer <- paste0("inc_POIs_", VPU) # layer of POIs used to derive first cut of segments
+
 split_meters <- 10000
 combine_meters <- 1000
 min_da_km <- 20
+
+# parallel factor for catchment reconciliation.
+para_reconcile <- 4
+para_split_flines <- 4
+
 nhd_outlet <- "nhd_outlet"
 nhd_catchment <- "nhd_catchment"
+
 refactored_layer <- "collapsed"
 reconciled_layer <- "reconciled"
-gf_outlets_layer <- "gf_outlets"
 split_layer <- "split_events"
 divide_layer <- "divides"
-refactor_table <- "refactor_lookup"
 agg_fline_layer <- "agg_fline"
 agg_cats_layer <- "agg_cats"
 outlets_layer <- "outlets"
 mapped_outlets_layer <- "mapped_outlets"
-lookup_VPU <- paste0(VPU, "_lookup")
+lookup_VPU_layer <- paste0(VPU, "_lookup")
 
+# Not used?
 # Making an "events" script to address "events" information
 #     i.e. gages, NID, TE Plants, waterbodies
-gage_gpkg <- "cache/Gage_info.gpkg"
-gage_cand <- "Potential_Gages"
-VPU_gage_Table <- file.path("cache", paste0("R", VPU, "_Gages.csv"))
-CONUS_gage_Table <- "data/gages_MDA.rds"
+# gage_gpkg <- "cache/Gage_info.gpkg"
+# gage_cand <- "Potential_Gages"
+# VPU_gage_Table <- file.path("cache", paste0("R", VPU, "_Gages.csv"))
+# CONUS_gage_Table <- "data/gages_MDA.rds"
 
 # Defined during NonDend.Rmd
 ND_gpkg <- file.path("cache", paste0("ND_", VPU,".gpkg"))
-lookup_miss <- paste0("lookup_missing_", VPU)
-divides_xwalk <- paste0("divides_nhd_", VPU)
-missing_cats <- paste0("miss_cats_", VPU)
-missing_terms <- paste0("miss_terms_", VPU)
-HRU_layer <- paste0("poi_cats_", VPU)
-cat_boundary <- paste0("boundaries_", VPU)
 
-# parallel factor for catchment reconciliation.
-para_reconcile <- 4
-para_split_flines <- 4
-
-data_paths <- jsonlite::read_json(file.path("cache", "data_paths.json"))
-nhdplus_path(data_paths$nhdplus_gdb)
-staged_nhd <- stage_national_data(nhdplus_data = data_paths$nhdplus_gdb,
-                                  output_path = data_paths$nhdplus_dir)
-
-# Manually verified
-crs_out <- st_crs('PROJCS["NAD83 / Conus Albers",GEOGCS["NAD83",DATUM["North American Datum 1983",SPHEROID["GRS 1980",6378137.0,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[1.0,1.0,-1.0,0.0,0.0,0.0,0.0],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519943295],AXIS["Geodetic longitude",EAST],AXIS["Geodetic latitude",NORTH],AUTHORITY["EPSG","4269"]],PROJECTION["Albers Equal Area",AUTHORITY["EPSG","9822"]],PARAMETER["central_meridian",-96.0],PARAMETER["latitude_of_origin",23.0],PARAMETER["standard_parallel_1",29.5],PARAMETER["false_easting",0.0],PARAMETER["false_northing",0.0],PARAMETER["standard_parallel_2",45.5],UNIT["m",1.0],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","5070"]]')
-
-# Gages checked after 01_Gage_Selection.Rmd
-gages_add <- c("02492360", "05537000", "05590000", "06244500", "06285400", 
-"06632400", "06659500", "06752000", "08302500", "09089500", "09235600", 
-"09279150", "09307500", "09343300", "09345200", "09537200", "10145400", 
-"10173450", "10194200", "10270830", "10293500", "11055501", "11055801", 
-"11065801", "11113001", "11118501", "11172945", "11230530", "11231600", 
-"11315030", "11358700", "11376050", "11376120", "11376460", "11403530", 
-"11404300", "11413300", "11427750", "11441900", "13176000", "02153700", 
-"04085108", "04086265", "04235600", "04237575", "05367102", "05382267", 
-"05430095", "05536890", "05575830", "06127060", "06230190", "06301850", 
-"06309460", "6349590", "06618300", "06630300", "06630330", "07083200", 
-"07279800", "07282075", "07282090", "07282097", "08260500", "08330800", 
-"09271550", "09274900", "09297600", "09297700", "09299600", "09363050", 
-"09407200", "10102300", "10157000", "10270885", "10270940", "10270970", 
-"10341950", "10347620", "11045050", "11075610", "11108134", "11201456", 
-"11201850", "11295100", "11351946", "11351948", "11351950", "11362900", 
-"11364200", "11364300", "11372350", "11374305", "11376040", "11376140", 
-"11376150", "11376160", "11376420", "11376440", "11389720", "11389740", 
-"11389780", "11396395", "11474780", "12051900", "12052200", "12099100", 
-"12209490", "12320700", "12447387", "12514095", "13012490", "13084400", 
-"13203510", "13292280", "13293350", "13294600", "13294640", "13303500", 
-"13316530", "13329770", "13342250", "50126150")
-
-comid_changes <- data.frame(gages = c("06076560", "09077200", "09520700", "09527900", "11023450"),
-                                    new_comid = c(12420373, 1327427, 20379283, 945030101, 20331218))
-
-
-gages_remove <- c("01484668", "255656081285000", "262007080321500", "02291715", 
-"02287400", "255754081314000", "02291717", "261017080565100", "255601081265300", 
-"02259627", "255738081300100", "02237733", "02274505", "02291670", "261004081082000", 
-"260823080524100", "02293263", "255843081333000", "02298928", "02294747", "02310286", 
-"02244690", "02321900", "02244651", "02246522", "02244601", "02240988", "02326280", 
-"02240982", "022409424", "02326993", "04087101", "04087113", "04085492", "04253294", 
-"04253296", "03352988", "03075650", "07010035", "07010034", "07010094", "07010030", 
-"054279449", "431655089393801", "05403042", "05412056", "073676607", "07380238", 
-"07369654", "05075720", "07032244", "07030295", "05075930", "05075694", "05056636", 
-"06720280", "06720255", "06720285", "05019000", "08075605", "08330940", "08329918", 
-"08329838", "08331105", "08061551", "08329939", "08480594", "08331130", "08330200", 
-"08055700", "08226600", "08329930", "08331140", "09319001", "090955285", "09095526", 
-"09152650", "09095528", "09106104", "09077300", "09095529", "09152600", "09429070", 
-"09522680", "09424150", "09512165", "09534550", "09479502", "09522660", "09423553", 
-"09531850", "09483010", "09533000", "09531900", "09512407", "09419745", "10310402", 
-"10293048", "10311250", "10243630", "10311260", "10302145", "10172791", "10312277", 
-"10336039", "10243640", "10312270", "10172220", "10141400", "10172650", "411403112200801", 
-"14066500", "14349500", "10127109", "13091733", "13157325", "13154510", "13152940", 
-"13155074", "13090999", "13093500", "13082060", "1315377299", "13091000", "13134556", 
-"13134520", "13157293", "131610556", "13095000", "13090998", "13069532", "13170350", 
-"13203510", "13172840", "12509612", "12473502", "12473880", "12509614", "12473560", 
-"12435840", "12471724", "12466101", "12473900", "12472950", "12472300", "12472380", 
-"12200730", "12212430", "11262890", "11060901", "10260776", "10260780", "11262895", 
-"10259920", "11108135", "11062411", "11097260", "10264555", "11062200", "10260778", 
-"11063682", "11073200", "11070263", "11208818", "11062401", "11236080", "11109396", 
-"10264675", "11208800", "11208565", "380626121383803", "11465690", "380629121385203", 
-"11452800", "11181330", "11465700", "380629121385208", "380626121383808", "11181335", 
-"11439950", "11452900", "11439300", "11452901", "11404901", "02084557", "02267000", 
-"04216200", "05018500", "08073000", "08061550", "08017400", "09512100", "09512200", 
-"09406300", "10143500", "10172200", "10142000", "11162800")
\ No newline at end of file
+# not used?
+# lookup_miss <- paste0("lookup_missing_", VPU)
+# divides_xwalk <- paste0("divides_nhd_", VPU)
+# missing_cats <- paste0("miss_cats_", VPU)
+# missing_terms <- paste0("miss_terms_", VPU)
+# HRU_layer <- paste0("poi_cats_", VPU)
+# cat_boundary <- paste0("boundaries_", VPU)