1 #include "libnav.h"

2 #include "nav.h"

19  // something wrong with fsol? Need to check

20  // check jd

21 double esdist_(int *iyr, int *iday, int *msec) {

26  double t = jd(*iyr, 1, *iday) - 2451545.0e0 + (*msec / 1000.e0 - 43200.e0) / 86400.e0;

27  // c Compute mean anomaly

28  double gs = 357.52772e0 + 0.9856002831e0 * t;

29  // Compute solar distance (AU)

30  const static double radeg_to_match_fortran = 57.295780181884766e0;

31  return 1.00014e0 - 0.01671e0 * cos(gs / radeg_to_match_fortran) - 0.00014e0 * cos(2.0e0 * gs / radeg_to_match_fortran);

32 }