Abstract
In this study, the Eulerian-Lagrangian method was extended to a three dimensional mass transport problem for a density dependent saturated-unsaturated flow. Validity of this method was proved in comparison with experiments and anther 3-dimensional program (3DFEMFAT). A new algorithm was also proposed here. It decreased a computation time sharply compared with a usual algorithm.