A modelling procedure of the two-dimensional advection-dispersion equation based on a finite volume method for the spatial discretization is presented in this paper. Control volumes are made up over a standard triangulation, which give the grid nodes where we obtain the numerical solutions. A first order accurate interpolation is used over each triangle that contains the control volume. Some idealized initial value problems are studied showing the practical advantages and disadvantages of the method in comparison with analytical solutions.