Description: Geodetic to geocentric transformation for ellipsoid given a,f.
The source code for the Fortran routine GD2GCE is available: