p_itca_01_PrepData_TIGER_DEM_NLCD_NHD.py was used to prepare spatial inputs for these test cases, including:
dem.asc
landcover.asc
imperviouscover.asc(optional, only needed when Flag_Buffer_Algorithm == 2016)
nhdfdr.asc (optional, only needed when Flag_FlowDirection_Source == NHD)
nhdwaters10.asc (optional, only needed when Flag_FlowDirection_Source == NHD)

p_itca_01_PrepData_TIGER_DEM_NLCD_NHD.py is available from https://code.itreetools.org/ESF/Hydro/CoolAir_NationalRuns


HydroPlus_Buffer_SSURGO_gen.py was used to prepare SSURGO data, including:
s_depth.asc (optional, only needed when Flag_Buffer_Algorithm == 2016)
s_ksat.asc (optional, only needed when Flag_Buffer_Algorithm == 2016)

HydroPlus_Buffer_SSURGO_gen.py is available from https://code.itreetools.org/ESF/Hydro/Buffer_git

SpatialBuffer Test Cases:
SaddleRiver_D8_2003 Saddle River watershed spanning New Jersey and New York, with HUC8 02030103 
SaddleRiver_D8_2003 uses nhdfdr.asc for flow direction, from the National Hydrography Database.
SaddleRiver_D8_2003 uses nhdwaters10.asc for identifying receiving waters, from the National Hydrography Database.

Note: Only Algorithm_FlowDirection = D8 is able to use the nhdfdr.asc data, when the Flag_FlowDirection_Source set as NHD. 
Note: The Algorithm_FlowDirection = MFD and DInf do not follow a single cell flow direction, and diverge from nhdfdr.asc.
Note: When Algorithm_FlowDirection is set to a non-D8 algorith when Flag_FlowDirection_Source = NHD, the program will force the Algorithm_FlowDirection to D8. 

All Test Cases:
StartDate_YYYYMMDD: All cases are simulating 1 year of loading, but are not using the MMDD values. 


USGS gauged basin ID 01390450 boundary polygon was delineated by USGS, evidently to capture stormwater drainage.
A map analysis of this USGS boundary polygon reveals it does not agree with the NHD boundary for this river system.
NHD river basin boundaries are shown in the Hydrodem files, and are used to generate the flow direction product.
Flow accumulation values generated using the NHD flow direction product may be erroneous, and not capture artificial drainage.

Note that NHD Hydrodem boundaries are aligned with the HUC 12 boundaries. 

Basin boundaries download:

1. HUC 12 boundaries can be download from 
option1: https://apps.nationalmap.gov/downloader/#/. 
Select  Hydrography (NHDPlus HR, NHD, WBD), then select Watershed Boundary Dataset (WBD)+ National +  FileGDB, search products, then download. GDB file name: WBD_National_GDB.gdb
option2: https://www.epa.gov/waterdata/nhdplus-national-data
Select NHDPlusV21_NationalData_WBDSnapshot_FileGDB_08.7z. GDB file name: NationalWBDSnapshot.gdb

2. USGS Streamgage boundaries can be download from 
option 1:  The HydroPlus team maintains a dataset located within the Buffer external drive under the path Constants > BasinGrid > USBasinsShape > USBasins.shp. 
To obtain this file, kindly inquire with the i-Tree team for assistance.
option 2: You can also access the USGS Streamgage boundaries dataset by visiting the following link https://www.sciencebase.gov/catalog/item/631405ccd34e36012efa32dc
This dataset is comprised of 18 zipped files, each containing boundary data for a total of 19,031 basins. This dataset was released in 2011.

