from skyfield.api import Topos, load

ts = load.timescale()
t = ts.now()

planets = load('de421.bsp')
mars = planets['mars']

# From the center of the Solar System (Barycentric)

barycentric = mars.at(t)
print(barycentric)

# From the center of the Sun (Heliocentric)

sun = planets['sun']
heliocentric = sun.at(t).observe(mars)
print(heliocentric)

ra, dec, distance = heliocentric.radec()
hours = ra.hours
degrees = ra._degrees

print(ra)
print(dec)
print(distance)

print("Hours: ", hours)
print("Degrees: ", degrees)

# From the center of the Earth (Geocentric)

earth = planets['earth']
astrometric = earth.at(t).observe(mars)
apparent = earth.at(t).observe(mars).apparent()

# From a place on Earth (Topocentric)

boston = earth + Topos('42.3583 N', '71.0603 W')
astrometric = boston.at(t).observe(mars)
apparent = boston.at(t).observe(mars).apparent()