The quantity sum_xh0_int was not being calculated correctly for the non ALLFRAC
case. It should not be 1.0 - sum_xh1_int, but real(mesh(1)*mesh(2)*mesh(3)) - sum_xh1_int.
The global average fractions are not normalized during this calculation hence
this factor instead of 1.0
This is now fixed.