shypFun.04110 {spsh}

Unimodal Brooks-Corey Model

Description

Calculates the soil hydraulic property function values based on given pressure heads (Kosugi 1996).

Usage

shypFun.04110(p, h)

Arguments

p

vector of the 6 Brooks-Corey model parameters, order is sensitve and has to be given as:

thr	residual water water content [cm cm-3]
ths	saturated water water content [cm cm-3]
alf	air entry pressure head [cm^-1]
bet	effective model parameter [ - ]
Ks	saturated conductivity [cm d-1]
tau	exponent of Se in the capillary conductivity model, sometimes denoted in the literature as l [-]

h

pressure heads [cm] for which the corresponding retention and conductivity values are calculated.

Details

The function solves analytically the spec. water capacity function and integral to the capillary bundle model.

Value

returns a list with calculations at specified h:

theta

calculated volumetric moisture content

Se

calculated saturation

cap

returns NA; not supported

psd

returns NA; not supported

Kh

Hydraulic conductivity values

Note

The Muealm integral is solved numerically

References

Brooks RH, Corey AT (1964). Hydraulic Properties of Porous Media. Hydrology Papers No 3, Colorado State University, Fort Collings, Colorado.

Examples

p      <- c(0.1, 0.4, .01, .3, 100, .5)
h      <- 10^seq(-2, 6.8, length = 197)
shyp.L <- shypFun.04110(p, h)