A new method to solve the quasi-3D shallow water equations is proposed. This method combines a suitable mass-preserving finite element approach in the horizontal plane with a conventional conforming finite element (or finite difference) scheme along the vertical