korba 0.7.1 → 0.9.0

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (15) hide show
  1. checksums.yaml +4 -4
  2. data/CHANGELOG.md +9 -0
  3. data/lib/korba/car.rb +82 -0
  4. data/lib/korba/constant.rb +2 -1
  5. data/lib/korba/kep.rb +8 -0
  6. data/lib/korba/orbit_utils.rb +6 -8
  7. data/lib/korba/propagator/kepler.rb +54 -0
  8. data/lib/korba/propagator/rk4.rb +63 -0
  9. data/lib/korba/propagator/sgp4.rb +2029 -0
  10. data/lib/korba/sgp4/elset_rec.rb +421 -419
  11. data/lib/korba/tle.rb +4 -4
  12. data/lib/korba/version.rb +1 -1
  13. data/lib/korba.rb +3 -1
  14. metadata +4 -2
  15. data/lib/korba/sgp4/sgp4.rb +0 -2027
data/lib/korba/sgp4/elset_rec.rb CHANGED
@@ -1,432 +1,434 @@
1
1
  module Korba
2
- #
3
- # This class implements the elsetrec data type from Vallado's SGP4 code.
4
- #
5
- # From SGP4.h
6
- # #define SGP4Version "SGP4 Version 2016-03-09"
7
- #
8
- # @author aholinch
9
- #
10
- #
11
- class ElsetRec
12
- attr_accessor :whichconst
13
- attr_accessor :satnum
14
- attr_accessor :epochyr
15
- attr_accessor :epochtynumrev
16
- attr_accessor :error
17
- attr_accessor :operationmode
18
- attr_accessor :init
19
- attr_accessor :method
20
- attr_accessor :a
21
- attr_accessor :altp
22
- attr_accessor :alta
23
- attr_accessor :epochdays
24
- attr_accessor :jdsatepoch
25
- attr_accessor :jdsatepochF
26
- attr_accessor :nddot
27
- attr_accessor :ndot
28
- attr_accessor :bstar
29
- attr_accessor :rcse
30
- attr_accessor :inclo
31
- attr_accessor :nodeo
32
- attr_accessor :ecco
33
- attr_accessor :argpo
34
- attr_accessor :mo
35
- attr_accessor :no_kozai
36
- attr_accessor :classification
37
- attr_accessor :intldesg
38
- attr_accessor :ephtype
39
- attr_accessor :elnum
40
- attr_accessor :revnum
41
- attr_accessor :no_unkozai
42
- attr_accessor :am
43
- attr_accessor :em
44
- attr_accessor :im
45
- attr_accessor :Om
46
- attr_accessor :om
47
- attr_accessor :mm
48
- attr_accessor :nm
49
- attr_accessor :t
50
- attr_accessor :tumin
51
- attr_accessor :mu
52
- attr_accessor :radiusearthkm
53
- attr_accessor :xke
54
- attr_accessor :j2
55
- attr_accessor :j3
56
- attr_accessor :j4
57
- attr_accessor :j3oj2
58
- attr_accessor :dia_mm
59
- attr_accessor :period_sec
60
- attr_accessor :active
61
- attr_accessor :not_orbital
62
- attr_accessor :rcs_m2
63
- attr_accessor :ep
64
- attr_accessor :inclp
65
- attr_accessor :nodep
66
- attr_accessor :argpp
67
- attr_accessor :mp
68
- attr_accessor :isimp
69
- attr_accessor :aycof
70
- attr_accessor :con41
71
- attr_accessor :cc1
72
- attr_accessor :cc4
73
- attr_accessor :cc5
74
- attr_accessor :d2
75
- attr_accessor :d3
76
- attr_accessor :d4
77
- attr_accessor :delmo
78
- attr_accessor :eta
79
- attr_accessor :argpdot
80
- attr_accessor :omgcof
81
- attr_accessor :sinmao
82
- attr_accessor :t2cof
83
- attr_accessor :t3cof
84
- attr_accessor :t4cof
85
- attr_accessor :t5cof
86
- attr_accessor :x1mth2
87
- attr_accessor :x7thm1
88
- attr_accessor :mdot
89
- attr_accessor :nodedot
90
- attr_accessor :xlcof
91
- attr_accessor :xmcof
92
- attr_accessor :nodecf
93
- attr_accessor :irez
94
- attr_accessor :d2201
95
- attr_accessor :d2211
96
- attr_accessor :d3210
97
- attr_accessor :d3222
98
- attr_accessor :d4410
99
- attr_accessor :d4422
100
- attr_accessor :d5220
101
- attr_accessor :d5232
102
- attr_accessor :d5421
103
- attr_accessor :d5433
104
- attr_accessor :dedt
105
- attr_accessor :del1
106
- attr_accessor :del2
107
- attr_accessor :del3
108
- attr_accessor :didt
109
- attr_accessor :dmdt
110
- attr_accessor :dnodt
111
- attr_accessor :domdt
112
- attr_accessor :e3
113
- attr_accessor :ee2
114
- attr_accessor :peo
115
- attr_accessor :pgho
116
- attr_accessor :pho
117
- attr_accessor :pinco
118
- attr_accessor :plo
119
- attr_accessor :se2
120
- attr_accessor :se3
121
- attr_accessor :sgh2
122
- attr_accessor :sgh3
123
- attr_accessor :sgh4
124
- attr_accessor :sh2
125
- attr_accessor :sh3
126
- attr_accessor :si2
127
- attr_accessor :si3
128
- attr_accessor :sl2
129
- attr_accessor :sl3
130
- attr_accessor :sl4
131
- attr_accessor :gsto
132
- attr_accessor :xfact
133
- attr_accessor :xgh2
134
- attr_accessor :xgh3
135
- attr_accessor :xgh4
136
- attr_accessor :xh2
137
- attr_accessor :xh3
138
- attr_accessor :xi2
139
- attr_accessor :xi3
140
- attr_accessor :xl2
141
- attr_accessor :xl3
142
- attr_accessor :xl4
143
- attr_accessor :xlamo
144
- attr_accessor :zmol
145
- attr_accessor :zmos
146
- attr_accessor :atime
147
- attr_accessor :xli
148
- attr_accessor :xni
149
- attr_accessor :snodm
150
- attr_accessor :cnodm
151
- attr_accessor :sinim
152
- attr_accessor :cosim
153
- attr_accessor :sinomm
154
- attr_accessor :cosomm
155
- attr_accessor :day
156
- attr_accessor :emsq
157
- attr_accessor :gam
158
- attr_accessor :rtemsq
159
- attr_accessor :s1
160
- attr_accessor :s2
161
- attr_accessor :s3
162
- attr_accessor :s4
163
- attr_accessor :s5
164
- attr_accessor :s6
165
- attr_accessor :s7
166
- attr_accessor :ss1
167
- attr_accessor :ss2
168
- attr_accessor :ss3
169
- attr_accessor :ss4
170
- attr_accessor :ss5
171
- attr_accessor :ss6
172
- attr_accessor :ss7
173
- attr_accessor :sz1
174
- attr_accessor :sz2
175
- attr_accessor :sz3
176
- attr_accessor :sz11
177
- attr_accessor :sz12
178
- attr_accessor :sz13
179
- attr_accessor :sz21
180
- attr_accessor :sz22
181
- attr_accessor :sz23
182
- attr_accessor :sz31
183
- attr_accessor :sz32
184
- attr_accessor :sz33
185
- attr_accessor :z1
186
- attr_accessor :z2
187
- attr_accessor :z3
188
- attr_accessor :z11
189
- attr_accessor :z12
190
- attr_accessor :z13
191
- attr_accessor :z21
192
- attr_accessor :z22
193
- attr_accessor :z23
194
- attr_accessor :z31
195
- attr_accessor :z32
196
- attr_accessor :z33
197
- attr_accessor :argpm
198
- attr_accessor :inclm
199
- attr_accessor :nodem
200
- attr_accessor :dndt
201
- attr_accessor :eccsq
202
- attr_accessor :ainv
203
- attr_accessor :ao
204
- attr_accessor :con42
205
- attr_accessor :cosio
206
- attr_accessor :cosio2
207
- attr_accessor :omeosq
208
- attr_accessor :posq
209
- attr_accessor :rp
210
- attr_accessor :rteosq
211
- attr_accessor :sinio
2
+ module Sgp4
3
+ #
4
+ # This class implements the elsetrec data type from Vallado's SGP4 code.
5
+ #
6
+ # From SGP4.h
7
+ # #define SGP4Version "SGP4 Version 2016-03-09"
8
+ #
9
+ # @author aholinch
10
+ #
11
+ #
12
+ class ElsetRec
13
+ attr_accessor :whichconst
14
+ attr_accessor :satnum
15
+ attr_accessor :epochyr
16
+ attr_accessor :epochtynumrev
17
+ attr_accessor :error
18
+ attr_accessor :operationmode
19
+ attr_accessor :init
20
+ attr_accessor :method
21
+ attr_accessor :a
22
+ attr_accessor :altp
23
+ attr_accessor :alta
24
+ attr_accessor :epochdays
25
+ attr_accessor :jdsatepoch
26
+ attr_accessor :jdsatepochF
27
+ attr_accessor :nddot
28
+ attr_accessor :ndot
29
+ attr_accessor :bstar
30
+ attr_accessor :rcse
31
+ attr_accessor :inclo
32
+ attr_accessor :nodeo
33
+ attr_accessor :ecco
34
+ attr_accessor :argpo
35
+ attr_accessor :mo
36
+ attr_accessor :no_kozai
37
+ attr_accessor :classification
38
+ attr_accessor :intldesg
39
+ attr_accessor :ephtype
40
+ attr_accessor :elnum
41
+ attr_accessor :revnum
42
+ attr_accessor :no_unkozai
43
+ attr_accessor :am
44
+ attr_accessor :em
45
+ attr_accessor :im
46
+ attr_accessor :Om
47
+ attr_accessor :om
48
+ attr_accessor :mm
49
+ attr_accessor :nm
50
+ attr_accessor :t
51
+ attr_accessor :tumin
52
+ attr_accessor :mu
53
+ attr_accessor :radiusearthkm
54
+ attr_accessor :xke
55
+ attr_accessor :j2
56
+ attr_accessor :j3
57
+ attr_accessor :j4
58
+ attr_accessor :j3oj2
59
+ attr_accessor :dia_mm
60
+ attr_accessor :period_sec
61
+ attr_accessor :active
62
+ attr_accessor :not_orbital
63
+ attr_accessor :rcs_m2
64
+ attr_accessor :ep
65
+ attr_accessor :inclp
66
+ attr_accessor :nodep
67
+ attr_accessor :argpp
68
+ attr_accessor :mp
69
+ attr_accessor :isimp
70
+ attr_accessor :aycof
71
+ attr_accessor :con41
72
+ attr_accessor :cc1
73
+ attr_accessor :cc4
74
+ attr_accessor :cc5
75
+ attr_accessor :d2
76
+ attr_accessor :d3
77
+ attr_accessor :d4
78
+ attr_accessor :delmo
79
+ attr_accessor :eta
80
+ attr_accessor :argpdot
81
+ attr_accessor :omgcof
82
+ attr_accessor :sinmao
83
+ attr_accessor :t2cof
84
+ attr_accessor :t3cof
85
+ attr_accessor :t4cof
86
+ attr_accessor :t5cof
87
+ attr_accessor :x1mth2
88
+ attr_accessor :x7thm1
89
+ attr_accessor :mdot
90
+ attr_accessor :nodedot
91
+ attr_accessor :xlcof
92
+ attr_accessor :xmcof
93
+ attr_accessor :nodecf
94
+ attr_accessor :irez
95
+ attr_accessor :d2201
96
+ attr_accessor :d2211
97
+ attr_accessor :d3210
98
+ attr_accessor :d3222
99
+ attr_accessor :d4410
100
+ attr_accessor :d4422
101
+ attr_accessor :d5220
102
+ attr_accessor :d5232
103
+ attr_accessor :d5421
104
+ attr_accessor :d5433
105
+ attr_accessor :dedt
106
+ attr_accessor :del1
107
+ attr_accessor :del2
108
+ attr_accessor :del3
109
+ attr_accessor :didt
110
+ attr_accessor :dmdt
111
+ attr_accessor :dnodt
112
+ attr_accessor :domdt
113
+ attr_accessor :e3
114
+ attr_accessor :ee2
115
+ attr_accessor :peo
116
+ attr_accessor :pgho
117
+ attr_accessor :pho
118
+ attr_accessor :pinco
119
+ attr_accessor :plo
120
+ attr_accessor :se2
121
+ attr_accessor :se3
122
+ attr_accessor :sgh2
123
+ attr_accessor :sgh3
124
+ attr_accessor :sgh4
125
+ attr_accessor :sh2
126
+ attr_accessor :sh3
127
+ attr_accessor :si2
128
+ attr_accessor :si3
129
+ attr_accessor :sl2
130
+ attr_accessor :sl3
131
+ attr_accessor :sl4
132
+ attr_accessor :gsto
133
+ attr_accessor :xfact
134
+ attr_accessor :xgh2
135
+ attr_accessor :xgh3
136
+ attr_accessor :xgh4
137
+ attr_accessor :xh2
138
+ attr_accessor :xh3
139
+ attr_accessor :xi2
140
+ attr_accessor :xi3
141
+ attr_accessor :xl2
142
+ attr_accessor :xl3
143
+ attr_accessor :xl4
144
+ attr_accessor :xlamo
145
+ attr_accessor :zmol
146
+ attr_accessor :zmos
147
+ attr_accessor :atime
148
+ attr_accessor :xli
149
+ attr_accessor :xni
150
+ attr_accessor :snodm
151
+ attr_accessor :cnodm
152
+ attr_accessor :sinim
153
+ attr_accessor :cosim
154
+ attr_accessor :sinomm
155
+ attr_accessor :cosomm
156
+ attr_accessor :day
157
+ attr_accessor :emsq
158
+ attr_accessor :gam
159
+ attr_accessor :rtemsq
160
+ attr_accessor :s1
161
+ attr_accessor :s2
162
+ attr_accessor :s3
163
+ attr_accessor :s4
164
+ attr_accessor :s5
165
+ attr_accessor :s6
166
+ attr_accessor :s7
167
+ attr_accessor :ss1
168
+ attr_accessor :ss2
169
+ attr_accessor :ss3
170
+ attr_accessor :ss4
171
+ attr_accessor :ss5
172
+ attr_accessor :ss6
173
+ attr_accessor :ss7
174
+ attr_accessor :sz1
175
+ attr_accessor :sz2
176
+ attr_accessor :sz3
177
+ attr_accessor :sz11
178
+ attr_accessor :sz12
179
+ attr_accessor :sz13
180
+ attr_accessor :sz21
181
+ attr_accessor :sz22
182
+ attr_accessor :sz23
183
+ attr_accessor :sz31
184
+ attr_accessor :sz32
185
+ attr_accessor :sz33
186
+ attr_accessor :z1
187
+ attr_accessor :z2
188
+ attr_accessor :z3
189
+ attr_accessor :z11
190
+ attr_accessor :z12
191
+ attr_accessor :z13
192
+ attr_accessor :z21
193
+ attr_accessor :z22
194
+ attr_accessor :z23
195
+ attr_accessor :z31
196
+ attr_accessor :z32
197
+ attr_accessor :z33
198
+ attr_accessor :argpm
199
+ attr_accessor :inclm
200
+ attr_accessor :nodem
201
+ attr_accessor :dndt
202
+ attr_accessor :eccsq
203
+ attr_accessor :ainv
204
+ attr_accessor :ao
205
+ attr_accessor :con42
206
+ attr_accessor :cosio
207
+ attr_accessor :cosio2
208
+ attr_accessor :omeosq
209
+ attr_accessor :posq
210
+ attr_accessor :rp
211
+ attr_accessor :rteosq
212
+ attr_accessor :sinio
212
213
 
213
- def initialize
214
- @whichconst = 2 #SGP4.wgs72
215
- @satnum = 0
216
- @epochyr = 0
217
- @epochtynumrev = 0
218
- @error = 0
219
- @operationmode = 0
220
- @init = 0
221
- @method = 0
222
- @a = 0.0
223
- @altp = 0.0
224
- @alta = 0.0
225
- @epochdays = 0.0
226
- @jdsatepoch = 0.0
227
- @jdsatepochF = 0.0
228
- @nddot = 0.0
229
- @ndot = 0.0
230
- @bstar = 0.0
231
- @rcse = 0.0
232
- @inclo = 0.0
233
- @nodeo = 0.0
234
- @ecco = 0.0
235
- @argpo = 0.0
236
- @mo = 0.0
237
- @no_kozai = 0.0
214
+ def initialize
215
+ @whichconst = 2 #SGP4.wgs72
216
+ @satnum = 0
217
+ @epochyr = 0
218
+ @epochtynumrev = 0
219
+ @error = 0
220
+ @operationmode = 0
221
+ @init = 0
222
+ @method = 0
223
+ @a = 0.0
224
+ @altp = 0.0
225
+ @alta = 0.0
226
+ @epochdays = 0.0
227
+ @jdsatepoch = 0.0
228
+ @jdsatepochF = 0.0
229
+ @nddot = 0.0
230
+ @ndot = 0.0
231
+ @bstar = 0.0
232
+ @rcse = 0.0
233
+ @inclo = 0.0
234
+ @nodeo = 0.0
235
+ @ecco = 0.0
236
+ @argpo = 0.0
237
+ @mo = 0.0
238
+ @no_kozai = 0.0
238
239
 
239
- # sgp4fix add new variables from tle
240
- @classification = "U"
241
- @intldesg = ""
242
- @ephtype = 0
243
- @elnum = 0
244
- @revnum = 0
240
+ # sgp4fix add new variables from tle
241
+ @classification = "U"
242
+ @intldesg = ""
243
+ @ephtype = 0
244
+ @elnum = 0
245
+ @revnum = 0
245
246
 
246
- # sgp4fix add unkozai'd variable
247
- @no_unkozai = 0.0
248
- # sgp4fix add singly averaged variables
249
- @am = 0.0
250
- @em = 0.0
251
- @im = 0.0
252
- @Om = 0.0
253
- @om = 0.0
254
- @mm = 0.0
255
- @nm = 0.0
256
- @t = 0.0
247
+ # sgp4fix add unkozai'd variable
248
+ @no_unkozai = 0.0
249
+ # sgp4fix add singly averaged variables
250
+ @am = 0.0
251
+ @em = 0.0
252
+ @im = 0.0
253
+ @Om = 0.0
254
+ @om = 0.0
255
+ @mm = 0.0
256
+ @nm = 0.0
257
+ @t = 0.0
257
258
 
258
- # sgp4fix add constant parameters to eliminate mutliple calls during execution
259
- @tumin = 0.0
260
- @mu = 0.0
261
- @radiusearthkm = 0.0
262
- @xke = 0.0
263
- @j2 = 0.0
264
- @j3 = 0.0
265
- @j4 = 0.0
266
- @j3oj2 = 0.0
259
+ # sgp4fix add constant parameters to eliminate mutliple calls during execution
260
+ @tumin = 0.0
261
+ @mu = 0.0
262
+ @radiusearthkm = 0.0
263
+ @xke = 0.0
264
+ @j2 = 0.0
265
+ @j3 = 0.0
266
+ @j4 = 0.0
267
+ @j3oj2 = 0.0
267
268
 
268
- # Additional elements to capture relevant TLE and object information:
269
- @dia_mm = 0 # RSO dia in mm
270
- @period_sec = 0.0 # Period in seconds
271
- @active = 0 # "Active S/C" flag (0=n, 1=y)
272
- @not_orbital = 0 # "Orbiting S/C" flag (0=n, 1=y)
273
- @rcs_m2 = 0.0 # "RCS (m^2)" storage
269
+ # Additional elements to capture relevant TLE and object information:
270
+ @dia_mm = 0 # RSO dia in mm
271
+ @period_sec = 0.0 # Period in seconds
272
+ @active = 0 # "Active S/C" flag (0=n, 1=y)
273
+ @not_orbital = 0 # "Orbiting S/C" flag (0=n, 1=y)
274
+ @rcs_m2 = 0.0 # "RCS (m^2)" storage
274
275
 
275
- # temporary variables because the original authors call the same method with different variables
276
- @ep = 0.0
277
- @inclp = 0.0
278
- @nodep = 0.0
279
- @argpp = 0.0
280
- @mp = 0.0
276
+ # temporary variables because the original authors call the same method with different variables
277
+ @ep = 0.0
278
+ @inclp = 0.0
279
+ @nodep = 0.0
280
+ @argpp = 0.0
281
+ @mp = 0.0
281
282
 
282
- @isimp = 0
283
- @aycof = 0.0
284
- @con41 = 0.0
285
- @cc1 = 0.0
286
- @cc4 = 0.0
287
- @cc5 = 0.0
288
- @d2 = 0.0
289
- @d3 = 0.0
290
- @d4 = 0.0
291
- @delmo = 0.0
292
- @eta = 0.0
293
- @argpdot = 0.0
294
- @omgcof = 0.0
295
- @sinmao = 0.0
296
- @t2cof = 0.0
297
- @t3cof = 0.0
298
- @t4cof = 0.0
299
- @t5cof = 0.0
300
- @x1mth2 = 0.0
301
- @x7thm1 = 0.0
302
- @mdot = 0.0
303
- @nodedot = 0.0
304
- @xlcof = 0.0
305
- @xmcof = 0.0
306
- @nodecf = 0.0
283
+ @isimp = 0
284
+ @aycof = 0.0
285
+ @con41 = 0.0
286
+ @cc1 = 0.0
287
+ @cc4 = 0.0
288
+ @cc5 = 0.0
289
+ @d2 = 0.0
290
+ @d3 = 0.0
291
+ @d4 = 0.0
292
+ @delmo = 0.0
293
+ @eta = 0.0
294
+ @argpdot = 0.0
295
+ @omgcof = 0.0
296
+ @sinmao = 0.0
297
+ @t2cof = 0.0
298
+ @t3cof = 0.0
299
+ @t4cof = 0.0
300
+ @t5cof = 0.0
301
+ @x1mth2 = 0.0
302
+ @x7thm1 = 0.0
303
+ @mdot = 0.0
304
+ @nodedot = 0.0
305
+ @xlcof = 0.0
306
+ @xmcof = 0.0
307
+ @nodecf = 0.0
307
308
 
308
- # deep space
309
- @irez = 0
310
- @d2201 = 0.0
311
- @d2211 = 0.0
312
- @d3210 = 0.0
313
- @d3222 = 0.0
314
- @d4410 = 0.0
315
- @d4422 = 0.0
316
- @d5220 = 0.0
317
- @d5232 = 0.0
318
- @d5421 = 0.0
319
- @d5433 = 0.0
320
- @dedt = 0.0
321
- @del1 = 0.0
322
- @del2 = 0.0
323
- @del3 = 0.0
324
- @didt = 0.0
325
- @dmdt = 0.0
326
- @dnodt = 0.0
327
- @domdt = 0.0
328
- @e3 = 0.0
329
- @ee2 = 0.0
330
- @peo = 0.0
331
- @pgho = 0.0
332
- @pho = 0.0
333
- @pinco = 0.0
334
- @plo = 0.0
335
- @se2 = 0.0
336
- @se3 = 0.0
337
- @sgh2 = 0.0
338
- @sgh3 = 0.0
339
- @sgh4 = 0.0
340
- @sh2 = 0.0
341
- @sh3 = 0.0
342
- @si2 = 0.0
343
- @si3 = 0.0
344
- @sl2 = 0.0
345
- @sl3 = 0.0
346
- @sl4 = 0.0
347
- @gsto = 0.0
348
- @xfact = 0.0
349
- @xgh2 = 0.0
350
- @xgh3 = 0.0
351
- @xgh4 = 0.0
352
- @xh2 = 0.0
353
- @xh3 = 0.0
354
- @xi2 = 0.0
355
- @xi3 = 0.0
356
- @xl2 = 0.0
357
- @xl3 = 0.0
358
- @xl4 = 0.0
359
- @xlamo = 0.0
360
- @zmol = 0.0
361
- @zmos = 0.0
362
- @atime = 0.0
363
- @xli = 0.0
364
- @xni = 0.0
365
- @snodm = 0.0
366
- @cnodm = 0.0
367
- @sinim = 0.0
368
- @cosim = 0.0
369
- @sinomm = 0.0
370
- @cosomm = 0.0
371
- @day = 0.0
372
- @emsq = 0.0
373
- @gam = 0.0
374
- @rtemsq = 0.0
375
- @s1 = 0.0
376
- @s2 = 0.0
377
- @s3 = 0.0
378
- @s4 = 0.0
379
- @s5 = 0.0
380
- @s6 = 0.0
381
- @s7 = 0.0
382
- @ss1 = 0.0
383
- @ss2 = 0.0
384
- @ss3 = 0.0
385
- @ss4 = 0.0
386
- @ss5 = 0.0
387
- @ss6 = 0.0
388
- @ss7 = 0.0
389
- @sz1 = 0.0
390
- @sz2 = 0.0
391
- @sz3 = 0.0
392
- @sz11 = 0.0
393
- @sz12 = 0.0
394
- @sz13 = 0.0
395
- @sz21 = 0.0
396
- @sz22 = 0.0
397
- @sz23 = 0.0
398
- @sz31 = 0.0
399
- @sz32 = 0.0
400
- @sz33 = 0.0
401
- @z1 = 0.0
402
- @z2 = 0.0
403
- @z3 = 0.0
404
- @z11 = 0.0
405
- @z12 = 0.0
406
- @z13 = 0.0
407
- @z21 = 0.0
408
- @z22 = 0.0
409
- @z23 = 0.0
410
- @z31 = 0.0
411
- @z32 = 0.0
412
- @z33 = 0.0
413
- @argpm = 0.0
414
- @inclm = 0.0
415
- @nodem = 0.0
416
- @dndt = 0.0
417
- @eccsq = 0.0
309
+ # deep space
310
+ @irez = 0
311
+ @d2201 = 0.0
312
+ @d2211 = 0.0
313
+ @d3210 = 0.0
314
+ @d3222 = 0.0
315
+ @d4410 = 0.0
316
+ @d4422 = 0.0
317
+ @d5220 = 0.0
318
+ @d5232 = 0.0
319
+ @d5421 = 0.0
320
+ @d5433 = 0.0
321
+ @dedt = 0.0
322
+ @del1 = 0.0
323
+ @del2 = 0.0
324
+ @del3 = 0.0
325
+ @didt = 0.0
326
+ @dmdt = 0.0
327
+ @dnodt = 0.0
328
+ @domdt = 0.0
329
+ @e3 = 0.0
330
+ @ee2 = 0.0
331
+ @peo = 0.0
332
+ @pgho = 0.0
333
+ @pho = 0.0
334
+ @pinco = 0.0
335
+ @plo = 0.0
336
+ @se2 = 0.0
337
+ @se3 = 0.0
338
+ @sgh2 = 0.0
339
+ @sgh3 = 0.0
340
+ @sgh4 = 0.0
341
+ @sh2 = 0.0
342
+ @sh3 = 0.0
343
+ @si2 = 0.0
344
+ @si3 = 0.0
345
+ @sl2 = 0.0
346
+ @sl3 = 0.0
347
+ @sl4 = 0.0
348
+ @gsto = 0.0
349
+ @xfact = 0.0
350
+ @xgh2 = 0.0
351
+ @xgh3 = 0.0
352
+ @xgh4 = 0.0
353
+ @xh2 = 0.0
354
+ @xh3 = 0.0
355
+ @xi2 = 0.0
356
+ @xi3 = 0.0
357
+ @xl2 = 0.0
358
+ @xl3 = 0.0
359
+ @xl4 = 0.0
360
+ @xlamo = 0.0
361
+ @zmol = 0.0
362
+ @zmos = 0.0
363
+ @atime = 0.0
364
+ @xli = 0.0
365
+ @xni = 0.0
366
+ @snodm = 0.0
367
+ @cnodm = 0.0
368
+ @sinim = 0.0
369
+ @cosim = 0.0
370
+ @sinomm = 0.0
371
+ @cosomm = 0.0
372
+ @day = 0.0
373
+ @emsq = 0.0
374
+ @gam = 0.0
375
+ @rtemsq = 0.0
376
+ @s1 = 0.0
377
+ @s2 = 0.0
378
+ @s3 = 0.0
379
+ @s4 = 0.0
380
+ @s5 = 0.0
381
+ @s6 = 0.0
382
+ @s7 = 0.0
383
+ @ss1 = 0.0
384
+ @ss2 = 0.0
385
+ @ss3 = 0.0
386
+ @ss4 = 0.0
387
+ @ss5 = 0.0
388
+ @ss6 = 0.0
389
+ @ss7 = 0.0
390
+ @sz1 = 0.0
391
+ @sz2 = 0.0
392
+ @sz3 = 0.0
393
+ @sz11 = 0.0
394
+ @sz12 = 0.0
395
+ @sz13 = 0.0
396
+ @sz21 = 0.0
397
+ @sz22 = 0.0
398
+ @sz23 = 0.0
399
+ @sz31 = 0.0
400
+ @sz32 = 0.0
401
+ @sz33 = 0.0
402
+ @z1 = 0.0
403
+ @z2 = 0.0
404
+ @z3 = 0.0
405
+ @z11 = 0.0
406
+ @z12 = 0.0
407
+ @z13 = 0.0
408
+ @z21 = 0.0
409
+ @z22 = 0.0
410
+ @z23 = 0.0
411
+ @z31 = 0.0
412
+ @z32 = 0.0
413
+ @z33 = 0.0
414
+ @argpm = 0.0
415
+ @inclm = 0.0
416
+ @nodem = 0.0
417
+ @dndt = 0.0
418
+ @eccsq = 0.0
418
419
 
419
- # for initl
420
- @ainv = 0.0
421
- @ao = 0.0
422
- @con42 = 0.0
423
- @cosio = 0.0
424
- @cosio2 = 0.0
425
- @omeosq = 0.0
426
- @posq = 0.0
427
- @rp = 0.0
428
- @rteosq = 0.0
429
- @sinio = 0.0
420
+ # for initl
421
+ @ainv = 0.0
422
+ @ao = 0.0
423
+ @con42 = 0.0
424
+ @cosio = 0.0
425
+ @cosio2 = 0.0
426
+ @omeosq = 0.0
427
+ @posq = 0.0
428
+ @rp = 0.0
429
+ @rteosq = 0.0
430
+ @sinio = 0.0
431
+ end
430
432
  end
431
433
  end
432
434
  end