Provides outputs of the ReFH model from catchment descriptors or user defined inputs
ReFH(
CDs = NULL,
Depth = NULL,
duration = NULL,
timestep = NULL,
scaled = NULL,
PlotTitle = NULL,
RP = NULL,
alpha = TRUE,
season = NULL,
AREA = NULL,
TP = NULL,
BR = NULL,
BL = NULL,
Cmax = NULL,
Cini = NULL,
BFini = NULL,
Rain = NULL
)catchment descriptors derived from either GetCDs or ImportCD
a numeric value. The depth of rainfall used as input in the estimation of a design hydrograph. The default is a two year rainfall
a numeric value. A duration for the design rainfall
a numeric value. A user defined data interval. The default changes depending on the estimated time to peak to formulate a sensible looking result
a numeric value of peak flow in m3/s
a character string. A user defined title for the ReFH plot
return period (Default = 2)
a logical argument with default TRUE. If TRUE the alpha adjustment is applied based on RP. If FALSE, no alpha adjustment is made
a choice of "summer" or "winter". The default is "summer" in urban catchments (URBEXT2000 > 0.03) and "winter" in rural catchments
numeric. Catchment area in km2.
numeric. Time to peak parameter (hours)
numeric. Baseflow recharge parameter
numeric. Baseflow lag parameter (hours)
numeric. Maximum soil moisture capacity parameter (mm)
numeric. Initial soil moisture content (mm)
numeric. Initial baseflow (m3/s)
numeric. User input rainfall. A numeric vector
A print out of parameters, a results data.frame, and a plot. First is a print of the parameters, initial conditions and the catchment area. The second is a data.frame with columns Rain, NetRain, Runoff, Baseflow, and TotalFlow. If the scale argument is used a numeric vector containing the scaled hydrograph is returned. The plot is of the ReFH output, with rainfall, net-rainfall, baseflow, runoff and total flow. If the scaled argument is used, a scaled hydrograph is plotted.
The ReFH is described in the Flood Estimation Handbook Supplementary Report No.1 (2007). The method to derive design rainfall profiles is described in the Flood Estimation Handbook (1999), volume 2. Users can also input their own rainfall with the 'Rain' argument. As a default, when catchment descriptors (CDs) are provided the ReFH function uses catchment descriptors to estimate the parameters of the ReFH model and the two year rainfall for the critical duration. The latter is based on a quadratic interpolation of the catchment descriptors RMED1H, RMED1D, and RMED2D. Parameters and initial conditions can also be inidvidually input by the user. If a parameter argument is used for one or more of the parameters, then these overwrite the CD derived parameters. If a value for the scaled argument is provided (m3/s), a scaled hydrograph is returned. The formulation of this function avoids the need for users to select a storm duration that is an odd integer multiple of the data interval. This is achieved as follows: for the design rainfall profile, the duration is split into a sequence by multiplying duration by the reciprocal of the data interval (timestep) and rounding up to the nearest odd integer. This forms the end of a sequence starting from one and progressing at an interval of two. This sequence is divided by the duration to form the vector to which the exponent b is raised in the calculation of z for the storm profile (see Flood Estimation Handbook (1999), volume 2)
# NOT RUN {
#Get CDs and apply the ReFH function
CDs.203018 <- GetCDs(203018)
ReFH(CDs.203018)
#Apply the ReFH function, scale to a 100-year flow estimate and change the plot title accordingly
ReFH(CDs.203018, scaled = 182, PlotTitle = "100-Year Design Hydrograph - Site 203018")
#Apply the ReFH function with a user defined initial baseflow
ReFH(CDs.203018, BFini = 6)
# }
Run the code above in your browser using DataLab