volume_dat2.c reading cube1.dat 
status=0, zmax=0.5, points=8, polys=6 

xmin=-0.500000, xmax=0.500000, ymin=-0.500000, ymax=0.500000 
zmin=-0.500000, zmax=0.500000 
enclosing area= 6, enclosing volume= 1 

final total area = 6, total volume = 1