palGeoc palGeoc - Convert geodetic position to geocentric

Description:

Convert geodetic position to geocentric.

Invocation:

void palGeoc( double p, double h, double $*$ r, double $*$z );

Arguments:

p = double (Given)

latitude (radians)

h = double (Given)

height above reference spheroid (geodetic, metres)

r = double $*$ (Returned)

distance from Earth axis (AU)

z = double $*$ (Returned)

distance from plane of Earth equator (AU)

Notes:

• Geocentric latitude can be obtained by evaluating atan2(z,r)

• Uses WGS84 reference ellipsoid and calls iauGd2gc