Change the way the CV mass matrix is found to use the hard coded 
sub cv element volume fraction * the element volume.

This adds a routine in femtools/FEFields.F90 called compute_cv_mass
which is called from from the added routine get_cv_mass in 
femtools/State_Fields.F90

These two routines are then called from all the places I could 
see that wanted a cv mass matrix. Then where these calls happen variables
that are called something like lumped_mass are changed to cv_mass.

This needs carefully checking. The compute_cv_mass will work 
for all linear and quadatic line and triangle. For tet the 
subcv fraction values need entering.

All cv test cases test this.