An underrelaxed-modified Picard iteration scheme for simulation of 3D wetting pattern under drip irrigation using Richards’ equation on non-orthogonal grids