Reach Mode

[dvalters]

Implements reach mode in HAIL-CAESAR

See the issue at #14

Summary

Reach model requires a DEM in the same format as the catchment model/mode. Water will exit from any edge of the DEM, but as per catchment mode, not be routed in NODATA cells.

Reach mode requires an additional file that contains the water and sediment inputs into the catchment. The format of the file is a standard plain text/ascii file.

[dvalters]

Closes #14

[dvalters]

Reach mode merge

• No sediment mode has been tested yet but hydro model only mode was working (serial only) with the idealised plane model and the Caersws test case.
• Simulation quite slow?!
• The addGS call in the first bit of sediment code caused segfaults, but this can be fixed at later date.
• Water depth DEMs and Output discarge files looking reasonable.
• Doesn't quite match outputs from the CL run with the idealised plane slope test case, results diverge quite a bit over time but then this is to be expected as discussed in other issues and via email exchange with Tom C. et al.

To do later

• Add in OpenMP parallelisation for quick speed up
• Test sediment routines
• (Separate sediment routing from hydro routine?)
• Further testing with Maggie et al.