Hydrologic model of the saturated zone beneath the Pajarito Plateau using the DPFEHM simulator. The model is differentiable and gradients can be efficiently computed with Zygote.
To run from in Julia, simply do
include("ex.jl")
The first time it is run it will download a mesh, boundary conditions, etc. from the waffle2017 model. Then it runs a steady-state groundwater flow simulation.