Added wb events to split events

workspace/05_hyRefactor_flines.Rmd

@@ -72,7 +72,14 @@ events <- read_sf("cache/gages_info.gpkg", "Gages") %>%
   filter(reach_meas - FromMeas > 10 & AreaSqKM > 2 & 
            ToMeas - reach_meas > 10 & LENGTHKM > (combine_meters / 1000)) %>%
   select(COMID = comid, REACHCODE = reachcode, REACH_meas = reach_meas, site_no) %>%
-  filter(!COMID %in% avoid$COMID)
+  filter(!COMID %in% avoid$COMID) %>%
+  mutate(event_type = "streamgage")
+
+if(!needs_layer(nav_gpkg, "WB_events")){
+  wb_events <- read_sf(nav_gpkg, "WB_events") %>%
+    filter(COMID %in% nhd$COMID) 
+  events <- rbind(events, select(wb_events, COMID, REACHCODE, REACH_meas, site_no = WBAREACOMI, event_type))
+}
 
 # Write out events and outelts
 write_sf(events, out_refac_gpkg, split_layer)