Made it into a scalar problem. Easier to check against an analytical solution and going to a 2D system is pretty straightforward.