Hello kamensky. Thanks for the tip on the improved flux calculation, it’s actually quite useful going forward. It didn’t solve the specific issue I have, and I’ve written up a new post describing the problem more thoroughly here (I couldn’t edit this post anymore) Gradient with a boundary point source