Nutation In Obliquity Of The Ecliptic - IAU 1980 Series
Date:   ±Y  M  D Time UT:   H M S ΔT:   ± H M S
Calendar Mode:   Julian            Gregorian
 
 Date     YMD = 2017 09 26 - Tuesday (Gregorian) 
 Time  UT HMS =  00 00 00.000 = +0 day = UTFrac
 ±Delta T HMS = +00 00 00.000 = +0 day = dTFrac

 JD12 = 2458023
 JD00 = 2458022.5 = (JD12 - 0.5)

 JDTT = 2458022.5 = (JD00 + UTFrac + dTFrac)
    T = 0.1773442847364818 Julian centuries from J2000.0

 -------------------------------------------------------
 Nutation In Obliquity of the Ecliptic (IAU 1980 Series)

 Delta epsilon = -0.001876811296° = -06.757"


To Compute The Nutation In Obliquity For Any Given JD Argument

The following algorithm computes the nutation in obliquity of the ecliptic, in degrees, for a given JDTT argument.   It is based on the 1980 IAU Theory of Nutation and Reduction.

This value is the correction to apply to the MEAN value of the obliquity to obtain the TRUE obliquity, with a theoretical precision of 0.001 arc sec.

As in most computer programming languages, the trigonometric functions (sin, cos, tan, et al.) in the pseudo-code here, assume radian arguments unless otherwise indicated.

In this case, the final result is converted into degrees for convenience.

True Obliquity = Mean Obliquity + Δε

Ref:
Astronomical Algorithms - 2nd Edition (1998)
Jean Meeus
Pub. Willmann-Bell, Inc.
ISBN 0-943396-61-1

Ref:
Explanatory Supplement to the Astronomical Almanac (1992)
p111 to p114
ISBN 0-935702-68-7

WikiRef:
http://en.wikipedia.org/wiki/Nutation


Elements
JDTTAstronomical JD number for given date/time (TT)
TTime factor in Julian centuries, corresponding to JDTT, reckoned from J2000.0
DegToRadThis is a conversion factor used to convert degrees into equivalent radians.  We multiply degrees by this factor to obtain the equivalent radians.
wNutation series accumulator variable
w1Mean elongation of the moon from the sun.
w2Mean solar anomaly.
w3Mean lunar anomaly.
w4Lunar argument of latitude.
w5Longitude of ascending node of lunar orbit on the ecliptic as measured from the mean equinox of date.
dEpsDegNutation in ecliptical longitude expressed in degrees.



Pseudocode Algorithm To Compute The Nutation In Obliquity of the Ecliptic
for any given JDTT Argument.


The following PHP function returns the nutation in obliquity in decimal degrees
for given JD argument.


The following C++/CPP function returns the nutation in obliquity in decimal degrees
for given JD argument.


The following VB.NET function returns the nutation in obliquity in decimal degrees.

VB.NET function to compute the nutation in obliquity in decimal degrees for given JDTT argument.


© Jay Tanner - PHP Science Labs - 2017