SGP4.h

Go to the documentation of this file.

#ifndef _h_
#define _h_
/*     ----------------------------------------------------------------
*
*                                 SGP4.h
*
*    this file contains the sgp4 procedures for analytical propagation
*    of a satellite. the code was originally released in the 1980 and 1986
*    spacetrack papers. a detailed discussion of the theory and history
*    may be found in the 2006 aiaa paper by vallado, crawford, hujsak,
*    and kelso.
*
*    current :
*               7 dec 15  david vallado
*                           fix jd, jdfrac
*    changes :
*               3 nov 14  david vallado
*                           update to msvs2013 c++
*              30 Dec 11  david vallado
*                           consolidate updated code version
*              30 Aug 10  david vallado
*                           delete unused variables in initl
*                           replace pow inetger 2, 3 with multiplies for speed
*               3 Nov 08  david vallado
*                           put returns in for error codes
*              29 sep 08  david vallado
*                           fix atime for faster operation in dspace
*                           add operationmode for afspc (a) or improved (i)
*                           performance mode
*              20 apr 07  david vallado
*                           misc fixes for constants
*              11 aug 06  david vallado
*                           chg lyddane choice back to strn3, constants, misc doc
*              15 dec 05  david vallado
*                           misc fixes
*              26 jul 05  david vallado
*                           fixes for paper
*                           note that each fix is preceded by a
*                           comment with "sgp4fix" and an explanation of
*                           what was changed
*              10 aug 04  david vallado
*                           2nd printing baseline working
*              14 may 01  david vallado
*                           2nd edition baseline
*                     80  norad
*                           original baseline
*       ----------------------------------------------------------------      */

#pragma once

//using namespace System;

#include <math.h>
#include <stdio.h>
#include <string.h>
#include <iostream>
#include "vrfutil/vrfutilDefines.h"


#define SGP4Version  "SGP4 Version 2016-03-09"

typedef enum
{
  wgs72old,
  wgs72,
  wgs84
} gravconsttype;

typedef struct elsetrec
{
  long int  satnum;
  int       epochyr, epochtynumrev;
  int       error;
  char      operationmode;
  char      init, method;

  /* Near Earth */
  int    isimp;
  double aycof  , con41  , cc1    , cc4      , cc5    , d2      , d3   , d4    ,
         delmo  , eta    , argpdot, omgcof   , sinmao , t       , t2cof, t3cof ,
         t4cof  , t5cof  , x1mth2 , x7thm1   , mdot   , nodedot, xlcof , xmcof ,
         nodecf;

  /* Deep Space */
  int    irez;
  double d2201  , d2211  , d3210  , d3222    , d4410  , d4422   , d5220 , d5232 ,
         d5421  , d5433  , dedt   , del1     , del2   , del3    , didt  , dmdt  ,
         dnodt  , domdt  , e3     , ee2      , peo    , pgho    , pho   , pinco ,
         plo    , se2    , se3    , sgh2     , sgh3   , sgh4    , sh2   , sh3   ,
         si2    , si3    , sl2    , sl3      , sl4    , gsto    , xfact , xgh2  ,
         xgh3   , xgh4   , xh2    , xh3      , xi2    , xi3     , xl2   , xl3   ,
         xl4    , xlamo  , zmol   , zmos     , atime  , xli     , xni;

  double a, altp, alta, epochdays, jdsatepoch, jdsatepochF, nddot, ndot,
        bstar, rcse, inclo, nodeo, ecco, argpo, mo, no_kozai;
  char  classification, intldesg[11];
  int   ephtype;
  long  elnum    , revnum;
  double no_unkozai;
  double am     , em     , im     , Om       , om     , mm      , nm;
  double tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2;

  long dia_mm; 
  double period_sec; 
  unsigned char active; 
  unsigned char not_orbital; 
  double rcs_m2; 

} elsetrec;


namespace SGP4Funcs 
{

   // public class SGP4Class
   // {

   DT_DLL_vrfutil bool sgp4init
      (
      gravconsttype whichconst, char opsmode, const int satn, const double epoch,
      const double xbstar, const double xndot, const double xnddot, const double xecco, const double xargpo,
      const double xinclo, const double xmo, const double xno,
      const double xnodeo, elsetrec& satrec
      );

   DT_DLL_vrfutil bool sgp4
      (
      elsetrec& satrec, double tsince,
      double r[3], double v[3]
      );

   DT_DLL_vrfutil void getgravconst
      (
      gravconsttype whichconst,
      double& tumin,
      double& mu,
      double& radiusearthkm,
      double& xke,
      double& j2,
      double& j3,
      double& j4,
      double& j3oj2
      );

   DT_DLL_vrfutil void twoline2rv
      (
      char      longstr1[130], char longstr2[130],
      char      typerun, char typeinput, char opsmode,
      gravconsttype       whichconst,
      double& startmfe, double& stopmfe, double& deltamin,
      elsetrec& satrec
      );

   DT_DLL_vrfutil double  gstime
      (
      double jdut1
      );

   DT_DLL_vrfutil double  sgn
      (
      double x
      );

   DT_DLL_vrfutil double  mag
      (
      double x[3]
      );

   DT_DLL_vrfutil void    cross
      (
      double vec1[3], double vec2[3], double outvec[3]
      );

   DT_DLL_vrfutil double  dot
      (
      double x[3], double y[3]
      );

   DT_DLL_vrfutil double  angle
      (
      double vec1[3],
      double vec2[3]
      );

   DT_DLL_vrfutil void    newtonnu
      (
      double ecc, double nu,
      double& e0, double& m
      );

   DT_DLL_vrfutil double  asinh
      (
      double xval
      );

   DT_DLL_vrfutil void    rv2coe
      (
      double r[3], double v[3], const double mu,
      double& p, double& a, double& ecc, double& incl, double& omega, double& argp,
      double& nu, double& m, double& arglat, double& truelon, double& lonper
      );

   DT_DLL_vrfutil void    jday
      (
      int year, int mon, int day, int hr, int minute, double sec,
      double& jd, double& jdFrac
      );

   DT_DLL_vrfutil void    days2mdhms
      (
      int year, double days,
      int& mon, int& day, int& hr, int& minute, double& sec
      );

   DT_DLL_vrfutil void    invjday
      (
      double jd, double jdFrac,
      int& year, int& mon, int& day,
      int& hr, int& minute, double& sec
      );

   DT_DLL_vrfutil void teme2ecef(double rteme[3], double vteme[3], double jdut1, double recef[3], double vecef[3]);

}  

#endif