hL.setup: Make DNAPL model layer heights appropriate to a MODFLOW data...
In CJBarry/DNAPL: DNAPL Source Term Framework

Description Usage Arguments Value Examples

Make DNAPL model layer heights appropriate to a MODFLOW data set

1 2	hL.setup(mfdata, x, y = NULL, ndlpmfl = 1L, bph = 0, z0 = "base", bph.warning = TRUE)

`mfdata`	NetCDF; a MODFLOW data set (see `GW.nc`)
`x, y`	for `xy.coords`; location of the source zone, in the same co-ordinate system as `mfdata`; one location only
`ndlpmfl`	integer `[1]` or `[mfNLAY]`, where `mfNLAY` is the number of layers in the MODFLOW model; number of DNAPL model layers per MODFLOW layer (including any which are always dry at the location, although the values for these layers won't make any difference)
`bph`	numeric [1] or [mfNLAY]; height of basal pools in each layer, or layer-by-layer: note that length-1 input for `bph` will put a basal pool in every MODFLOW layer, not just the bottom; put `c(0, 0, ..., 0, h)` if you want just the bottom layer to have a basal pool
`z0`	numeric `[1]` or `"base"`; elevation of the bottom of the DNAPL source zone, with the same datum as `mfdata`; `"base"` instructs to read the bottom of the MODFLOW model at this location
`bph.warning`	logical [1]; `hL.setup` will give a warning with length-1 input for `bph`, because the result may be unexpected. If you really do want a basal pool in each layer, you can set `bph.warning = FALSE` to suppress the warning. Do this if you know what you are doing.

numeric []; heights of layers of the DNAPL model, from top to bottom

mfdata <- RNetCDF::open.nc(system.file("rflow_mf_demo.nc",
                                       package = "Rflow"))

hL.setup(mfdata, 625, 825, 1L)
hL.setup(mfdata, 625, 825, 3L)
hL.setup(mfdata, 625, 825, 3L, z0 = 10)

# include a height-1 basal pool
hL.setup(mfdata, 625, 825, 3L, 1)