add center_of_mass function

[loriab]

We needed this yesterday on models.Molecule. Posting it in case anyone wants a starter project. Should test on example that does and doesn't have fix_com set. (True should give origin.)