Infiltration routine used in TOPMODEL, based on the Green-Ampt model.
infiltration(rain, parameters)The function returns a vector with the same length as the input vector rain representing infiltration.
A vector of rain data (m per timestep)
A vector containing 3 parameters (see below for the exact structure)
Wouter Buytaert, Imperial College London, based on original FORTRAN code from Lancaster University
This function gives direct access to the infiltration routine implemented in topmodel().
The function needs three parameters: c(dt, CD, K0, m), where:
| dt | The timestep (hours) |
| CD | capillary drive, see Morel-Seytoux and Khanji (1974) |
| K0 | Surface hydraulic conductivity (m/h) |
| m | Model parameter controlling the rate of decline of transmissivity in the soil profile, see Beven, 1984 |
Morel-Seytoux, H.J., Khanji, J., 1974. Derivation of an Equation of Infiltration. Water Resources Research, 10, 795-800. Beven, K., 1984. Infiltration into a Class of Vertically Non-Uniform Soils. Hydrological Sciences Journal 29, 425-434.
See also https://github.com/ICHydro/topmodel for examples.
topmodel