I'm aware of the NASA orbit viewer but it limits you to years > 1600 AD.
Is there an orbit viewer for any date?
Astronomy Stack Exchange is a question and answer site for astronomers and astrophysicists. It only takes a minute to sign up.
Sign up to join this communityI'm aware of the NASA orbit viewer but it limits you to years > 1600 AD.
Is there an orbit viewer for any date?
From inspecting the JPL ephemeris website, the ephemeris with the longest time span as of this writing is DE441, which spans from August 15th, -13200 to March 15th, 17191. It is less accurate than DE440 (which covers the time span 1550-2650 AD), but is more stable over long time periods. For reference, between the years 1970-2020 the difference between DE441 and DE440 is less than 2 meters.
You can access this ephemeris using astropy. For example, I prepared the following script from astropy's tutorial, evaluating the Moon's position for Julian year -3000. The script will take a while to run the first time (since it needs to download the ephemeris), but subsequent runs will be much faster.
import numpy as np
import matplotlib.pyplot as plt
import astropy.units as u
import astropy.time
import astropy.coordinates
import astropy.visualization
astropy.visualization.quantity_support()
astropy.coordinates.solar_system_ephemeris.set('de441_part-1')
time_stop = astropy.time.Time('-3000', format='jyear')
time_start = time_stop - 27 * u.day
time = np.linspace(time_start, time_stop, num=1000)
location = astropy.coordinates.EarthLocation.from_geocentric(0 * u.m, 0 * u.m, 0 * u.m)
coords = astropy.coordinates.get_moon(time, location).cartesian
plt.figure(figsize=(9, 9))
plt.gca().set_aspect('equal')
plt.plot(coords.x, coords.y)
plt.scatter(coords[~0].x, coords[~0].y)
plt.savefig('moon.png')