De-integrate the pressure gradient term by parts.
Change over to a pressure boundary condition on the right
solution.xml.gz -> ns-solution.xml.gz
Change Re = 40. Newton fails on Re=50, but Newton+NGMRES converges (try it out with Krylov later)