FIG. 9. Flow domain showing boundary conditions. The surface boundary condition is a specified flux, J_w⁰, whereas the bottom and side boundary conditions are specified as no flow conditions by setting the gradient in matric potential, , in the x and z directions equal to zero.