1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
//! Parameter assignments for SGP4 algorithm as expressed
//! by Vallado \[2].
//!
//! ```text
//! C$Procedure ZZSGP4 ( SGP4 parameters )
//! C
//! C$ Abstract
//! C
//! C Parameter assignments for SGP4 algorithm as expressed
//! C by Vallado [2].
//! C
//! C$ Disclaimer
//! C
//! C THIS SOFTWARE AND ANY RELATED MATERIALS WERE CREATED BY THE
//! C CALIFORNIA INSTITUTE OF TECHNOLOGY (CALTECH) UNDER A U.S.
//! C GOVERNMENT CONTRACT WITH THE NATIONAL AERONAUTICS AND SPACE
//! C ADMINISTRATION (NASA). THE SOFTWARE IS TECHNOLOGY AND SOFTWARE
//! C PUBLICLY AVAILABLE UNDER U.S. EXPORT LAWS AND IS PROVIDED "AS-IS"
//! C TO THE RECIPIENT WITHOUT WARRANTY OF ANY KIND, INCLUDING ANY
//! C WARRANTIES OF PERFORMANCE OR MERCHANTABILITY OR FITNESS FOR A
//! C PARTICULAR USE OR PURPOSE (AS SET FORTH IN UNITED STATES UCC
//! C SECTIONS 2312-2313) OR FOR ANY PURPOSE WHATSOEVER, FOR THE
//! C SOFTWARE AND RELATED MATERIALS, HOWEVER USED.
//! C
//! C IN NO EVENT SHALL CALTECH, ITS JET PROPULSION LABORATORY, OR NASA
//! C BE LIABLE FOR ANY DAMAGES AND/OR COSTS, INCLUDING, BUT NOT
//! C LIMITED TO, INCIDENTAL OR CONSEQUENTIAL DAMAGES OF ANY KIND,
//! C INCLUDING ECONOMIC DAMAGE OR INJURY TO PROPERTY AND LOST PROFITS,
//! C REGARDLESS OF WHETHER CALTECH, JPL, OR NASA BE ADVISED, HAVE
//! C REASON TO KNOW, OR, IN FACT, SHALL KNOW OF THE POSSIBILITY.
//! C
//! C RECIPIENT BEARS ALL RISK RELATING TO QUALITY AND PERFORMANCE OF
//! C THE SOFTWARE AND ANY RELATED MATERIALS, AND AGREES TO INDEMNIFY
//! C CALTECH AND NASA FOR ALL THIRD-PARTY CLAIMS RESULTING FROM THE
//! C ACTIONS OF RECIPIENT IN THE USE OF THE SOFTWARE.
//! C
//! C$ Required_Reading
//! C
//! C None.
//! C
//! C$ Keywords
//! C
//! C None.
//! C
//! C$ Declarations
//! C
//! C None.
//! C
//! C$ Brief_I/O
//! C
//! C None.
//! C
//! C$ Detailed_Input
//! C
//! C None.
//! C
//! C$ Detailed_Output
//! C
//! C None.
//! C
//! C$ Parameters
//! C
//! C None.
//! C
//! C$ Exceptions
//! C
//! C None.
//! C
//! C$ Files
//! C
//! C None.
//! C
//! C$ Particulars
//! C
//! C None.
//! C
//! C$ Examples
//! C
//! C J2 = GEOPHS(K_J2)
//! C J3 = GEOPHS(K_J3)
//! C J4 = GEOPHS(K_J4)
//! C ER = GEOPHS(K_ER)
//! C XKE = GEOPHS(K_KE)
//! C
//! C TUMIN = 1.D0/XKE
//! C J3OJ2 = J3/J2
//! C
//! C$ Restrictions
//! C
//! C None.
//! C
//! C$ Literature_References
//! C
//! C [1] Hoots, F. R., and Roehrich, R. L. 1980. "Models for
//! C Propagation of the NORAD Element Sets." Spacetrack Report #3.
//! C U.S. Air Force: Aerospace Defense Command.
//! C
//! C [2] Vallado, David, Crawford, Paul, Hujsak, Richard, and Kelso, T.S.
//! C 2006. Revisiting Spacetrack Report #3. Paper AIAA 2006-6753
//! C presented at the AIAA/AAS Astrodynamics Specialist Conference,
//! C August 21-24, 2006. Keystone, CO.
//! C
//! C$ Author_and_Institution
//! C
//! C E. D. Wright (JPL)
//! C
//! C$ Version
//! C
//! C- SPICELIB Version 1.0.0, MAY-27-2020 (EDW)
//! C
//! C Updated descriptions of GEOPHS constants to be consistent
//! C with what's used in other routines.
//! C
//! C- SPICELIB Version 1.0.0 22-JUL-2014 (EDW)
//! C
//! C-&
//!
//! C$ Index_Entries
//! C
//! C SGP4
//! C
//! C-&
//!
//!
//! C
//! C WGS gravitational constants IDs.
//! C
//!
//! INTEGER NGRAVS
//! PARAMETER ( NGRAVS = 3 )
//!
//! INTEGER NGRAVC
//! PARAMETER ( NGRAVC = 8 )
//!
//!
//! INTEGER WGS721
//! PARAMETER ( WGS721 = 1 )
//!
//! INTEGER WGS72
//! PARAMETER ( WGS72 = 2 )
//!
//! INTEGER WGS84
//! PARAMETER ( WGS84 = 3 )
//!
//!
//! C
//! C Gravitational constant indices.
//! C
//!
//! INTEGER P_RAD
//! PARAMETER ( P_RAD = 1 )
//!
//! INTEGER P_XKE
//! PARAMETER ( P_XKE = 2 )
//!
//! INTEGER P_MU
//! PARAMETER ( P_MU = 3 )
//!
//! INTEGER P_TUMN
//! PARAMETER ( P_TUMN = 4 )
//!
//! INTEGER P_J2
//! PARAMETER ( P_J2 = 5 )
//!
//! INTEGER P_J3
//! PARAMETER ( P_J3 = 6 )
//!
//! INTEGER P_J4
//! PARAMETER ( P_J4 = 7 )
//!
//! INTEGER P_J3J2
//! PARAMETER ( P_J3J2 = 8 )
//!
//!
//! C
//! C The following parameters give the indices in the GEOPHS
//! C array of the various geophysical parameters needed for
//! C the two line element sets.
//! C
//! C K_J2 --- index of J2 gravitational harmonic for earth
//! C K_J3 --- index of J3 gravitational harmonic for earth
//! C K_J4 --- index of J4 gravitational harmonic for earth
//! C K_KE --- index of KE = sqrt(GM) in earth-radii**1.5/MIN
//! C K_QO --- index of high altitude bound for atmospheric
//! C model in km
//! C K_SO --- index of low altitude bound for atmospheric
//! C model in km
//! C K_ER --- index of earth equatorial radius in km
//! C K_AE --- index of distance units/earth radius
//! C
//!
//! INTEGER K_J2
//! PARAMETER ( K_J2 = 1 )
//!
//! INTEGER K_J3
//! PARAMETER ( K_J3 = 2 )
//!
//! INTEGER K_J4
//! PARAMETER ( K_J4 = 3 )
//!
//! INTEGER K_KE
//! PARAMETER ( K_KE = 4 )
//!
//! INTEGER K_QO
//! PARAMETER ( K_QO = 5 )
//!
//! INTEGER K_SO
//! PARAMETER ( K_SO = 6 )
//!
//! INTEGER K_ER
//! PARAMETER ( K_ER = 7 )
//!
//! INTEGER K_AE
//! PARAMETER ( K_AE = 8 )
//!
//! INTEGER NGEO
//! PARAMETER ( NGEO = K_AE )
//!
//! C
//! C Operation mode values, OPMODE.
//! C
//!
//! INTEGER AFSPC
//! PARAMETER ( AFSPC = 1 )
//!
//! INTEGER IMPRVD
//! PARAMETER ( IMPRVD = 2 )
//!
//! C
//! C An enumeration of the various components of the
//! C elements array---ELEMS
//! C
//! C KNDT20 --- location of NDT20
//! C KNDD60 --- location of NDD60
//! C KBSTAR --- location of BSTAR
//! C KINCL --- location of INCL
//! C KNODE0 --- location of NODE0
//! C KECC --- location of ECC
//! C KOMEGA --- location of OMEGA
//! C KMO --- location of MO
//! C KNO --- location of NO
//! C
//!
//! INTEGER KNDT20
//! PARAMETER ( KNDT20 = 1 )
//!
//! INTEGER KNDD60
//! PARAMETER ( KNDD60 = 2 )
//!
//! INTEGER KBSTAR
//! PARAMETER ( KBSTAR = 3 )
//!
//! INTEGER KINCL
//! PARAMETER ( KINCL = 4 )
//!
//! INTEGER KNODE0
//! PARAMETER ( KNODE0 = 5 )
//!
//! INTEGER KECC
//! PARAMETER ( KECC = 6 )
//!
//! INTEGER KOMEGA
//! PARAMETER ( KOMEGA = 7 )
//!
//! INTEGER KMO
//! PARAMETER ( KMO = 8 )
//!
//! INTEGER KNO
//! PARAMETER ( KNO = 9 )
//!
//! INTEGER KEPOCH
//! PARAMETER ( KEPOCH = 10 )
//!
//! INTEGER NELEMS
//! PARAMETER ( NELEMS = KEPOCH )
//!
//!
//! ```
pub const NGRAVS: i32 = 3;
pub const NGRAVC: i32 = 8;
pub const WGS721: i32 = 1;
pub const WGS72: i32 = 2;
pub const WGS84: i32 = 3;
pub const P_RAD: i32 = 1;
pub const P_XKE: i32 = 2;
pub const P_MU: i32 = 3;
pub const P_TUMN: i32 = 4;
pub const P_J2: i32 = 5;
pub const P_J3: i32 = 6;
pub const P_J4: i32 = 7;
pub const P_J3J2: i32 = 8;
pub const K_J2: i32 = 1;
pub const K_J3: i32 = 2;
pub const K_J4: i32 = 3;
pub const K_KE: i32 = 4;
pub const K_QO: i32 = 5;
pub const K_SO: i32 = 6;
pub const K_ER: i32 = 7;
pub const K_AE: i32 = 8;
pub const NGEO: i32 = K_AE;
pub const AFSPC: i32 = 1;
pub const IMPRVD: i32 = 2;
pub const KNDT20: i32 = 1;
pub const KNDD60: i32 = 2;
pub const KBSTAR: i32 = 3;
pub const KINCL: i32 = 4;
pub const KNODE0: i32 = 5;
pub const KECC: i32 = 6;
pub const KOMEGA: i32 = 7;
pub const KMO: i32 = 8;
pub const KNO: i32 = 9;
pub const KEPOCH: i32 = 10;
pub const NELEMS: i32 = KEPOCH;