Description: Geocentric to geodetic transformation for ellipsoid given a,f.
The source code for the Fortran routine GC2GDE is available: