Created wcs_util.py, added add_celestial_axis. Tests to follow.