Hi all,<br><br>I&#39;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&#39;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.<br>
<br>Thanks,<br>Brian<br>