Remove the two darcy test cases that were there in 1d and 
replace them with versions that do 1,2 and 3d.

Change the sign of the ct_rhs add to in Divergence_CV.F90 
for incompressible flow.

Adjust the darcy test that was there to test this so that 
it runs and the 1,2 and 3d version use weak velocity BC.

Add a similar test that uses two phase but one with a large 
absorption such as to be the solid porous media. This tests
the multiphase ct_rhs add to sign change.

The mms_ns case still has issues.