[Spice_discussion] discrepancy SPICE vs horizons

Hans-Peter Doerr doerr at kis.uni-freiburg.de
Tue Feb 5 15:13:54 PST 2013


Dear Navigators,

while testing a SPICE-based code, I saw a significant discrepancy
between values obtained via horizons and those computed via CSPICE. I
boiled it down to the following code snippet:

   SpiceDouble et, lt, r, rr, relstate[6], los[3];
   furnsh_c ("kernels/pck00010.tpc");
   furnsh_c ("kernels/naif0010.tls");
   furnsh_c ("kernels/de405.bsp");
   furnsh_c ("kernels/earth_latest_high_prec.bpc");
   furnsh_c ("kernels//earthstns_fx_050714.bsp");
   furnsh_c ("kernels/earth_fixed_itrf93.tf");

   str2et_c ("2011-Aug-01 06:00:00.0000", &et);
   spkezr_c ("Sun", et, "J2000", "NONE", "DSS-12", relstate, &lt);
   unorm_c (relstate, los, &r);
   rr = vdot_c (los, &relstate[3]);

   printf ("dist: %.16E  v_los: % .16E\n", r, rr);

this gives:
   dist: 1.5185605413773715E+08  v_los: -4.6934776127220079E-02

the 'earth_fixex_itrf93' text kernels binds EARTH_FIXED to ITRF93 as 
suggested.

whith these horizons settings:

   Target body name: Sun (10)                 {source: DE405}
   Center body name: Earth (399)              {source: DE405}
   Center-site name: DSS 12
   EOP file        : eop.130204.p130428
   EOP coverage    : DATA-BASED 1962-JAN-20 TO 2013-FEB-04. PREDICTS-> 
2013-APR-27
   Output type     : GEOMETRIC cartesian states
   Coordinate systm: Earth Mean Equator and Equinox of Reference Epoch


I get (jd, range, range-rate):

   2455774.750000000, 1.518560571902532E+08, -4.540414938447372E-02

these values are off by about 3 km and 1.5 m/s. I don't see what could 
be the issue here.

Any Ideas?


Thanks!

  Hans-Peter


More information about the Spice_discussion mailing list