A mass-conservative method for the integration of the two-dimensional groundwater (Boussinesq) equation