SoilWaterStorage: Soil Moisture Mapping
In ecor/AnalyzeMoistureinSuperslab: TO DO

Description Usage Arguments Examples

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

1
2
3

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

`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`
`...`	further aguments for `stackApply`

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)