The convergence difficulty associated with solving the non-linearity of the Richards's equation is often encountered in numerical simulation of infiltration into dry soils driven by gravity. In previous work we proposed a simple approach to solve the gravity term for one-dimensional infiltration, and test on a wide range of problems shows considerable improvement over the existing Picard and Newton methods. The purpose of this note is to explain how to extend the method to multi-dimensional finite element solution and demonstrate its superiority to the standard Picard and Newton methods.