EarthOrientation.jl

Calculate Earth orientation parameters from IERS tables in Julia.

Installation

The package can be installed through Julia's package manager:

Pkg.add("EarthOrientation")

Quickstart

Fetch the latest [IERS][iers-link] tables:

using EarthOrientation
EarthOrientation.update()

Get the current Earth orientation parameters, e.g. for polar motion:

xp, yp = polarmotion(now()) # arcseconds

Or the current difference between UT1 and UTC and the associated prediction error:

ΔUT1 = getΔUT1(now()) # seconds
ΔUT1_err = getΔUT1_err(now()) # seconds

Available data

• Polar motion:
  • x-coordinate of Earth's north pole: getxp
  • y-coordinate of Earth's north pole: getyp
  • both: polarmotion
• Earth rotation
  • Difference between UT1 and UTC: getΔUT1
  • Difference between UT1 and TAI: getΔUT1_TAI
  • Excess length of day: getlod
• Precession and nutation based on the 1980 IAU conventions
  • Correction to the nutation of the ecliptic: getdψ
  • Correction to the obliquity of the ecliptic: getdϵ
  • both: precession_nutation80
• Precession and nutation based on the 2000 IAU conventions
  • Correction to the celestial pole's x-coordinate: getdx
  • Correction to the celestial pole's y-coordinate: getdy
  • both: precession_nutation00

There is an associated function that returns the prediction error for each data type, e.g. getxp_err.

Manual Data Management

By default the files downloaded by EarthOrientation.update() will be used. It is also possible to pass different finals.all and finals2000A.all files in CSV format.

using EarthOrientation

push!(EOP_DATA, "finals.csv", "finals2000A.csv")