[CIG-LONG] Conservation of Mass

[Walter Landry]

Brian Wilson <bjwilson83 at gmail.com> wrote:
> Hi all,
> 
> I'm working on a model of lithospheric drip and despite having zero flux on
> the boundaries, I am getting downward velocities with no compensating
> upwelling, which would seem to violate the divergence free condition. Has
> anybody else had a similar problem and/or figured out a solution? I've
> attached a copy of my input file if anyone wants to look at it; everything
> is scaled as follows: distance: km, mass: 10^15 kg, time: m.y.

To solve this with the current elements, you have to add a
HydrostaticTerm.  Depending on your setup, that may solve your
problems.

Cheers,
Walter Landry
walter at geodynamics.org