shypFun.01310: van Genuchten-Mualem trimodal Soil Hydraulic Propterty Model

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/shypFun.01310.R

Description

trimodal van Genuchten-Mualem functions for the retention and hydraulic conductivity curves \insertCiteDurner.1994spsh.

Usage

1

Arguments

p

vector of the 9 van Genuchten-Mualem model parameters, order is sensitve and has to be given as:

thrresidual water water content [cm cm-3]
thssaturated water water content [cm cm-3]
alf1van Genuchten alpha [cm-3]
n1van Genuchten n [-]
w1fraction of the first modality [-], w2 is internally computed as w2 = 1-w1
alf2van Genuchten alpha of the second modality[cm-3]
n2van Genuchten n of the second modality [-]
w2fraction of the second modality [-], w3 is internally computed as w3 = 1-w1-w2, in resFun ensures w3 >=0
alf3van Genuchten alpha of the third modality[cm-3]
n3van Genuchten n of the third modality [-]
Kssaturated conductivity [cm d-1]
tauexponent 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

specific water capacity function

psd

pore size distribution

Kh

Hydraulic conductivity values

Author(s)

Tobias KD Weber , tobias.weber@uni-hohenheim.de

References

\insertRef

Durner.1994spsh

See Also

\insertRef

Weber.2019spsh
\insertRefWeber.2017aspsh
\insertRefWeber.2017bspsh

Examples

1
2
3
4
5
6
p <- c("thr" = 0.1, "ths" = 0.4, alf1 = .5, "n1" = 3,
       "w1" = .5, "alf2" = 0.01, "n2" = 2, 
       "w2" = .3, "alf3" = 0.01, "n3" = 1.6, 
       "Ks" = 100, "tau" = .5)
h <- 10^seq(-2, 6.8, length = 197)
shyp.L <- shypFun.01310(p, h)

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