spsh (version 1.1.0)

sncFun: Non-capillary Saturation Function to Extend Other Functions

Description

The general purpose method to calculate numerically the effective non-capillary saturation is directly obtained from any arbritrary expression for the rescaled capillary saturation function as described by Weber.2019spsh. Examples of capillary saturation functions are the well known vanGenuchten.1980spsh, Kosugi.1996spsh, Fredlund.1994spsh functions.

Usage

sncFun(h, scap)

Arguments

h

A vector of n pressure head values for which scap was calculated

scap

vector of n monotonically decreasing capillary saturation function values calculated by shypFun, rescaled between 0 and 1.

Value

A vector of n elements with calculated saturation content of the non-capillary part.

References

vanGenuchten.1980spsh Kosugi.1996spsh Fredlund.1994spsh

Examples

# NOT RUN {
# set variables
p <- c(0.1, 0.4, 0.01, 2, 100, .5)
h <- 10^seq(-2, 6.8, length = 197)

# Calculate the capillary and non-capillary saturation function.
Se <- shypFun(p, h, shpmodel = "01110")$Se     
Snc <- sncFun(Se)

# }