DDpg.DNmodel: Dual-domain Pool and Ganglia model
In CJBarry/DNAPL: DNAPL Source Term Framework

Description Usage Arguments Value

This model has separate domains for ganglia and pools. Ganglia are flushed efficiently (effectively Gamma = Inf) and pools are flushed inefficiently, with a poolGamma depletion parameter that defaults to 1.

1 2	DDpg.DNmodel(wg, wpm, hp, Srn, Srw, phi, rho, Cs, hL, NLAY = length(hL), poolGamma = 1)

`wg`	numeric `[1]` or `[NLAY]`; ganglia width
`wpm`	numeric `[1]` or `[NLAY]`; maximum pool width (may be `Inf` in the lowest layer, in which case no DNAPL is allowed to escape)
`hp`	numeric `[1]` or `[NLAY]`; pool height, or total height of pools in a layer
`Srn`	numeric `[1]`; residual saturations of NAPL and water
`Srw`	numeric `[1]`; residual saturations of NAPL and water
`phi`	numeric `[1]` or `[NLAY]`; bulk porosity of each layer
`rho`	numeric `[1]`; solvent density (ensure consistent units, probably kg/m^3)
`Cs`	numeric `[1]`; solvent solubility in water (ensure consistent units, probably kg/m^3 which is the same as g/l)
`hL`	numeric `[1]` or `[NLAY]`; the height of each layer; note that higher layers will contain more NAPL mass
`NLAY`	integer `[1]`; number of layers in the DNAPL model
`poolGamma`	numeric `[1]` or `[NLAY]`; the empirical depletion power to be applied to pools, analogous to `Gamma` in `cstG.DNmodel` which applies to the bulk NAPL