Comment out this surface integral (that probably ought to be here).
Makes convergence of the Krylov methods very slow.