#include <stdio.h>
#include "calceph.h"

/*-----------------------------------------------------------------*/
/* main program */
/*-----------------------------------------------------------------*/
int main()
{
 int res;
 double AU, EMRAT, GM_Mer;
 double jd0=2451624;
 double dt=0.5E0;
 double PV[6];

 /* open the ephemeris file */
 res = calceph_sopen("example1.dat");
 if (res)
 {
    printf("The ephemeris is already opened\n");
   /* print the values of AU, EMRAT and GM_Mer */
   if (calceph_sgetconstant("AU", &AU))
       printf("AU=%23.16E\n", AU);

   if (calceph_sgetconstant("EMRAT", &EMRAT))
     printf("EMRAT=%23.16E\n", EMRAT);

   if (calceph_sgetconstant("GM_Mer", &GM_Mer))
     printf("GM_Mer=%23.16E\n", GM_Mer);

   /* compute and print the coordinates */
   /* the geocentric moon coordinates in AU and AU/day */
   calceph_scompute(jd0, dt, 10, 3, PV);
   printcoord(PV,"geocentric coordinates of the Moon in AU and AU/day");

   /* the value TT-TDB */
   calceph_scompute(jd0, dt, 16, 0, PV);
   printf("TT-TDB = %23.16E\n", PV[0]);

   /* the heliocentric coordinates of Mars */
   calceph_scompute(jd0, dt, 4, 11, PV);
   printcoord(PV,"heliocentric coordinates of Mars");

   /* close the ephemeris file */
   calceph_sclose();
   printf("The ephemeris is already closed\n");
 }
 else
 {
   printf("The ephemeris can't be opened\n");
 }
 return res;
}