rmsd.py and molecule.py working, but support for coordinate subsets incomplete