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: