This is the repository for the paper "Error Assessment for Height Above the Nearest Drainage Inundation Mapping."
Authors: Lukas Godbout, Jeff Y. Zheng, Sayan Dey, Damilola Eyelade, David Maidment, and Paola Passalacqua.
Please direct all communication to LukasGodbout@utexas.edu.

This repository includes four base folders:

1. Scripts
2. Inputs
3. Results

The scripts folder contains the four scripts designed for this study in addition to an initial scripts folder. All scripts are designed to run on Python 2.7.13 (32 bit).
The initial scripts folder contains scripts and a tutorial designed and provided by Dr. Xing Zheng, following the method presented in Zheng, X., D.G. Tarboton, D.R. Maidment, Y.Y. Liu & P. Passalacqua. (2018). River Channel Geometry and Rating Curve Estimation using Height Above the Nearest Drainage. Journal of the American Water Resources Association (JAWRA), in press.
Please refer to the tutorial provided on how to create the input files used for our study.

The four scripts designed for this study are to be run individually for each river. The scripts shown use Guadalupe River as the example. Guadalupe river is in HUC-6 121002.
The Rating_Curve_Comparison script takes the input data and transforms it to allow for comparison between SRCs and HEC-RAS reach averaged rating curves.
The Rating_Curve_Evaluation script calculates statistics to evaluate SRC performance against HEC-RAS reach averaged rating curves.
The HAND_Moving_Window_Analysis script calculates new slope values for short reaches over a specified moving window length.
The Averaged_NRMSE script aggregates the new SRCs to calculate mean normalized RMSE for the set of moving window lengths.

The input folder contains the files output from the initial scripts for each river and information available online.
Each river has a NHD cross section spreadsheet containing the list of cross sections within each reach.
Each river has a COMID Max Station spreadsheet containing the most upstream cross section of each reach.
Each river has a Final Rating Curves spreadsheet containing all the stage height discharge relationships for each reach.
There are two general files for each HUC-6, where the Hydroprop spreadsheet contains reach geometry information and the NHD spreadsheet contains relevant NHDPlus attributes.
For other rivers new required input files can be found categorized by HUC-6 at https://web.corral.tacc.utexas.edu/nfiedata/HAND/.

The results folder contains the main result spreadsheets and the depth difference spreadsheets.
The main result spreadsheets display for each reach the RMSE, normalized RMSE, max error, range, percent bias and R squared.
The depth difference spreadsheets display for each reach the median and mean absolute depth difference, reach length, reach slope, bankfull width and bankfull depth values.
The scaling analysis spreadsheets show the mean normalized RMSE values for a set of moving window lengths.