DNST: DNAPL source term model

Description Usage Arguments Value Examples

View source: R/DNST_solver.R

Description

DNAPL source term model

Usage

1
2
3
4
DNST(result.file, description = "", DNAPLmodel, uhist,
  uhist.J_units = "kg/day", fsphist = data.frame(year = c(1974, 1990), f =
  c(1, 0)), fw = 1, fi = 1, start.t = 9100, end.t, dt = 20, x = NULL,
  y = NULL, z0 = "base", mfdata, qh = NULL, check.qh = TRUE)

Arguments

result.file

character string; file to save results to (should end with extension ".rds")

description

character string; information about the model and what it represents; this is saved with the results for your record

DNAPLmodel

DNAPLmodel S4 object; see DNAPLmodel

uhist

data.frame:
$year (num)
$cons (num): consumption rate
the transient import rate of solvent to the site (usage history)

uhist.J_units

character string; units of flux (or usage rate) in the $cons column of uhist: "kg/day", (default), "kg/year" or "t/year" (metric tonnes per year)

fsphist

data.frame:
$year (num)
$f (num): fraction of imported solvent that is spilt, rather than disposed of safely
the transient safe disposal fraction of solvent by the site

fw

numeric [1]; proportion of imported solvent that becomes waste (between 0 and 1, 1 is suggested)

fi

numeric [1]; proportion of spilt solvent that infiltrates to water table (between 0 1, 1 is suggested)

start.t, end.t

numeric [1]; start and end time of the model, in days since 30/12/1899

dt

numeric [1]; time step length, in days, usually in the order of tens of days

x, y

x and y co-ordinates of the spill with the same origin and units as mfdata; used to find the grid cell from which to read MODFLOW results; alternatively, if mfdata is not used, x and y can simply be given as information to be stored with the model

z0

numeric [1] or "base"; elevation of the base of the DNAPL model, with the same datum as mfdata; "base" instructs to use the bottom elevation of the lowest layer at x,y, read from mfdata

mfdata

NetCDF object; MODFLOW data in NetCDF format (see GW.nc), from which the transient horizontal Darcy velocity through the source zone is determined

qh

1-layer: function(t) or numeric [1];
multi-layer: list [NLAY] of function(t), one per layer or numeric [NLAY], where NLAY is the number of layers in the DNAPLmodel ([email protected]);
if mfdata is missing, the horizontal Darcy velocity at the source zone may be given explicity as a function of time or a single constant value (or list of functions, one per layer, or numeric vector of values per layer if the model is multi-layer)

check.qh

logical [1]; whether to check the final set of qh functions for NA results (could be time-consuming for lots of time steps, but makes error catching easier)

Value

A DNAPLSourceTerm S4 object

Examples

1
# to follow

CJBarry/DNAPL documentation built on Aug. 2, 2017, 6:28 p.m.