New version of the incompressible adjoint solver

The test cases are going to fail because the change in the Jacobians