palEl2ue palEl2ue - Transform conventional elements into "universal" form

Description:
Transform conventional osculating elements into "universal" form.

Invocation:
void palEl2ue ( double date, int jform, double epoch, double orbinc, double anode, double perih, double aorq, double e, double aorl, double dm, double u[13], int $*$jstat );

Arguments:

date = double (Given)
Epoch (TT MJD) of osculation (Note 3)
jform = int (Given)
Element set actually returned (1-3; Note 6)
epoch = double (Given)
Epoch of elements (TT MJD)
orbinc = double (Given)
inclination (radians)
anode = double (Given)
longitude of the ascending node (radians)
perih = double (Given)
longitude or argument of perihelion (radians)
aorq = double (Given)
mean distance or perihelion distance (AU)
e = double (Given)
eccentricity
aorl = double (Given)
mean anomaly or longitude (radians, JFORM=1,2 only)
dm = double (Given)
daily motion (radians, JFORM=1 only)
u = double [13] (Returned)
Universal orbital elements (Note 1) (0) combined mass (M$+$m) (1) total energy of the orbit (alpha) (2) reference (osculating) epoch (t0) (3-5) position at reference epoch (r0) (6-8) velocity at reference epoch (v0) (9) heliocentric distance at reference epoch (10) r0.v0 (11) date (t) (12) universal eccentric anomaly (psi) of date, approx
jstat = int $*$ (Returned)
status: 0 = OK
  • 1 = illegal JFORM

  • 2 = illegal E

  • 3 = illegal AORQ

  • 4 = illegal DM

  • 5 = numerical error


Notes:

Option JFORM=1, suitable for the major planets:

EPOCH = epoch of elements (TT MJD) ORBINC = inclination i (radians) ANODE = longitude of the ascending node, big omega (radians) PERIH = longitude of perihelion, curly pi (radians) AORQ = mean distance, a (AU) E = eccentricity, e (range 0 to $<$1) AORL = mean longitude L (radians) DM = daily motion (radians)

Option JFORM=2, suitable for minor planets:

EPOCH = epoch of elements (TT MJD) ORBINC = inclination i (radians) ANODE = longitude of the ascending node, big omega (radians) PERIH = argument of perihelion, little omega (radians) AORQ = mean distance, a (AU) E = eccentricity, e (range 0 to $<$1) AORL = mean anomaly M (radians)

Option JFORM=3, suitable for comets:

EPOCH = epoch of perihelion (TT MJD) ORBINC = inclination i (radians) ANODE = longitude of the ascending node, big omega (radians) PERIH = argument of perihelion, little omega (radians) AORQ = perihelion distance, q (AU) E = eccentricity, e (range 0 to 10)


See Also
Everhart & Pitkin, Am.J.Phys. 51, 712 (1983).