Bug fix to new shrink-sphere-centre algorithm
Add enclosed mean density property to profile
Minor optimization to virial_radius
Update _com .c file to match .pyx