sncFun.01110: Unimodal van Genuchten Non-Capillary Saturation Model

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

Description

Analytical implementation of the non-capillary saturation function \insertCitevanGenuchten.1980spsh.

Usage

sncFun.01110(p_snc, h)

Arguments

p_snc	vector of the 2 van Genuchten Mualem model and h0, the order is sensitve and has to be given as: alf1 van Genuchten alpha [cm-3] n1 van Genuchten n [-] h0 pressure head representing oven dryness given in pF, i.e. log[10](|pressure head| [cm])
h	pressure heads [cm] for which the corresponding retention and conductivity values are calculated.

Details

The function is Eq. Table 1-C1 in insertRefStreck.2020spsh using eq 21 which can be used to accelerate the convergence of the sum. The analytical solution presented in sncFun.01110 only requires the Brooks-Corey model parameters

Value

returns a list with calculations at specified h:

snc	non-capillary saturation

References

\insertRef

vanGenuchten.1980spsh \insertRefWeber.2019spsh \insertRefStreck.2020spsh

Examples

p      <- c(0.1, 0.4, .01, 2, 100, .5)
# add h0
p_snc  <- c(p[3:4], 6.8)
h      <- 10^seq(-2, 6.8, length = 197)
Se     <- shypFun.01110(p, h)$Se
snc    <- sncFun.01110(p_snc, h)

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