A two-constituent solute-transport model for ground water having variable density

A numerical model has been developed to simulate solute transport and dispersion of either one or two constituents in groundwater where there is two-dimensional, density-dependent flow. The model is a modified version of the one documented by Konikow and Bredehoeft (1978), which uses finite-difference methods and the method of characteristics to solve the flow and transport equations. The model was tested on an idealized seawater intrusion problem for which an analytical solution has been developed. The results were nearly identical to those of other numerical models tested on the same problem. A description of the formats for the input data, a sample of input and output for a two-constituent example problem, and a listing of the Fortran program are presented. (Author 's abstract)