spsh (version 1.1.0)

numMualem: Function to Numerically Compute the Mualem Integral

Description

This function will calculate numerically Mualems Integral Mualem.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.

Value

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

Details

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

References

Mualem.1976spsh

Examples

# NOT RUN {
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)
# }