spsh (version 1.1.0)

sncFun.04110: Unimodal Brooks-Corey Non-Capillary Saturation Model

Description

Analytical implementation of the non-capillary saturation function from the Brooks-Corey function Brooks.1964spsh.

Usage

sncFun.04110(p_snc, h)

Arguments

p_snc

vector of the 2 Brooks-Corey model parameters, order is sensitve and has to be given as:

alf	air entry pressure head [cm^-1]	bet
effective model parameter [ - ]	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-A in insertRefStreck.2020spsh The analytical solution presented in sncFun.04110 only requires the Brooks-Corey model specific parameters and h0.

References

Brooks.1964spsh Streck.2020spsh

Examples

# NOT RUN {
p     <- c(0.1, 0.4, .02, 2, 100, .5)
p_snc <- c(p[3:4], 6.8)
h     <- 10^seq(-2, 6.8, length = 197)
Se    <- shypFun.04110(p, h)$Se
snc   <- sncFun.04110(p_snc, h)
# }