On Fri, Feb 22, 2013 at 06:49:46PM +0100, Ian Hinder wrote:
Are you using sufficient ghost zones?
Dissipation checks for that and aborts if it is not the case. The only reason I can think of right now why random things could happen is if the primary variables (not the rhs) are not updated in the ghost- or boundary zones before Dissipation is called in post_rhs.
Frank