spsh (version 1.1.0)

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

Description

Analytical implementation of the non-capillary saturation function vanGenuchten.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.

Value

returns a list with calculations at specified h:

snc

non-capillary saturation

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

References

vanGenuchten.1980spsh Weber.2019spsh Streck.2020spsh

Examples

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