root='/home/bcarry/data/mining/muse/'
  list='sci.csv'

  dirSB = root+'fovs/'

  cutDIT = 300.
  fov = 1.5 / 60.  ;-1.5 arcmin
  MJDtoJD = 2400000.5



  readcol, root+list, object, ra, dec, filter, exptime, mjd, progId, $
           survey, id, cat, type, tech, release, instId, delim=',', $
           format='(A,F,F,A,F,D,A,A,A,A,A,A,A,A)', /silent
  nbMuse=n_elements(ra)






  ecLon = fltarr(nbMuse)
  ecLat = fltarr(nbMuse)

  for kM=0, nbMuse-1 do begin
     if dec[kM] eq 0 then dec[kM]+=1e-6
    ec = frameCoord_eq2ec( ra[kM], dec[kM] )
    ecLon[kM]=ec[0]
    ecLat[kM]=ec[1]
;    print, ra[kM], dec[kM], ecLon[kM], ecLat[kM]
  endfor

  sel20 = where( abs(ecLat) le 20 and exptime le cutDIT, nb20 )
  sel10 = where( abs(ecLat) le 10 and exptime le cutDIT, nb10 )
  print, nbMuse, nb20, nb10

;  cgWindow
;  cgPlot, ecLon, ecLat, psym='Open circle', /AddCmd, $
;          xStyle=1, xRange=[0,360],    $
;          yStyle=1, yRange=[-90, 90]
;
;  cgPlot, /over, ecLon[sel20], ecLat[sel20], symSize=1.3, psym='Filled circle', /AddCmd, color='Cornflower blue'
;  cgPlot, /over, ecLon[sel10], ecLat[sel10], symSize=1.3, psym='Filled circle', /AddCmd, color='Red'


  sel  = sel20
  nbSel=nb20

  sso=0
  nbSSo=0
  for kSel=0, nbSel-1 do begin


    if kSel mod 20 eq 0 then print, kSel

   ;---- Get skybot
    sb = voSkyBoT_ConeSearch( ra[sel[kSel]], dec[sel[kSel]], fov, mjd[sel[kSel]]+MJDtoJD, $
                              observer='309', from='BCarry-MUSE')


    dimS=size(sb)
    if dimS[dimS[0]+1] eq 8 then begin
      sso=[sso,kSel]
      nbSSo++
    endif

  endfor
  sso=sso[1:nbSSo]

  forprint, object, ra, dec, exptime, mjd, mjd+MJDtoJD, progId, $
            id, cat, type, release, $
            format='(A-20,1x,F9.5,1x,F9.5,1x,F9.4,1x,D12.4,1x,D12.4,1x,A-14,1x,'+$
            'A-30,1x,A-10,1x,A-10,1x,A-12)', textout=dirSB+'sso.dat', $
            subset=sel[sso]


end