/* Most of these variables are used only in calculations
 * for the outer planets:
 */
double sinV = 0.0;
double cosV = 0.0;
double sin2V = 0.0;
double cos2V = 0.0;
double sinW = 0.0;
double sinze = 0.0;
double sin2ze = 0.0;
double sin3ze = 0.0;
double sin4ze = 0.0;
double cosze = 0.0;
double cos2ze = 0.0;
double cos3ze = 0.0;
double cos4ze = 0.0;
double sinQ = 0.0;
double sin2Q = 0.0;
double sin3Q = 0.0;
double cosQ = 0.0;
double cos2Q = 0.0;
double cos3Q = 0.0;
double sin5ze = 0.0;
double cos5ze = 0.0;
double sin1H = 0.0;
double sin2H = 0.0;
double cos1H = 0.0;
double cos2H = 0.0;
double sinth = 0.0;
double sin2th = 0.0;
double costh = 0.0;
double cos2th = 0.0;
double M1 = 0.0;
double M2 = 0.0;
double M4 = 0.0;
double M5 = 0.0;
double M6 = 0.0;
double M7 = 0.0;
double M8 = 0.0;
double nu = 0.0;
double P = 0.0;
double Q = 0.0;
double S = 0.0;
double V = 0.0;
double W = 0.0;
double ze = 0.0;
double psi = 0.0;
double H = 0.0;
double G = 0.0;
double eta = 0.0;
double th = 0.0;
double f = 0.0;

/* mean anomaly of the earth
 */
double M = 0.0;
extern double M;

/* reusable time variable,
 * centuries since 1900.0
 */
double T = 0.0;
extern double T;


#include "planet.h"
#include "kep.h"


/* This subroutine calculates mean anomalies of the planets.
 * The mean anomalies of earth, uranus and neptune are in
 * oearth.c, ouranus.c, and oneptune.c.
 */
manoms(J)
double J;
{

T = (J - 2415020.0)/36525.0;

#if BandS
/* mean anomaly of the earth (and sun) */
f = (( -0.0000033*T - 0.000150)*T + 35999.04975)*T + 358.47583;
M = mod360(f);
#endif

/* mean anomaly of Mercury */
M1 = (0.000007*T + 149472.51529)*T + 102.27938;
M1 = mod360(M1);
/* Venus */
M2 = (0.001286*T + 58517.80387)*T + 212.60322;
M2 = mod360(M2);
/* Mars */
M4 = (0.000181*T + 19139.85475)*T + 319.51913;
M4 = mod360(M4);
/* Jupiter */
M5 = (-0.000722*T + 3034.69202)*T + 225.32833;
M5 = mod360(M5);
/* Saturn */
M6 = (-0.000502*T + 1221.55147)*T + 175.46622;
M6 = mod360(M6);
}