| 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: 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.
|
Parent Directory
| 
Revision Log