SoilWaterStorage: Soil Moisture Mapping

In ecor/geotopsim: Analysis of GEOtop Output Data

Description

Vertical or Slope-Normal Integration/Aggregation of soil water content

Usage

SoilWaterStorage(theta, psi, layers = 1, fun = sum, na.rm = TRUE,
  psi_thres = 0, comparison = c("<", ">", "<=", ">=", "==", "!="),
  indices = 1, bedrock.depth = NULL, ...)

Arguments

theta	soil water content RasterBrick-class object
psi	soil pressure head content RasterBrick-class object
layers	layer thickness. It is a scalar or a vector whose length is equal to the number of theta layers.
fun	function aggregation
na.rm	a logical value for fun indicating whether NA values should be stripped before the computation proceeds.
psi_thres	threshold value. It is 0 for separating satureted and unsatureted zones, or NA if the two zones are not distinguished.
comparison	comparison operator for psi .See Comparison. Default is the first element of c("<",">","<=",">=","==","!=")
indices	see stackApply
bedrock.depth	map of deprock depth. If it is not NULL, theta cells under bedrock are not considered.
...	further aguments for stackApply

Examples

data(PsiTheta)

 t <- 5
 dz <- 50 ## 50 millimiters layer depth
 unsatWaterVolume <- SoilWaterStorage(theta[[t]],psi[[t]],layer=dz,comparison="<",psi_thres=0)
 satWaterVolume <- SoilWaterStorage(theta[[t]],psi[[t]],layer=dz,comparison=">=",psi_thres=0)

ecor/geotopsim documentation built on May 15, 2019, 10:05 p.m.