This matlab program solves the 1D Non Equilibrium Richards Equation (NERE) in initially dry snow considering a water flux at the top boundary (can be modified) and afree-drainage boundary condition at the bottom.
In this program, snow is assumed homoegeneous and the snow density and grain size can be modified in the file, as well as somw model parameters.
If using this program and a modified version of this program, please cite: Leroux, N.R. and J.W. Pomeroy (2018), Simulation of Capillary Overshoot in Snow Combining Trapping of the Wetting Phase with a Non-Equilibrium Richards Equation Model, Water Resources Research, doi: 10.1029/2018WR022969.