verbose=1


;------------------------------------------------------------------------------------------;
;------------------------------------------------------------------------------------------;
;--- TAG ---    I --- General User Choices                         ------------------------;
;------------------------------------------------------------------------------------------;
;------------------------------------------------------------------------------------------;
  if verbose eq 1 then begin
     print, ''
     print, '--------------------------------------------------------------------------------'
     print, '  I. Load Configuration'
  endif


 ;--I.1-- Directory  ----------------------------------------------------------------------;
 ;--I.1.1-- Root Directory
  spawn, 'hostname', host
  case strTrim(host) of
     'hyperion': root = '/data/binary/summary/'
     'endymion': root = '/home/bcarry/work/data/binary/summary/'
;     'gnanzi':   root = 'PON_TU_PATH'
     'gnanzi':   root = '/home/mpajuelo/W/referencias/aster_data/EAR_A_COMPIL_5_BINMP_V8_0/data/'
     else: stop
  endcase
  
  



;------------------------------------------------------------------------------------------;
;------------------------------------------------------------------------------------------;
;--- TAG ---   II --- Retrieve the Binary Systems                  ------------------------;
;------------------------------------------------------------------------------------------;
;------------------------------------------------------------------------------------------;
  if verbose eq 1 then begin
     print, ''
     print, '--------------------------------------------------------------------------------'
     print, '  II. Load Data'
  endif
  bin = pds_binMP_getCatalog()
  nbBin = n_elements(bin)

  bin.dyn=strTrim(bin.dyn,2)
;  uDyn = uniq( dyn, sort(dyn) )
;  print, dyn[uDyn]
;
;  ast = ['AMO','APO','ATE','JL4','JL5','MB','MC','ML5']



  ;-select binaries with DETERMINED density and NOT ASSUMED DENSITY
  sel = where( strCmp(bin.dens.code,'-') and $
               ( strCmp(bin.dyn,'AMO') or $
                 strCmp(bin.dyn,'APO') or $
                 strCmp(bin.dyn,'ATE') or $
                 strCmp(bin.dyn,'JL4') or $
                 strCmp(bin.dyn,'JL5') or $
                 strCmp(bin.dyn,'MB') or $
                 strCmp(bin.dyn,'MC') or $
                 strCmp(bin.dyn,'ML5') ), nbSel )

 ;--- select only asteroids
  sel = where( strCmp(bin.dyn,'AMO') or $
               strCmp(bin.dyn,'APO') or $
               strCmp(bin.dyn,'ATE') or $
               strCmp(bin.dyn,'JL4') or $
               strCmp(bin.dyn,'JL5') or $
               strCmp(bin.dyn,'MB')  or $
               strCmp(bin.dyn,'MC')  or $
               strCmp(bin.dyn,'ML5') , nbSel )

 ;--- select everything
;  sel = indgen(nbBin)
;  nbSel=nbBin

  print, 'select ',nbSel,' among ',nbBin,' binaries'
;  cgHistoPlot, 100*bin[sel].dens.unc/bin[sel].dens.val, minInput=0, maxInput=300


  readcol, '/astrodata/Catalog/PDS/articles/2015-Icarus-247-Lindsay.csv', $
           delimiter=',', lNum, lName, lTax, format='(L,A,A)',/Silent
  lName = strTrim(lName)

  
;  print, bin[sel].num
  samp = bin[sel]

  taxo   = strArr(nbSel)
  scheme = strArr(nbSel)

  withNum = where( samp.num ne 0, nbNum, com=noNum, nCom=nbNoNum )

;;;;;;;;;;;;;
 ;-con numero
  for i=0, nbNum-1 do begin
     k=withNum[i]

     pL = where( samp[k].num eq lNum)
     if pL[0] ne -1 then begin
;       print, samp[k].num, lNum[pL],strTrim(lTax[pL]), format='(I6,1x,I6,1x,A-12)'
       taxo[k] = strTrim(lTax[pL])
       scheme[k] = 'Lindsay'
     endif

     info = pds_taxonomy_getCatalog(samp[k].num )
     if size(info,/type) eq 8 then begin
       taxo[k] = strTrim(info.taxo,2)
       scheme[k] = strTrim(info.scheme,2)
     endif
  endfor



  ;---sin numero
  for i=0, nbNoNum-1 do begin
     k=noNum[i]

     pL = where( samp[k].des eq lName )
     if pL[0] ne -1 then begin
;       print, samp[k].des, lName[pL],strTrim(lTax[pL]), format='(A-10,1x,A-10,1x,A-12)'
       taxo[k] = strTrim(lTax[pL])
       scheme[k] = 'Bus-DeMeo'
     endif

;     info = pds_taxonomy_getCatalog(samp[k].des )
;     if size(info,/type) eq 8 then begin
;       taxo[k] = strTrim(info.taxo,2)
;       scheme[k] = strTrim(info.scheme,2)
;     endif
  endfor
;;;;;;;;;;;;;;;;;

  taxo =strtrim(taxo,2)
  noTax = where( strCmp(taxo,'-') or strCmp(taxo,''), nbNoTax, comp=tax, ncomp=nbTax )

  print, nbTax, nbSel, nbBin

  forprint, samp.num, samp.des, scheme, taxo, format='(I6,1x,A-15,1x,A-7,1x,A-9)'



end