Richard’s equation was approximated by finite-difference numerical scheme to model water infiltration profile in variably unsaturated soil. The published data of Philip’s semi-analytical solution was used to validate the simulated results from the numerical scheme. A discrepancy was found between the simulated and the published semi-analytical results. Morris method as a global sensitivity tool was used as an alternative to local sensitivity analysis to assess the results discrepancy. Morris method with different sampling strategies were tested, of which Manhattan distance method have resulted a better sensitivity measures and also a better scan of input space than Euclidean method. Moreover, Morris method at and Manhattan distance sampling strategy, with only 2 extra simulation runs than local sensitivity analysis, was able to produce reliable sensitivity measures ( , ). The sensitivity analysis results were cross-validated by Sobol’ variance-based method with 150,000 simulation runs. The global sensitivity tool has identified three important parameters, of which spatial discretization size was the sole reason of the discrepancy observed. In addition, a high proportion of total output variance contributed by parameters and is suggesting a greater significant digits is required to reduce its input uncertainty range.