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 2026/10/31
0h TT to 2026/12/02 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 "ECL" first for the ecliptic coordinates, with at least
SIZE 031
-Then "EQ" for the equatorial
coordinates ( SIZE 039 )
-And then "AZ" 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 R30 = coordinates of the Sun, the Moon, Mercury, Venus, Mars, Jupiter, Saturn, Uranus, Neptune & Pluto.
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 "AZ"
• 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 |
S -EPH2026NOV V ECL EQ AZ |
Subroutine
that is called by "V" Section Header Ecliptic Coordinates of the Sun, the Moon & the Planets Takes day of month & time and calls "V" Ecliptic -> Equatorial Coordinates Equatorial -> Azimuthal Coordinates |
-"ECL" "EQ" &
"AZ" calculate & store the
coordinates in registers R01 thru R27 as follows:
>>> h0 is the height, corrected for refraction
| Celestial Body | Registers | "ECL" | "EQ" | "AZ" |
| 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 2026/11/24 at 16h41m TT
-Enter the day of the month and the
time expressed in Terrestrial
Time ( TT )
24
ENTER^
16.41
XEQ "ECL"
>>>>
L0 = 242°410754
= R01
RDN B0 = 0°000205
= R02
RDN R0 = 0.98727071
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 [ 2026/10/31 0h TT , 2026/12/02
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 "ECL", use "EQ" 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 2026/11/24 at 16h41m
TT
After executing "ECLI"
XEQ "EQ" or simply R/S if you've just executed
"ECL"
>>>>
RA0
= 16h01m20s84
= R01
RDN
Decl 0 = -20°38'29"71
= R02
RDN
eps = 23°437749
= R31
-The distances in R03-R06-.....-R27 are unchanged.
-Cf paragraph 4°) for
the other results
3°) Azimuthal
Topocentric Coordinates
-AFTER executing "ECL" & "EQ" use "AZ" 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 2026/11/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
= 13h07m00s33
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 | "ECL" | "EQ" | "AZ" |
| SUN |
R01 R02 R03 |
242.410754 0.000205 0.98727071 |
16.012084 -20.382971 unchanged |
135.573453 21.513058 21.535162 |
| MOON |
R04 R05 R06 |
63.446461 4.943753 0.0024104855 |
4.013139 25.413278 unchanged |
-40.465138 -18.480788 -18.480788 |
| MERCURY |
R07 R08 R09 |
223.414039 2.159460 1.09025096 |
14.463062 -13.481320 unchanged |
149.083160 37.121857 37.133356 |
| VENUS |
R10 R11 R12 |
205.003483 0.495264 0.38552732 |
13.332391 -9.130294 unchanged |
170.255783 46.574135 46.583457 |
| MARS |
R13 R14 R15 |
149.480784 2.147275 1.23248469 |
10.092780 13.394975 unchanged |
-105.111648 45.133626 45.143279 |
| JUPITER |
R16 R17 R18 |
146.482361 0.777432 5.13098523 |
9.555649 13.250724 unchanged |
-102.571325 42.211244 42.221494 |
| SATURN |
R19 R20 R21 |
8.170168 -2.611029 8.80761354 |
0.340805 0.502616 unchanged |
14.242291 -54.561913 -54.561913 |
| URANUS |
R22 R23 R24 |
63.727893 -0.148235 18.44432847 |
4.065975 20.445988 unchanged |
-45.042385 -20.565604 -20.565604 |
| NEPTUNE |
R25 R26 R27 |
1.707374 -1.402496 29.38149227 |
0.082981 -0.362828 unchanged |
25.455994 -54.291352 -54.291352 |
| True obliquity of the ecliptic |
R31 |
/ |
23.437749 |
unchanged |
|
Local Sidereal Time |
R32 |
/ |
/ |
13.070033 |
-This subroutine may be used for itself to calculate the geocentric
ecliptic coordinates
-First initialize R00 before
executing "V".
-With the example above,
R00 = 0.5434461806
WARNING !!!
6°) Refraction
-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 about 0"12 if -0°32'58"0 <= h <=
90°
References:
[1] Aldo Vitagliano SOLEX http://www.solexorb.it/
[2] ftp://ssd.jpl.nasa.gov/pub/eph/planets/ascii/
[3] Jean Meeus -
"Astronomical Algorithms" - Willmann-Bell
- ISBN 0-943396-61-1