palPv2ue
palPv2ue
- Universal elements to position and velocity
- Description:
- Construct a universal element set based on an instantaneous position
and velocity.
- Invocation:
- void palPv2ue( const double pv[6], double date, double pmass,
double u[13], int
jstat );
- Arguments:
-
-
pv = double [6] (Given)
-
Heliocentric x,y,z,xdot,ydot,zdot of date, (AU,AU/s; Note 1)
-
date = double (Given)
-
Date (TT modified Julian Date = JD-2400000.5)
-
pmass = double (Given)
-
Mass of the planet (Sun=1; note 2)
-
u = double [13] (Returned)
-
Universal orbital elements (Note 3)
- (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 PMASS
- -2 = too close to Sun
- -3 = too slow
- Notes:
- The PV 6-vector can be with respect to any chosen inertial frame,
and the resulting universal-element set will be with respect to
the same frame. A common choice will be mean equator and ecliptic
of epoch J2000.
- The mass, PMASS, is important only for the larger planets. For
most purposes (e.g. asteroids) use 0D0. Values less than zero
are illegal.
- The "universal" elements are those which define the orbit for the
purposes of the method of universal variables (see reference).
They consist of the combined mass of the two bodies, an epoch,
and the position and velocity vectors (arbitrary reference frame)
at that epoch. The parameter set used here includes also various
quantities that can, in fact, be derived from the other
information. This approach is taken to avoiding unnecessary
computation and loss of accuracy. The supplementary quantities
are (i) alpha, which is proportional to the total energy of the
orbit, (ii) the heliocentric distance at epoch, (iii) the
outwards component of the velocity at the given epoch, (iv) an
estimate of psi, the "universal eccentric anomaly" at a given
date and (v) that date.
- Reference: Everhart, E. & Pitkin, E.T., Am.J.Phys. 51, 712, 1983.