Overview
1°)
Ecliptic Geocentric
Coordinates
2°)
Equatorial Geocentric
Coordinates
3°)
Azimuthal Topocentric Coordinates
4°) Numerical
Results
-These programs compute
accurate positions of the Sun, the Moon and the major
planets ( This month, not enough room for Pluto )
for a short time-span
of 32 days, i-e 2022/07/31 0h TT to 2022/09/01
0h TT
-The longitudes & latitudes and the right-ascensions & declinations
are geocentric apparent,
referred to the
true equator & equinox of the date, corrected
for aberration and light-time.
-The precision is about 0"01 for the longitudes & latitudes and of
the order of 3 E-8 AU for the distances ( 5 E-11
AU for the Moon ).
-The distances are true
distances.
-The azimuthal ( topocentric ) coordinates are also given, corrected for parallax & diurnal aberration.
-These coordinates are calculated by polynomials fitted to the JPL Ephemerides
DE441
Notes:
-Always execute "ECLI" first for the ecliptic coordinates, with at least
SIZE 031
-Then "EQUA" for the
equatorial coordinates ( SIZE 039 )
-And then "AZIM" for
the azimuthal coordinates with at least SIZE 041.
-The azimuths are reckoned clockwise from North.
-Longitudes are positive
East.
Data Registers
R00 = ( DOM - 16 ) / 16 ( from -1 to +1 ) Terrestrial Time ( TT )
R01 thru R27 = coordinates of the Sun, the Moon, Mercury, Venus, Mars, Jupiter, Saturn, Uranus & Neptune.
R31 = True obliquity of the ecliptic ( deg )
R32 = Local Sidereal
Time ( hh.mnss )
• R33 = Longitude of the observer ( ° ' " )
positive East
• R34 = Latitude
of the observer ( ° ' " )
Registers R33-R34-R35 are to be initialized
before executing "AZIM"
• R35 = Observer
altitude in meters
( R36 to R40: temporary data storage )
XROM | Function | Desciption |
24,00 24,01 24,02 24,03 24,04 24,05 24,06 24,07 24,08 24,09 24,10 24,11 24,12 |
-EPH2022AUG SUN MOON MER VEN MAR JUP SAT URA NEP ECLI EQUA AZIM |
Section Header Ecliptic Coordinates of the Sun Ecliptic Coordinates of the Moon Ecliptic Coordinates of Mercury Ecliptic Coordinates of Venus Ecliptic Coordinates of Mars Ecliptic Coordinates of Jupiter Ecliptic Coordinates of Saturn Ecliptic Coordinates of Uranus Ecliptic Coordinates of Neptune Ecliptic Coordinates of the Sun, Planets & Moon Ecliptic -> Equatorial Coordinates Equatorial -> Azimuthal Coordinates |
-"ECLI" "EQUA"
& "AZIM" calculate & store the
coordinates in registers R01 thru R27 as follows:
>>> h0 is the height, corrected for refraction
Celestial Body | Registers | "ECLI" | "EQUA" | "AZIM" |
SUN |
R01 R02 R03 |
Eclipt Longitude ( deg ) Eclipt Latitude ( deg ) Dist from Earth ( AU ) |
Right-Ascens(hh;mnss) Declination ( ° ' " ) Dist from Earth ( AU ) |
Azimuth ( ° ' " ) height ( ° ' " ) h0 ( ° ' " ) |
MOON |
R04 R05 R06 |
Eclipt Longitude ( deg ) Eclipt Latitude ( deg ) Dist from Earth ( AU ) |
Right-Ascens(hh;mnss) Declination ( ° ' " ) Dist from Earth ( AU ) |
Azimuth ( ° ' " ) height ( ° ' " ) h0 ( ° ' " ) |
MERCURY |
R07 R08 R09 |
Eclipt Longitude ( deg ) Eclipt Latitude ( deg ) Dist from Earth ( AU ) |
Right-Ascens(hh;mnss) Declination ( ° ' " ) Dist from Earth ( AU ) |
Azimuth ( ° ' " ) height ( ° ' " ) h0 ( ° ' " ) |
VENUS |
R10 R11 R12 |
Eclipt Longitude ( deg ) Eclipt Latitude ( deg ) Dist from Earth ( AU ) |
Right-Ascens(hh;mnss) Declination ( ° ' " ) Dist from Earth ( AU ) |
Azimuth ( ° ' " ) height ( ° ' " ) h0 ( ° ' " ) |
MARS |
R13 R14 R15 |
Eclipt Longitude ( deg ) Eclipt Latitude ( deg ) Dist from Earth ( AU ) |
Right-Ascens(hh;mnss) Declination ( ° ' " ) Dist from Earth ( AU ) |
Azimuth ( ° ' " ) height ( ° ' " ) h0 ( ° ' " ) |
JUPITER |
R16 R17 R18 |
Eclipt Longitude ( deg ) Eclipt Latitude ( deg ) Dist from Earth ( AU ) |
Right-Ascens(hh;mnss) Declination ( ° ' " ) Dist from Earth ( AU ) |
Azimuth ( ° ' " ) height ( ° ' " ) h0 ( ° ' " ) |
SATURN |
R19 R20 R21 |
Eclipt Longitude ( deg ) Eclipt Latitude ( deg ) Dist from Earth ( AU ) |
Right-Ascens(hh;mnss) Declination ( ° ' " ) Dist from Earth ( AU ) |
Azimuth ( ° ' " ) height ( ° ' " ) h0 ( ° ' " ) |
URANUS |
R22 R23 R24 |
Eclipt Longitude ( deg ) Eclipt Latitude ( deg ) Dist from Earth ( AU ) |
Right-Ascens(hh;mnss) Declination ( ° ' " ) Dist from Earth ( AU ) |
Azimuth ( ° ' " ) height ( ° ' " ) h0 ( ° ' " ) |
NEPTUNE |
R25 R26 R27 |
Eclipt Longitude ( deg ) Eclipt Latitude ( deg ) Dist from Earth ( AU ) |
Right-Ascens(hh;mnss) Declination ( ° ' " ) Dist from Earth ( AU ) |
Azimuth ( ° ' " ) height ( ° ' " ) h0 ( ° ' " ) |
1°) Ecliptic Geocentric Coordinates of the Sun, the
Moon & the major Planets
STACK | INPUTS | OUTPUTS |
Z | / | R0 ( AU ) |
Y | Day of the Month | B0 ( deg ) |
X | HH.MNSS(TT) | L0 ( deg ) |
Where L = Longitude B = Latitude R = radius vector
Example: Calculate the apparent geocentric ecliptic coordinates of the Sun, the Moon and the planets on 2022/08/24 at 16h41m TT
-Enter the day of the month and the time
expressed in Terrestrial Time ( TT )
24
ENTER^
16.41
XEQ "ECLI"
>>>> L0 =
151°502272 = R01
---Execution time = 81s---
RDN B0
= 0°000223
= R02
RDN R0
= 1.01102590 AU = R03
Notes:
-All the angles are expressed in decimal degrees.
-Cf paragraph 4°)
for the other results.
-If you key in a date outside the interval [ 2022/07/31 0h TT , 2022/09/01
0h TT ] you'll get a DATA ERROR message.
-However, this program
may probably be used a few hours outside the prescribed
interval: set F25 and R/S
-But the precision is
less guaranteed and the results may even become completely
meaningless several days before 00 or after 32, especially
for the Moon.
2°) Equatorial Geocentric Coordinates
-AFTER executing "ECLI", use "EQUA" to get the equatorial coordinates
-The right-ascensions
are expressed in hh.mnss and the declinations in °
' "
-They replace the ecliptic
longitudes & latitudes ( cf the tableau in the
paragraph above )
-"EQUA" also calculates the true obliquity of the ecliptic which is returned
in Z-register
-A polynomial is also
used for that.
STACK | INPUTS | OUTPUTS |
Z | / | eps ( deg ) |
Y | / | Decl0 ( ° ' " ) |
X | / | RA0 ( hh.mnss ) |
Where RA = Right-Ascension Decl = declination eps = true obliquity of the ecliptic
Example: Calculate the apparent geocentric equatorial
coordinates of the Sun, the Moon and the planets on 2022/08/24 at 16h41m
TT
After executing "ECLI"
XEQ "EQUA"
or simply R/S if you've just executed "ECLI"
>>>> RA0 =
10h14m05s22 = R01
---Execution time = 45s---
RDN
Decl 0 = 10°56'24"61
= R02
RDN
eps = 23°4381483
= R31
-The distances in R03-R06-.....-R27 are unchanged.
-Cf paragraph 4°) for
the other results
3°) Azimuthal
Topocentric Coordinates
-AFTER executing "ECLI" & "EQUA" use "AZIM" to get the horizontal coordinates
-The azimuths & heights
are expressed in ° ' "
-The heights corrected for refraction are also computed and replace the
distances in R03 R06 ..... R27
STACK | INPUTS | OUTPUTS |
Z | / | h0 ( ° ' " ) |
Y | / | h ( ° ' " ) |
X | / | Az ( ° ' " ) |
Long = longitude ( positive East )
Az = Azimuth ( clockwise from North )
|
Where
Lat = latitude
h = height
> of the
Sun
Alt = altitude in meters
h0 = height ( corrected for refraction
) |
Example: Calculate the apparent topocentric
azimuthal coordinates of the Sun, the Moon
and the planets on 2022/08/24 at 16h41m TT
at the Palomar Observatory, Longitude
= 116°51'50"4 W Latitude = 33°21'22"4
N Altitude = 1706 m
>>> After executing "ECLI" & "EQUA"
-116.51504 STO 33
which are
the topocentric coordinates of the Sun.
>>> We also have the local sidereal time in R32 = LST
= 7h04m08s72
Notes:
-Cf paragraph 4°) for the other results.
-The difference TT - UTC
= 69.184 seconds.
-> h0 is often meaningless
when h < 0
Celestial Body | Registers | "ECLI" | "EQUA" | "AZIM" |
SUN |
R01 R02 R03 |
151.502272 0.000223 1.01102590 |
10.140522 10.562461 unchanged |
105.543785 41.111464 41.121972 |
MOON |
R04 R05 R06 |
121.752813 4.847109 0.0026999530 |
8.204435 24.292539 unchanged |
112.565219 70.471551 70.473536 |
MERCURY |
R07 R08 R09 |
178.580643 -1.749010 0.97153433 |
11.520052 -1.022496 unchanged |
101.015663 14.232361 14.270136 |
VENUS |
R10 R11 R12 |
135.845877 0.889567 1.63817328 |
9.141963 16.560587 unchanged |
111.144741 56.285886 56.293661 |
MARS |
R13 R14 R15 |
62.519502 -1.420701 1.00051658 |
4.025901 19.161257 unchanged |
-97.355364 47.241337 47.250575 |
JUPITER |
R16 R17 R18 |
7.558742 -1.549336 4.10442679 |
0.301280 1.343223 unchanged |
-83.593282 -6.121269 -6.121269 |
SATURN |
R19 R20 R21 |
-38.841995 -1.313495 8.87117168 |
21.355517 -15.412353 unchanged |
-72.040171 -51.312335 -51.312335 |
URANUS |
R22 R23 R24 |
48.920827 -0.368178 19.44395328 |
3.061705 17.053612 unchanged |
-91.294625 34.332757 34.345016 |
NEPTUNE |
R25 R26 R27 |
-5.352613 -1.217783 28.98325705 |
23.421683 -3.143993 unchanged |
-81.122784 -18.493638 -18.493638 |
True obliquity of the ecliptic |
R31 |
/ |
23.4381483 |
unchanged |
Local Sidereal Time |
R32 |
/ |
/ |
7.040872 |
5°) Sun-Moon-Mercury-Venus-Mars-Jupiter-Saturn-Uranus-Neptune
-All these subroutines may be used for themselves to calculate the geocentric
ecliptic coordinates
-First initialize R00
before executing them.
STACK | INPUTS | OUTPUTS |
Z | / | R ( AU ) |
Y | / | B ( deg ) |
X | / | L ( deg ) |
Where L = Longitude B = Latitude R = radius vector
Example: The same one, which corresponds
to R00 = 0.5434461806
>>> Likewise with the Moon, Mercury, ........... , & Neptune
( see above the numerical values
)
WARNING !!!
-Unlike "ECLI" , these routines do not check that R00 is between
-1 and +1
Remark:
-The apparent heights are calculated by a refraction formula which approximates
the Pulkovo refraction tables
for standard conditions
of temperature & pressure ( T = 15°C , P
= 1013.25 mbar, humidity = 0 , wave length = 0.59µ
)
-The precision is better than 0"06 over the whole range [ -0°32'58"0
, 90° ]
h0 ~ h + 1° / 62.95929 /
Tan ( h + 4°8043 / ( h + 7°0822 / ( h +11°1187
/ ( h + 38°2290 / ( h + 9°9098 ) ) ) ) )
References:
[1] Aldo Vitagliano SOLEX http://www.solexorb.it/