The Advection-Dispersion Equation (p. 342) is solved numerically using
an implicit finite difference scheme which, in principle, is unconditionally
stable and has negligible numerical dispersion. A correction term has been
introduced in order to reduce the third order truncation error. This correction
term makes it possible to simulate advection-dispersion of concentration
profiles with very steep fronts.