@@ -36,6 +36,10 @@ discard block |
||
36 | 36 | /* the year 1999 assuming that two-digit years in the range 00-56 */ |
37 | 37 | /* correspond to 2000-2056. Until the two-line element set format */ |
38 | 38 | /* is changed, it is only valid for dates through 2056 December 31. */ |
39 | + |
|
40 | + /** |
|
41 | + * @param double $epoch |
|
42 | + */ |
|
39 | 43 | public static function Julian_Date_of_Epoch($epoch) |
40 | 44 | { |
41 | 45 | $year = 0; |
@@ -54,6 +58,11 @@ discard block |
||
54 | 58 | } |
55 | 59 | |
56 | 60 | /* Equivalent to the C modf function */ |
61 | + |
|
62 | + /** |
|
63 | + * @param double $x |
|
64 | + * @param integer $ipart |
|
65 | + */ |
|
57 | 66 | public static function modf($x, &$ipart) { |
58 | 67 | $ipart = (int)$x; |
59 | 68 | return $x - $ipart; |
@@ -63,6 +72,10 @@ discard block |
||
63 | 72 | /* of Day 0.0 of {year}. This function is used to calculate the */ |
64 | 73 | /* Julian Date of any date by using Julian_Date_of_Year, DOY, */ |
65 | 74 | /* and Fraction_of_Day. */ |
75 | + |
|
76 | + /** |
|
77 | + * @param integer $year |
|
78 | + */ |
|
66 | 79 | public static function Julian_Date_of_Year($year) |
67 | 80 | { |
68 | 81 | /* Astronomical Formulae for Calculators, Jean Meeus, */ |
@@ -85,6 +98,10 @@ discard block |
||
85 | 98 | /* 1999, as described above. The function ThetaG_JD provides the */ |
86 | 99 | /* same calculation except that it is based on an input in the */ |
87 | 100 | /* form of a Julian Date. */ |
101 | + |
|
102 | + /** |
|
103 | + * @param double $epoch |
|
104 | + */ |
|
88 | 105 | public static function ThetaG($epoch, Predict_DeepArg $deep_arg) |
89 | 106 | { |
90 | 107 | /* Reference: The 1992 Astronomical Almanac, page B6. */ |
@@ -157,6 +174,10 @@ discard block |
||
157 | 174 | /* (approximately the same as UTC) and ET (now referred to as TDT).*/ |
158 | 175 | /* This function is based on a least squares fit of data from 1950 */ |
159 | 176 | /* to 1991 and will need to be updated periodically. */ |
177 | + |
|
178 | + /** |
|
179 | + * @param double $year |
|
180 | + */ |
|
160 | 181 | public static function Delta_ET($year) |
161 | 182 | { |
162 | 183 | /* Values determined using data from 1950-1991 in the 1990 |
@@ -204,7 +225,7 @@ discard block |
||
204 | 225 | * |
205 | 226 | * @param Predict_TLE $tle The TLE object |
206 | 227 | * |
207 | - * @return int |
|
228 | + * @return double |
|
208 | 229 | */ |
209 | 230 | public static function getEpochTimeStamp(Predict_TLE $tle) |
210 | 231 | { |
@@ -30,194 +30,194 @@ |
||
30 | 30 | */ |
31 | 31 | class Predict_Time |
32 | 32 | { |
33 | - /* The function Julian_Date_of_Epoch returns the Julian Date of */ |
|
34 | - /* an epoch specified in the format used in the NORAD two-line */ |
|
35 | - /* element sets. It has been modified to support dates beyond */ |
|
36 | - /* the year 1999 assuming that two-digit years in the range 00-56 */ |
|
37 | - /* correspond to 2000-2056. Until the two-line element set format */ |
|
38 | - /* is changed, it is only valid for dates through 2056 December 31. */ |
|
39 | - public static function Julian_Date_of_Epoch($epoch) |
|
40 | - { |
|
41 | - $year = 0; |
|
42 | - |
|
43 | - /* Modification to support Y2K */ |
|
44 | - /* Valid 1957 through 2056 */ |
|
45 | - $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
46 | - if ($year < 57) { |
|
47 | - $year = $year + 2000; |
|
48 | - } else { |
|
49 | - $year = $year + 1900; |
|
50 | - } |
|
51 | - /* End modification */ |
|
52 | - |
|
53 | - return self::Julian_Date_of_Year($year) + $day; |
|
54 | - } |
|
55 | - |
|
56 | - /* Equivalent to the C modf function */ |
|
57 | - public static function modf($x, &$ipart) { |
|
58 | - $ipart = (int)$x; |
|
59 | - return $x - $ipart; |
|
60 | - } |
|
61 | - |
|
62 | - /* The function Julian_Date_of_Year calculates the Julian Date */ |
|
63 | - /* of Day 0.0 of {year}. This function is used to calculate the */ |
|
64 | - /* Julian Date of any date by using Julian_Date_of_Year, DOY, */ |
|
65 | - /* and Fraction_of_Day. */ |
|
66 | - public static function Julian_Date_of_Year($year) |
|
67 | - { |
|
68 | - /* Astronomical Formulae for Calculators, Jean Meeus, */ |
|
69 | - /* pages 23-25. Calculate Julian Date of 0.0 Jan year */ |
|
70 | - $year = $year - 1; |
|
71 | - $i = (int) ($year / 100); |
|
72 | - $A = $i; |
|
73 | - $i = (int) ($A / 4); |
|
74 | - $B = (int) (2 - $A + $i); |
|
75 | - $i = (int) (365.25 * $year); |
|
76 | - $i += (int) (30.6001 * 14); |
|
77 | - $jdoy = $i + 1720994.5 + $B; |
|
78 | - |
|
79 | - return $jdoy; |
|
80 | - } |
|
81 | - |
|
82 | - /* The function ThetaG calculates the Greenwich Mean Sidereal Time */ |
|
83 | - /* for an epoch specified in the format used in the NORAD two-line */ |
|
84 | - /* element sets. It has now been adapted for dates beyond the year */ |
|
85 | - /* 1999, as described above. The function ThetaG_JD provides the */ |
|
86 | - /* same calculation except that it is based on an input in the */ |
|
87 | - /* form of a Julian Date. */ |
|
88 | - public static function ThetaG($epoch, Predict_DeepArg $deep_arg) |
|
89 | - { |
|
90 | - /* Reference: The 1992 Astronomical Almanac, page B6. */ |
|
91 | - // double year,day,UT,jd,TU,GMST,_ThetaG; |
|
92 | - |
|
93 | - /* Modification to support Y2K */ |
|
94 | - /* Valid 1957 through 2056 */ |
|
95 | - $year = 0; |
|
96 | - $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
97 | - |
|
98 | - if ($year < 57) { |
|
99 | - $year += 2000; |
|
100 | - } else { |
|
101 | - $year += 1900; |
|
102 | - } |
|
103 | - /* End modification */ |
|
104 | - |
|
105 | - $UT = fmod($day, $day); |
|
106 | - $jd = self::Julian_Date_of_Year($year) + $day; |
|
107 | - $TU = ($jd - 2451545.0) / 36525; |
|
108 | - $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
109 | - $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
110 | - $deep_arg->ds50 = $jd - 2433281.5 + $UT; |
|
111 | - |
|
112 | - return Predict_Math::FMod2p(6.3003880987 * $deep_arg->ds50 + 1.72944494); |
|
113 | - } |
|
114 | - |
|
115 | - /* See the ThetaG doc block above */ |
|
116 | - public static function ThetaG_JD($jd) |
|
117 | - { |
|
118 | - /* Reference: The 1992 Astronomical Almanac, page B6. */ |
|
119 | - $UT = Predict_Math::Frac($jd + 0.5); |
|
120 | - $jd = $jd - $UT; |
|
121 | - $TU = ($jd - 2451545.0) / 36525; |
|
122 | - $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
123 | - $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
124 | - |
|
125 | - return Predict::twopi * $GMST / Predict::secday; |
|
126 | - } |
|
127 | - |
|
128 | - /** |
|
129 | - * Read the system clock and return the current Julian day. From phpPredict |
|
130 | - * |
|
131 | - * @return float |
|
132 | - */ |
|
133 | - public static function get_current_daynum() { |
|
134 | - // Gets the current decimal day number from microtime |
|
135 | - |
|
136 | - list($usec, $sec) = explode(' ', microtime()); |
|
137 | - return self::unix2daynum($sec, $usec); |
|
138 | - } |
|
139 | - |
|
140 | - /** |
|
141 | - * Converts a standard unix timestamp and optional |
|
142 | - * milliseconds to a daynum |
|
143 | - * |
|
144 | - * @param int $sec Seconds from the unix epoch |
|
145 | - * @param int $usec Optional milliseconds |
|
146 | - * |
|
147 | - * @return float |
|
148 | - */ |
|
149 | - public static function unix2daynum($sec, $usec = 0) |
|
150 | - { |
|
151 | - $time = ((($sec + $usec) / 86400.0) - 3651.0); |
|
152 | - return $time + 2444238.5; |
|
153 | - } |
|
154 | - |
|
155 | - /* The function Delta_ET has been added to allow calculations on */ |
|
156 | - /* the position of the sun. It provides the difference between UT */ |
|
157 | - /* (approximately the same as UTC) and ET (now referred to as TDT).*/ |
|
158 | - /* This function is based on a least squares fit of data from 1950 */ |
|
159 | - /* to 1991 and will need to be updated periodically. */ |
|
160 | - public static function Delta_ET($year) |
|
161 | - { |
|
162 | - /* Values determined using data from 1950-1991 in the 1990 |
|
33 | + /* The function Julian_Date_of_Epoch returns the Julian Date of */ |
|
34 | + /* an epoch specified in the format used in the NORAD two-line */ |
|
35 | + /* element sets. It has been modified to support dates beyond */ |
|
36 | + /* the year 1999 assuming that two-digit years in the range 00-56 */ |
|
37 | + /* correspond to 2000-2056. Until the two-line element set format */ |
|
38 | + /* is changed, it is only valid for dates through 2056 December 31. */ |
|
39 | + public static function Julian_Date_of_Epoch($epoch) |
|
40 | + { |
|
41 | + $year = 0; |
|
42 | + |
|
43 | + /* Modification to support Y2K */ |
|
44 | + /* Valid 1957 through 2056 */ |
|
45 | + $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
46 | + if ($year < 57) { |
|
47 | + $year = $year + 2000; |
|
48 | + } else { |
|
49 | + $year = $year + 1900; |
|
50 | + } |
|
51 | + /* End modification */ |
|
52 | + |
|
53 | + return self::Julian_Date_of_Year($year) + $day; |
|
54 | + } |
|
55 | + |
|
56 | + /* Equivalent to the C modf function */ |
|
57 | + public static function modf($x, &$ipart) { |
|
58 | + $ipart = (int)$x; |
|
59 | + return $x - $ipart; |
|
60 | + } |
|
61 | + |
|
62 | + /* The function Julian_Date_of_Year calculates the Julian Date */ |
|
63 | + /* of Day 0.0 of {year}. This function is used to calculate the */ |
|
64 | + /* Julian Date of any date by using Julian_Date_of_Year, DOY, */ |
|
65 | + /* and Fraction_of_Day. */ |
|
66 | + public static function Julian_Date_of_Year($year) |
|
67 | + { |
|
68 | + /* Astronomical Formulae for Calculators, Jean Meeus, */ |
|
69 | + /* pages 23-25. Calculate Julian Date of 0.0 Jan year */ |
|
70 | + $year = $year - 1; |
|
71 | + $i = (int) ($year / 100); |
|
72 | + $A = $i; |
|
73 | + $i = (int) ($A / 4); |
|
74 | + $B = (int) (2 - $A + $i); |
|
75 | + $i = (int) (365.25 * $year); |
|
76 | + $i += (int) (30.6001 * 14); |
|
77 | + $jdoy = $i + 1720994.5 + $B; |
|
78 | + |
|
79 | + return $jdoy; |
|
80 | + } |
|
81 | + |
|
82 | + /* The function ThetaG calculates the Greenwich Mean Sidereal Time */ |
|
83 | + /* for an epoch specified in the format used in the NORAD two-line */ |
|
84 | + /* element sets. It has now been adapted for dates beyond the year */ |
|
85 | + /* 1999, as described above. The function ThetaG_JD provides the */ |
|
86 | + /* same calculation except that it is based on an input in the */ |
|
87 | + /* form of a Julian Date. */ |
|
88 | + public static function ThetaG($epoch, Predict_DeepArg $deep_arg) |
|
89 | + { |
|
90 | + /* Reference: The 1992 Astronomical Almanac, page B6. */ |
|
91 | + // double year,day,UT,jd,TU,GMST,_ThetaG; |
|
92 | + |
|
93 | + /* Modification to support Y2K */ |
|
94 | + /* Valid 1957 through 2056 */ |
|
95 | + $year = 0; |
|
96 | + $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
97 | + |
|
98 | + if ($year < 57) { |
|
99 | + $year += 2000; |
|
100 | + } else { |
|
101 | + $year += 1900; |
|
102 | + } |
|
103 | + /* End modification */ |
|
104 | + |
|
105 | + $UT = fmod($day, $day); |
|
106 | + $jd = self::Julian_Date_of_Year($year) + $day; |
|
107 | + $TU = ($jd - 2451545.0) / 36525; |
|
108 | + $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
109 | + $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
110 | + $deep_arg->ds50 = $jd - 2433281.5 + $UT; |
|
111 | + |
|
112 | + return Predict_Math::FMod2p(6.3003880987 * $deep_arg->ds50 + 1.72944494); |
|
113 | + } |
|
114 | + |
|
115 | + /* See the ThetaG doc block above */ |
|
116 | + public static function ThetaG_JD($jd) |
|
117 | + { |
|
118 | + /* Reference: The 1992 Astronomical Almanac, page B6. */ |
|
119 | + $UT = Predict_Math::Frac($jd + 0.5); |
|
120 | + $jd = $jd - $UT; |
|
121 | + $TU = ($jd - 2451545.0) / 36525; |
|
122 | + $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
123 | + $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
124 | + |
|
125 | + return Predict::twopi * $GMST / Predict::secday; |
|
126 | + } |
|
127 | + |
|
128 | + /** |
|
129 | + * Read the system clock and return the current Julian day. From phpPredict |
|
130 | + * |
|
131 | + * @return float |
|
132 | + */ |
|
133 | + public static function get_current_daynum() { |
|
134 | + // Gets the current decimal day number from microtime |
|
135 | + |
|
136 | + list($usec, $sec) = explode(' ', microtime()); |
|
137 | + return self::unix2daynum($sec, $usec); |
|
138 | + } |
|
139 | + |
|
140 | + /** |
|
141 | + * Converts a standard unix timestamp and optional |
|
142 | + * milliseconds to a daynum |
|
143 | + * |
|
144 | + * @param int $sec Seconds from the unix epoch |
|
145 | + * @param int $usec Optional milliseconds |
|
146 | + * |
|
147 | + * @return float |
|
148 | + */ |
|
149 | + public static function unix2daynum($sec, $usec = 0) |
|
150 | + { |
|
151 | + $time = ((($sec + $usec) / 86400.0) - 3651.0); |
|
152 | + return $time + 2444238.5; |
|
153 | + } |
|
154 | + |
|
155 | + /* The function Delta_ET has been added to allow calculations on */ |
|
156 | + /* the position of the sun. It provides the difference between UT */ |
|
157 | + /* (approximately the same as UTC) and ET (now referred to as TDT).*/ |
|
158 | + /* This function is based on a least squares fit of data from 1950 */ |
|
159 | + /* to 1991 and will need to be updated periodically. */ |
|
160 | + public static function Delta_ET($year) |
|
161 | + { |
|
162 | + /* Values determined using data from 1950-1991 in the 1990 |
|
163 | 163 | Astronomical Almanac. See DELTA_ET.WQ1 for details. */ |
164 | 164 | |
165 | - $delta_et = 26.465 + 0.747622 * ($year - 1950) + |
|
166 | - 1.886913 * sin(Predict::twopi * ($year - 1975) / 33); |
|
167 | - |
|
168 | - return $delta_et; |
|
169 | - } |
|
170 | - |
|
171 | - /** |
|
172 | - * Converts a daynum to a unix timestamp. From phpPredict. |
|
173 | - * |
|
174 | - * @param float $dn Julian Daynum |
|
175 | - * |
|
176 | - * @return float |
|
177 | - */ |
|
178 | - public static function daynum2unix($dn) { |
|
179 | - // Converts a daynum to a UNIX timestamp |
|
180 | - |
|
181 | - return (86400.0 * ($dn - 2444238.5 + 3651.0)); |
|
182 | - } |
|
183 | - |
|
184 | - /** |
|
185 | - * Converts a daynum to a readable time format. |
|
186 | - * |
|
187 | - * @param float $dn The julian date |
|
188 | - * @param string $zone The zone string, defaults to America/Los_Angeles |
|
189 | - * @param string $format The date() function's format string. Defaults to m-d-Y H:i:s |
|
190 | - * |
|
191 | - * @return string |
|
192 | - */ |
|
193 | - public static function daynum2readable($dn, $zone = 'America/Los_Angeles', $format = 'm-d-Y H:i:s') |
|
194 | - { |
|
195 | - $unix = self::daynum2unix($dn); |
|
196 | - $date = new DateTime("@" . round($unix)); |
|
197 | - $dateTimezone = new DateTimezone($zone); |
|
198 | - $date->setTimezone($dateTimezone); |
|
199 | - return $date->format($format); |
|
200 | - } |
|
201 | - |
|
202 | - /** |
|
203 | - * Returns the unix timestamp of a TLE's epoch |
|
204 | - * |
|
205 | - * @param Predict_TLE $tle The TLE object |
|
206 | - * |
|
207 | - * @return int |
|
208 | - */ |
|
209 | - public static function getEpochTimeStamp(Predict_TLE $tle) |
|
210 | - { |
|
211 | - $year = $tle->epoch_year; |
|
212 | - $day = $tle->epoch_day; |
|
213 | - $sec = round(86400 * $tle->epoch_fod); |
|
214 | - |
|
215 | - $zone = new DateTimeZone('GMT'); |
|
216 | - $date = new DateTime(); |
|
217 | - $date->setTimezone($zone); |
|
218 | - $date->setDate($year, 1, 1); |
|
219 | - $date->setTime(0, 0, 0); |
|
220 | - |
|
221 | - return $date->format('U') + (86400 * $day) + $sec - 86400; |
|
222 | - } |
|
165 | + $delta_et = 26.465 + 0.747622 * ($year - 1950) + |
|
166 | + 1.886913 * sin(Predict::twopi * ($year - 1975) / 33); |
|
167 | + |
|
168 | + return $delta_et; |
|
169 | + } |
|
170 | + |
|
171 | + /** |
|
172 | + * Converts a daynum to a unix timestamp. From phpPredict. |
|
173 | + * |
|
174 | + * @param float $dn Julian Daynum |
|
175 | + * |
|
176 | + * @return float |
|
177 | + */ |
|
178 | + public static function daynum2unix($dn) { |
|
179 | + // Converts a daynum to a UNIX timestamp |
|
180 | + |
|
181 | + return (86400.0 * ($dn - 2444238.5 + 3651.0)); |
|
182 | + } |
|
183 | + |
|
184 | + /** |
|
185 | + * Converts a daynum to a readable time format. |
|
186 | + * |
|
187 | + * @param float $dn The julian date |
|
188 | + * @param string $zone The zone string, defaults to America/Los_Angeles |
|
189 | + * @param string $format The date() function's format string. Defaults to m-d-Y H:i:s |
|
190 | + * |
|
191 | + * @return string |
|
192 | + */ |
|
193 | + public static function daynum2readable($dn, $zone = 'America/Los_Angeles', $format = 'm-d-Y H:i:s') |
|
194 | + { |
|
195 | + $unix = self::daynum2unix($dn); |
|
196 | + $date = new DateTime("@" . round($unix)); |
|
197 | + $dateTimezone = new DateTimezone($zone); |
|
198 | + $date->setTimezone($dateTimezone); |
|
199 | + return $date->format($format); |
|
200 | + } |
|
201 | + |
|
202 | + /** |
|
203 | + * Returns the unix timestamp of a TLE's epoch |
|
204 | + * |
|
205 | + * @param Predict_TLE $tle The TLE object |
|
206 | + * |
|
207 | + * @return int |
|
208 | + */ |
|
209 | + public static function getEpochTimeStamp(Predict_TLE $tle) |
|
210 | + { |
|
211 | + $year = $tle->epoch_year; |
|
212 | + $day = $tle->epoch_day; |
|
213 | + $sec = round(86400 * $tle->epoch_fod); |
|
214 | + |
|
215 | + $zone = new DateTimeZone('GMT'); |
|
216 | + $date = new DateTime(); |
|
217 | + $date->setTimezone($zone); |
|
218 | + $date->setDate($year, 1, 1); |
|
219 | + $date->setTime(0, 0, 0); |
|
220 | + |
|
221 | + return $date->format('U') + (86400 * $day) + $sec - 86400; |
|
222 | + } |
|
223 | 223 | } |
@@ -42,7 +42,7 @@ discard block |
||
42 | 42 | |
43 | 43 | /* Modification to support Y2K */ |
44 | 44 | /* Valid 1957 through 2056 */ |
45 | - $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
45 | + $day = self::modf($epoch*1E-3, $year)*1E3; |
|
46 | 46 | if ($year < 57) { |
47 | 47 | $year = $year + 2000; |
48 | 48 | } else { |
@@ -55,7 +55,7 @@ discard block |
||
55 | 55 | |
56 | 56 | /* Equivalent to the C modf function */ |
57 | 57 | public static function modf($x, &$ipart) { |
58 | - $ipart = (int)$x; |
|
58 | + $ipart = (int) $x; |
|
59 | 59 | return $x - $ipart; |
60 | 60 | } |
61 | 61 | |
@@ -68,12 +68,12 @@ discard block |
||
68 | 68 | /* Astronomical Formulae for Calculators, Jean Meeus, */ |
69 | 69 | /* pages 23-25. Calculate Julian Date of 0.0 Jan year */ |
70 | 70 | $year = $year - 1; |
71 | - $i = (int) ($year / 100); |
|
71 | + $i = (int) ($year/100); |
|
72 | 72 | $A = $i; |
73 | - $i = (int) ($A / 4); |
|
73 | + $i = (int) ($A/4); |
|
74 | 74 | $B = (int) (2 - $A + $i); |
75 | - $i = (int) (365.25 * $year); |
|
76 | - $i += (int) (30.6001 * 14); |
|
75 | + $i = (int) (365.25*$year); |
|
76 | + $i += (int) (30.6001*14); |
|
77 | 77 | $jdoy = $i + 1720994.5 + $B; |
78 | 78 | |
79 | 79 | return $jdoy; |
@@ -93,7 +93,7 @@ discard block |
||
93 | 93 | /* Modification to support Y2K */ |
94 | 94 | /* Valid 1957 through 2056 */ |
95 | 95 | $year = 0; |
96 | - $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
96 | + $day = self::modf($epoch*1E-3, $year)*1E3; |
|
97 | 97 | |
98 | 98 | if ($year < 57) { |
99 | 99 | $year += 2000; |
@@ -104,12 +104,12 @@ discard block |
||
104 | 104 | |
105 | 105 | $UT = fmod($day, $day); |
106 | 106 | $jd = self::Julian_Date_of_Year($year) + $day; |
107 | - $TU = ($jd - 2451545.0) / 36525; |
|
108 | - $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
109 | - $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
107 | + $TU = ($jd - 2451545.0)/36525; |
|
108 | + $GMST = 24110.54841 + $TU*(8640184.812866 + $TU*(0.093104 - $TU*6.2E-6)); |
|
109 | + $GMST = Predict_Math::Modulus($GMST + Predict::secday*Predict::omega_E*$UT, Predict::secday); |
|
110 | 110 | $deep_arg->ds50 = $jd - 2433281.5 + $UT; |
111 | 111 | |
112 | - return Predict_Math::FMod2p(6.3003880987 * $deep_arg->ds50 + 1.72944494); |
|
112 | + return Predict_Math::FMod2p(6.3003880987*$deep_arg->ds50 + 1.72944494); |
|
113 | 113 | } |
114 | 114 | |
115 | 115 | /* See the ThetaG doc block above */ |
@@ -118,11 +118,11 @@ discard block |
||
118 | 118 | /* Reference: The 1992 Astronomical Almanac, page B6. */ |
119 | 119 | $UT = Predict_Math::Frac($jd + 0.5); |
120 | 120 | $jd = $jd - $UT; |
121 | - $TU = ($jd - 2451545.0) / 36525; |
|
122 | - $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
123 | - $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
121 | + $TU = ($jd - 2451545.0)/36525; |
|
122 | + $GMST = 24110.54841 + $TU*(8640184.812866 + $TU*(0.093104 - $TU*6.2E-6)); |
|
123 | + $GMST = Predict_Math::Modulus($GMST + Predict::secday*Predict::omega_E*$UT, Predict::secday); |
|
124 | 124 | |
125 | - return Predict::twopi * $GMST / Predict::secday; |
|
125 | + return Predict::twopi*$GMST/Predict::secday; |
|
126 | 126 | } |
127 | 127 | |
128 | 128 | /** |
@@ -148,7 +148,7 @@ discard block |
||
148 | 148 | */ |
149 | 149 | public static function unix2daynum($sec, $usec = 0) |
150 | 150 | { |
151 | - $time = ((($sec + $usec) / 86400.0) - 3651.0); |
|
151 | + $time = ((($sec + $usec)/86400.0) - 3651.0); |
|
152 | 152 | return $time + 2444238.5; |
153 | 153 | } |
154 | 154 | |
@@ -162,8 +162,8 @@ discard block |
||
162 | 162 | /* Values determined using data from 1950-1991 in the 1990 |
163 | 163 | Astronomical Almanac. See DELTA_ET.WQ1 for details. */ |
164 | 164 | |
165 | - $delta_et = 26.465 + 0.747622 * ($year - 1950) + |
|
166 | - 1.886913 * sin(Predict::twopi * ($year - 1975) / 33); |
|
165 | + $delta_et = 26.465 + 0.747622*($year - 1950) + |
|
166 | + 1.886913*sin(Predict::twopi*($year - 1975)/33); |
|
167 | 167 | |
168 | 168 | return $delta_et; |
169 | 169 | } |
@@ -178,7 +178,7 @@ discard block |
||
178 | 178 | public static function daynum2unix($dn) { |
179 | 179 | // Converts a daynum to a UNIX timestamp |
180 | 180 | |
181 | - return (86400.0 * ($dn - 2444238.5 + 3651.0)); |
|
181 | + return (86400.0*($dn - 2444238.5 + 3651.0)); |
|
182 | 182 | } |
183 | 183 | |
184 | 184 | /** |
@@ -193,7 +193,7 @@ discard block |
||
193 | 193 | public static function daynum2readable($dn, $zone = 'America/Los_Angeles', $format = 'm-d-Y H:i:s') |
194 | 194 | { |
195 | 195 | $unix = self::daynum2unix($dn); |
196 | - $date = new DateTime("@" . round($unix)); |
|
196 | + $date = new DateTime("@".round($unix)); |
|
197 | 197 | $dateTimezone = new DateTimezone($zone); |
198 | 198 | $date->setTimezone($dateTimezone); |
199 | 199 | return $date->format($format); |
@@ -210,7 +210,7 @@ discard block |
||
210 | 210 | { |
211 | 211 | $year = $tle->epoch_year; |
212 | 212 | $day = $tle->epoch_day; |
213 | - $sec = round(86400 * $tle->epoch_fod); |
|
213 | + $sec = round(86400*$tle->epoch_fod); |
|
214 | 214 | |
215 | 215 | $zone = new DateTimeZone('GMT'); |
216 | 216 | $date = new DateTime(); |
@@ -218,6 +218,6 @@ discard block |
||
218 | 218 | $date->setDate($year, 1, 1); |
219 | 219 | $date->setTime(0, 0, 0); |
220 | 220 | |
221 | - return $date->format('U') + (86400 * $day) + $sec - 86400; |
|
221 | + return $date->format('U') + (86400*$day) + $sec - 86400; |
|
222 | 222 | } |
223 | 223 | } |
@@ -6,30 +6,30 @@ |
||
6 | 6 | /* Common arguments between deep-space functions */ |
7 | 7 | class Predict_DeepArg |
8 | 8 | { |
9 | - /* Used by dpinit part of Deep() */ |
|
10 | - public $eosq; |
|
11 | - public $sinio; |
|
12 | - public $cosio; |
|
13 | - public $betao; |
|
14 | - public $aodp; |
|
15 | - public $theta2; |
|
16 | - public $sing; |
|
17 | - public $cosg; |
|
18 | - public $betao2; |
|
19 | - public $xmdot; |
|
20 | - public $omgdot; |
|
21 | - public $xnodot; |
|
22 | - public $xnodp; |
|
9 | + /* Used by dpinit part of Deep() */ |
|
10 | + public $eosq; |
|
11 | + public $sinio; |
|
12 | + public $cosio; |
|
13 | + public $betao; |
|
14 | + public $aodp; |
|
15 | + public $theta2; |
|
16 | + public $sing; |
|
17 | + public $cosg; |
|
18 | + public $betao2; |
|
19 | + public $xmdot; |
|
20 | + public $omgdot; |
|
21 | + public $xnodot; |
|
22 | + public $xnodp; |
|
23 | 23 | |
24 | - /* Used by dpsec and dpper parts of Deep() */ |
|
25 | - public $xll; |
|
26 | - public $omgadf; |
|
27 | - public $xnode; |
|
28 | - public $em; |
|
29 | - public $xinc; |
|
30 | - public $xn; |
|
31 | - public $t; |
|
24 | + /* Used by dpsec and dpper parts of Deep() */ |
|
25 | + public $xll; |
|
26 | + public $omgadf; |
|
27 | + public $xnode; |
|
28 | + public $em; |
|
29 | + public $xinc; |
|
30 | + public $xn; |
|
31 | + public $t; |
|
32 | 32 | |
33 | - /* Used by thetg and Deep() */ |
|
34 | - public $ds50; |
|
33 | + /* Used by thetg and Deep() */ |
|
34 | + public $ds50; |
|
35 | 35 | } |
@@ -3,16 +3,16 @@ |
||
3 | 3 | /** Brief satellite pass info. */ |
4 | 4 | class Predict_Pass |
5 | 5 | { |
6 | - public $satname; /*!< satellite name */ |
|
7 | - public $aos; /*!< AOS time in "jul_utc" */ |
|
8 | - public $tca; /*!< TCA time in "jul_utc" */ |
|
9 | - public $los; /*!< LOS time in "jul_utc" */ |
|
10 | - public $max_el; /*!< Maximum elevation during pass */ |
|
11 | - public $aos_az; /*!< Azimuth at AOS */ |
|
12 | - public $los_az; /*!< Azimuth at LOS */ |
|
13 | - public $orbit; /*!< Orbit number */ |
|
14 | - public $maxel_az; /*!< Azimuth at maximum elevation */ |
|
15 | - public $vis; /*!< Visibility string, e.g. VSE, -S-, V-- */ |
|
16 | - public $details = array(); /*!< List of pass_detail_t entries */ |
|
17 | - public $max_apparent_magnitude = null; /* maximum apparent magnitude, experimental */ |
|
6 | + public $satname; /*!< satellite name */ |
|
7 | + public $aos; /*!< AOS time in "jul_utc" */ |
|
8 | + public $tca; /*!< TCA time in "jul_utc" */ |
|
9 | + public $los; /*!< LOS time in "jul_utc" */ |
|
10 | + public $max_el; /*!< Maximum elevation during pass */ |
|
11 | + public $aos_az; /*!< Azimuth at AOS */ |
|
12 | + public $los_az; /*!< Azimuth at LOS */ |
|
13 | + public $orbit; /*!< Orbit number */ |
|
14 | + public $maxel_az; /*!< Azimuth at maximum elevation */ |
|
15 | + public $vis; /*!< Visibility string, e.g. VSE, -S-, V-- */ |
|
16 | + public $details = array(); /*!< List of pass_detail_t entries */ |
|
17 | + public $max_apparent_magnitude = null; /* maximum apparent magnitude, experimental */ |
|
18 | 18 | } |
@@ -3,16 +3,16 @@ |
||
3 | 3 | /** Brief satellite pass info. */ |
4 | 4 | class Predict_Pass |
5 | 5 | { |
6 | - public $satname; /*!< satellite name */ |
|
7 | - public $aos; /*!< AOS time in "jul_utc" */ |
|
8 | - public $tca; /*!< TCA time in "jul_utc" */ |
|
9 | - public $los; /*!< LOS time in "jul_utc" */ |
|
10 | - public $max_el; /*!< Maximum elevation during pass */ |
|
11 | - public $aos_az; /*!< Azimuth at AOS */ |
|
12 | - public $los_az; /*!< Azimuth at LOS */ |
|
13 | - public $orbit; /*!< Orbit number */ |
|
6 | + public $satname; /*!< satellite name */ |
|
7 | + public $aos; /*!< AOS time in "jul_utc" */ |
|
8 | + public $tca; /*!< TCA time in "jul_utc" */ |
|
9 | + public $los; /*!< LOS time in "jul_utc" */ |
|
10 | + public $max_el; /*!< Maximum elevation during pass */ |
|
11 | + public $aos_az; /*!< Azimuth at AOS */ |
|
12 | + public $los_az; /*!< Azimuth at LOS */ |
|
13 | + public $orbit; /*!< Orbit number */ |
|
14 | 14 | public $maxel_az; /*!< Azimuth at maximum elevation */ |
15 | - public $vis; /*!< Visibility string, e.g. VSE, -S-, V-- */ |
|
16 | - public $details = array(); /*!< List of pass_detail_t entries */ |
|
15 | + public $vis; /*!< Visibility string, e.g. VSE, -S-, V-- */ |
|
16 | + public $details = array(); /*!< List of pass_detail_t entries */ |
|
17 | 17 | public $max_apparent_magnitude = null; /* maximum apparent magnitude, experimental */ |
18 | 18 | } |
@@ -6,78 +6,78 @@ |
||
6 | 6 | /* static data for DEEP */ |
7 | 7 | class Predict_DeepStatic |
8 | 8 | { |
9 | - public $thgr; |
|
10 | - public $xnq; |
|
11 | - public $xqncl; |
|
12 | - public $omegaq; |
|
13 | - public $zmol; |
|
14 | - public $zmos; |
|
15 | - public $savtsn; |
|
16 | - public $ee2; |
|
17 | - public $e3; |
|
18 | - public $xi2; |
|
19 | - public $xl2; |
|
20 | - public $xl3; |
|
21 | - public $xl4; |
|
22 | - public $xgh2; |
|
23 | - public $xgh3; |
|
24 | - public $xgh4; |
|
25 | - public $xh2; |
|
26 | - public $xh3; |
|
27 | - public $sse; |
|
28 | - public $ssi; |
|
29 | - public $ssg; |
|
30 | - public $xi3; |
|
31 | - public $se2; |
|
32 | - public $si2; |
|
33 | - public $sl2; |
|
34 | - public $sgh2; |
|
35 | - public $sh2; |
|
36 | - public $se3; |
|
37 | - public $si3; |
|
38 | - public $sl3; |
|
39 | - public $sgh3; |
|
40 | - public $sh3; |
|
41 | - public $sl4; |
|
42 | - public $sgh4; |
|
43 | - public $ssl; |
|
44 | - public $ssh; |
|
45 | - public $d3210; |
|
46 | - public $d3222; |
|
47 | - public $d4410; |
|
48 | - public $d4422; |
|
49 | - public $d5220; |
|
50 | - public $d5232; |
|
51 | - public $d5421; |
|
52 | - public $d5433; |
|
53 | - public $del1; |
|
54 | - public $del2; |
|
55 | - public $del3; |
|
56 | - public $fasx2; |
|
57 | - public $fasx4; |
|
58 | - public $fasx6; |
|
59 | - public $xlamo; |
|
60 | - public $xfact; |
|
61 | - public $xni; |
|
62 | - public $atime; |
|
63 | - public $stepp; |
|
64 | - public $stepn; |
|
65 | - public $step2; |
|
66 | - public $preep; |
|
67 | - public $pl; |
|
68 | - public $sghs; |
|
69 | - public $xli; |
|
70 | - public $d2201; |
|
71 | - public $d2211; |
|
72 | - public $sghl; |
|
73 | - public $sh1; |
|
74 | - public $pinc; |
|
75 | - public $pe; |
|
76 | - public $shs; |
|
77 | - public $zsingl; |
|
78 | - public $zcosgl; |
|
79 | - public $zsinhl; |
|
80 | - public $zcoshl; |
|
81 | - public $zsinil; |
|
82 | - public $zcosil; |
|
9 | + public $thgr; |
|
10 | + public $xnq; |
|
11 | + public $xqncl; |
|
12 | + public $omegaq; |
|
13 | + public $zmol; |
|
14 | + public $zmos; |
|
15 | + public $savtsn; |
|
16 | + public $ee2; |
|
17 | + public $e3; |
|
18 | + public $xi2; |
|
19 | + public $xl2; |
|
20 | + public $xl3; |
|
21 | + public $xl4; |
|
22 | + public $xgh2; |
|
23 | + public $xgh3; |
|
24 | + public $xgh4; |
|
25 | + public $xh2; |
|
26 | + public $xh3; |
|
27 | + public $sse; |
|
28 | + public $ssi; |
|
29 | + public $ssg; |
|
30 | + public $xi3; |
|
31 | + public $se2; |
|
32 | + public $si2; |
|
33 | + public $sl2; |
|
34 | + public $sgh2; |
|
35 | + public $sh2; |
|
36 | + public $se3; |
|
37 | + public $si3; |
|
38 | + public $sl3; |
|
39 | + public $sgh3; |
|
40 | + public $sh3; |
|
41 | + public $sl4; |
|
42 | + public $sgh4; |
|
43 | + public $ssl; |
|
44 | + public $ssh; |
|
45 | + public $d3210; |
|
46 | + public $d3222; |
|
47 | + public $d4410; |
|
48 | + public $d4422; |
|
49 | + public $d5220; |
|
50 | + public $d5232; |
|
51 | + public $d5421; |
|
52 | + public $d5433; |
|
53 | + public $del1; |
|
54 | + public $del2; |
|
55 | + public $del3; |
|
56 | + public $fasx2; |
|
57 | + public $fasx4; |
|
58 | + public $fasx6; |
|
59 | + public $xlamo; |
|
60 | + public $xfact; |
|
61 | + public $xni; |
|
62 | + public $atime; |
|
63 | + public $stepp; |
|
64 | + public $stepn; |
|
65 | + public $step2; |
|
66 | + public $preep; |
|
67 | + public $pl; |
|
68 | + public $sghs; |
|
69 | + public $xli; |
|
70 | + public $d2201; |
|
71 | + public $d2211; |
|
72 | + public $sghl; |
|
73 | + public $sh1; |
|
74 | + public $pinc; |
|
75 | + public $pe; |
|
76 | + public $shs; |
|
77 | + public $zsingl; |
|
78 | + public $zcosgl; |
|
79 | + public $zsinhl; |
|
80 | + public $zcoshl; |
|
81 | + public $zsinil; |
|
82 | + public $zcosil; |
|
83 | 83 | } |
@@ -12,26 +12,26 @@ |
||
12 | 12 | */ |
13 | 13 | class Predict_PassDetail |
14 | 14 | { |
15 | - public $time; /*!< time in "jul_utc" */ |
|
16 | - public $pos; /*!< Raw unprocessed position at time */ |
|
17 | - public $vel; /*!< Raw unprocessed velocity at time */ |
|
18 | - public $velo; |
|
19 | - public $az; |
|
20 | - public $el; |
|
21 | - public $range; |
|
22 | - public $range_rate; |
|
23 | - public $lat; |
|
24 | - public $lon; |
|
25 | - public $alt; |
|
26 | - public $ma; |
|
27 | - public $phase; |
|
28 | - public $footprint; |
|
29 | - public $vis; |
|
30 | - public $orbit; |
|
15 | + public $time; /*!< time in "jul_utc" */ |
|
16 | + public $pos; /*!< Raw unprocessed position at time */ |
|
17 | + public $vel; /*!< Raw unprocessed velocity at time */ |
|
18 | + public $velo; |
|
19 | + public $az; |
|
20 | + public $el; |
|
21 | + public $range; |
|
22 | + public $range_rate; |
|
23 | + public $lat; |
|
24 | + public $lon; |
|
25 | + public $alt; |
|
26 | + public $ma; |
|
27 | + public $phase; |
|
28 | + public $footprint; |
|
29 | + public $vis; |
|
30 | + public $orbit; |
|
31 | 31 | |
32 | - public function __construct() |
|
33 | - { |
|
34 | - $this->pos = new Predict_Vector(); |
|
35 | - $this->vel = new Predict_Vector(); |
|
36 | - } |
|
32 | + public function __construct() |
|
33 | + { |
|
34 | + $this->pos = new Predict_Vector(); |
|
35 | + $this->vel = new Predict_Vector(); |
|
36 | + } |
|
37 | 37 | } |
@@ -12,9 +12,9 @@ |
||
12 | 12 | */ |
13 | 13 | class Predict_PassDetail |
14 | 14 | { |
15 | - public $time; /*!< time in "jul_utc" */ |
|
16 | - public $pos; /*!< Raw unprocessed position at time */ |
|
17 | - public $vel; /*!< Raw unprocessed velocity at time */ |
|
15 | + public $time; /*!< time in "jul_utc" */ |
|
16 | + public $pos; /*!< Raw unprocessed position at time */ |
|
17 | + public $vel; /*!< Raw unprocessed velocity at time */ |
|
18 | 18 | public $velo; |
19 | 19 | public $az; |
20 | 20 | public $el; |
@@ -5,8 +5,8 @@ |
||
5 | 5 | */ |
6 | 6 | class Predict_ObsSet |
7 | 7 | { |
8 | - public $az = 0.0; /*!< Azimuth [deg] */ |
|
9 | - public $el = 0.0; /*!< Elevation [deg] */ |
|
10 | - public $range = 0.0; /*!< Range [km] */ |
|
11 | - public $range_rate = 0.0; /*!< Velocity [km/sec] */ |
|
8 | + public $az = 0.0; /*!< Azimuth [deg] */ |
|
9 | + public $el = 0.0; /*!< Elevation [deg] */ |
|
10 | + public $range = 0.0; /*!< Range [km] */ |
|
11 | + public $range_rate = 0.0; /*!< Velocity [km/sec] */ |
|
12 | 12 | } |
@@ -5,8 +5,8 @@ |
||
5 | 5 | */ |
6 | 6 | class Predict_ObsSet |
7 | 7 | { |
8 | - public $az = 0.0; /*!< Azimuth [deg] */ |
|
9 | - public $el = 0.0; /*!< Elevation [deg] */ |
|
10 | - public $range = 0.0; /*!< Range [km] */ |
|
11 | - public $range_rate = 0.0; /*!< Velocity [km/sec] */ |
|
8 | + public $az = 0.0; /*!< Azimuth [deg] */ |
|
9 | + public $el = 0.0; /*!< Elevation [deg] */ |
|
10 | + public $range = 0.0; /*!< Range [km] */ |
|
11 | + public $range_rate = 0.0; /*!< Velocity [km/sec] */ |
|
12 | 12 | } |
@@ -12,221 +12,221 @@ |
||
12 | 12 | */ |
13 | 13 | class Predict_TLE |
14 | 14 | { |
15 | - public $header; /* Header line of TLE file */ |
|
16 | - public $line1; /* Line 1 of TLE */ |
|
17 | - public $line2; /* Line 2 of TLE */ |
|
18 | - public $epoch; /*!< Epoch Time in NORAD TLE format YYDDD.FFFFFFFF */ |
|
19 | - public $epoch_year; /*!< Epoch: year */ |
|
20 | - public $epoch_day; /*!< Epoch: day of year */ |
|
21 | - public $epoch_fod; /*!< Epoch: Fraction of day. */ |
|
22 | - public $xndt2o; /*!< 1. time derivative of mean motion */ |
|
23 | - public $xndd6o; /*!< 2. time derivative of mean motion */ |
|
24 | - public $bstar; /*!< Bstar drag coefficient. */ |
|
25 | - public $xincl; /*!< Inclination */ |
|
26 | - public $xnodeo; /*!< R.A.A.N. */ |
|
27 | - public $eo; /*!< Eccentricity */ |
|
28 | - public $omegao; /*!< argument of perigee */ |
|
29 | - public $xmo; /*!< mean anomaly */ |
|
30 | - public $xno; /*!< mean motion */ |
|
31 | - |
|
32 | - public $catnr; /*!< Catalogue Number. */ |
|
33 | - public $elset; /*!< Element Set number. */ |
|
34 | - public $revnum; /*!< Revolution Number at epoch. */ |
|
35 | - |
|
36 | - public $sat_name; /*!< Satellite name string. */ |
|
37 | - public $idesg; /*!< International Designator. */ |
|
38 | - public $status; /*!< Operational status. */ |
|
39 | - |
|
40 | - /* values needed for squint calculations */ |
|
41 | - public $xincl1; |
|
42 | - public $xnodeo1; |
|
43 | - public $omegao1; |
|
44 | - |
|
45 | - |
|
46 | - /* Converts the strings in a raw two-line element set */ |
|
47 | - /* to their intended numerical values. No processing */ |
|
48 | - /* of these values is done, e.g. from deg to rads etc. */ |
|
49 | - /* This is done in the select_ephemeris() function. */ |
|
50 | - public function __construct($header, $line1, $line2) |
|
51 | - { |
|
52 | - if (!$this->Good_Elements($line1, $line2)) { |
|
53 | - throw new Predict_Exception('Invalid TLE contents'); |
|
54 | - } |
|
55 | - |
|
56 | - $this->header = $header; |
|
57 | - $this->line1 = $line1; |
|
58 | - $this->line2 = $line2; |
|
59 | - |
|
60 | - /** Decode Card 1 **/ |
|
61 | - /* Satellite's catalogue number */ |
|
62 | - $this->catnr = (int) substr($line1, 2, 5); |
|
63 | - |
|
64 | - /* International Designator for satellite */ |
|
65 | - $this->idesg = substr($line1, 9, 8); |
|
66 | - |
|
67 | - /* Epoch time; this is the complete, unconverted epoch. */ |
|
68 | - /* Replace spaces with 0 before casting, as leading spaces are allowed */ |
|
69 | - $this->epoch = (float) str_replace(' ', '0', substr($line1, 18, 14)); |
|
70 | - |
|
71 | - /* Now, convert the epoch time into year, day |
|
15 | + public $header; /* Header line of TLE file */ |
|
16 | + public $line1; /* Line 1 of TLE */ |
|
17 | + public $line2; /* Line 2 of TLE */ |
|
18 | + public $epoch; /*!< Epoch Time in NORAD TLE format YYDDD.FFFFFFFF */ |
|
19 | + public $epoch_year; /*!< Epoch: year */ |
|
20 | + public $epoch_day; /*!< Epoch: day of year */ |
|
21 | + public $epoch_fod; /*!< Epoch: Fraction of day. */ |
|
22 | + public $xndt2o; /*!< 1. time derivative of mean motion */ |
|
23 | + public $xndd6o; /*!< 2. time derivative of mean motion */ |
|
24 | + public $bstar; /*!< Bstar drag coefficient. */ |
|
25 | + public $xincl; /*!< Inclination */ |
|
26 | + public $xnodeo; /*!< R.A.A.N. */ |
|
27 | + public $eo; /*!< Eccentricity */ |
|
28 | + public $omegao; /*!< argument of perigee */ |
|
29 | + public $xmo; /*!< mean anomaly */ |
|
30 | + public $xno; /*!< mean motion */ |
|
31 | + |
|
32 | + public $catnr; /*!< Catalogue Number. */ |
|
33 | + public $elset; /*!< Element Set number. */ |
|
34 | + public $revnum; /*!< Revolution Number at epoch. */ |
|
35 | + |
|
36 | + public $sat_name; /*!< Satellite name string. */ |
|
37 | + public $idesg; /*!< International Designator. */ |
|
38 | + public $status; /*!< Operational status. */ |
|
39 | + |
|
40 | + /* values needed for squint calculations */ |
|
41 | + public $xincl1; |
|
42 | + public $xnodeo1; |
|
43 | + public $omegao1; |
|
44 | + |
|
45 | + |
|
46 | + /* Converts the strings in a raw two-line element set */ |
|
47 | + /* to their intended numerical values. No processing */ |
|
48 | + /* of these values is done, e.g. from deg to rads etc. */ |
|
49 | + /* This is done in the select_ephemeris() function. */ |
|
50 | + public function __construct($header, $line1, $line2) |
|
51 | + { |
|
52 | + if (!$this->Good_Elements($line1, $line2)) { |
|
53 | + throw new Predict_Exception('Invalid TLE contents'); |
|
54 | + } |
|
55 | + |
|
56 | + $this->header = $header; |
|
57 | + $this->line1 = $line1; |
|
58 | + $this->line2 = $line2; |
|
59 | + |
|
60 | + /** Decode Card 1 **/ |
|
61 | + /* Satellite's catalogue number */ |
|
62 | + $this->catnr = (int) substr($line1, 2, 5); |
|
63 | + |
|
64 | + /* International Designator for satellite */ |
|
65 | + $this->idesg = substr($line1, 9, 8); |
|
66 | + |
|
67 | + /* Epoch time; this is the complete, unconverted epoch. */ |
|
68 | + /* Replace spaces with 0 before casting, as leading spaces are allowed */ |
|
69 | + $this->epoch = (float) str_replace(' ', '0', substr($line1, 18, 14)); |
|
70 | + |
|
71 | + /* Now, convert the epoch time into year, day |
|
72 | 72 | and fraction of day, according to: |
73 | 73 | |
74 | 74 | YYDDD.FFFFFFFF |
75 | 75 | */ |
76 | 76 | |
77 | - // Adjust for 2 digit year through 2056 |
|
78 | - $this->epoch_year = (int) substr($line1, 18, 2); |
|
79 | - if ($this->epoch_year > 56) { |
|
80 | - $this->epoch_year = $this->epoch_year + 1900; |
|
81 | - } else { |
|
82 | - $this->epoch_year = $this->epoch_year + 2000; |
|
83 | - } |
|
77 | + // Adjust for 2 digit year through 2056 |
|
78 | + $this->epoch_year = (int) substr($line1, 18, 2); |
|
79 | + if ($this->epoch_year > 56) { |
|
80 | + $this->epoch_year = $this->epoch_year + 1900; |
|
81 | + } else { |
|
82 | + $this->epoch_year = $this->epoch_year + 2000; |
|
83 | + } |
|
84 | 84 | |
85 | - /* Epoch day */ |
|
86 | - $this->epoch_day = (int) substr($line1, 20, 3); |
|
85 | + /* Epoch day */ |
|
86 | + $this->epoch_day = (int) substr($line1, 20, 3); |
|
87 | 87 | |
88 | - /* Epoch fraction of day */ |
|
89 | - $this->epoch_fod = (float) substr($line1, 23, 9); |
|
88 | + /* Epoch fraction of day */ |
|
89 | + $this->epoch_fod = (float) substr($line1, 23, 9); |
|
90 | 90 | |
91 | 91 | |
92 | - /* Satellite's First Time Derivative */ |
|
93 | - $this->xndt2o = (float) substr($line1, 33, 10); |
|
92 | + /* Satellite's First Time Derivative */ |
|
93 | + $this->xndt2o = (float) substr($line1, 33, 10); |
|
94 | 94 | |
95 | - /* Satellite's Second Time Derivative */ |
|
96 | - $this->xndd6o = (float) (substr($line1, 44, 1) . '.' . substr($line1, 45, 5) . 'E' . substr($line1, 50, 2)); |
|
95 | + /* Satellite's Second Time Derivative */ |
|
96 | + $this->xndd6o = (float) (substr($line1, 44, 1) . '.' . substr($line1, 45, 5) . 'E' . substr($line1, 50, 2)); |
|
97 | 97 | |
98 | - /* Satellite's bstar drag term |
|
98 | + /* Satellite's bstar drag term |
|
99 | 99 | FIXME: How about buff[0] ???? |
100 | 100 | */ |
101 | - $this->bstar = (float) (substr($line1, 53, 1) . '.' . substr($line1, 54, 5) . 'E' . substr($line1, 59, 2)); |
|
101 | + $this->bstar = (float) (substr($line1, 53, 1) . '.' . substr($line1, 54, 5) . 'E' . substr($line1, 59, 2)); |
|
102 | 102 | |
103 | - /* Element Number */ |
|
104 | - $this->elset = (int) substr($line1, 64, 4); |
|
103 | + /* Element Number */ |
|
104 | + $this->elset = (int) substr($line1, 64, 4); |
|
105 | 105 | |
106 | - /** Decode Card 2 **/ |
|
107 | - /* Satellite's Orbital Inclination (degrees) */ |
|
108 | - $this->xincl = (float) substr($line2, 8, 8); |
|
106 | + /** Decode Card 2 **/ |
|
107 | + /* Satellite's Orbital Inclination (degrees) */ |
|
108 | + $this->xincl = (float) substr($line2, 8, 8); |
|
109 | 109 | |
110 | - /* Satellite's RAAN (degrees) */ |
|
111 | - $this->xnodeo = (float) substr($line2, 17, 8); |
|
110 | + /* Satellite's RAAN (degrees) */ |
|
111 | + $this->xnodeo = (float) substr($line2, 17, 8); |
|
112 | 112 | |
113 | - /* Satellite's Orbital Eccentricity */ |
|
114 | - $this->eo = (float) ('.' . substr($line2, 26, 7)); |
|
113 | + /* Satellite's Orbital Eccentricity */ |
|
114 | + $this->eo = (float) ('.' . substr($line2, 26, 7)); |
|
115 | 115 | |
116 | - /* Satellite's Argument of Perigee (degrees) */ |
|
117 | - $this->omegao = (float) substr($line2, 34, 8); |
|
116 | + /* Satellite's Argument of Perigee (degrees) */ |
|
117 | + $this->omegao = (float) substr($line2, 34, 8); |
|
118 | 118 | |
119 | - /* Satellite's Mean Anomaly of Orbit (degrees) */ |
|
120 | - $this->xmo = (float) substr($line2, 43, 8); |
|
119 | + /* Satellite's Mean Anomaly of Orbit (degrees) */ |
|
120 | + $this->xmo = (float) substr($line2, 43, 8); |
|
121 | 121 | |
122 | - /* Satellite's Mean Motion (rev/day) */ |
|
123 | - $this->xno = (float) substr($line2, 52, 11); |
|
122 | + /* Satellite's Mean Motion (rev/day) */ |
|
123 | + $this->xno = (float) substr($line2, 52, 11); |
|
124 | 124 | |
125 | - /* Satellite's Revolution number at epoch */ |
|
126 | - $this->revnum = (float) substr($line2, 63, 5); |
|
127 | - } |
|
125 | + /* Satellite's Revolution number at epoch */ |
|
126 | + $this->revnum = (float) substr($line2, 63, 5); |
|
127 | + } |
|
128 | 128 | |
129 | - /* Calculates the checksum mod 10 of a line from a TLE set and */ |
|
130 | - /* returns true if it compares with checksum in column 68, else false.*/ |
|
131 | - /* tle_set is a character string holding the two lines read */ |
|
132 | - /* from a text file containing NASA format Keplerian elements. */ |
|
133 | - /* NOTE!!! The stuff about two lines is not quite true. |
|
129 | + /* Calculates the checksum mod 10 of a line from a TLE set and */ |
|
130 | + /* returns true if it compares with checksum in column 68, else false.*/ |
|
131 | + /* tle_set is a character string holding the two lines read */ |
|
132 | + /* from a text file containing NASA format Keplerian elements. */ |
|
133 | + /* NOTE!!! The stuff about two lines is not quite true. |
|
134 | 134 | The function assumes that tle_set[0] is the begining |
135 | 135 | of the line and that there are 68 elements - see the consumer |
136 | 136 | */ |
137 | - public function Checksum_Good($tle_set) |
|
138 | - { |
|
139 | - if (strlen($tle_set) < 69) { |
|
140 | - return false; |
|
141 | - } |
|
142 | - |
|
143 | - $checksum = 0; |
|
144 | - |
|
145 | - for ($i = 0; $i < 68; $i++) { |
|
146 | - if (($tle_set[$i] >= '0') && ($tle_set[$i] <= '9')) { |
|
147 | - $value = $tle_set[$i] - '0'; |
|
148 | - } else if ($tle_set[$i] == '-' ) { |
|
149 | - $value = 1; |
|
150 | - } else { |
|
151 | - $value = 0; |
|
152 | - } |
|
153 | - |
|
154 | - $checksum += $value; |
|
155 | - } |
|
156 | - |
|
157 | - $checksum %= 10; |
|
158 | - $check_digit = $tle_set[68] - '0'; |
|
159 | - |
|
160 | - return $checksum == $check_digit; |
|
161 | - } |
|
162 | - |
|
163 | - /* Carries out various checks on a TLE set to verify its validity */ |
|
164 | - /* $line1 is the first line of the TLE, $line2 is the second line */ |
|
165 | - /* from a text file containing NASA format Keplerian elements. */ |
|
166 | - public function Good_Elements($line1, $line2) |
|
167 | - { |
|
168 | - /* Verify checksum of both lines of a TLE set */ |
|
169 | - if (!$this->Checksum_Good($line1) || !$this->Checksum_Good($line2)) { |
|
170 | - return false; |
|
171 | - } |
|
172 | - |
|
173 | - /* Check the line number of each line */ |
|
174 | - if (($line1[0] != '1') || ($line2[0] != '2')) { |
|
175 | - return false; |
|
176 | - } |
|
177 | - |
|
178 | - /* Verify that Satellite Number is same in both lines */ |
|
179 | - if (strncmp($line1[2], $line2[2], 5) != 0) { |
|
180 | - return false; |
|
181 | - } |
|
182 | - |
|
183 | - /* Check that various elements are in the right place */ |
|
184 | - if (($line1[23] != '.') || |
|
185 | - ($line1[34] != '.') || |
|
186 | - ($line2[11] != '.') || |
|
187 | - ($line2[20] != '.') || |
|
188 | - ($line2[37] != '.') || |
|
189 | - ($line2[46] != '.') || |
|
190 | - ($line2[54] != '.') || |
|
191 | - (strncmp(substr($line1, 61), ' 0 ', 3) != 0)) { |
|
192 | - |
|
193 | - return false; |
|
194 | - } |
|
195 | - |
|
196 | - return true; |
|
197 | - } |
|
198 | - |
|
199 | - /** |
|
200 | - * A function to allow checksum creation of a line. This is driven by |
|
201 | - * the fact that some TLEs from SpaceTrack are missing checksum numbers. |
|
202 | - * You can use this to create a checksum for a line, but you should |
|
203 | - * probably have confidence that the TLE data itself is good. YMMV. |
|
204 | - * |
|
205 | - * @throws Predict_Exception if the line is not exactly 68 chars |
|
206 | - * @return string |
|
207 | - */ |
|
208 | - static public function createChecksum($line) |
|
209 | - { |
|
210 | - if (strlen($line) != 68) { |
|
211 | - throw Predict_Exception('Invalid line, needs to e 68 chars'); |
|
212 | - } |
|
213 | - |
|
214 | - $checksum = 0; |
|
215 | - |
|
216 | - for ($i = 0; $i < 68; $i++) { |
|
217 | - if (($line[$i] >= '0') && ($line[$i] <= '9')) { |
|
218 | - $value = (int) $line[$i]; |
|
219 | - } else if ($line[$i] == '-' ) { |
|
220 | - $value = 1; |
|
221 | - } else { |
|
222 | - $value = 0; |
|
223 | - } |
|
224 | - |
|
225 | - $checksum += $value; |
|
226 | - } |
|
227 | - |
|
228 | - $checksum %= 10; |
|
229 | - |
|
230 | - return $checksum; |
|
231 | - } |
|
137 | + public function Checksum_Good($tle_set) |
|
138 | + { |
|
139 | + if (strlen($tle_set) < 69) { |
|
140 | + return false; |
|
141 | + } |
|
142 | + |
|
143 | + $checksum = 0; |
|
144 | + |
|
145 | + for ($i = 0; $i < 68; $i++) { |
|
146 | + if (($tle_set[$i] >= '0') && ($tle_set[$i] <= '9')) { |
|
147 | + $value = $tle_set[$i] - '0'; |
|
148 | + } else if ($tle_set[$i] == '-' ) { |
|
149 | + $value = 1; |
|
150 | + } else { |
|
151 | + $value = 0; |
|
152 | + } |
|
153 | + |
|
154 | + $checksum += $value; |
|
155 | + } |
|
156 | + |
|
157 | + $checksum %= 10; |
|
158 | + $check_digit = $tle_set[68] - '0'; |
|
159 | + |
|
160 | + return $checksum == $check_digit; |
|
161 | + } |
|
162 | + |
|
163 | + /* Carries out various checks on a TLE set to verify its validity */ |
|
164 | + /* $line1 is the first line of the TLE, $line2 is the second line */ |
|
165 | + /* from a text file containing NASA format Keplerian elements. */ |
|
166 | + public function Good_Elements($line1, $line2) |
|
167 | + { |
|
168 | + /* Verify checksum of both lines of a TLE set */ |
|
169 | + if (!$this->Checksum_Good($line1) || !$this->Checksum_Good($line2)) { |
|
170 | + return false; |
|
171 | + } |
|
172 | + |
|
173 | + /* Check the line number of each line */ |
|
174 | + if (($line1[0] != '1') || ($line2[0] != '2')) { |
|
175 | + return false; |
|
176 | + } |
|
177 | + |
|
178 | + /* Verify that Satellite Number is same in both lines */ |
|
179 | + if (strncmp($line1[2], $line2[2], 5) != 0) { |
|
180 | + return false; |
|
181 | + } |
|
182 | + |
|
183 | + /* Check that various elements are in the right place */ |
|
184 | + if (($line1[23] != '.') || |
|
185 | + ($line1[34] != '.') || |
|
186 | + ($line2[11] != '.') || |
|
187 | + ($line2[20] != '.') || |
|
188 | + ($line2[37] != '.') || |
|
189 | + ($line2[46] != '.') || |
|
190 | + ($line2[54] != '.') || |
|
191 | + (strncmp(substr($line1, 61), ' 0 ', 3) != 0)) { |
|
192 | + |
|
193 | + return false; |
|
194 | + } |
|
195 | + |
|
196 | + return true; |
|
197 | + } |
|
198 | + |
|
199 | + /** |
|
200 | + * A function to allow checksum creation of a line. This is driven by |
|
201 | + * the fact that some TLEs from SpaceTrack are missing checksum numbers. |
|
202 | + * You can use this to create a checksum for a line, but you should |
|
203 | + * probably have confidence that the TLE data itself is good. YMMV. |
|
204 | + * |
|
205 | + * @throws Predict_Exception if the line is not exactly 68 chars |
|
206 | + * @return string |
|
207 | + */ |
|
208 | + static public function createChecksum($line) |
|
209 | + { |
|
210 | + if (strlen($line) != 68) { |
|
211 | + throw Predict_Exception('Invalid line, needs to e 68 chars'); |
|
212 | + } |
|
213 | + |
|
214 | + $checksum = 0; |
|
215 | + |
|
216 | + for ($i = 0; $i < 68; $i++) { |
|
217 | + if (($line[$i] >= '0') && ($line[$i] <= '9')) { |
|
218 | + $value = (int) $line[$i]; |
|
219 | + } else if ($line[$i] == '-' ) { |
|
220 | + $value = 1; |
|
221 | + } else { |
|
222 | + $value = 0; |
|
223 | + } |
|
224 | + |
|
225 | + $checksum += $value; |
|
226 | + } |
|
227 | + |
|
228 | + $checksum %= 10; |
|
229 | + |
|
230 | + return $checksum; |
|
231 | + } |
|
232 | 232 | } |
@@ -12,30 +12,30 @@ discard block |
||
12 | 12 | */ |
13 | 13 | class Predict_TLE |
14 | 14 | { |
15 | - public $header; /* Header line of TLE file */ |
|
16 | - public $line1; /* Line 1 of TLE */ |
|
17 | - public $line2; /* Line 2 of TLE */ |
|
18 | - public $epoch; /*!< Epoch Time in NORAD TLE format YYDDD.FFFFFFFF */ |
|
15 | + public $header; /* Header line of TLE file */ |
|
16 | + public $line1; /* Line 1 of TLE */ |
|
17 | + public $line2; /* Line 2 of TLE */ |
|
18 | + public $epoch; /*!< Epoch Time in NORAD TLE format YYDDD.FFFFFFFF */ |
|
19 | 19 | public $epoch_year; /*!< Epoch: year */ |
20 | - public $epoch_day; /*!< Epoch: day of year */ |
|
21 | - public $epoch_fod; /*!< Epoch: Fraction of day. */ |
|
22 | - public $xndt2o; /*!< 1. time derivative of mean motion */ |
|
23 | - public $xndd6o; /*!< 2. time derivative of mean motion */ |
|
24 | - public $bstar; /*!< Bstar drag coefficient. */ |
|
25 | - public $xincl; /*!< Inclination */ |
|
26 | - public $xnodeo; /*!< R.A.A.N. */ |
|
27 | - public $eo; /*!< Eccentricity */ |
|
28 | - public $omegao; /*!< argument of perigee */ |
|
29 | - public $xmo; /*!< mean anomaly */ |
|
30 | - public $xno; /*!< mean motion */ |
|
31 | - |
|
32 | - public $catnr; /*!< Catalogue Number. */ |
|
33 | - public $elset; /*!< Element Set number. */ |
|
34 | - public $revnum; /*!< Revolution Number at epoch. */ |
|
35 | - |
|
36 | - public $sat_name; /*!< Satellite name string. */ |
|
37 | - public $idesg; /*!< International Designator. */ |
|
38 | - public $status; /*!< Operational status. */ |
|
20 | + public $epoch_day; /*!< Epoch: day of year */ |
|
21 | + public $epoch_fod; /*!< Epoch: Fraction of day. */ |
|
22 | + public $xndt2o; /*!< 1. time derivative of mean motion */ |
|
23 | + public $xndd6o; /*!< 2. time derivative of mean motion */ |
|
24 | + public $bstar; /*!< Bstar drag coefficient. */ |
|
25 | + public $xincl; /*!< Inclination */ |
|
26 | + public $xnodeo; /*!< R.A.A.N. */ |
|
27 | + public $eo; /*!< Eccentricity */ |
|
28 | + public $omegao; /*!< argument of perigee */ |
|
29 | + public $xmo; /*!< mean anomaly */ |
|
30 | + public $xno; /*!< mean motion */ |
|
31 | + |
|
32 | + public $catnr; /*!< Catalogue Number. */ |
|
33 | + public $elset; /*!< Element Set number. */ |
|
34 | + public $revnum; /*!< Revolution Number at epoch. */ |
|
35 | + |
|
36 | + public $sat_name; /*!< Satellite name string. */ |
|
37 | + public $idesg; /*!< International Designator. */ |
|
38 | + public $status; /*!< Operational status. */ |
|
39 | 39 | |
40 | 40 | /* values needed for squint calculations */ |
41 | 41 | public $xincl1; |
@@ -93,12 +93,12 @@ discard block |
||
93 | 93 | $this->xndt2o = (float) substr($line1, 33, 10); |
94 | 94 | |
95 | 95 | /* Satellite's Second Time Derivative */ |
96 | - $this->xndd6o = (float) (substr($line1, 44, 1) . '.' . substr($line1, 45, 5) . 'E' . substr($line1, 50, 2)); |
|
96 | + $this->xndd6o = (float) (substr($line1, 44, 1).'.'.substr($line1, 45, 5).'E'.substr($line1, 50, 2)); |
|
97 | 97 | |
98 | 98 | /* Satellite's bstar drag term |
99 | 99 | FIXME: How about buff[0] ???? |
100 | 100 | */ |
101 | - $this->bstar = (float) (substr($line1, 53, 1) . '.' . substr($line1, 54, 5) . 'E' . substr($line1, 59, 2)); |
|
101 | + $this->bstar = (float) (substr($line1, 53, 1).'.'.substr($line1, 54, 5).'E'.substr($line1, 59, 2)); |
|
102 | 102 | |
103 | 103 | /* Element Number */ |
104 | 104 | $this->elset = (int) substr($line1, 64, 4); |
@@ -111,7 +111,7 @@ discard block |
||
111 | 111 | $this->xnodeo = (float) substr($line2, 17, 8); |
112 | 112 | |
113 | 113 | /* Satellite's Orbital Eccentricity */ |
114 | - $this->eo = (float) ('.' . substr($line2, 26, 7)); |
|
114 | + $this->eo = (float) ('.'.substr($line2, 26, 7)); |
|
115 | 115 | |
116 | 116 | /* Satellite's Argument of Perigee (degrees) */ |
117 | 117 | $this->omegao = (float) substr($line2, 34, 8); |
@@ -145,7 +145,7 @@ discard block |
||
145 | 145 | for ($i = 0; $i < 68; $i++) { |
146 | 146 | if (($tle_set[$i] >= '0') && ($tle_set[$i] <= '9')) { |
147 | 147 | $value = $tle_set[$i] - '0'; |
148 | - } else if ($tle_set[$i] == '-' ) { |
|
148 | + } else if ($tle_set[$i] == '-') { |
|
149 | 149 | $value = 1; |
150 | 150 | } else { |
151 | 151 | $value = 0; |
@@ -216,7 +216,7 @@ discard block |
||
216 | 216 | for ($i = 0; $i < 68; $i++) { |
217 | 217 | if (($line[$i] >= '0') && ($line[$i] <= '9')) { |
218 | 218 | $value = (int) $line[$i]; |
219 | - } else if ($line[$i] == '-' ) { |
|
219 | + } else if ($line[$i] == '-') { |
|
220 | 220 | $value = 1; |
221 | 221 | } else { |
222 | 222 | $value = 0; |
@@ -6,8 +6,8 @@ |
||
6 | 6 | */ |
7 | 7 | class Predict_Geodetic |
8 | 8 | { |
9 | - public $lat; /*!< Lattitude [rad] */ |
|
10 | - public $lon; /*!< Longitude [rad] */ |
|
11 | - public $alt; /*!< Altitude [km] */ |
|
12 | - public $theta; |
|
9 | + public $lat; /*!< Lattitude [rad] */ |
|
10 | + public $lon; /*!< Longitude [rad] */ |
|
11 | + public $alt; /*!< Altitude [km] */ |
|
12 | + public $theta; |
|
13 | 13 | } |
@@ -8,32 +8,32 @@ |
||
8 | 8 | */ |
9 | 9 | class Predict_SGSDPStatic |
10 | 10 | { |
11 | - public $aodp; |
|
12 | - public $aycof; |
|
13 | - public $c1; |
|
14 | - public $c4; |
|
15 | - public $c5; |
|
16 | - public $cosio; |
|
17 | - public $d2; |
|
18 | - public $d3; |
|
19 | - public $d4; |
|
20 | - public $delmo; |
|
21 | - public $omgcof; |
|
22 | - public $eta; |
|
23 | - public $omgdot; |
|
24 | - public $sinio; |
|
25 | - public $xnodp; |
|
26 | - public $sinmo; |
|
27 | - public $t2cof; |
|
28 | - public $t3cof; |
|
29 | - public $t4cof; |
|
30 | - public $t5cof; |
|
31 | - public $x1mth2; |
|
32 | - public $x3thm1; |
|
33 | - public $x7thm1; |
|
34 | - public $xmcof; |
|
35 | - public $xmdot; |
|
36 | - public $xnodcf; |
|
37 | - public $xnodot; |
|
38 | - public $xlcof; |
|
11 | + public $aodp; |
|
12 | + public $aycof; |
|
13 | + public $c1; |
|
14 | + public $c4; |
|
15 | + public $c5; |
|
16 | + public $cosio; |
|
17 | + public $d2; |
|
18 | + public $d3; |
|
19 | + public $d4; |
|
20 | + public $delmo; |
|
21 | + public $omgcof; |
|
22 | + public $eta; |
|
23 | + public $omgdot; |
|
24 | + public $sinio; |
|
25 | + public $xnodp; |
|
26 | + public $sinmo; |
|
27 | + public $t2cof; |
|
28 | + public $t3cof; |
|
29 | + public $t4cof; |
|
30 | + public $t5cof; |
|
31 | + public $x1mth2; |
|
32 | + public $x3thm1; |
|
33 | + public $x7thm1; |
|
34 | + public $xmcof; |
|
35 | + public $xmdot; |
|
36 | + public $xnodcf; |
|
37 | + public $xnodot; |
|
38 | + public $xlcof; |
|
39 | 39 | } |