cG: Constant and Convering power mass-flux relationship model
In CJBarry/DNAPL: DNAPL Source Term Framework

Usage Arguments Value

cstG.DNmodel(wg, wpm, hp, Gamma, Srn, Srw, phi, rho, Cs, hL,
  NLAY = length(hL))

cnvG.DNmodel(wg, wpm, hp, Gamma0, Srn, Srw, phi, rho, Cs, hL,
  NLAY = length(hL))

`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
`Gamma, Gamma0`	numeric `[1]` or `[NLAY]`; empirical source depletion parameter, positive; small values (<1) imply a more persistent source term for `cstG.DNmodel`, \Gamma is constant throughout the model, but for `cnvG.DNmodel`, \Gamma converges linearly from `Gamma0` (at peak mass) to 1 as the NAPL mass depletes.
`Srn, 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