code: def _parse_call_args(self,*args,**kwargs):
          
          interp= kwargs.get('interp',self._useInterp)
          if len(args) == 5:
              raise IOError("Must specify phi for streamdf")
          elif len(args) == 6:
              if kwargs.get('aAInput',False):
                  if isinstance(args[0],(int,float,numpy.float32,numpy.float64)):
                      out= numpy.empty((6,1))
                  else:
                      out= numpy.empty((6,len(args[0])))
                  for ii in range(6):
                      out[ii,:]= args[ii]
                  return out
              else:
                  return self._approxaA(*args,interp=interp)
          elif isinstance(args[0],Orbit):
              o= args[0]
              return self._approxaA(o.R(),o.vR(),o.vT(),o.z(),o.vz(),o.phi(),
                                    interp=interp)
          elif isinstance(args[0],list) and isinstance(args[0][0],Orbit):
              R, vR, vT, z, vz, phi= [], [], [], [], [], []
              for o in args[0]:
                  R.append(o.R())
                  vR.append(o.vR())
                  vT.append(o.vT())
                  z.append(o.z())
                  vz.append(o.vz())
                  phi.append(o.phi())
              return self._approxaA(numpy.array(R),numpy.array(vR),
                                    numpy.array(vT),numpy.array(z),
                                    numpy.array(vz),numpy.array(phi),
                                    interp=interp)
docs: Helper function to parse the arguments to the __call__ and related functions,
          return [6,nobj] array of frequencies (:3) and angles (3:)