first try to integrate astropy.units into sun.py