code: def spkgeo(targ, et, ref, obs):
      
      targ = ctypes.c_int(targ)
      et = ctypes.c_double(et)
      ref = stypes.stringToCharP(ref)
      obs = ctypes.c_int(obs)
      state = stypes.emptyDoubleVector(6)
      lt = ctypes.c_double()
      libspice.spkgeo_c(targ, et, ref, obs, state, ctypes.byref(lt))
      return stypes.cVectorToPython(state), lt.value
docs: Compute the geometric state (position and velocity) of a target
      body relative to an observing body.
  
      http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/spkgeo_c.html
  
      :param targ: Target body.
      :type targ: int
      :param et: Target epoch.
      :type et: float
      :param ref: Target reference frame.
      :type ref: str
      :param obs: Observing body.
      :type obs: int
      :return: State of target, Light time.
      :rtype: tuple