Read_GMSH assumed that there were only two types of elements in the msh file and worked out the number of volume elements by subtracting the number of facets from the total number of elements. In 3D if gmsh produces interval elements then this didn't actually give the number of tets/hexes. This merge fixes this by working out the number of volume elements properly (from the number of tets/hexes actually in the file). This in effect drops any interval elements, which is the same behaviour as produced by going through gmsh2triangle.