code: def _position_and_velocity_TEME_km(self, t):
          
          sat = self.model
          minutes_past_epoch = (t._utc_float() - sat.jdsatepoch) * 1440.0
          if getattr(minutes_past_epoch, 'shape', None):
              position = []
              velocity = []
              error = []
              for m in minutes_past_epoch:
                  p, v = sgp4(sat, m)
                  position.append(p)
                  velocity.append(v)
                  error.append(sat.error_message)
              return array(position).T, array(velocity).T, error
          else:
              position, velocity = sgp4(sat, minutes_past_epoch)
              return array(position), array(velocity), sat.error_message
docs: Return the raw true equator mean equinox (TEME) vectors from SGP4.
  
          Returns a tuple of NumPy arrays ``([x y z], [xdot ydot zdot])``
          expressed in kilometers and kilometers per second.  Note that we
          assume the TLE epoch to be a UTC date, per AIAA 2006-6753.