Abstract: A two-dimensional transient model for flow through saturated-unsaturated porous media is developed, The model numerically solves the pressure head dependent or moisture content dependent form of Richard's equation. The model code uses isoparametric quadratic triangular and/or quadrilateral finite elements for the geometric representation and for the weak Galerkin spacial integrations. An implicit, unconditionally stable single-step numerical time integration scheme with an oscilliatory noise reduction option is utilized for the temporal discretization. The highly efficient...