function geocentric2geodetic,geocentric
;geodetic=geodetic latitude in radians

f=1/298.257222101d
esquared=2*f-f^2
;GRS80 ellipsoid, use 1/298.257223563 for WGS84

geodetic=geocentric
geodetic[0]=atan(tan(geocentric[0])/(1-esquared))
return,geodetic

end