~ ~ ~ PURPOSE ~ ~ ~
This subroutine routes flow at any required time step through the reach
using a constant storage coefficient
Routing method: Variable Storage routing
~ ~ ~ INCOMING VARIABLES ~ ~ ~
name |units |definition
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
~ ~ ~ OUTGOING VARIABLES ~ ~ ~
name |units |definition
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
~ ~ ~ LOCAL DEFINITIONS ~ ~ ~
name |units |definition
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
ii |none |counter (hour)
inhyd |none |inflow hydrograph storage location number
jrch |none |reach number
scoef |none |storage coefficient
topw |m |width of channel at water level
vol |m^3 H2O |volume of water in reach
wtrin |m^3 H2O |water entering reach during hour
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
~ ~ ~ SUBROUTINES/FUNCTIONS CALLED ~ ~ ~
Intrinsic: Sum, Min, Sqrt
~ ~ ~ ~ ~ ~ END SPECIFICATIONS ~ ~ ~ ~ ~ ~
subroutine developed by A. Van Griensven,
Hydrology-Vrije Universiteit Brussel, Belgium
Modified by Jeahak Jeong, Blackland Research, Temple, USA
volume at start of day subdaily time step flow rate during time step m3/s = m3 / (24 / (1/day) * 3600) update volume of water in reach - m3 find where flow fits in rating curve (0.1,1.0,1.5 * bankfull flow rate) keep max travel time at 1.5 bankfull flood routing using variable storage coeffiecient interpolated travel time add incoming flow to channel and flood plain storage based on ratio from rating curve - floodplain/channel volumes Variable Storage Coefficient - calculate volume of water leaving reach during timestep subtract routed outflow from from channel storage
calculate flow from flood plain wetland storage below principal spillway is wetland storage - above is flood plain storage
readjust channel and flood plain storage?
calculate transmission losses calculate evaporation calculate width of channel at water level set volume of water in channel at end of hour
rcurv%flo_rate > 0.
end of sub-daily loop
save storage volume for next day and set outflow for day