For the linear equation solution a preconditioned conjugate gradient method is used.
Solute transport is handled by an upstream weighting scheme.
The two-dimensional flow can be either in a horizontal or in a vertical plane. Furthermore, the model can handle radial symmetric simulations. The code contains a restart feature for simulation of changing boundary conditions. Information on the root zone includes wilting pressure, maximum transpiration rate, root effectiveness function, and root zone geometric data. Analytical fuctions must be used for relative hydraulic relationship and moisture characteristic curve fuctions.