Download the Richards' equation implementation for MODFLOW 6 and the complete input files for the three benchmark problems presented in the paper.
This page accompanies the paper Evaluation and Implementation of Fully 3-D Richards’ Equation in MODFLOW 6, by Sebastián Vázquez-Gasty, Sorab Panday, Christian Langevin, and Tirthankar Roy — a collaboration between GSI Environmental, S.S. Papadopulos & Associates and the University of Nebraska-Lincoln.
The work implements a fully three-dimensional Richards’ equation solver within the MODFLOW 6 framework, so that the unsaturated zone and the groundwater system are represented together in a single model. The formulation uses the mixed form of Richards’ equation with van Genuchten constitutive relations and Newton-Raphson linearization.