How to set gravitation potential

[kongdd]

The Richards Equation is:

$$ \frac{\partial \theta}{\partial t}=\frac{\partial}{\partial z} [D \frac{\partial \theta}{\partial z}+K ] $$

It seems that Fronts.jl ignored the gravitation potential term $\frac{\partial K}{\partial z}$.
Any suggestion how to solve this problem?

[gerlero]

Indeed Fronts implements the horizontal case of the Richards equation, where such a term is neglected (which is the case in pretty much all of our problems in microfluidics).

I'm not against adding support for such a term, but it'd need to dispatch to the FiniteDifference() algorithm because the default BoltzmannODE() algorithm is not compatible with that term.