> This is all dodgy, and there are several edge effects that will > probably result in you not getting what you want. Daniel, do you see > an easier way to get this d(phi)/dx=constant*d2(phi)/dy2 boundary > condition that Jason wants?
I would say that the surface flux that Jason wants is,
surfaceFlux = mask * constant * phi.getFaceGrad().dot((0,
1)).getDivergence().getArithmeticFaceValue().dot((1,0)
Then add surfaceFlux.getDivergence() to the equation.
Cheers
--
Daniel Wheeler
