calculating volume in surface agrees with NKS
have partly working code to determine voxels on atom of interest, need to find why all the voxel and atom coordinates are off by half the cube size...