I've found the gradient sign error and committed the fix to trunk/, branches/version-1_2, and branches/anisotropy[*]. Unfortunately, it doesn't seem to do anything to reduce the magnitude of your divergences. I suspect that may be a result of mesh non- orthogonality, rather than normal misorientation. I'll run it by Daniel when he gets back next week(?).
On Feb 4, 2008, at 4:29 PM, asreeve wrote:
I think this is why my previous scripts involving a convection term where not working (I was getting negative concentrations in some cases).
Yes, it seems to fix that problem (which sort of surprises me, because the error in sign seemed to impact very, very little unless you explicitly asked for a gradient).
[*] As Daniel said, you should switch over to trunk, as your anisotropic diffusion is merged there now and then we can terminate anisotropy branch to remove any confusion.
