Multilayer, sharp-interface, saltwater intrusion model. The code is delivered with complete documentation and user's manual.
Technical support is available.
The equations include flux terms that account for vertical leakage in intervening aquitard units. A fully implicit solution technique for a multiphase flow problem is used to obtain the numerical solution. Spatial discretization is performed using rectangular finite elements. Element matrices are computed using highly efficient influence-coefficient formulas. Nonlinearity resulting from water table and interface boundary conditions are handled via a modified Newton-Raphson algorithm that incorporates upstream weighting of relative permeabilities.