Enabling reach routing at headwater HRUs

ESCOMP / mizuRoute

Reach-based river routing model

http://escomp.github.io/mizuRoute/

GNU General Public License v3.0

42 stars 52 forks source link

Enabling reach routing at headwater HRUs #467

Closed nmizukami closed 3 months ago

nmizukami commented 4 months ago

Currently, reach routing in a headwater HRU is not performed because lateral flow (BASIN_QR) is always added at the bottom of reach. Note nothing is coming from the top of the headwater reach!

This PR provides an option to add lateral flow at the top of reach in the headwater, so the headwater reach has volume of water as well as discharge from the reach.

New control variable is hw_drain_point (integer)-- 1-> add top of reach, 2-> add bottom of reach. Default is 2 as currently done.

nmizukami commented 3 months ago

The answers were not expected to change with this pull request when hw_drain_point is set to 2 (i.e., adding lateral flow at the bottom as currently done) for all the routing methods except MC routing. MC routing volume check is changed in this PR, resulting in bit-to-bit change in discharge values (nccmp detect this difference)