Consegui obter os vetores orbitais cartesianos de todos os principais órgãos de HORIZON apenas na época J2000. Eu poderia estender a cobertura adiante através do tempo. É fácil obter sobrecarga de dados fazendo isso. Minha simulação é modelada usando apenas as Leis da Gravitação e Movimento. Isso fornece resultados surpreendentemente próximos aos publicados. Executar o sistema solar para trás (invertendo os vetores de velocidade) me deu os vetores iniciais de volta a 1900. Isso é tudo que eu precisava e os resultados foram próximos o suficiente para meus propósitos. Eu ainda tenho os arquivos CSV.
Eu também tive todos os tipos de problemas com a interface dos horizontes. Por exemplo, alterar a data não teve efeito no valor dos vetores. ou seja: todas as datas de início especificadas têm os mesmos valores. Ultimamente, não consegui duplicar essa façanha. Obviamente, existem alguns problemas sérios nessa interface, principalmente nos últimos tempos.
Eu sei que os dados que obtive estavam corretos porque se correlacionam perfeitamente com eventos publicados, por exemplo: o trânsito recente de Mercúrio.
Eu também ainda estou procurando esse tipo de dados.
Você provavelmente já se mudou há muito tempo, mas, apenas para referência, as condições iniciais que o HORIZONS usa são mencionadas ("header.431_572") em ftp://ssd.jpl.nasa.gov/pub/eph/planets/Linux/ README.txt, mas o único local em que eu poderia encontrá-los no "formulário da tabela" é no meu próprio repositório git:
https://github.com/barrycarter/bcapps/blob/master/ASTRO/header.431_572
Os valores são explicados em http://ilrs.gsfc.nasa.gov/docs/2014/196C.pdf a partir da página 39, "VI. Condições e constantes iniciais", especialmente nas tabelas que começam com a Tabela 4 na página 47 e terminando com Tabela 13 na página 74.
Escrevi scripts para configurar as condições iniciais e resolver numericamente as equações diferenciais usando o Mathematica; portanto, o seguinte pode 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
Quando reduzo o tamanho da etapa o suficiente (o tamanho padrão da etapa do Mathematica é muito grande), meus resultados correspondem aos de HORIZONS:
https://github.com/barrycarter/bcapps/blob/master/ASTRO/bc-integrate-compare.m
fonte