kwb.hantush (version 0.3.0)

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

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

# NOT RUN {
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
)
# }