Pude obtener los vectores orbitales cartesianos para todos los cuerpos principales de HORIZON solo en la época J2000. Podría extender la cobertura a través del tiempo. Es fácil obtener una sobrecarga de datos haciendo esto. Mi simulación está modelada usando solo las Leyes de Gravitación y Movimiento. Esto da resultados que son sorprendentemente cercanos a los publicados. Ejecutar el sistema solar hacia atrás (invirtiendo los vectores de velocidad) me ha devuelto los vectores iniciales a 1900. Esto es todo lo que necesitaba y los resultados fueron lo suficientemente cercanos para mis propósitos. Todavía tengo los archivos CSV.
También he tenido todo tipo de problemas con la interfaz de horizontes. Por ejemplo, cambiar la fecha no tuvo efecto en el valor de los vectores. es decir: todas las fechas de inicio especificadas tienen los mismos valores. Últimamente, no he podido duplicar esta hazaña. Obviamente, hay algunos problemas graves con esta interfaz, especialmente últimamente.
Sé que los datos que obtuve eran correctos porque se correlaciona perfectamente con los eventos publicados, por ejemplo, el tránsito reciente de Mercurio.
Yo también sigo buscando este tipo de datos.
Probablemente haya avanzado desde hace mucho tiempo, pero, solo como referencia, se mencionan las condiciones iniciales que utiliza HORIZONTES ("header.431_572") en ftp://ssd.jpl.nasa.gov/pub/eph/planets/Linux/ README.txt pero el único lugar donde podría encontrarlos en "forma de tabla" es en mi propio repositorio git:
https://github.com/barrycarter/bcapps/blob/master/ASTRO/header.431_572
Los valores se explican en http://ilrs.gsfc.nasa.gov/docs/2014/196C.pdf a partir de la página 39, "VI. Condiciones iniciales y constantes", especialmente en las tablas que comienzan con la Tabla 4 en la página 47 y terminando con la Tabla 13 en la página 74.
He escrito guiones para configurar las condiciones iniciales y resolver numéricamente las ecuaciones diferenciales usando Mathematica, por lo que lo siguiente puede ser útil:
https://github.com/barrycarter/bcapps/blob/master/ASTRO/README
https://github.com/barrycarter/bcapps/blob/master/ASTRO/bc-header-values.pl
https://github.com/barrycarter/bcapps/blob/master/ASTRO/bc-integrate.m
Cuando reduzco el tamaño de paso lo suficiente (el tamaño de paso predeterminado de Mathematica es demasiado grande), mis resultados coinciden estrechamente con los de HORIZONTES:
https://github.com/barrycarter/bcapps/blob/master/ASTRO/bc-integrate-compare.m
fuente