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	kusanagi	1.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);