modelConfiguration: Model: prepare model configuration

In KWB-R/kwb.demeau: Heat tracer study SVH

Description

Model: prepare model configuration

Usage

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

## Not run: 
### 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)

## End(Not run)

KWB-R/kwb.demeau documentation built on Sept. 10, 2019, 12:19 p.m.