KJ135
/
sgp4
mirror of http://172.16.0.12:5000/LAMPSARToolSoftware/sgp4.git
4
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Comments

feature/19
Daniel Warner 2011-04-01 02:11:20 +01:00
parent bcefdb6b6f
commit 4a7d07535e
2 changed files with 6 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
SGDP4.cpp
View File

@ -924,7 +924,7 @@ void SGDP4::DeepSpaceCalculateLunarSolarTerms(const double t, double& pe, double
static const double ZEL = 0.05490; static const double ZEL = 0.05490;
/* /*
* update solar terms * calculate solar terms
*/ */
double zm = d_zmos_ + ZNS * t; double zm = d_zmos_ + ZNS * t;
if (first_run_) if (first_run_)
@ -940,7 +940,7 @@ void SGDP4::DeepSpaceCalculateLunarSolarTerms(const double t, double& pe, double
const double shs = d_sh2_ * f2 + d_sh3_ * f3; const double shs = d_sh2_ * f2 + d_sh3_ * f3;
/* /*
* update lunar terms * calculate lunar terms
*/ */
zm = d_zmol_ + ZNL * t; zm = d_zmol_ + ZNL * t;
if (first_run_) if (first_run_)
@ -956,7 +956,7 @@ void SGDP4::DeepSpaceCalculateLunarSolarTerms(const double t, double& pe, double
const double shl = d_xh2_ * f2 + d_xh3_ * f3; const double shl = d_xh2_ * f2 + d_xh3_ * f3;
/* /*
* merge computed values * merge calculated values
*/ */
pe = ses + sel; pe = ses + sel;
pinc = sis + sil; pinc = sis + sil;
@ -1154,10 +1154,6 @@ void SGDP4::DeepSpaceCalcDotTerms(double& xndot, double& xnddt, double& xldot) c
} else { } else {
/*
* TODO: check ArgumentPerigee() and i_omgdot_
* are correct to use
*/
const double xomi = ArgumentPerigee() + i_omgdot_ * d_atime_; const double xomi = ArgumentPerigee() + i_omgdot_ * d_atime_;
const double x2omi = xomi + xomi; const double x2omi = xomi + xomi;
const double x2li = d_xli_ + d_xli_; const double x2li = d_xli_ + d_xli_;
@ -1185,7 +1181,7 @@ void SGDP4::DeepSpaceCalcDotTerms(double& xndot, double& xnddt, double& xldot) c
} }
xldot = d_xni_ + d_xfact_; xldot = d_xni_ + d_xfact_;
xnddt = xnddt * xldot; xnddt *= xldot;
} }
/* /*

4
main.cpp
View File

@ -144,7 +144,7 @@ void RunTle(Tle tle, double start, double end, double inc) {
void RunTest() { void RunTest() {
//#if 0 //#if 0
/* /*
# ------------------ Verification test cases ---------------------- # ------------------ Verification test cases ----------------------
@ -216,7 +216,7 @@ void RunTest() {
RunTle(Tle("11801", RunTle(Tle("11801",
"1 11801U 80230.29629788 .01431103 00000-0 14311-1 13", "1 11801U 80230.29629788 .01431103 00000-0 14311-1 13",
"2 11801 46.7916 230.4354 7318036 47.4722 10.4117 2.28537848 13"), 0.0, 1440.0, 360.0); "2 11801 46.7916 230.4354 7318036 47.4722 10.4117 2.28537848 13"), 0.0, 1440.0, 360.0);
//#endif //#endif
/* /*
# EUTELSAT 1-F1 (ECS1)## fig lyddane choice in GSFC at 2080 min # EUTELSAT 1-F1 (ECS1)## fig lyddane choice in GSFC at 2080 min