This merge will add the capability for the enforce discrete property solenoidal projection when 
the lagrange multiplier is CG to test the divergence of the interpolated field with the CV dual. 
This is akin to what was added for the velocity-pressure solve previously via this branch. 
This adds an option to the schema, commenting that this will make the matrix non symmetric.

There is also a small correction to an addition of the known velocity divergence to the projection 
rhs to use ctp_m rather than ct_m (in this solenoidal projection routine).

This is tested via the solenoidal unittest, with an extra case added to check for CV divergence. 
Non used code was also removed from here.