shypFun: Wrapper Function for all Supported Soil Hydraulic Property...

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

Description

This function allows to select soil hydraulic property models.

Usage

shypFun(p, h, shpmodel = "01110", ivap.query = NULL)

Arguments

p	Vector of the model parameters, order is sensitve cf respective model documentation.
h	Pressure heads [cm] for which the corresponding retention and conductivity values are calculated.
shpmodel	character 01110 unimodel van Genuchten-Mualem model, with the contraint of m = 1-1/n \insertCitevanGenuchten.1980spsh 01210 bimodel van Genuchten-Mualem model, with the contraint of m_i = 1-1/n_i \insertCiteDurner.1994spsh 01310 trimodal van Genuchten-Mualem model, with the contraint of m_i = 1-1/n_i \insertCiteDurner.1994spsh 02110 unimodel Kosugi 2 parametric-Mualem model \insertCiteKosugi.1996spsh 03110 unimodel Fredlund-Xing-Mualem model, with the contraint of m = 1-1/n \insertCiteFredlund.1994spsh
ivap.query	NULL no isothermal vapour conductivity will be calculated with Kvap Model type for isothermal vapour conductivity, see Details of function KvapFun for model codes

Value

a list with calculated soil hydraulic properties at specified h:

theta	calculated volumetric moisture content
Se	calculated saturation
Scap	effective saturation (of the capillary part if FM is specified
cap	specific water capacity function (of the capillary part if FM is specified
psd	pore size distribution (of the capillary part if FM is specified
Kh	total hydraulic conductivity

if FM specified, additionally:

thetacap	calculated volumetric moisture content of the capillary part
thetanc	calculated volumetric moisture content of the non-capillary part
Snc	effective saturation of the non-capillary part
Kcap	hydraulic conductivity of the capillary
Knc	hydraulic conductivity of the non-capillary
Kvap	isothermal vapour conductivity
Krcap	relative hydraulic conductivity of the capillary
Krnc	relative hydraulic conductivity of the non-capillary

Note

the function is used to assign a new function variable with a function which calculates the soil hydraulic properties according to sepcified shpmodel and model specified by ivap.query

Author(s)

Tobias KD Weber , [e-mail]

References

\insertRef

Durner.1994spsh
\insertRefFredlund.1994spsh \insertRefKosugi.1996spsh \insertRefvanGenuchten.1980spsh \insertRefWeber.2019spsh

Examples

 data("shpdata1")
 retdata <- shpdata1$LFH1$wrc[!is.na(shpdata1$LFH1$wrc[,1]),]
 condata <- shpdata1$LFH1$hcc
 
 # assign auxiliary variables
 pF <- seq(-3, 6.8, length = 501)
 h <- 10^pF
 
 # assign list of parameters for the van Genuchten-Mualem model
 parL <- list("p" = c("thr"= 0.05, "ths" = 0.45, "alf1" = 0.01, "n" = 2, "Ks" = 100, "tau" = .5),
              "psel" = c(1, 1, 0, 1, 1, 1),
              "plo" = c(0.001 , 0.2, 0.001, 1.1, 1, -2),
              "pup" = c(0.3, 0.95, 1, 10, 1e4, 10))
 
 # calculate soil hydraulic property function values
 
 shyp.L <- shypFun(parL$p, h, shpmodel = "01110", ivap.query = NULL)
 wrc <- shyp.L$theta
 hcc <- log10(shyp.L$Kh)
 
 # # PLOT THE MEASURED WATER RETENTION CURVE
 ticksatmin <- -2
 tcllen <- 0.4
 ticksat <- seq(ticksatmin,5,1)
 pow <- ticksatmin:6
 
 par(mfrow = c(1,2), tcl = tcllen)
 plot(retdata, ylim = c(.0,.50), xlim = c(0, 6.8), ylab = "", xlab = "",
      col = "darkgrey", axes = FALSE, main = "Water Retention Curve", cex = 2)
 lines(log10(h), wrc, col = "darkblue", lwd = 2)
 legend("topright", c("observed", "simulated"),pch = c(1,NA),
        lty = c(NA,1), lwd = 2, bty = "n", cex = 1.3,col = c("darkgrey", "darkblue"))
 axis(1, at = pow, labels = parse(text = paste('10^',(pow), sep = "")), tcl = tcllen)
 axis(2, tcl = tcllen)
 axis(4, labels = NA)
 axis(3, labels = NA)
 mtext("pressure head |h| [cm]", 1, cex = 1.2, line = 2.8)
 mtext("vol. water content [ - ]", 2, cex = 1.2, line = 2.8)
 box()
 
 # PLOT THE MEASURED HYDRAULIC CONDUCTIVITY CURVE
 plot(condata, ylim = c(-8,2), xlim = c(0, 6.8), ylab = "", xlab = "", col = "darkgrey",
      axes = FALSE, main = "Hydraulic Conductivity Curve", cex = 2)
 lines(log10(h), hcc, col = "darkblue", lwd = 2)
 legend("topright", c("observed", "simulated"), pch = c(1,NA),
        lty = c(NA,1), lwd = 2, bty = "n", cex = 1.3, col = c("darkgrey","darkblue"))
 axis(1, at = pow, labels = parse(text = paste('10^',(pow), sep = "")), tcl = tcllen)
 axis(2)
 axis(4, labels = NA)
 axis(3, labels = NA)
 mtext("log10 K [cm/d]", 2, cex = 1.2, line = 2.8)
 mtext("pressure head |h| [cm]",1 , cex = 1.2, line = 2.8)
 box()
 par(mfrow = c(1,1))
 
 ## Not run: 
 # HOW TO WRITE A MATER.IN FOR HYDRUS-1D
 
 mater_out <- cbind(shyp.L[['theta']], h, shyp.L[['Kh']], abs(shyp.L[['cap']]))
 
 materWriteFun <- function(mater_out.L, path = getwd(), sample) {
       
       # Function to write a Mater.in
       
       # ARGUMENTS
       
       # mater_outdata frame of 4 columns of calculated SHP values at h and length n. 
       # 1. Column: THETA
       # 2. Column: H(negative pressure heads)
       # 3. Column: K
       # 4. Column: C(positive)
       # path character specifying the path where the MATER.IN should be saved
       # sample optional chr for sample name: NULL = no name given
       
       n <- dim(mater_out)[1]
       sink(file.path(path, paste(sample, "MATER.IN", sep = "")))
       cat(c("iCap", "\n", "1", "\n", "NTab", "\n", n, "\n"))
       cat(c("\t","THETA", "\t\t\t","H","\t\t\t","K","\t\t\t","C"))
       cat("\n")
       
       write.table(format(mater_out, justify = "right"),
                   row.names = FALSE, col.names = FALSE, quote = FALSE)
       sink()
 }
 
## End(Not run)

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