Kernel: Python 3 (Ubuntu Linux)
https://rhodesmill.org/skyfield/
from skyfield.api import load planets = load('de421.bsp') earth, mars = planets['earth'], planets['mars'] ts = load.timescale() t = ts.now() astrometric = earth.at(t).observe(mars) ra, dec, distance = astrometric.radec() print(ra) print(dec) print(distance)
from skyfield.api import Topos boston = earth + Topos('42.3583 N', '71.0636 W') astrometric = boston.at(t).observe(mars) alt, az, d = astrometric.apparent().altaz() print(alt) print(az)