Apparently working strong boundary conditions test. To get this working it was necessary to impose the boundary conditions on the initial guess passed to the nonlinear solver.