unit.hydrograph.storage.cascade: Unit hydrograph linear storage cascade
In RHydro: Classes and methods for hydrological modelling and analysis

Description Usage Arguments Details Value Author(s) Examples

Compute the shape of a Unit hydrograph linear storage cascade with the given parameters and apply it to given rainfall series.

1 2	unit.hydrograph.storage.cascade(n, k, delta_t) apply_uhg(rainfall, uhg_params, delta_t)

`n`	number of linear storages [-]
`k`	recession constant (same unit as `delta_t`)
`delta_t`	time resolution of return vector [time].
`rainfall`	vector of rainfall series in time resolution of `delta_t`.
`uhg_params`	vector with named components "n" and "k" (parameters of UHG) OR matrix with named columns "n", "k" and "a". In the latter case, a multiple UHG is applied with the weights in `"a"`.

The resulting vector h would, in principle, be of infinite length because of the asymptotic behaviour of the storage outflow. For practical reasons, it is trimmed so that at least 99% of the volume is reproduced.

vector h with ordinates of unit hydrograph or modelled hydrograph, corresponding to temporal resolution of delta_t starting with t=0

Till Francke

#simple unit hydrograph ordinates
delta_t = 0.5
uh = unit.hydrograph.storage.cascade(n=2, k=3, delta_t = delta_t)
plot(seq(from=delta_t, by=delta_t, along.with=uh) , uh, type="l")

#compute hydrograph from multiple UHGs:
rainfall <- fuse.DATA[,"P"]
uhg_params=rbind( c(n=2,  k=3, a=0.4),
                  c(n=20, k=3, a=0.6))	#two-element unit hydrograph, weighted 40 and 60 %, resp.
plot(apply_uhg(rainfall, uhg_params, delta_t), type="l")