A 2-D ADI model was set up from the depth-averaged hydrodynamic continuity equation and momentum equation and the depth integrated conservation equation of solute transport.