aggregate to NHD catchments

[jsadler2]

For the PGDL water temperature sub-quest to account for GW influence, the gridmet data will need to be aggregated to the nhd network. The DO project might also use this more fine-scale data.

I think the only things we'll need to change are:

1. make available a new geopackage with the nhd catchment polygons for the DRB. I believe it makes most sense for this to come from "https://github.com/USGS-R/drb-network-prep/"
2. Revise the Snakemake code to accept the layer name as variable instead of hard coded

   drb_gridmet_tools/Snakefile

   Line 34 in d2a4daf

   gdf = gpd.read_file(input[0], layer="GFv1_catchments_edited")

3. Probably change the "prefix" value to say something about "nhd" so we can distinguish the output

   drb_gridmet_tools/Snakefile

   Line 73 in d2a4daf

   g2s_file_prefix='drb_',

4. Change the value on this line to the nhd catchment geopackage

   drb_gridmet_tools/Snakefile

   Line 9 in d2a4daf

   drb_hru_file = "data/GFv1_catchments_edited.gpkg"

cc @msleckman @lekoenig