FACT++  1.0
double eraEqeq94 ( double  date1,
double  date2 
)

Definition at line 3 of file eqeq94.c.

References eraAnpm(), eraNut80(), eraObl80(), ERFA_D2PI, ERFA_DAS2R, ERFA_DJ00, ERFA_DJC, and t.

Referenced by eraGst94(), and t_eqeq94().

57 {
58  double t, om, dpsi, deps, eps0, ee;
59 
60 /* Interval between fundamental epoch J2000.0 and given date (JC). */
61  t = ((date1 - ERFA_DJ00) + date2) / ERFA_DJC;
62 
63 /* Longitude of the mean ascending node of the lunar orbit on the */
64 /* ecliptic, measured from the mean equinox of date. */
65  om = eraAnpm((450160.280 + (-482890.539
66  + (7.455 + 0.008 * t) * t) * t) * ERFA_DAS2R
67  + fmod(-5.0 * t, 1.0) * ERFA_D2PI);
68 
69 /* Nutation components and mean obliquity. */
70  eraNut80(date1, date2, &dpsi, &deps);
71  eps0 = eraObl80(date1, date2);
72 
73 /* Equation of the equinoxes. */
74  ee = dpsi*cos(eps0) + ERFA_DAS2R*(0.00264*sin(om) + 0.000063*sin(om+om));
75 
76  return ee;
77 
78 }
#define ERFA_DJ00
Definition: erfam.h:87
void eraNut80(double date1, double date2, double *dpsi, double *deps)
Definition: nut80.c:3
#define ERFA_DAS2R
Definition: erfam.h:60
#define ERFA_D2PI
Definition: erfam.h:48
#define ERFA_DJC
Definition: erfam.h:81
double eraAnpm(double a)
Definition: anpm.c:3
double eraObl80(double date1, double date2)
Definition: obl80.c:3
TT t
Definition: test_client.c:26

+ Here is the call graph for this function:

+ Here is the caller graph for this function: