Model: prepare model configuration

Source: R/model.R

Model: prepare model configuration

modelConfiguration(modelStructure, pondTemp = 20, gwTempIni = 12,
  infRate = 0.03, depthToWaterTable = 6, hydraulicGradient = 0.001,
  bnd = list(tmp = gwTempIni, ntx = 4, ntc = 1),
  hk = kwb.vs2dh::vs2dh.ConfigureGenuchten(ratioKzKh = 1, ss = 0, satKh =
  750, porosity = 0.2, alpha = 2.3, rmc = 0, beta = 5.8),
  ht = kwb.vs2dh::vs2dh.ConfigureTrans(), iniOutputTime = 1/(3600 *
  24), minSimTime = 0.5, maxSimTime = 31, outputTimeStep = 1,
  solver = kwb.vs2dh::vs2dh.ConfigureBasicSolver(),
  rSolver = kwb.vs2dh::vs2dh.ConfigureRechargePeriodSolver())

Arguments

modelStructure

as retrieved by convRealCoordinatesToNodes()
pondTemp

constant pond temperature (default: 20)
gwTempIni

initial groundwater temperature (default: 12)
infRate

infiltration rate per unit area (Default: 0.9514151 m/d)
depthToWaterTable

water table below ground level (default: 6 m)
hydraulicGradient

hydraulic gradient between left & right model boundary (Default: 0), if positive flow is from left to right, if negative from right to left
bnd

list of structure list(temp=VALUE, ntx=VALUE, ntx=VALUE) passed to function defineHeadBoundary (i.e. boundary with seepage face), if bnd=NULL left/right boundaries are no-flow boundaries
hk

hydraulic properties of soil as retrieved by kwb.vs2dh::vs2dh.ConfigureGenuchten()
ht

transport properties of soil as retrieved by kwb.vs2dh::vs2dh.ConfigureTrans()
iniOutputTime

automatically output results after 1 second of simulation
minSimTime

minimum simulation time in days (default: 0.5)
maxSimTime

maximum simulation time in days (default: 31)
outputTimeStep

at which timestep are the results printed (default: 1 ), i.e. each day
solver

general solver (Default: kwb.vs2dh::vs2dh.ConfigureBasicSolver())
rSolver

recharge period solver (Default: kwb.vs2dh::vs2dh.ConfigureRechargePeriodSolver())

Value

SVH model configuration

See also

defineHeadBoundary for valid additional arguments

Examples

if (FALSE) {
### Importing GIS features
shp.dir <- system.file("extdata", "qgis", package="kwb.demeau")
shp.files <- dir(path = shp.dir, pattern = ".shp", full.names = TRUE)
gisData <- importShapefiles(shp.files)

### Optionally remove some features 
gisData <- removeFeatures(gisData = gisData, ignoreFeatureIDs = c(3,20))
modelStructure <- convRealCoordinatesToNodes(gisData = gisData)

### Model config
conf <- modelConfiguration(modelStructure = modelStructure)

### Running the configuration in VS2DH
res <- kwb.vs2dh::vs2di.runConfig(conf = conf, openTargetDir = TRUE)

### Plotting results
kwb.vs2dh::vs2dh.plotObservationPoints(
  paras = "TEMP", paraLabel = "Temperature", data = res$obsPoints
)
kwb.vs2dh::vs2dh.plotVariables(para = "Temp", data = res$variables)
}