For PISMDistributedHydrology only:  Realized you can't fail on checking P bounds if the geometry can change and thus the overburden pressure moves so the bounds fail for that reason.