added mercury_reference.py with observational constraints