code: def sphlat(r, colat, lons):
      
      r = ctypes.c_double(r)
      colat = ctypes.c_double(colat)
      lons = ctypes.c_double(lons)
      radius = ctypes.c_double()
      lon = ctypes.c_double()
      lat = ctypes.c_double()
      libspice.sphcyl_c(r, colat, lons, ctypes.byref(radius), ctypes.byref(lon),
                        ctypes.byref(lat))
      return radius.value, lon.value, lat.value
docs: Convert from spherical coordinates to latitudinal coordinates.
  
      http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/sphlat_c.html
  
      :param r: Distance of the point from the origin.
      :type r: float
      :param colat: Angle of the point from positive z axis (radians).
      :type colat: float
      :param lons: Angle of the point from the XZ plane (radians).
      :type lons: float
      :return:
              Distance of a point from the origin,
              Angle of the point from the XZ plane in radians,
              Angle of the point from the XY plane in radians.
      :rtype: tuple