ch_rthr Subroutine

subroutine ch_rthr()

Uses

  • proc~~ch_rthr~~UsesGraph proc~ch_rthr ch_rthr module~basin_module basin_module proc~ch_rthr->module~basin_module module~channel_data_module channel_data_module proc~ch_rthr->module~channel_data_module module~channel_module channel_module proc~ch_rthr->module~channel_module module~climate_module climate_module proc~ch_rthr->module~climate_module module~hydrograph_module hydrograph_module proc~ch_rthr->module~hydrograph_module module~sd_channel_module sd_channel_module proc~ch_rthr->module~sd_channel_module module~time_module time_module proc~ch_rthr->module~time_module module~hydrograph_module->module~basin_module module~hydrograph_module->module~time_module

~ ~ ~ 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

Arguments

None

Calls

proc~~ch_rthr~~CallsGraph proc~ch_rthr ch_rthr proc~chrc_interp chrc_interp proc~ch_rthr->proc~chrc_interp