-
Notifications
You must be signed in to change notification settings - Fork 3
Expand file tree
/
Copy pathelp2000orb.js
More file actions
360 lines (321 loc) · 27.9 KB
/
elp2000orb.js
File metadata and controls
360 lines (321 loc) · 27.9 KB
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
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
// from http://www.mathworks.com/matlabcentral/fileexchange/39130-orbital-elements-of-the-moon/content/elp2000.m
// http://www.mathworks.com/matlabcentral/mlc-downloads/downloads/submissions/39130/versions/2/previews/elp2000.m/index.html
// Terms also found in https://github.com/bragboy/horoscope/blob/master/lib/horoscope/planets/nod.rb
// This function calculates the osculating classical orbital elements of the
// Moon in the mean ecliptic and mean equinox of date coordinate system. It is
// based on the book Lunar Tables and Programs From 4000 B.C. TO A.D. 8000 by
// Michelle Chapront-Touze and Jean Chapront. This book and its optional
// companion software are available from Willmann-Bell (www.willbell.com).
/*
% osculating orbital elements of the moon
% mean ecliptic and mean equinox of date
% input
% tjd = tdt julian day
% output
% oev(1) = semimajor axis (kilometers)
% oev(2) = orbital eccentricity (non-dimensional)
% (0 <= eccentricity < 1)
% oev(3) = orbital inclination (radians)
% (0 <= inclination <= pi)
% oev(4) = argument of perigee (radians)
% (0 <= argument of perigee <= 2 pi)
% oev(5) = right ascension of ascending node (radians)
% (0 <= raan <= 2 pi)
% oev(6) = true anomaly (radians)
% (0 <= true anomaly <= 2 pi)
*/
function elp2000orb(theory, jy2k, elems, addGM, addEpoch, off)
{
off = off || 0;
jy2k = jy2k || 0;
elems = elems || [];
var t = jy2k / 100;
var t2 = t * t;
var t3 = t * t2;
var t4 = t * t3;
var t5 = t * t4;
var t2e4 = t2 * 1e-4;
var t3e6 = t3 * 1e-6;
var t4e8 = t4 * 1e-8;
// % conversion factor - degrees to radians
var dtr = Math.PI / 180;
var sa = 3400.4 * Math.cos(dtr * (235.7004 + 890534.2230 * t - 32.601 * t2e4 + 3.664 * t3e6 - 1.769 * t4e8))
- 635.6 * Math.cos(dtr * (100.7370 + 413335.3554 * t - 122.571 * t2e4 - 10.684 * t3e6 + 5.028 * t4e8))
- 235.6 * Math.cos(dtr * (134.9634 + 477198.8676 * t + 89.970 * t2e4 + 14.348 * t3e6 - 6.797 * t4e8))
+ 218.1 * Math.cos(dtr * (238.1713 + 854535.1727 * t - 31.065 * t2e4 + 3.623 * t3e6 - 1.769 * t4e8))
+ 181.0 * Math.cos(dtr * (10.6638 + 1367733.0907 * t + 57.370 * t2e4 + 18.011 * t3e6 - 8.566 * t4e8))
- 39.9 * Math.cos(dtr * (103.2079 + 377336.3051 * t - 121.035 * t2e4 - 10.724 * t3e6 + 5.028 * t4e8))
- 38.4 * Math.cos(dtr * (233.2295 + 926533.2733 * t - 34.136 * t2e4 + 3.705 * t3e6 - 1.769 * t4e8))
+ 33.8 * Math.cos(dtr * (336.4374 + 1303869.5784 * t - 155.171 * t2e4 - 7.020 * t3e6 + 3.259 * t4e8))
+ 28.8 * Math.cos(dtr * (111.4008 + 1781068.4461 * t - 65.201 * t2e4 + 7.328 * t3e6 - 3.538 * t4e8))
+ 12.6 * Math.cos(dtr * (13.1347 + 1331734.0404 * t + 58.906 * t2e4 + 17.971 * t3e6 - 8.566 * t4e8))
+ 11.4 * Math.cos(dtr * (186.5442 + 966404.0351 * t - 68.058 * t2e4 - 0.567 * t3e6 + 0.232 * t4e8))
- 11.1 * Math.cos(dtr * (222.5657 - 441199.8173 * t - 91.506 * t2e4 - 14.307 * t3e6 + 6.797 * t4e8))
- 10.2 * Math.cos(dtr * (269.9268 + 954397.7353 * t + 179.941 * t2e4 + 28.695 * t3e6 - 13.594 * t4e8))
+ 9.7 * Math.cos(dtr * (145.6272 + 1844931.9583 * t + 147.340 * t2e4 + 32.359 * t3e6 - 15.363 * t4e8))
+ 9.6 * Math.cos(dtr * (240.6422 + 818536.1225 * t - 29.529 * t2e4 + 3.582 * t3e6 - 1.769 * t4e8))
+ 8.0 * Math.cos(dtr * (297.8502 + 445267.1115 * t - 16.300 * t2e4 + 1.832 * t3e6 - 0.884 * t4e8))
- 6.2 * Math.cos(dtr * (132.4925 + 513197.9179 * t + 88.434 * t2e4 + 14.388 * t3e6 - 6.797 * t4e8))
+ 6.0 * Math.cos(dtr * (173.5506 + 1335801.3346 * t - 48.901 * t2e4 + 5.496 * t3e6 - 2.653 * t4e8))
+ 3.7 * Math.cos(dtr * (113.8717 + 1745069.3958 * t - 63.665 * t2e4 + 7.287 * t3e6 - 3.538 * t4e8))
+ 3.6 * Math.cos(dtr * (338.9083 + 1267870.5281 * t - 153.636 * t2e4 - 7.061 * t3e6 + 3.259 * t4e8))
+ 3.2 * Math.cos(dtr * (246.3642 + 2258267.3137 * t + 24.769 * t2e4 + 21.675 * t3e6 - 10.335 * t4e8))
- 3.0 * Math.cos(dtr * (8.1929 + 1403732.1410 * t + 55.834 * t2e4 + 18.052 * t3e6 - 8.566 * t4e8))
+ 2.3 * Math.cos(dtr * (98.2661 + 449334.4057 * t - 124.107 * t2e4 - 10.643 * t3e6 + 5.028 * t4e8))
- 2.2 * Math.cos(dtr * (357.5291 + 35999.0503 * t - 1.536 * t2e4 + 0.041 * t3e6 + 0.000 * t4e8))
- 2.0 * Math.cos(dtr * (38.5872 + 858602.4669 * t - 138.871 * t2e4 - 8.852 * t3e6 + 4.144 * t4e8))
- 1.8 * Math.cos(dtr * (105.6788 + 341337.2548 * t - 119.499 * t2e4 - 10.765 * t3e6 + 5.028 * t4e8))
- 1.7 * Math.cos(dtr * (201.4740 + 826670.7108 * t - 245.142 * t2e4 - 21.367 * t3e6 + 10.057 * t4e8))
+ 1.6 * Math.cos(dtr * (184.1196 + 401329.0556 * t + 125.428 * t2e4 + 18.579 * t3e6 - 8.798 * t4e8))
- 1.4 * Math.cos(dtr * (308.4192 - 489205.1674 * t + 158.029 * t2e4 + 14.915 * t3e6 - 7.029 * t4e8))
+ 1.3 * Math.cos(dtr * (325.7736 - 63863.5122 * t - 212.541 * t2e4 - 25.031 * t3e6 + 11.826 * t4e8));
// % semimajor axis
var sapp = - 0.55 * Math.cos(dtr * (238.2 + 854535.2 * t))
+ 0.10 * Math.cos(dtr * (103.2 + 377336.3 * t))
+ 0.10 * Math.cos(dtr * (233.2 + 926533.3 * t));
var sma = 383397.6 + sa + sapp * t;
// % orbital eccentricity
var se = 0.014217 * Math.cos(dtr * (100.7370 + 413335.3554 * t - 122.571 * t2e4 - 10.684 * t3e6 + 5.028 * t4e8))
+ 0.008551 * Math.cos(dtr * (325.7736 - 63863.5122 * t - 212.541 * t2e4 - 25.031 * t3e6 + 11.826 * t4e8))
- 0.001383 * Math.cos(dtr * (134.9634 + 477198.8676 * t + 89.970 * t2e4 + 14.348 * t3e6 - 6.797 * t4e8))
+ 0.001353 * Math.cos(dtr * (10.6638 + 1367733.0907 * t + 57.370 * t2e4 + 18.011 * t3e6 - 8.566 * t4e8))
- 0.001146 * Math.cos(dtr * (66.5106 + 349471.8432 * t - 335.112 * t2e4 - 35.715 * t3e6 + 16.854 * t4e8))
- 0.000915 * Math.cos(dtr * (201.4740 + 826670.7108 * t - 245.142 * t2e4 - 21.367 * t3e6 + 10.057 * t4e8))
+ 0.000869 * Math.cos(dtr * (103.2079 + 377336.3051 * t - 121.035 * t2e4 - 10.724 * t3e6 + 5.028 * t4e8))
- 0.000628 * Math.cos(dtr * (235.7004 + 890534.2230 * t - 32.601 * t2e4 + 3.664 * t3e6 - 1.769 * t4e8))
- 0.000393 * Math.cos(dtr * (291.5472 - 127727.0245 * t - 425.082 * t2e4 - 50.062 * t3e6 + 23.651 * t4e8))
+ 0.000284 * Math.cos(dtr * (328.2445 - 99862.5625 * t - 211.005 * t2e4 - 25.072 * t3e6 + 11.826 * t4e8))
- 0.000278 * Math.cos(dtr * (162.8868 - 31931.7561 * t - 106.271 * t2e4 - 12.516 * t3e6 + 5.913 * t4e8))
- 0.000240 * Math.cos(dtr * (269.9268 + 954397.7353 * t + 179.941 * t2e4 + 28.695 * t3e6 - 13.594 * t4e8))
+ 0.000230 * Math.cos(dtr * (111.4008 + 1781068.4461 * t - 65.201 * t2e4 + 7.328 * t3e6 - 3.538 * t4e8))
+ 0.000229 * Math.cos(dtr * (167.2476 + 762807.1986 * t - 457.683 * t2e4 - 46.398 * t3e6 + 21.882 * t4e8))
- 0.000202 * Math.cos(dtr * ( 83.3826 - 12006.2998 * t + 247.999 * t2e4 + 29.262 * t3e6 - 13.826 * t4e8))
+ 0.000190 * Math.cos(dtr * (190.8102 - 541062.3799 * t - 302.511 * t2e4 - 39.379 * t3e6 + 18.623 * t4e8))
+ 0.000177 * Math.cos(dtr * (357.5291 + 35999.0503 * t - 1.536 * t2e4 + 0.041 * t3e6 + 0.000 * t4e8))
+ 0.000153 * Math.cos(dtr * (32.2842 + 285608.3309 * t - 547.653 * t2e4 - 60.746 * t3e6 + 28.679 * t4e8))
- 0.000137 * Math.cos(dtr * (44.8902 + 1431596.6029 * t + 269.911 * t2e4 + 43.043 * t3e6 - 20.392 * t4e8))
+ 0.000122 * Math.cos(dtr * (145.6272 + 1844931.9583 * t + 147.340 * t2e4 + 32.359 * t3e6 - 15.363 * t4e8))
+ 0.000116 * Math.cos(dtr * (302.2110 + 1240006.0662 * t - 367.713 * t2e4 - 32.051 * t3e6 + 15.085 * t4e8))
- 0.000111 * Math.cos(dtr * (203.9449 + 790671.6605 * t - 243.606 * t2e4 - 21.408 * t3e6 + 10.057 * t4e8))
- 0.000108 * Math.cos(dtr * (68.9815 + 313472.7929 * t - 333.576 * t2e4 - 35.756 * t3e6 + 16.854 * t4e8))
+ 0.000096 * Math.cos(dtr * (336.4374 + 1303869.5784 * t - 155.171 * t2e4 - 7.020 * t3e6 + 3.259 * t4e8))
- 0.000090 * Math.cos(dtr * (98.2661 + 449334.4057 * t - 124.107 * t2e4 - 10.643 * t3e6 + 5.028 * t4e8))
+ 0.000090 * Math.cos(dtr * (13.1347 + 1331734.0404 * t + 58.906 * t2e4 + 17.971 * t3e6 - 8.566 * t4e8))
+ 0.000056 * Math.cos(dtr * (55.8468 - 1018261.2475 * t - 392.482 * t2e4 - 53.726 * t3e6 + 25.420 * t4e8))
- 0.000056 * Math.cos(dtr * (238.1713 + 854535.1727 * t - 31.065 * t2e4 + 3.623 * t3e6 - 1.769 * t4e8))
+ 0.000052 * Math.cos(dtr * (308.4192 - 489205.1674 * t + 158.029 * t2e4 + 14.915 * t3e6 - 7.029 * t4e8))
- 0.000050 * Math.cos(dtr * (133.0212 + 698943.6863 * t - 670.224 * t2e4 - 71.429 * t3e6 + 33.708 * t4e8))
- 0.000049 * Math.cos(dtr * (267.9846 + 1176142.5540 * t - 580.254 * t2e4 - 57.082 * t3e6 + 26.911 * t4e8))
- 0.000049 * Math.cos(dtr * (184.1196 + 401329.0556 * t + 125.428 * t2e4 + 18.579 * t3e6 - 8.798 * t4e8))
- 0.000045 * Math.cos(dtr * (49.1562 - 75869.8120 * t + 35.458 * t2e4 + 4.231 * t3e6 - 2.001 * t4e8))
+ 0.000044 * Math.cos(dtr * (257.3208 - 191590.5367 * t - 637.623 * t2e4 - 75.093 * t3e6 + 35.477 * t4e8))
+ 0.000042 * Math.cos(dtr * (105.6788 + 341337.2548 * t - 119.499 * t2e4 - 10.765 * t3e6 + 5.028 * t4e8))
+ 0.000042 * Math.cos(dtr * (160.4159 + 4067.2942 * t - 107.806 * t2e4 - 12.475 * t3e6 + 5.913 * t4e8))
+ 0.000040 * Math.cos(dtr * (246.3642 + 2258267.3137 * t + 24.769 * t2e4 + 21.675 * t3e6 - 10.335 * t4e8))
- 0.000040 * Math.cos(dtr * (156.5838 - 604925.8921 * t - 515.053 * t2e4 - 64.410 * t3e6 + 30.448 * t4e8))
+ 0.000036 * Math.cos(dtr * (169.7185 + 726808.1483 * t - 456.147 * t2e4 - 46.439 * t3e6 + 21.882 * t4e8))
+ 0.000029 * Math.cos(dtr * (113.8717 + 1745069.3958 * t - 63.665 * t2e4 + 7.287 * t3e6 - 3.538 * t4e8))
- 0.000029 * Math.cos(dtr * (297.8502 + 445267.1115 * t - 16.300 * t2e4 + 1.832 * t3e6 - 0.884 * t4e8))
- 0.000028 * Math.cos(dtr * (294.0181 - 163726.0747 * t - 423.546 * t2e4 - 50.103 * t3e6 + 23.651 * t4e8))
+ 0.000027 * Math.cos(dtr * (263.6238 + 381403.5993 * t - 228.841 * t2e4 - 23.199 * t3e6 + 10.941 * t4e8))
- 0.000026 * Math.cos(dtr * (358.0578 + 221744.8187 * t - 760.194 * t2e4 - 85.777 * t3e6 + 40.505 * t4e8))
- 0.000026 * Math.cos(dtr * (8.1929 + 1403732.1410 * t + 55.834 * t2e4 + 18.052 * t3e6 - 8.566 * t4e8));
var sedp = -0.0022 * Math.cos(dtr * (103.2 + 377336.3 * t));
var ecc = 0.055544 + se + 1e-3 * t * sedp;
// % sine of half the inclination
var sg = 0.0011776 * Math.cos(dtr * (49.1562 - 75869.8120 * t + 35.458 * t2e4 + 4.231 * t3e6 - 2.001 * t4e8))
- 0.0000971 * Math.cos(dtr * (235.7004 + 890534.2230 * t - 32.601 * t2e4 + 3.664 * t3e6 - 1.769 * t4e8))
+ 0.0000908 * Math.cos(dtr * (186.5442 + 966404.0351 * t - 68.058 * t2e4 - 0.567 * t3e6 + 0.232 * t4e8))
+ 0.0000623 * Math.cos(dtr * (83.3826 - 12006.2998 * t + 247.999 * t2e4 + 29.262 * t3e6 - 13.826 * t4e8))
+ 0.0000483 * Math.cos(dtr * (51.6271 - 111868.8623 * t + 36.994 * t2e4 + 4.190 * t3e6 - 2.001 * t4e8))
+ 0.0000348 * Math.cos(dtr * (100.7370 + 413335.3554 * t - 122.571 * t2e4 - 10.684 * t3e6 + 5.028 * t4e8))
- 0.0000316 * Math.cos(dtr * (308.4192 - 489205.1674 * t + 158.029 * t2e4 + 14.915 * t3e6 - 7.029 * t4e8))
- 0.0000253 * Math.cos(dtr * (46.6853 - 39870.7617 * t + 33.922 * t2e4 + 4.272 * t3e6 - 2.001 * t4e8))
- 0.0000141 * Math.cos(dtr * (274.1928 - 553068.6797 * t - 54.513 * t2e4 - 10.116 * t3e6 + 4.797 * t4e8))
+ 0.0000127 * Math.cos(dtr * (325.7736 - 63863.5122 * t - 212.541 * t2e4 - 25.031 * t3e6 + 11.826 * t4e8))
+ 0.0000117 * Math.cos(dtr * (184.1196 + 401329.0556 * t + 125.428 * t2e4 + 18.579 * t3e6 - 8.798 * t4e8))
- 0.0000078 * Math.cos(dtr * (98.3124 - 151739.6240 * t + 70.916 * t2e4 + 8.462 * t3e6 - 4.001 * t4e8))
- 0.0000063 * Math.cos(dtr * (238.1713 + 854535.1727 * t - 31.065 * t2e4 + 3.623 * t3e6 - 1.769 * t4e8))
+ 0.0000063 * Math.cos(dtr * (134.9634 + 477198.8676 * t + 89.970 * t2e4 + 14.348 * t3e6 - 6.797 * t4e8))
+ 0.0000036 * Math.cos(dtr * (321.5076 + 1443602.9027 * t + 21.912 * t2e4 + 13.780 * t3e6 - 6.566 * t4e8))
- 0.0000035 * Math.cos(dtr * (10.6638 + 1367733.0907 * t + 57.370 * t2e4 + 18.011 * t3e6 - 8.566 * t4e8))
+ 0.0000024 * Math.cos(dtr * (149.8932 + 337465.5434 * t - 87.113 * t2e4 - 6.453 * t3e6 + 3.028 * t4e8))
+ 0.0000024 * Math.cos(dtr * (170.9849 - 930404.9848 * t + 66.523 * t2e4 + 0.608 * t3e6 - 0.232 * t4e8));
var sgp = - 0.0203 * Math.cos(dtr * (125.0 - 1934.1 * t))
+ 0.0034 * Math.cos(dtr * (220.2 - 1935.5 * t));
var gamma = 0.0449858 + sg + 1e-3 * sgp;
// % longitude of perigee
var sp = - 15.448 * Math.sin(dtr * (100.7370 + 413335.3554 * t - 122.571 * t2e4 - 10.684 * t3e6 + 5.028 * t4e8))
- 9.642 * Math.sin(dtr * (325.7736 - 63863.5122 * t - 212.541 * t2e4 - 25.031 * t3e6 + 11.826 * t4e8))
- 2.721 * Math.sin(dtr * (134.9634 + 477198.8676 * t + 89.970 * t2e4 + 14.348 * t3e6 - 6.797 * t4e8))
+ 2.607 * Math.sin(dtr * (66.5106 + 349471.8432 * t - 335.112 * t2e4 - 35.715 * t3e6 + 16.854 * t4e8))
+ 2.085 * Math.sin(dtr * (201.4740 + 826670.7108 * t - 245.142 * t2e4 - 21.367 * t3e6 + 10.057 * t4e8))
+ 1.477 * Math.sin(dtr * (10.6638 + 1367733.0907 * t + 57.370 * t2e4 + 18.011 * t3e6 - 8.566 * t4e8))
+ 0.968 * Math.sin(dtr * (291.5472 - 127727.0245 * t - 425.082 * t2e4 - 50.062 * t3e6 + 23.651 * t4e8))
- 0.949 * Math.sin(dtr * (103.2079 + 377336.3051 * t - 121.035 * t2e4 - 10.724 * t3e6 + 5.028 * t4e8))
- 0.703 * Math.sin(dtr * (167.2476 + 762807.1986 * t - 457.683 * t2e4 - 46.398 * t3e6 + 21.882 * t4e8))
- 0.660 * Math.sin(dtr * (235.7004 + 890534.2230 * t - 32.601 * t2e4 + 3.664 * t3e6 - 1.769 * t4e8))
- 0.577 * Math.sin(dtr * (190.8102 - 541062.3799 * t - 302.511 * t2e4 - 39.379 * t3e6 + 18.623 * t4e8))
- 0.524 * Math.sin(dtr * (269.9268 + 954397.7353 * t + 179.941 * t2e4 + 28.695 * t3e6 - 13.594 * t4e8))
- 0.482 * Math.sin(dtr * (32.2842 + 285608.3309 * t - 547.653 * t2e4 - 60.746 * t3e6 + 28.679 * t4e8))
+ 0.452 * Math.sin(dtr * (357.5291 + 35999.0503 * t - 1.536 * t2e4 + 0.041 * t3e6 + 0.000 * t4e8))
- 0.381 * Math.sin(dtr * (302.2110 + 1240006.0662 * t - 367.713 * t2e4 - 32.051 * t3e6 + 15.085 * t4e8))
- 0.342 * Math.sin(dtr * (328.2445 - 99862.5625 * t - 211.005 * t2e4 - 25.072 * t3e6 + 11.826 * t4e8))
- 0.312 * Math.sin(dtr * (44.8902 + 1431596.6029 * t + 269.911 * t2e4 + 43.043 * t3e6 - 20.392 * t4e8))
+ 0.282 * Math.sin(dtr * (162.8868 - 31931.7561 * t - 106.271 * t2e4 - 12.516 * t3e6 + 5.913 * t4e8))
+ 0.255 * Math.sin(dtr * (203.9449 + 790671.6605 * t - 243.606 * t2e4 - 21.408 * t3e6 + 10.057 * t4e8))
+ 0.252 * Math.sin(dtr * (68.9815 + 313472.7929 * t - 333.576 * t2e4 - 35.756 * t3e6 + 16.854 * t4e8))
- 0.211 * Math.sin(dtr * (83.3826 - 12006.2998 * t + 247.999 * t2e4 + 29.262 * t3e6 - 13.826 * t4e8))
+ 0.193 * Math.sin(dtr * (267.9846 + 1176142.5540 * t - 580.254 * t2e4 - 57.082 * t3e6 + 26.911 * t4e8))
+ 0.191 * Math.sin(dtr * (133.0212 + 698943.6863 * t - 670.224 * t2e4 - 71.429 * t3e6 + 33.708 * t4e8))
- 0.184 * Math.sin(dtr * (55.8468 - 1018261.2475 * t - 392.482 * t2e4 - 53.726 * t3e6 + 25.420 * t4e8))
+ 0.182 * Math.sin(dtr * (145.6272 + 1844931.9583 * t + 147.340 * t2e4 + 32.359 * t3e6 - 15.363 * t4e8))
- 0.158 * Math.sin(dtr * (257.3208 - 191590.5367 * t - 637.623 * t2e4 - 75.093 * t3e6 + 35.477 * t4e8))
+ 0.148 * Math.sin(dtr * (156.5838 - 604925.8921 * t - 515.053 * t2e4 - 64.410 * t3e6 + 30.448 * t4e8))
- 0.111 * Math.sin(dtr * (169.7185 + 726808.1483 * t - 456.147 * t2e4 - 46.439 * t3e6 + 21.882 * t4e8))
+ 0.101 * Math.sin(dtr * (13.1347 + 1331734.0404 * t + 58.906 * t2e4 + 17.971 * t3e6 - 8.566 * t4e8))
+ 0.100 * Math.sin(dtr * (358.0578 + 221744.8187 * t - 760.194 * t2e4 - 85.777 * t3e6 + 40.505 * t4e8))
+ 0.087 * Math.sin(dtr * (98.2661 + 449334.4057 * t - 124.107 * t2e4 - 10.643 * t3e6 + 5.028 * t4e8))
+ 0.080 * Math.sin(dtr * (42.9480 + 1653341.4216 * t - 490.283 * t2e4 - 42.734 * t3e6 + 20.113 * t4e8))
+ 0.080 * Math.sin(dtr * (222.5657 - 441199.8173 * t - 91.506 * t2e4 - 14.307 * t3e6 + 6.797 * t4e8))
+ 0.077 * Math.sin(dtr * (294.0181 - 163726.0747 * t - 423.546 * t2e4 - 50.103 * t3e6 + 23.651 * t4e8))
- 0.073 * Math.sin(dtr * (280.8834 - 1495460.1151 * t - 482.452 * t2e4 - 68.074 * t3e6 + 32.217 * t4e8))
- 0.071 * Math.sin(dtr * (304.6819 + 1204007.0159 * t - 366.177 * t2e4 - 32.092 * t3e6 + 15.085 * t4e8))
- 0.069 * Math.sin(dtr * (233.7582 + 1112279.0417 * t - 792.795 * t2e4 - 82.113 * t3e6 + 38.736 * t4e8))
- 0.067 * Math.sin(dtr * (34.7551 + 249609.2807 * t - 546.117 * t2e4 - 60.787 * t3e6 + 28.679 * t4e8))
- 0.067 * Math.sin(dtr * (263.6238 + 381403.5993 * t - 228.841 * t2e4 - 23.199 * t3e6 + 10.941 * t4e8))
+ 0.055 * Math.sin(dtr * (21.6203 - 1082124.7597 * t - 605.023 * t2e4 - 78.757 * t3e6 + 37.246 * t4e8))
+ 0.055 * Math.sin(dtr * (308.4192 - 489205.1674 * t + 158.029 * t2e4 + 14.915 * t3e6 -7.029 * t4e8))
- 0.054 * Math.sin(dtr * (8.7216 + 1589477.9094 * t - 702.824 * t2e4 - 67.766 * t3e6 + 31.939 * t4e8))
- 0.052 * Math.sin(dtr * (179.8536 + 1908795.4705 * t + 359.881 * t2e4 + 57.390 * t3e6 - 27.189 * t4e8))
- 0.050 * Math.sin(dtr * (98.7948 + 635080.1741 * t - 882.765 * t2e4 - 96.461 * t3e6 + 45.533 * t4e8))
- 0.049 * Math.sin(dtr * (128.6604 - 95795.2683 * t - 318.812 * t2e4 - 37.547 * t3e6 + 17.738 * t4e8))
- 0.047 * Math.sin(dtr * (17.3544 + 425341.6552 * t - 370.570 * t2e4 - 39.946 * t3e6 + 18.854 * t4e8))
- 0.044 * Math.sin(dtr * (160.4159 + 4067.2942 * t - 107.806 * t2e4 - 12.475 * t3e6 + 5.913 * t4e8))
- 0.043 * Math.sin(dtr * (238.1713 + 854535.1727 * t - 31.065 * t2e4 + 3.623 * t3e6 - 1.769 * t4e8))
+ 0.042 * Math.sin(dtr * (270.4555 + 1140143.5037 * t - 578.718 * t2e4 - 57.123 * t3e6 + 26.911 * t4e8))
- 0.042 * Math.sin(dtr * (132.4925 + 513197.9179 * t + 88.434 * t2e4 + 14.388 * t3e6 - 6.797 * t4e8))
- 0.041 * Math.sin(dtr * (122.3573 - 668789.4043 * t - 727.594 * t2e4 - 89.441 * t3e6 + 42.274 * t4e8))
- 0.040 * Math.sin(dtr * (105.6788 + 341337.2548 * t - 119.499 * t2e4 - 10.765 * t3e6 + 5.028 * t4e8))
+ 0.038 * Math.sin(dtr * (135.4921 + 662944.6361 * t - 668.688 * t2e4 - 71.470 * t3e6 + 33.708 * t4e8))
- 0.037 * Math.sin(dtr * (242.3910 - 51857.2124 * t - 460.540 * t2e4 - 54.293 * t3e6 + 25.652 * t4e8))
+ 0.036 * Math.sin(dtr * (336.4374 + 1303869.5784 * t - 155.171 * t2e4 - 7.020 * t3e6 + 3.259 * t4e8))
+ 0.035 * Math.sin(dtr * (223.0943 - 255454.0489 * t - 850.164 * t2e4 - 100.124 * t3e6 + 47.302 * t4e8))
- 0.034 * Math.sin(dtr * (193.2811 - 577061.4302 * t - 300.976 * t2e4 - 39.419 * t3e6 + 18.623 * t4e8))
+ 0.031 * Math.sin(dtr * (87.6023 - 918398.6850 * t - 181.476 * t2e4 - 28.654 * t3e6 + 13.594 * t4e8));
var spp = 2.4 * Math.sin(dtr * (103.2 + 377336.3 * t));
var lp = 83.353 + 4069.0137 * t - 103.238 * t2e4
- 12.492 * t3e6 + 5.263 * t4e8 + sp + 1e-3 * t * spp;
// % longitude of the ascending node
var sr = - 1.4979 * Math.sin(dtr * (49.1562 - 75869.8120 * t + 35.458 * t2e4 + 4.231 * t3e6 - 2.001 * t4e8))
- 0.1500 * Math.sin(dtr * (357.5291 + 35999.0503 * t - 1.536 * t2e4 + 0.041 * t3e6 + 0.000 * t4e8))
- 0.1226 * Math.sin(dtr * (235.7004 + 890534.2230 * t - 32.601 * t2e4 + 3.664 * t3e6 - 1.769 * t4e8))
+ 0.1176 * Math.sin(dtr * (186.5442 + 966404.0351 * t - 68.058 * t2e4 - 0.567 * t3e6 + 0.232 * t4e8))
- 0.0801 * Math.sin(dtr * (83.3826 - 12006.2998 * t + 247.999 * t2e4 + 29.262 * t3e6 - 13.826 * t4e8))
- 0.0616 * Math.sin(dtr * (51.6271 - 111868.8623 * t + 36.994 * t2e4 + 4.190 * t3e6 - 2.001 * t4e8))
+ 0.0490 * Math.sin(dtr * (100.7370 + 413335.3554 * t - 122.571 * t2e4 - 10.684 * t3e6 + 5.028 * t4e8))
+ 0.0409 * Math.sin(dtr * (308.4192 - 489205.1674 * t + 158.029 * t2e4 + 14.915 * t3e6 - 7.029 * t4e8))
+ 0.0327 * Math.sin(dtr * (134.9634 + 477198.8676 * t + 89.970 * t2e4 + 14.348 * t3e6 - 6.797 * t4e8))
+ 0.0324 * Math.sin(dtr * (46.6853 - 39870.7617 * t + 33.922 * t2e4 + 4.272 * t3e6 - 2.001 * t4e8))
+ 0.0196 * Math.sin(dtr * (98.3124 - 151739.6240 * t + 70.916 * t2e4 + 8.462 * t3e6 - 4.001 * t4e8))
+ 0.0180 * Math.sin(dtr * (274.1928 - 553068.6797 * t - 54.513 * t2e4 - 10.116 * t3e6 + 4.797 * t4e8))
+ 0.0150 * Math.sin(dtr * (325.7736 - 63863.5122 * t - 212.541 * t2e4 - 25.031 * t3e6 + 11.826 * t4e8))
- 0.0150 * Math.sin(dtr * (184.1196 + 401329.0556 * t + 125.428 * t2e4 + 18.579 * t3e6 - 8.798 * t4e8))
- 0.0078 * Math.sin(dtr * (238.1713 + 854535.1727 * t - 31.065 * t2e4 + 3.623 * t3e6 - 1.769 * t4e8))
- 0.0045 * Math.sin(dtr * (10.6638 + 1367733.0907 * t + 57.370 * t2e4 + 18.011 * t3e6 - 8.566 * t4e8))
+ 0.0044 * Math.sin(dtr * (321.5076 + 1443602.9027 * t + 21.912 * t2e4 + 13.780 * t3e6 - 6.566 * t4e8))
- 0.0042 * Math.sin(dtr * (162.8868 - 31931.7561 * t - 106.271 * t2e4 - 12.516 * t3e6 + 5.913 * t4e8))
- 0.0031 * Math.sin(dtr * (170.9849 - 930404.9848 * t + 66.523 * t2e4 + 0.608 * t3e6 - 0.232 * t4e8))
+ 0.0031 * Math.sin(dtr * (103.2079 + 377336.3051 * t - 121.035 * t2e4 - 10.724 * t3e6 + 5.028 * t4e8))
+ 0.0029 * Math.sin(dtr * (222.6120 - 1042273.8471 * t + 103.516 * t2e4 + 4.798 * t3e6 - 2.232 * t4e8))
+ 0.0028 * Math.sin(dtr * (184.0733 + 1002403.0853 * t - 69.594 * t2e4 - 0.526 * t3e6 + 0.232 * t4e8));
var srp = 25.9 * Math.sin(dtr * (125.0 - 1934.1 * t))
- 4.3 * Math.sin(dtr * (220.2 - 1935.5 * t));
var srpp = 0.38 * Math.sin(dtr * (357.5 + 35999.1 * t));
var raan = 125.0446 - 1934.13618 * t + 20.762 * t2e4
+ 2.139 * t3e6 - 1.650 * t4e8 + sr
+ 1e-3 * (srp + srpp * t);
// % mean longitude
var sl = - 0.92581 * Math.sin(dtr * (235.7004 + 890534.2230 * t - 32.601 * t2e4 + 3.664 * t3e6 - 1.769 * t4e8))
+ 0.33262 * Math.sin(dtr * (100.7370 + 413335.3554 * t - 122.571 * t2e4 - 10.684 * t3e6 + 5.028 * t4e8))
- 0.18402 * Math.sin(dtr * (357.5291 + 35999.0503 * t - 1.536 * t2e4 + 0.041 * t3e6 + 0.000 * t4e8))
+ 0.11007 * Math.sin(dtr * (134.9634 + 477198.8676 * t + 89.970 * t2e4 + 14.348 * t3e6 - 6.797 * t4e8))
- 0.06055 * Math.sin(dtr * (238.1713 + 854535.1727 * t - 31.065 * t2e4 + 3.623 * t3e6 - 1.769 * t4e8))
+ 0.04741 * Math.sin(dtr * (325.7736 - 63863.5122 * t - 212.541 * t2e4 - 25.031 * t3e6 + 11.826 * t4e8))
- 0.03086 * Math.sin(dtr * (10.6638 + 1367733.0907 * t + 57.370 * t2e4 + 18.011 * t3e6 - 8.566 * t4e8))
+ 0.02184 * Math.sin(dtr * (103.2079 + 377336.3051 * t - 121.035 * t2e4 - 10.724 * t3e6 + 5.028 * t4e8))
+ 0.01645 * Math.sin(dtr * (49.1562 - 75869.8120 * t + 35.458 * t2e4 + 4.231 * t3e6 - 2.001 * t4e8))
+ 0.01022 * Math.sin(dtr * (233.2295 + 926533.2733 * t - 34.136 * t2e4 + 3.705 * t3e6 - 1.769 * t4e8))
- 0.00756 * Math.sin(dtr * (336.4374 + 1303869.5784 * t - 155.171 * t2e4 - 7.020 * t3e6 + 3.259 * t4e8))
- 0.00530 * Math.sin(dtr * (222.5657 - 441199.8173 * t - 91.506 * t2e4 - 14.307 * t3e6 + 6.797 * t4e8))
- 0.00496 * Math.sin(dtr * (162.8868 - 31931.7561 * t - 106.271 * t2e4 - 12.516 * t3e6 + 5.913 * t4e8))
- 0.00472 * Math.sin(dtr * (297.8502 + 445267.1115 * t - 16.300 * t2e4 + 1.832 * t3e6 - 0.884 * t4e8))
- 0.00271 * Math.sin(dtr * (240.6422 + 818536.1225 * t - 29.529 * t2e4 + 3.582 * t3e6 - 1.769 * t4e8))
+ 0.00264 * Math.sin(dtr * (132.4925 + 513197.9179 * t + 88.434 * t2e4 + 14.388 * t3e6 - 6.797 * t4e8))
- 0.00254 * Math.sin(dtr * (186.5442 + 966404.0351 * t - 68.058 * t2e4 - 0.567 * t3e6 + 0.232 * t4e8))
+ 0.00234 * Math.sin(dtr * (269.9268 + 954397.7353 * t + 179.941 * t2e4 + 28.695 * t3e6 - 13.594 * t4e8))
- 0.00220 * Math.sin(dtr * (13.1347 + 1331734.0404 * t + 58.906 * t2e4 + 17.971 * t3e6 - 8.566 * t4e8))
- 0.00202 * Math.sin(dtr * (355.0582 + 71998.1006 * t - 3.072 * t2e4 + 0.082 * t3e6 + 0.000 * t4e8))
+ 0.00167 * Math.sin(dtr * (328.2445 - 99862.5625 * t - 211.005 * t2e4 - 25.072 * t3e6 + 11.826 * t4e8))
- 0.00143 * Math.sin(dtr * (173.5506 + 1335801.3346 * t - 48.901 * t2e4 + 5.496 * t3e6 - 2.653 * t4e8))
- 0.00121 * Math.sin(dtr * (98.2661 + 449334.4057 * t - 124.107 * t2e4 - 10.643 * t3e6 + 5.028 * t4e8))
- 0.00116 * Math.sin(dtr * (145.6272 + 1844931.9583 * t + 147.340 * t2e4 + 32.359 * t3e6 - 15.363 * t4e8))
+ 0.00102 * Math.sin(dtr * (105.6788 + 341337.2548 * t - 119.499 * t2e4 - 10.765 * t3e6 + 5.028 * t4e8))
- 0.00090 * Math.sin(dtr * (184.1196 + 401329.0556 * t + 125.428 * t2e4 + 18.579 * t3e6 - 8.798 * t4e8))
- 0.00086 * Math.sin(dtr * (338.9083 + 1267870.5281 * t - 153.636 * t2e4 - 7.061 * t3e6 + 3.259 * t4e8))
- 0.00078 * Math.sin(dtr * (111.4008 + 1781068.4461 * t - 65.201 * t2e4 + 7.328 * t3e6 - 3.538 * t4e8))
+ 0.00069 * Math.sin(dtr * (323.3027 - 27864.4619 * t - 214.077 * t2e4 - 24.990 * t3e6 + 11.826 * t4e8))
+ 0.00066 * Math.sin(dtr * (51.6271 - 111868.8623 * t + 36.994 * t2e4 + 4.190 * t3e6 - 2.001 * t4e8))
+ 0.00065 * Math.sin(dtr * (38.5872 + 858602.4669 * t - 138.871 * t2e4 - 8.852 * t3e6 + 4.144 * t4e8))
- 0.00060 * Math.sin(dtr * (83.3826 - 12006.2998 * t + 247.999 * t2e4 + 29.262 * t3e6 - 13.826 * t4e8))
+ 0.00054 * Math.sin(dtr * (201.4740 + 826670.7108 * t - 245.142 * t2e4 - 21.367 * t3e6 + 10.057 * t4e8))
- 0.00052 * Math.sin(dtr * (308.4192 - 489205.1674 * t + 158.029 * t2e4 + 14.915 * t3e6 - 7.029 * t4e8))
+ 0.00048 * Math.sin(dtr * (8.1929 + 1403732.1410 * t + 55.834 * t2e4 + 18.052 * t3e6 - 8.566 * t4e8))
- 0.00041 * Math.sin(dtr * (46.6853 - 39870.7617 * t + 33.922 * t2e4 + 4.272 * t3e6 - 2.001 * t4e8))
- 0.00033 * Math.sin(dtr * (274.1928 - 553068.6797 * t - 54.513 * t2e4 - 10.116 * t3e6 + 4.797 * t4e8))
+ 0.00030 * Math.sin(dtr * (160.4159 + 4067.2942 * t - 107.806 * t2e4 - 12.475 * t3e6 + 5.913 * t4e8));
var slp = 3.96 * Math.sin(dtr * (119.7 + 131.8 * t))
+ 1.96 * Math.sin(dtr * (125.0 - 1934.1 * t));
var slpp = 0.463 * Math.sin(dtr * (357.5 + 35999.1 * t))
+ 0.152 * Math.sin(dtr * (238.2 + 854535.2 * t))
- 0.071 * Math.sin(dtr * (27.8 + 131.8 * t))
- 0.055 * Math.sin(dtr * (103.2 + 377336.3 * t))
- 0.026 * Math.sin(dtr * (233.2 + 926533.3 * t));
var slppp = 14 * Math.sin(dtr * (357.5 + 35999.1 * t))
+ 5 * Math.sin(dtr * (238.2 + 854535.2 * t));
var lambda = 218.31665 + 481267.88134 * t - 13.268 * t2e4
+ 1.856 * t3e6 - 1.534 * t4e8 + sl
+ 1e-3 * (slp + slpp * t + slppp * t2e4);
/*
% output
% oev(1) = semimajor axis (kilometers)
% oev(2) = orbital eccentricity (non-dimensional)
% (0 <= eccentricity < 1)
% oev(3) = orbital inclination (radians)
% (0 <= inclination <= pi)
% oev(4) = argument of perigee (radians)
% (0 <= argument of perigee <= 2 pi)
% oev(5) = right ascension of ascending node (radians)
% (0 <= raan <= 2 pi)
% oev(6) = true anomaly (radians)
% (0 <= true anomaly <= 2 pi)
*/
var a = sma / 1.495978707e+8; // KM2AU
var e = ecc;
var i = 2.0 * Math.asin(gamma);
var w = dtr * (lp - raan);
var O = dtr * raan;
var M = dtr * (lambda - lp);
var W = w + O;
elems[off++] = a;
elems[off++] = (W + M) % (2.0 * Math.PI);
elems[off++] = e * Math.cos(W);
elems[off++] = e * Math.sin(W);
elems[off++] = Math.cos(O) * Math.sin(i / 2);
elems[off++] = Math.sin(O) * Math.sin(i / 2);
// update optional elements
if (addGM) elems[off++] = theory.GM;
if (addEpoch) elems[off++] = jy2k;
return elems;
}
elp2000orb = VSOP(elp2000orb, GMJY.ear + GMJY.moon)