numMualem: Function to Numerically Compute the Mualem Integral

In spsh: Estimation and Prediction of Parameters of Various Soil Hydraulic Property Models

Description

This function will calculate numerically Mualems Integral \insertCiteMualem.1976spsh and return Hydraulic Conductivity Values.

Usage

numMualem(h, scap, pcon = NA)

Arguments

h	vector of length l pressure head values.
scap	Capillary saturation [cm3 cm-3].
pcon	vector of soil hydraulic conductivity model parameters, the first argument is q and the second r.

Details

The numerical solution of Mualems integral relies on the trapezoidal rul of integration.

Value

returns a vector of length l of calulcated conductivity values at h.

Author(s)

Tobias KD Weber , [e-mail]

References

\insertRef

Mualem.1976spsh

Examples

h <- 10^seq(-3, 6.8, length = 501)
p = c(.05, .5, .01, 1.8, 100, .5)
shyp.L <- shypFun.01110(p, h)

Ks <- p[5]
tau <- p[6]
Se <- shyp.L[['Se']]
Khrnum <- numMualem(h, pcon = tau, scap = Se) 

Khnum <- Ks * Se^tau * Khrnum

plot(log10(h), log10(shyp.L[['Kh']]), ylim = c(-10, 2.3), 
     xlim = c(-1,6), ylab = "log10 Kunsat [ cm/d ]", xlab = "pF [ - ]", type = "l", lwd = 8)
lines(log10(h), log10(Khnum), col = "red", lwd = 2)

spsh documentation built on April 14, 2020, 6:37 p.m.