Allow multiple QM metal centers

New feature. Need to replace Spec1.sh and Spec2.sh with calls to
so_metal1.py and so_metal2.py (note to self: fix this tomorrow)