1 |
// Upgraded to Delphi 2009: Sebastian Zierer
|
2 |
|
3 |
(* ***** BEGIN LICENSE BLOCK *****
|
4 |
* Version: MPL 1.1
|
5 |
*
|
6 |
* The contents of this file are subject to the Mozilla Public License Version
|
7 |
* 1.1 (the "License"); you may not use this file except in compliance with
|
8 |
* the License. You may obtain a copy of the License at
|
9 |
* http://www.mozilla.org/MPL/
|
10 |
*
|
11 |
* Software distributed under the License is distributed on an "AS IS" basis,
|
12 |
* WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
|
13 |
* for the specific language governing rights and limitations under the
|
14 |
* License.
|
15 |
*
|
16 |
* The Original Code is TurboPower SysTools
|
17 |
*
|
18 |
* The Initial Developer of the Original Code is
|
19 |
* TurboPower Software
|
20 |
*
|
21 |
* Portions created by the Initial Developer are Copyright (C) 1996-2002
|
22 |
* the Initial Developer. All Rights Reserved.
|
23 |
*
|
24 |
* Contributor(s):
|
25 |
*
|
26 |
* ***** END LICENSE BLOCK ***** *)
|
27 |
|
28 |
{*********************************************************}
|
29 |
{* SysTools: StAstroP.pas 4.04 *}
|
30 |
{*********************************************************}
|
31 |
{* SysTools: Astronomical Routines (general Planetary) *}
|
32 |
{*********************************************************}
|
33 |
|
34 |
{$I StDefine.inc}
|
35 |
|
36 |
{ ************************************************************** }
|
37 |
{ Sources: }
|
38 |
{ 1. Astronomical Algorithms, Jean Meeus, Willmann-Bell, 1991. }
|
39 |
{ }
|
40 |
{ 2. Planetary and Lunar Coordinates (1984-2000), U.S. Govt, }
|
41 |
{ 1983. }
|
42 |
{ }
|
43 |
{ 3. Supplement to the American Ephemeris and Nautical Almanac,}
|
44 |
{ U.S. Govt, 1964. }
|
45 |
{ }
|
46 |
{ 4. MPO96-MPO98 source files, Brian D. Warner, 1995-2000. }
|
47 |
{ }
|
48 |
{ ************************************************************** }
|
49 |
|
50 |
unit StAstroP;
|
51 |
|
52 |
interface
|
53 |
|
54 |
const
|
55 |
StdDate = 2451545.0; {Ast. Julian Date for J2000 Epoch}
|
56 |
OB2000 = 0.409092804; {J2000 obliquity of the ecliptic (radians)}
|
57 |
|
58 |
type
|
59 |
TStEclipticalCord = packed record
|
60 |
L0,
|
61 |
B0,
|
62 |
R0 : Double;
|
63 |
end;
|
64 |
|
65 |
TStRectangularCord = packed record
|
66 |
X,
|
67 |
Y,
|
68 |
Z : Double;
|
69 |
end;
|
70 |
|
71 |
TStPlanetsRec = packed record
|
72 |
RA,
|
73 |
DC,
|
74 |
Elong : Double;
|
75 |
end;
|
76 |
TStPlanetsArray = array[1..8] of TStPlanetsRec;
|
77 |
|
78 |
|
79 |
procedure PlanetsPos(JD : Double; var PA : TStPlanetsArray);
|
80 |
|
81 |
|
82 |
implementation
|
83 |
|
84 |
uses
|
85 |
Windows,
|
86 |
StDate, StMerc, StVenus, StMars, StJup, StSaturn, StUranus, StNeptun,
|
87 |
StPluto, StMath;
|
88 |
|
89 |
var
|
90 |
PlanEC : TStEclipticalCord;
|
91 |
PlanRC,
|
92 |
SunRC : TStRectangularCord;
|
93 |
SunEQ : TStPlanetsRec;
|
94 |
|
95 |
|
96 |
{--------------------------------------------------------------------------}
|
97 |
|
98 |
function RealAngle(Value2, Value1, Start : Double) : Double;
|
99 |
begin
|
100 |
Result := Start;
|
101 |
|
102 |
if (Value1 = 0) then begin
|
103 |
if Value2 > 0 then
|
104 |
Result := Pi / 2.0
|
105 |
else
|
106 |
Result := 3.0 * Pi / 2.0;
|
107 |
end else begin
|
108 |
if (Value2 > 0.0) then begin
|
109 |
if (Value1 < 0.0) then
|
110 |
Result := Start + Pi
|
111 |
else
|
112 |
Result := Start;
|
113 |
end else begin
|
114 |
if (Value2 = 0) then begin
|
115 |
if Value1 > 0 then
|
116 |
Result := 0
|
117 |
else
|
118 |
Result := Pi;
|
119 |
end else begin
|
120 |
if (Value2 < 0) then begin
|
121 |
if (Value1 < 0) then
|
122 |
Result := Start + Pi
|
123 |
else
|
124 |
Result := Start + (2.0 * Pi)
|
125 |
end;
|
126 |
end;
|
127 |
end;
|
128 |
end;
|
129 |
end;
|
130 |
|
131 |
{--------------------------------------------------------------------------}
|
132 |
|
133 |
function SunOfDate(JD : Double) : TStRectangularCord;
|
134 |
{-compute J2000 XYZ coordinates of the Sun}
|
135 |
var
|
136 |
T0,
|
137 |
A,
|
138 |
L,
|
139 |
B,
|
140 |
RV,
|
141 |
TX,
|
142 |
TY,
|
143 |
TZ : Double;
|
144 |
|
145 |
begin
|
146 |
T0 := (JD - StdDate) / 365250;
|
147 |
|
148 |
{solar longitude}
|
149 |
L := 175347046
|
150 |
+ 3341656 * cos(4.6692568 + 6283.07585*T0)
|
151 |
+ 34894 * cos(4.6261000 + 12566.1517*T0)
|
152 |
+ 3497 * cos(2.7441000 + 5753.3849*T0)
|
153 |
+ 3418 * cos(2.8289000 + 3.5231*T0)
|
154 |
+ 3136 * cos(3.6277000 + 77713.7715*T0)
|
155 |
+ 2676 * cos(4.4181000 + 7860.4194*T0)
|
156 |
+ 2343 * cos(6.1352000 + 3930.2097*T0)
|
157 |
+ 1324 * cos(0.7425000 + 11506.7698*T0)
|
158 |
+ 1273 * cos(2.0371000 + 529.6910*T0)
|
159 |
+ 1199 * cos(1.1096000 + 1577.3435*T0)
|
160 |
+ 990 * cos(5.2330000 + 5884.9270*T0)
|
161 |
+ 902 * cos(2.0450000 + 26.1490*T0)
|
162 |
+ 857 * cos(3.5080000 + 398.149*T0)
|
163 |
+ 780 * cos(1.1790000 + 5223.694*T0)
|
164 |
+ 753 * cos(2.5330000 + 5507.553*T0)
|
165 |
+ 505 * cos(4.5830000 + 18849.228*T0)
|
166 |
+ 492 * cos(4.2050000 + 775.523*T0)
|
167 |
+ 357 * cos(2.9200000 + 0.067*T0)
|
168 |
+ 317 * cos(5.8490000 + 11790.626*T0)
|
169 |
+ 284 * cos(1.8990000 + 796.298*T0)
|
170 |
+ 271 * cos(0.3150000 + 10977.079*T0)
|
171 |
+ 243 * cos(0.3450000 + 5486.778*T0)
|
172 |
+ 206 * cos(4.8060000 + 2544.314*T0)
|
173 |
+ 205 * cos(1.8690000 + 5573.143*T0)
|
174 |
+ 202 * cos(2.4580000 + 6069.777*T0)
|
175 |
+ 156 * cos(0.8330000 + 213.299*T0)
|
176 |
+ 132 * cos(3.4110000 + 2942.463*T0)
|
177 |
+ 126 * cos(1.0830000 + 20.775*T0)
|
178 |
+ 115 * cos(0.6450000 + 0.980*T0)
|
179 |
+ 103 * cos(0.6360000 + 4694.003*T0)
|
180 |
+ 102 * cos(0.9760000 + 15720.839*T0)
|
181 |
+ 102 * cos(4.2670000 + 7.114*T0)
|
182 |
+ 99 * cos(6.2100000 + 2146.170*T0)
|
183 |
+ 98 * cos(0.6800000 + 155.420*T0)
|
184 |
+ 86 * cos(5.9800000 +161000.690*T0)
|
185 |
+ 85 * cos(1.3000000 + 6275.960*T0)
|
186 |
+ 85 * cos(3.6700000 + 71430.700*T0)
|
187 |
+ 80 * cos(1.8100000 + 17260.150*T0);
|
188 |
|
189 |
A := 628307584999.0
|
190 |
+ 206059 * cos(2.678235 + 6283.07585*T0)
|
191 |
+ 4303 * cos(2.635100 + 12566.1517*T0)
|
192 |
+ 425 * cos(1.590000 + 3.523*T0)
|
193 |
+ 119 * cos(5.796000 + 26.298*T0)
|
194 |
+ 109 * cos(2.966000 + 1577.344*T0)
|
195 |
+ 93 * cos(2.590000 + 18849.23*T0)
|
196 |
+ 72 * cos(1.140000 + 529.69*T0)
|
197 |
+ 68 * cos(1.870000 + 398.15*T0)
|
198 |
+ 67 * cos(4.410000 + 5507.55*T0)
|
199 |
+ 59 * cos(2.890000 + 5223.69*T0)
|
200 |
+ 56 * cos(2.170000 + 155.42*T0)
|
201 |
+ 45 * cos(0.400000 + 796.30*T0)
|
202 |
+ 36 * cos(0.470000 + 775.52*T0)
|
203 |
+ 29 * cos(2.650000 + 7.11*T0)
|
204 |
+ 21 * cos(5.340000 + 0.98*T0)
|
205 |
+ 19 * cos(1.850000 + 5486.78*T0)
|
206 |
+ 19 * cos(4.970000 + 213.30*T0)
|
207 |
+ 17 * cos(2.990000 + 6275.96*T0)
|
208 |
+ 16 * cos(0.030000 + 2544.31*T0);
|
209 |
L := L + (A * T0);
|
210 |
|
211 |
A := 8722 * cos(1.0725 + 6283.0758*T0)
|
212 |
+ 991 * cos(3.1416)
|
213 |
+ 295 * cos(0.437 + 12566.1520*T0)
|
214 |
+ 27 * cos(0.050 + 3.52*T0)
|
215 |
+ 16 * cos(5.190 + 26.30*T0)
|
216 |
+ 16 * cos(3.69 + 155.42*T0)
|
217 |
+ 9 * cos(0.30 + 18849.23*T0)
|
218 |
+ 9 * cos(2.06 + 77713.77*T0);
|
219 |
L := L + (A * sqr(T0));
|
220 |
|
221 |
A := 289 * cos(5.842 + 6283.076*T0)
|
222 |
+ 21 * cos(6.05 + 12566.15*T0)
|
223 |
+ 3 * cos(5.20 + 155.42*T0)
|
224 |
+ 3 * cos(3.14);
|
225 |
L := L + (A * sqr(T0) * T0);
|
226 |
L := L / 1.0E+8;
|
227 |
|
228 |
|
229 |
{solar latitude}
|
230 |
B := 280 * cos(3.199 + 84334.662*T0)
|
231 |
+ 102 * cos(5.422 + 5507.553*T0)
|
232 |
+ 80 * cos(3.88 + 5223.69*T0)
|
233 |
+ 44 * cos(3.70 + 2352.87*T0)
|
234 |
+ 32 * cos(4.00 + 1577.34*T0);
|
235 |
B := B / 1.0E+8;
|
236 |
|
237 |
A := 227778 * cos(3.413766 + 6283.07585*T0)
|
238 |
+ 3806 * cos(3.3706 + 12566.1517*T0)
|
239 |
+ 3620
|
240 |
+ 72 * cos(3.33 + 18849.23*T0)
|
241 |
+ 8 * cos(3.89 + 5507.55*T0)
|
242 |
+ 8 * cos(1.79 + 5223.69*T0)
|
243 |
+ 6 * cos(5.20 + 2352.87*T0);
|
244 |
B := B + (A * T0 / 1.0E+8);
|
245 |
|
246 |
A := 9721 * cos(5.1519 + 6283.07585*T0)
|
247 |
+ 233 * cos(3.1416)
|
248 |
+ 134 * cos(0.644 + 12566.152*T0)
|
249 |
+ 7 * cos(1.07 + 18849.23*T0);
|
250 |
B := B + (A * sqr(T0) / 1.0E+8);
|
251 |
|
252 |
A := 276 * cos(0.595 + 6283.076*T0)
|
253 |
+ 17 * cos(3.14)
|
254 |
+ 4 * cos(0.12 + 12566.15*T0);
|
255 |
B := B + (A * sqr(T0) * T0 / 1.0E+8);
|
256 |
|
257 |
|
258 |
{solar radius vector (astronomical units)}
|
259 |
RV := 100013989
|
260 |
+ 1670700 * cos(3.0984635 + 6283.07585*T0)
|
261 |
+ 13956 * cos(3.05525 + 12566.15170*T0)
|
262 |
+ 3084 * cos(5.1985 + 77713.7715*T0)
|
263 |
+ 1628 * cos(1.1739 + 5753.3849*T0)
|
264 |
+ 1576 * cos(2.8649 + 7860.4194*T0)
|
265 |
+ 925 * cos(5.453 + 11506.770*T0)
|
266 |
+ 542 * cos(4.564 + 3930.210*T0)
|
267 |
+ 472 * cos(3.661 + 5884.927*T0)
|
268 |
+ 346 * cos(0.964 + 5507.553*T0)
|
269 |
+ 329 * cos(5.900 + 5223.694*T0)
|
270 |
+ 307 * cos(0.299 + 5573.143*T0)
|
271 |
+ 243 * cos(4.273 + 11790.629*T0)
|
272 |
+ 212 * cos(5.847 + 1577.344*T0)
|
273 |
+ 186 * cos(5.022 + 10977.079*T0)
|
274 |
+ 175 * cos(3.012 + 18849.228*T0)
|
275 |
+ 110 * cos(5.055 + 5486.778*T0)
|
276 |
+ 98 * cos(0.89 + 6069.78*T0)
|
277 |
+ 86 * cos(5.69 + 15720.84*T0)
|
278 |
+ 86 * cos(1.27 +161000.69*T0)
|
279 |
+ 65 * cos(0.27 + 17260.15*T0)
|
280 |
+ 63 * cos(0.92 + 529.69*T0)
|
281 |
+ 57 * cos(2.01 + 83996.85*T0)
|
282 |
+ 56 * cos(5.24 + 71430.70*T0)
|
283 |
+ 49 * cos(3.25 + 2544.31*T0)
|
284 |
+ 47 * cos(2.58 + 775.52*T0)
|
285 |
+ 45 * cos(5.54 + 9437.76*T0)
|
286 |
+ 43 * cos(6.01 + 6275.96*T0)
|
287 |
+ 39 * cos(5.36 + 4694.00*T0)
|
288 |
+ 38 * cos(2.39 + 8827.39*T0)
|
289 |
+ 37 * cos(0.83 + 19651.05*T0)
|
290 |
+ 37 * cos(4.90 + 12139.55*T0)
|
291 |
+ 36 * cos(1.67 + 12036.46*T0)
|
292 |
+ 35 * cos(1.84 + 2942.46*T0)
|
293 |
+ 33 * cos(0.24 + 7084.90*T0)
|
294 |
+ 32 * cos(0.18 + 5088.63*T0)
|
295 |
+ 32 * cos(1.78 + 398.15*T0)
|
296 |
+ 28 * cos(1.21 + 6286.60*T0)
|
297 |
+ 28 * cos(1.90 + 6279.55*T0)
|
298 |
+ 26 * cos(4.59 + 10447.39*T0);
|
299 |
RV := RV / 1.0E+8;
|
300 |
|
301 |
A := 103019 * cos(1.107490 + 6283.075850*T0)
|
302 |
+ 1721 * cos(1.0644 + 12566.1517*T0)
|
303 |
+ 702 * cos(3.142)
|
304 |
+ 32 * cos(1.02 + 18849.23*T0)
|
305 |
+ 31 * cos(2.84 + 5507.55*T0)
|
306 |
+ 25 * cos(1.32 + 5223.69*T0)
|
307 |
+ 18 * cos(1.42 + 1577.34*T0)
|
308 |
+ 10 * cos(5.91 + 10977.08*T0)
|
309 |
+ 9 * cos(1.42 + 6275.96*T0)
|
310 |
+ 9 * cos(0.27 + 5486.78*T0);
|
311 |
RV := RV + (A * T0 / 1.0E+8);
|
312 |
|
313 |
A := 4359 * cos(5.7846 + 6283.0758*T0)
|
314 |
+ 124 * cos(5.579 + 12566.152*T0)
|
315 |
+ 12 * cos(3.14)
|
316 |
+ 9 * cos(3.63 + 77713.77*T0)
|
317 |
+ 6 * cos(1.87 + 5573.14*T0)
|
318 |
+ 3 * cos(5.47 + 18849.23*T0);
|
319 |
RV := RV + (A * sqr(T0) / 1.0E+8);
|
320 |
|
321 |
L := (L + PI);
|
322 |
L := Frac(L / 2.0 / PI) * 2.0 * Pi;
|
323 |
if L < 0 then
|
324 |
L := L + (2.0*PI);
|
325 |
B := -B;
|
326 |
|
327 |
TX := RV * cos(B) * cos(L);
|
328 |
TY := RV * cos(B) * sin(L);
|
329 |
TZ := RV * sin(B);
|
330 |
|
331 |
Result.X := TX + 4.40360E-7 * TY - 1.90919E-7 * TZ;
|
332 |
Result.Y := -4.79966E-7 * TX + 0.917482137087 * TY - 0.397776982902 * TZ;
|
333 |
Result.Z := 0.397776982902 * TY + 0.917482137087 * TZ;
|
334 |
end;
|
335 |
|
336 |
{--------------------------------------------------------------------------}
|
337 |
|
338 |
function EclipticToRectangular(Longitude, Latitude,
|
339 |
RadiusVector : Double) : TStRectangularCord;
|
340 |
var
|
341 |
var1,
|
342 |
var2,
|
343 |
var3 : Double;
|
344 |
begin
|
345 |
var1 := RadiusVector * cos(Longitude) * cos(Latitude);
|
346 |
var2 := RadiusVector * sin(Longitude) * cos(Latitude);
|
347 |
var3 := RadiusVector * sin(Latitude);
|
348 |
|
349 |
Result.X := var1;
|
350 |
Result.Y := var2 * cos(OB2000) - var3 * sin(OB2000);
|
351 |
Result.Z := var2 * sin(OB2000) + var3 * cos(OB2000);
|
352 |
end;
|
353 |
|
354 |
{--------------------------------------------------------------------------}
|
355 |
|
356 |
function RADec(Planet, Sun : TStRectangularCord;
|
357 |
ComputeElong : Boolean) : TStPlanetsRec;
|
358 |
var
|
359 |
var1,
|
360 |
var2,
|
361 |
var3,
|
362 |
var4,
|
363 |
var5 : Double;
|
364 |
begin
|
365 |
FillChar(Result, SizeOf(TStPlanetsRec), #0);
|
366 |
|
367 |
var1 := Sun.X + Planet.X;
|
368 |
var2 := Sun.Y + Planet.Y;
|
369 |
var3 := Sun.Z + Planet.Z;
|
370 |
|
371 |
var4 := arctan(var2/var1);
|
372 |
var4 := RealAngle(var2, var1, var4) * radcor;
|
373 |
|
374 |
var5 := sqrt(sqr(var1) + sqr(var2) + sqr(var3));
|
375 |
var3 := StInvsin(var3/var5) * radcor;
|
376 |
|
377 |
Result.RA := var4;
|
378 |
Result.DC := var3;
|
379 |
|
380 |
var4 := Result.RA / radcor;
|
381 |
var3 := Result.DC / radcor;
|
382 |
|
383 |
if (ComputeElong) then begin
|
384 |
var1 := sin(SunEQ.DC/radcor) * sin(var3);
|
385 |
var2 := cos(SunEQ.DC/radcor) * cos(var3) * cos(SunEQ.RA/radcor - var4);
|
386 |
Result.Elong := StInvcos(var1+var2) * radcor;
|
387 |
end;
|
388 |
end;
|
389 |
|
390 |
{--------------------------------------------------------------------------}
|
391 |
|
392 |
function MercuryPosition(JD : Double) : TStPlanetsRec;
|
393 |
begin
|
394 |
PlanEC := ComputeMercury(JD);
|
395 |
PlanRC := EclipticToRectangular(PlanEC.L0, PlanEC.B0, PlanEC.R0);
|
396 |
Result := RADec(PlanRC, SunRC, True);
|
397 |
end;
|
398 |
|
399 |
{--------------------------------------------------------------------------}
|
400 |
|
401 |
function VenusPosition(JD : Double) : TStPlanetsRec;
|
402 |
begin
|
403 |
PlanEC := ComputeVenus(JD);
|
404 |
PlanRC := EclipticToRectangular(PlanEC.L0, PlanEC.B0, PlanEC.R0);
|
405 |
Result := RADec(PlanRC, SunRC, True);
|
406 |
end;
|
407 |
|
408 |
{--------------------------------------------------------------------------}
|
409 |
|
410 |
function MarsPosition(JD : Double) : TStPlanetsRec;
|
411 |
begin
|
412 |
PlanEC := ComputeMars(JD);
|
413 |
PlanRC := EclipticToRectangular(PlanEC.L0, PlanEC.B0, PlanEC.R0);
|
414 |
Result := RADec(PlanRC, SunRC, True);
|
415 |
end;
|
416 |
|
417 |
{--------------------------------------------------------------------------}
|
418 |
|
419 |
function JupiterPosition(JD : Double) : TStPlanetsRec;
|
420 |
begin
|
421 |
PlanEC := ComputeJupiter(JD);
|
422 |
PlanRC := EclipticToRectangular(PlanEC.L0, PlanEC.B0, PlanEC.R0);
|
423 |
Result := RADec(PlanRC, SunRC, True);
|
424 |
end;
|
425 |
|
426 |
{--------------------------------------------------------------------------}
|
427 |
|
428 |
function SaturnPosition(JD : Double) : TStPlanetsRec;
|
429 |
begin
|
430 |
PlanEC := ComputeSaturn(JD);
|
431 |
PlanRC := EclipticToRectangular(PlanEC.L0, PlanEC.B0, PlanEC.R0);
|
432 |
Result := RADec(PlanRC, SunRC, True);
|
433 |
end;
|
434 |
|
435 |
{--------------------------------------------------------------------------}
|
436 |
|
437 |
function UranusPosition(JD : Double) : TStPlanetsRec;
|
438 |
begin
|
439 |
PlanEC := ComputeUranus(JD);
|
440 |
PlanRC := EclipticToRectangular(PlanEC.L0, PlanEC.B0, PlanEC.R0);
|
441 |
Result := RADec(PlanRC, SunRC, True);
|
442 |
end;
|
443 |
|
444 |
{--------------------------------------------------------------------------}
|
445 |
|
446 |
function NeptunePosition(JD : Double) : TStPlanetsRec;
|
447 |
begin
|
448 |
PlanEC := ComputeNeptune(JD);
|
449 |
PlanRC := EclipticToRectangular(PlanEC.L0, PlanEC.B0, PlanEC.R0);
|
450 |
Result := RADec(PlanRC, SunRC, True);
|
451 |
end;
|
452 |
|
453 |
{--------------------------------------------------------------------------}
|
454 |
|
455 |
function PlutoPosition(JD : Double) : TStPlanetsRec;
|
456 |
begin
|
457 |
PlanEC := ComputePluto(JD);
|
458 |
PlanRC := EclipticToRectangular(PlanEC.L0, PlanEC.B0, PlanEC.R0);
|
459 |
Result := RADec(PlanRC, SunRC, True);
|
460 |
end;
|
461 |
|
462 |
{--------------------------------------------------------------------------}
|
463 |
|
464 |
procedure PlanetsPos(JD : Double; var PA : TStPlanetsArray);
|
465 |
var
|
466 |
I : Integer;
|
467 |
Sun : TStRectangularCord;
|
468 |
begin
|
469 |
{find Sun's Rectangular Coordinates}
|
470 |
SunRC := SunofDate(JD);
|
471 |
|
472 |
FillChar(SunEQ, SizeOf(TStPlanetsRec), #0);
|
473 |
FillChar(Sun, SizeOf(TStRectangularCord), #0);
|
474 |
|
475 |
{find Sun's RA/Dec}
|
476 |
SunEQ := RADec(SunRC, Sun, False);
|
477 |
PA[1] := PlutoPosition(JD);
|
478 |
|
479 |
{find RA/Dec of each planet}
|
480 |
for I := 1 to 8 do begin
|
481 |
case I of
|
482 |
1 : PA[I] := MercuryPosition(JD);
|
483 |
2 : PA[I] := VenusPosition(JD);
|
484 |
3 : PA[I] := MarsPosition(JD);
|
485 |
4 : PA[I] := JupiterPosition(JD);
|
486 |
5 : PA[I] := SaturnPosition(JD);
|
487 |
6 : PA[I] := UranusPosition(JD);
|
488 |
7 : PA[I] := NeptunePosition(JD);
|
489 |
8 : PA[I] := PlutoPosition(JD);
|
490 |
end;
|
491 |
end;
|
492 |
end;
|
493 |
|
494 |
|
495 |
end.
|