yodaUtility/sgp4/globals.h

//
// globals.h
//
#pragma once

#include "math.h"

const double PI           = 3.141592653589793;
const double TWOPI        = 2.0 * PI;
const double RADS_PER_DEG = PI / 180.0;

const double GM           = 398601.2;   // Earth gravitational constant, km^3/sec^2
const double GEOSYNC_ALT  = 42241.892;  // km
const double EARTH_DIA    = 12800.0;    // km
const double DAY_SIDERAL  = (23 * 3600) + (56 * 60) + 4.09;  // sec
const double DAY_24HR     = (24 * 3600);   // sec

const double AE           = 1.0;
const double AU           = 149597870.0;  // Astronomical unit (km) (IAU 76)
const double SR           = 696000.0;     // Solar radius (km)      (IAU 76)
const double TWOTHRD      = 2.0 / 3.0;
const double XKMPER_WGS72 = 6378.135;     // Earth equatorial radius - km (WGS '72)
const double F            = 1.0 / 298.26; // Earth flattening (WGS '72)
const double GE           = 398600.8;     // Earth gravitational constant (WGS '72)
const double J2           = 1.0826158E-3; // J2 harmonic (WGS '72)
const double J3           = -2.53881E-6;  // J3 harmonic (WGS '72)
const double J4           = -1.65597E-6;  // J4 harmonic (WGS '72)
const double CK2          = J2 / 2.0;
const double CK4          = -3.0 * J4 / 8.0;
const double XJ3          = J3;
const double E6A          = 1.0e-06;
const double QO           = AE + 120.0 / XKMPER_WGS72;
const double S            = AE + 78.0  / XKMPER_WGS72;
const double HR_PER_DAY   = 24.0;          // Hours per day   (solar)
const double MIN_PER_DAY  = 1440.0;        // Minutes per day (solar)
const double SEC_PER_DAY  = 86400.0;       // Seconds per day (solar)
const double OMEGA_E      = 1.00273790934; // earth rotation per sideral day
const double XKE          = sqrt(3600.0 * GE /           //sqrt(ge) ER^3/min^2
                                (XKMPER_WGS72 * XKMPER_WGS72 * XKMPER_WGS72)); 
const double QOMS2T       = pow((QO - S), 4);            //(QO - S)^4 ER^4

// Utility functions
double sqr   (const double x);
double Fmod2p(const double arg);
double AcTan (const double sinx, double cosx);

double rad2deg(const double);
double deg2rad(const double);
1	//
2	// globals.h
3	//
4	#pragma once
5
6	#include "math.h"
7
8	const double PI = 3.141592653589793;
9	const double TWOPI = 2.0 * PI;
10	const double RADS_PER_DEG = PI / 180.0;
11
12	const double GM = 398601.2; // Earth gravitational constant, km^3/sec^2
13	const double GEOSYNC_ALT = 42241.892; // km
14	const double EARTH_DIA = 12800.0; // km
15	const double DAY_SIDERAL = (23 * 3600) + (56 * 60) + 4.09; // sec
16	const double DAY_24HR = (24 * 3600); // sec
17
18	const double AE = 1.0;
19	const double AU = 149597870.0; // Astronomical unit (km) (IAU 76)
20	const double SR = 696000.0; // Solar radius (km) (IAU 76)
21	const double TWOTHRD = 2.0 / 3.0;
22	const double XKMPER_WGS72 = 6378.135; // Earth equatorial radius - km (WGS '72)
23	const double F = 1.0 / 298.26; // Earth flattening (WGS '72)
24	const double GE = 398600.8; // Earth gravitational constant (WGS '72)
25	const double J2 = 1.0826158E-3; // J2 harmonic (WGS '72)
26	const double J3 = -2.53881E-6; // J3 harmonic (WGS '72)
27	const double J4 = -1.65597E-6; // J4 harmonic (WGS '72)
28	const double CK2 = J2 / 2.0;
29	const double CK4 = -3.0 * J4 / 8.0;
30	const double XJ3 = J3;
31	const double E6A = 1.0e-06;
32	const double QO = AE + 120.0 / XKMPER_WGS72;
33	const double S = AE + 78.0 / XKMPER_WGS72;
34	const double HR_PER_DAY = 24.0; // Hours per day (solar)
35	const double MIN_PER_DAY = 1440.0; // Minutes per day (solar)
36	const double SEC_PER_DAY = 86400.0; // Seconds per day (solar)
37	const double OMEGA_E = 1.00273790934; // earth rotation per sideral day
38	const double XKE = sqrt(3600.0 * GE / //sqrt(ge) ER^3/min^2
39	(XKMPER_WGS72 * XKMPER_WGS72 * XKMPER_WGS72));
40	const double QOMS2T = pow((QO - S), 4); //(QO - S)^4 ER^4
41
42	// Utility functions
43	double sqr (const double x);
44	double Fmod2p(const double arg);
45	double AcTan (const double sinx, double cosx);
46
47	double rad2deg(const double);
48	double deg2rad(const double);