hantushDistancesBaseProps: Hantush distances & base properties: allows input of vector...

In KWB-R/kwb.hantush: Calculation of Groundwater Mounding Beneath an Infiltration Basin

Description

Hantush distances & base properties: allows input of vector of x,y coordinates and also a vector for one of the base properties

Usage

hantushDistancesBaseProps(x = seq(0, 200, 5), y = rep(0, length(x)),
  baseProps = baseProperties(time = 2^(0:6), infiltrationRate = 1,
  basinWidth = 10, basinLength = 50, horizConductivity = 10, iniHead = 10,
  specificYield = 0.2), dbg = FALSE)

Arguments

x	vector with distances from the center of the recharge basin in the x direction (L), (Default: every 5 meter between 0-200m)
y	vector with distances from the center of the recharge basin in the y direction (L), (Default: 0 times length of x)
baseProps	as retrieved by baseProperties(), but one property is allowed to be a vector (e.g. infiltrationRate = c(1,2,4))
dbg	If True additional debug messages are printed on screen

Value

List with sublists "dat" (x,y,head & WLincrease), "changedBaseProp.Name" (name of base property with multiple values) and "baseProps" (complete base parameterisation)

See Also

baseProperties for basic model properties

Examples

baseProps <- baseProperties(
  time = 2^(0:6),
  infiltrationRate = 1,
  basinWidth = 10,
  basinLength = 50,
  horizConductivity = 10,
  iniHead = 10,
  specificYield = 0.2,
  numberTimeSteps = 10
)
res <- hantushDistancesBaseProps(baseProps = baseProps)
cols <- length(unique(res$dat[[res$changedBaseProp.Name]]))
mainTxt <- sprintf("Changed baseProperty: %s", res$changedBaseProp.Name)
xyplot(WLincrease ~ x,
  groups = res$dat[[res$changedBaseProp.Name]],
  data = res$dat,
  type = "b",
  auto.key = list(columns = cols),
  main = mainTxt
)

KWB-R/kwb.hantush documentation built on Oct. 1, 2019, 3:32 p.m.