import ephem import math import time from datetime import datetime, timedelta greenwich = ephem.Observer() greenwich.lat = "0" greenwich.lon = "0" greenwich.date = datetime.utcnow() sun = ephem.Sun(greenwich) sun.compute(greenwich.date) sun_lon = math.degrees(sun.ra - greenwich.sidereal_time() ) if sun_lon < -180.0 : sun_lon = 360.0 + sun_lon elif sun_lon > 180.0 : sun_lon = sun_lon - 360.0 sun_lat = math.degrees(sun.dec) print('%s %s' % (sun_lon, sun_lat)) quit()