root='/home/bcarry/data/379_Huenna/'
  dirEphem=root+'ephem/'

  id=379

  cutDist = 0.01


  tS='2008-01-01T00:00:00'
  tE='2014-01-01T00:00:00'
  step=10
  stepS=strTrim(string(step,format='(I)'),2)+'d'
  
  readcol, root+'perturbers.csv', num, name, H, i,e,a, $
           format='(I,A,F,F,F,F)', delimiter=',', /Silent
  nbP = n_elements(num)

  jdS=date_conv(tS,'JULIAN')
  jdE=date_conv(tE,'JULIAN')
  nbd= (jdE-jdS)/step


  eBin=voMiriade_CallEphemCC(id, jdS, nbd=nbd, step=stepS, /Helio, /Julian, /Web, $
                             tcoor=2, rPlane=2, mime='text',dump=dirEphem+'379.dat' ) 

  empty={num:0, name:'', dist:0., when:''}
  enc = replicate(empty,nbP)

  for kP=0, nbP-1 do begin

    file=dirEphem+strTrim(string(num[kP],format='(I)'),2)+'.dat'
    if not file_test(file,/Read) then begin
      ePer=voMiriade_CallEphemCC(num[kP], jdS, nbd=nbd, step=stepS, /Helio, /Julian, /Web, $
                                 tcoor=2, rPlane=2, mime='text',dump=file ) 
    endif else begin
      ePer=voMiriade_ReadEphemCC(file,tcoor=2, rPlane=2, mime='text',date='julian')
    endelse

    distSq = (eBin.x-ePer.x)^2. + $
             (eBin.y-ePer.y)^2. + $
             (eBin.z-ePer.z)^2.

    minDist=sqrt(min(distSq,pMin))
    enc[kP].num  = num[kP]
    enc[kP].name = name[kP]
    enc[kP].dist = minDist
    enc[kP].when = ePer[pMin].iso


    if enc[kP].dist le cutDist then begin
       print, kp, enc[kP].num, enc[kP].name, enc[kP].dist, enc[kP].when, $
              format='(I5,2x,I5,2x,A-14,2x,F9.6,2x,A-15)'
    endif else begin
       print, kp, enc[kP].num, enc[kP].name, $
              format='(I5,2x,I5,2x,A-14)'
    endelse

  endfor




end