The work of Niccolò Tubini is going who already developed a very solid Richards1D (with ponding), code coupled with the energy budget is proceeding towards a 2D version on unstructured grids coupled, at present with a 1D de Saint-Venant equation. This is the summary of the work done so far given at the Italian Hydrological society meeting in Bologna.
Actually the de Saint-Venant coupling is not yet ready but it will be very soon. Stay tuned. Click on the figure for getting the presentation.