Learn R Programming

spsh (version 1.1.0)

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

Description

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

Usage

shypFun.01310(p, h)

Arguments

p

vector of the 9 van Genuchten-Mualem 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] alf1 van Genuchten alpha [cm-3]
n1 van Genuchten n [-] w1
fraction of the first modality [-], w2 is internally computed as w2 = 1-w1 alf2 van Genuchten alpha of the second modality[cm-3]
n2 van Genuchten n of the second modality [-] w2
fraction of the second modality [-], w3 is internally computed as w3 = 1-w1-w2, in resFun ensures w3 >=0 alf3 van Genuchten alpha of the third modality[cm-3]
n3 van Genuchten n of the third modality [-] Ks
saturated conductivity [cm d-1] thr residual water water content [cm cm-3]
h

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

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

Details

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

References

Durner.1994spsh

See Also

Weber.2019spsh Weber.2017aspsh Weber.2017bspsh

Examples

Run this code
# NOT RUN {
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)

# }

Run the code above in your browser using DataLab