DNST: DNAPL source term model
In CJBarry/DNAPL: DNAPL Source Term Framework

Description Usage Arguments Value Examples

DNAPL source term model

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)

`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` (`DNAPLmodel@NLAY`); 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)