Application of an improved P(m)-SOR iteration method for flow in partially saturated soils