one more file, might be unstable
sphere, grd_split: split populations by metallicity, velocity, and
position in Plummer sphere. change in overall structure, introduced
matrix representation for \hat{p}_MW calculations, use vector
calculations by numpy for speedup. Still buggy.
first debug round grd_split, deactivated exception for sanitize