A Douglas‐Jones Predictor‐Corrector Program for Simulating One‐Dimensional Unsaturated Flow in Soil

Abstract

Abstract. A fully documented program to represent one‐dimensional unsaturated flow in soil is described. The program is based on a Douglas‐Jones finite‐difference implicit method to solve the Richards equation. An implicit linearization scheme is used to estimate the hydraulic conductivity and specific moisture capacity functions. Predicted moisture content profiles compared with two Galerkin finite‐element solutions and field observations on a Panoche clay loam soil show very good agreement.