/[projects]/dao/DelphiScanner/Components/tpsystools_4.04/source/StPluto.pas
ViewVC logotype

Annotation of /dao/DelphiScanner/Components/tpsystools_4.04/source/StPluto.pas

Parent Directory Parent Directory | Revision Log Revision Log

Revision 2671 - (hide annotations) (download)
Tue Aug 25 18:15:15 2015 UTC (8 years, 10 months ago) by torben
File size: 9289 byte(s) 
Added tpsystools component
1 torben 2671 // 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: StPluto.pas 4.04 *}
30     {*********************************************************}
31     {* SysTools: Astronomical Routines (for Pluto) *}
32     {*********************************************************}
33    
34     {$I StDefine.inc}
35    
36     unit StPluto;
37    
38     interface
39    
40     uses
41     StAstroP, StMath;
42    
43     function ComputePluto(JD : Double) : TStEclipticalCord;
44    
45     implementation
46    
47     function ComputePluto(JD : Double) : TStEclipticalCord;
48     var
49     T,
50     J, S, P,
51     L,
52     B,
53     R : Double;
54     begin
55     T := (JD - 2451545.0) / 36525.0;
56     J := ( 34.35 + 3034.9057 * T) / radcor;
57     S := ( 50.08 + 1222.1138 * T) / radcor;
58     P := (238.96 + 144.9600 * T) / radcor;
59    
60     L :=
61     - 19798886 * sin(P) + 19848454 * cos(P)
62     + 897499 * sin(2*P) - 4955707 * cos(2*P)
63     + 610820 * sin(3*P) + 1210521 * cos(3*P)
64     - 341639 * sin(4*P) - 189719 * cos(4*P)
65     + 129027 * sin(5*P) - 34863 * cos(5*P)
66     - 38215 * sin(6*P) + 31061 * cos(6*P)
67     + 20349 * sin(S-P) - 9886 * cos(S-P)
68     - 4045 * sin(S) - 4904 * cos(S)
69     - 5885 * sin(S+P) - 3238 * cos(S+P)
70     - 3812 * sin(S+2*P) + 3011 * cos(s+2*P)
71     - 601 * sin(S+3*P) + 3468 * cos(S+3*P)
72     + 1237 * sin(2*(S-P)) + 463 * cos(2*(S-P))
73     + 1086 * sin(2*S-P) - 911 * cos(2*S-P)
74     + 595 * sin(2*S) - 1229 * cos(2*S)
75     + 2484 * sin(J-S) - 485 * cos(J-S)
76     + 839 * sin(J-S+P) - 1414 * cos(J-S+P)
77     - 964 * sin(J-3*P) + 1059 * cos(J-3*P)
78     - 2303 * sin(J-2*P) - 1038 * cos(J-2*P)
79     + 7049 * sin(J-P) + 747 * cos(J-P)
80     + 1179 * sin(J) - 358 * cos(J)
81     + 393 * sin(J+P) - 63 * cos(J+P)
82     + 111 * sin(J+2*P) - 268 * cos(J+2*P)
83     - 52 * sin(J+3*P) - 154 * cos(J+3*P)
84     - 78 * sin(J+4*P) - 30 * cos(J+4*P)
85     - 34 * sin(J+S-3*P) - 26 * cos(J+S-3*P)
86     - 43 * sin(J+S-2*P) + 1 * cos(J+S-2*P)
87     - 15 * sin(J+S-P) + 21 * cos(J+S-P)
88     - 1 * sin(J+S) + 15 * cos(J+S)
89     + 4 * sin(J+S+P) + 7 * cos(J+S+P)
90     + 1 * sin(J+S+3*P) + 5 * cos(J+S+3*P)
91     + 8 * sin(2*J-6*P) + 3 * cos(2*J-6*P)
92     - 3 * sin(2*J-5*P) + 6 * cos(2*J-5*P)
93     + 6 * sin(2*J-4*P) - 13 * cos(2*J-4*P)
94     + 10 * sin(2*J-3*P) + 22 * cos(2*J-3*P)
95     - 57 * sin(2*J-2*P) - 32 * cos(2*J-2*P)
96     + 157 * sin(2*J-P) - 46 * cos(2*J-P)
97     + 12 * sin(2*J) - 18 * cos(2*J)
98     - 4 * sin(2*J-P) + 8 * cos(2*J-P)
99     - 5 * sin(2*(J+P)) + 0 * sin(2*(J+P))
100     + 3 * sin(2*J+3*P) + 4 * cos(2*J+3*P)
101     - 1 * sin(3*J-2*P) - 1 * cos(3*J-2*P)
102     + 6 * sin(3*J-P) - 3 * cos(3*J-P)
103     - 1 * sin(3*J) - 2 * cos(3*J);
104     Result.L0 := (238.956785 + 144.96*T + (L/1000000)) / radcor;
105    
106     B := -5453098 * sin(P) - 14974876 * cos(P)
107     + 3527363 * sin(2*P) + 1672673 * cos(2*P)
108     - 1050939 * sin(3*P) + 327763 * cos(3*P)
109     + 178691 * sin(4*P) - 291925 * cos(4*P)
110     + 18763 * sin(5*P) + 100448 * cos(5*P)
111     - 30594 * sin(6*P) - 25838 * cos(6*P)
112     + 4965 * sin(S-P) + 11263 * cos(S-P)
113     + 310 * sin(S) - 132 * cos(S)
114     + 2036 * sin(S+P) - 947 * cos(S+P)
115     - 2 * sin(S+2*P) - 674 * cos(S+2*P)
116     - 329 * sin(S+3*P) - 563 * cos(S+3*P)
117     - 64 * sin(2*(S-P)) + 39 * cos(2*(S-P))
118     - 94 * sin(2*S-P) + 210 * cos(2*S-P)
119     - 8 * sin(2*S) - 160 * cos(2*S)
120     + 177 * sin(J-S) + 259 * cos(J-S)
121     + 17 * sin(J-S+P) + 234 * cos(J-S+P)
122     + 582 * sin(J-3*P) - 285 * cos(J-3*P)
123     - 298 * sin(J-2*P) + 692 * cos(J-2*P)
124     + 157 * sin(J-P) + 201 * cos(J-P)
125     + 304 * sin(J) + 825 * cos(J)
126     - 124 * sin(J+P) - 29 * cos(J+P)
127     + 15 * sin(J+2*P) + 8 * cos(J+2*P)
128     + 7 * sin(J+3*P) + 15 * cos(J+3*P)
129     + 2 * sin(J+4*P) + 2 * cos(J+4*P)
130     + 4 * sin(J+S-3*P) + 2 * cos(J+S-3*P)
131     + 3 * sin(J+S-2*P) + 0 * cos(J+S-2*P)
132     + 1 * sin(J+S-P) - 1 * cos(J+S-P)
133     + 0 * sin(J+S) - 2 * cos(J+S)
134     + 1 * sin(J+S+P) - 0 * cos(J+S+P)
135     + 1 * sin(J+S+3*P) - 1 * cos(J+S+3*P)
136     - 2 * sin(2*J-6*P) - 3 * cos(2*J-6*P)
137     + 1 * sin(2*J-5*P) + 2 * cos(2*J-5*P)
138     - 8 * sin(2*J-4*P) + 2 * cos(2*J-4*P)
139     + 10 * sin(2*J-3*P) - 7 * cos(2*J-3*P)
140     + 0 * sin(2*J-2*P) + 21 * cos(2*J-2*P)
141     + 8 * sin(2*J-P) + 5 * cos(2*J-P)
142     + 13 * sin(2*J) + 16 * cos(2*J)
143     - 2 * sin(2*J-P) - 3 * cos(2*J-P)
144     + 0 * sin(2*(J+P)) + 0 * cos(2*(J+P))
145     + 0 * sin(2*J+3*P) + 1 * cos(2*J+3*P)
146     + 0 * sin(3*J-2*P) + 1 * cos(3*J-2*P)
147     + 0 * sin(3*J-P) + 0 * cos(3*J-P)
148     + 0 * sin(3*J) + 1 * cos(3*J);
149     Result.B0 := (-3.908202 + B/1000000) / radcor;
150    
151     R := 66867334 * sin(P) + 68955876 * cos(P)
152     - 11826086 * sin(2*P) - 333765 * cos(2*P)
153     + 1593657 * sin(3*P) - 1439953 * cos(3*P)
154     - 18948 * sin(4*P) + 482443 * cos(4*P)
155     - 66634 * sin(5*P) - 85576 * cos(5*P)
156     + 30841 * sin(6*P) - 5765 * cos(6*P)
157     - 6140 * sin(S-P) + 22254 * cos(S-P)
158     + 4434 * sin(S) + 4443 * cos(S)
159     - 1518 * sin(S+P) + 641 * cos(S+P)
160     - 5 * sin(S+2*P) + 792 * cos(s+2*P)
161     + 518 * sin(S+3*P) + 518 * cos(S+3*P)
162     - 13 * sin(2*(S-P)) - 221 * cos(2*(S-P))
163     + 837 * sin(2*S-P) - 494 * cos(2*S-P)
164     - 281 * sin(2*S) + 616 * cos(2*S)
165     + 260 * sin(J-S) - 395 * cos(J-S)
166     - 191 * sin(J-S+P) - 396 * cos(J-S+P)
167     - 3218 * sin(J-3*P) + 370 * cos(J-3*P)
168     + 8019 * sin(J-2*P) - 7689 * cos(J-2*P)
169     + 105 * sin(J-P) + 45637 * cos(J-P)
170     + 8623 * sin(J) + 8444 * cos(J)
171     - 896 * sin(J+P) - 801 * cos(J+P)
172     + 208 * sin(J+2*P) - 122 * cos(J+2*P)
173     - 133 * sin(J+3*P) + 65 * cos(J+3*P)
174     - 16 * sin(J+4*P) + 1 * cos(J+4*P)
175     - 22 * sin(J+S-3*P) + 7 * cos(J+S-3*P)
176     - 8 * sin(J+S-2*P) + 16 * cos(J+S-2*P)
177     + 2 * sin(J+S-P) + 9 * cos(J+S-P)
178     + 12 * sin(J+S) + 5 * cos(J+S)
179     + 1 * sin(J+S+P) - 3 * cos(J+S+P)
180     + 1 * sin(J+S+3*P) + 0 * cos(J+S+3*P)
181     + 9 * sin(2*J-6*P) + 5 * cos(2*J-6*P)
182     + 2 * sin(2*J-5*P) - 1 * cos(2*J-5*P)
183     + 14 * sin(2*J-4*P) + 10 * cos(2*J-4*P)
184     - 65 * sin(2*J-3*P) + 12 * cos(2*J-3*P)
185     + 126 * sin(2*J-2*P) - 233 * cos(2*J-2*P)
186     + 270 * sin(2*J-P) + 1068 * cos(2*J-P)
187     + 254 * sin(2*J) + 155 * cos(2*J)
188     - 26 * sin(2*J-P) - 2 * cos(2*J-P)
189     + 7 * sin(2*(J+P)) + 0 * cos(2*(J+P))
190     - 11 * sin(2*J+3*P) + 4 * cos(2*J+3*P)
191     + 4 * sin(3*J-2*P) - 14 * cos(3*J-2*P)
192     + 18 * sin(3*J-P) + 35 * cos(3*J-P)
193     + 13 * sin(3*J) + 3 * cos(3*J);
194     Result.R0 := 40.7247248 + R / 10000000;
195     end;
196    
197    
198     end.
  ViewVC Help
Powered by ViewVC 1.1.20