A finite element dual mesh method to calculate Nodal Darcy velocities in nonhomogeneous and anisotropic aquifers