Description Usage Arguments Details Value Author(s) References Examples
This function will calculate numerically Mualems Integral \insertCiteMualem.1976spsh and return Hydraulic Conductivity Values.
1 |
h |
vector of length |
scap |
Capillary saturation [cm3 cm-3]. |
pcon |
vector of soil hydraulic conductivity model parameters, the first argument is |
The numerical solution of Mualems integral relies on the trapezoidal rul of integration.
returns a vector of length l
of calulcated conductivity values at h
.
Tobias KD Weber , tobias.weber@uni-hohenheim.de
Mualem.1976spsh
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.