Book contents
- Frontmatter
- Contents
- Preface
- Using your personal computer for astronomy
- DEFAULT: default value input routnine & YESNO: ‘Y’ or ‘N’ input routine
- MINSEC: converts between decimal hours/degrees and minutes/seconds form
- JULDAY: calendar date to Julian day number since 1900 January 0.5
- CALDAY: Julian day number since 1900 January 0.5 to calendar date
- TIME: converts between local civil and sidereal times
- EQHOR: converts between equatorial and horizon coordinates
- HRANG: converts between right ascension and hour angle
- OBLIQ: calculates the value of the obliquity of the ecliptic
- NUTAT: finds corrections for nutation in longitude and obliquity
- EQECL: converts between equatorial and ecliptic coordinates
- EQGAL: converts between equatorial and galactic coordinates
- GENCON: converts between any of the coordinate systems
- PRCESS1: approximate precession of equatorial coordinates & PRCESS2: rigorous precession of equatorial coordinates
- PARALLX: converts between geocentric and apparent position
- REFRACT: calculates the effect of atmospheric refraction
- RISET: finds the circumstances of rising and setting
- ANOMALY: solves Kepler's equation for elliptical motion
- SUN: finds the ecliptic coordinates of the Sun
- SUNRS: finds the circumstances of sunrise and sunset
- PELMENT: returns the orbital elements of the major planets
- PLANS: finds the position of a planet
- MOON: finds the position and parallax of the Moon
- MOONRS: finds the circumstances of moonrise and moonset
- MOONNF: finds the times of new and full moon
- ECLIPSE: finds the circumstances of lunar and solar eclipses
- DISPLAY: displays an eclipse in graphical form
- ELOSC: finds positions from osculating elliptical elements
- RELEM: converts elliptic orbital elements from one epoch to another
- PCOMET: finds the position of a comet from parabolic elements
- PFIT: finds parabolic elements from observations & EFIT: finds elliptical elements from observations
- List of variables
- Bibliography
- Index
- PROGRAMS AVAILABLE ON DISK
PRCESS1: approximate precession of equatorial coordinates & PRCESS2: rigorous precession of equatorial coordinates
Published online by Cambridge University Press: 17 February 2010
- Frontmatter
- Contents
- Preface
- Using your personal computer for astronomy
- DEFAULT: default value input routnine & YESNO: ‘Y’ or ‘N’ input routine
- MINSEC: converts between decimal hours/degrees and minutes/seconds form
- JULDAY: calendar date to Julian day number since 1900 January 0.5
- CALDAY: Julian day number since 1900 January 0.5 to calendar date
- TIME: converts between local civil and sidereal times
- EQHOR: converts between equatorial and horizon coordinates
- HRANG: converts between right ascension and hour angle
- OBLIQ: calculates the value of the obliquity of the ecliptic
- NUTAT: finds corrections for nutation in longitude and obliquity
- EQECL: converts between equatorial and ecliptic coordinates
- EQGAL: converts between equatorial and galactic coordinates
- GENCON: converts between any of the coordinate systems
- PRCESS1: approximate precession of equatorial coordinates & PRCESS2: rigorous precession of equatorial coordinates
- PARALLX: converts between geocentric and apparent position
- REFRACT: calculates the effect of atmospheric refraction
- RISET: finds the circumstances of rising and setting
- ANOMALY: solves Kepler's equation for elliptical motion
- SUN: finds the ecliptic coordinates of the Sun
- SUNRS: finds the circumstances of sunrise and sunset
- PELMENT: returns the orbital elements of the major planets
- PLANS: finds the position of a planet
- MOON: finds the position and parallax of the Moon
- MOONRS: finds the circumstances of moonrise and moonset
- MOONNF: finds the times of new and full moon
- ECLIPSE: finds the circumstances of lunar and solar eclipses
- DISPLAY: displays an eclipse in graphical form
- ELOSC: finds positions from osculating elliptical elements
- RELEM: converts elliptic orbital elements from one epoch to another
- PCOMET: finds the position of a comet from parabolic elements
- PFIT: finds parabolic elements from observations & EFIT: finds elliptical elements from observations
- List of variables
- Bibliography
- Index
- PROGRAMS AVAILABLE ON DISK
Summary
The right ascension and declination of a celestial body change slowly with time because of luni-solar precession, a gyrating motion of the Earth's axis caused by the gravitational effects of the Sun and the Moon on the equatorial bulge of the Earth. The Earth's axis describes a cone about a line perpendicular to the ecliptic through the centre of the Earth with half-angle approximately 23.5 degrees, completing one circuit in about 25 800 years. The change is thus about 50 arcseconds per year and this must be allowed for when comparing equatorial coordinates at two different epochs.
PRCESS1
Two subroutines are given here which will correct equatorial coordinates for the effects of precession between one epoch and another. The first of these, PRCESS1 (2500), uses an approximate method which is simple and compact, but which nevertheless gives good results. You should find it sufficient for most purposes. The corrections, XI and Yl, which have to be made to the right ascension, X, and declination, Y, measured at epoch A, to find the corresponding values at epoch B, are well approximated by the expressions:
X1 = (MP + (NP × SIN(X) × TAN(Y))/15) × NY
Y1 = NP × COS(X) × NY,
where NY is the number of years between the two epochs, and MP and NP are ‘constants’ which change slowly with time.
- Type
- Chapter
- Information
- Astronomy with your Personal Computer , pp. 84 - 92Publisher: Cambridge University PressPrint publication year: 1990