initial implementation of water fraction diagnostic routine
fix lagrange multiplier in CTDCast.water_fraction

...and add a simple test

also enforce output of kronecker product to be "csr" for speed