visDIR = '/observ/vista/'
  filterDir = visDIR+'filters/'
  demeoDir  = visDIR+'bus-demeo/'

  doplot=01


 ;---I -- Mean reflectance from Bus-DeMeo
  demeo  = taxoRead( demeoDIR+'meanspectra.tab' )
 ;-             Z      Y   J       H       K
  filtwave = [0.877,1.020,1.252  ,1.645  ,2.147]

  taxo = taxo2vista( demeo, filterDIR, /plot )


 ;---II-- Read Bus-DeMeo summary Table
  readfmt, demeoDir+'tableDeMeo', 'I6,4x,A20,22x,A3', bdNum, bdName, bdClass, /silent
    bdName = strtrim(bdName,2)
    bdClass= strtrim(bdClass,2)
  uniqClass = bdClass(uniq(bdClass,sort(bdClass)))

  cClass=['A','B','C','Cb','Cg','Cgh','Ch','D','K','L','O','Q','R','S','Sa','Sq','Sqw','Sr','Srw','Sv','Svw','Sw','T','V','Vw','X','Xc','Xe','Xk']
  cID   =[ 1 , 2 , 3 , 3  , 3  , 3   , 3  , 4 , 5 , 5 , 6 , 6 , 6 , 6 ,  6 ,  6 ,  6  ,  6 ,  6  ,  6 ,  6  ,  6 , 7 , 8 , 8  , 7  , 7  , 7 , 7  ]
  cLabel=['U','A','B','C','D','K','S','X','V']


 ;---III-- Reflectance of each Bus-DeMeo spectrum
  spawn, 'ls -1 '+demeoDir+'a*.txt', list
  N=n_elements(list)

  ast = replicate({num:0L, name:'', class:'', classId:0, $
                   cID:0, cName:'', $
	           ref:fltarr(5), unc:fltarr(5), col:fltarr(5), $
                   c:{z:0., y:0., j:0., h:0., k:0.}, $
                   u:{z:0., y:0., j:0., h:0., k:0.}  },N)


  for kA=0, N-1 do begin
    print, kA, ' ', list(kA)
    cols = spec2vista( list(kA), demeoDIR, plot=doPlot, dirFig=demeoDir+'figs/' )

    
    ast(kA).num    = bdNum(kA)
    ast(kA).name   = bdName(kA)
    ast(kA).class  = bdClass(kA)
    ast(kA).classId= where(bdClass(kA) eq uniqClass)

    ast(kA).cID    = cID( where( bdClass(kA) eq cClass ) )
    ast(kA).cName  = cLabel(ast(kA).cID)
 
    ast(kA).ref = cols.ref
    ast(kA).unc = cols.unc
    ast(kA).col = cols.col

    print, ast(kA).num,",", ast(kA).name,",", ast(kA).class,",", ast(kA).cName


  endfor



  forprint, ast.num, ast.name, ast.class, $
            ast.ref(0), ast.ref(1), ast.ref(2), ast.ref(3), ast.ref(4), $
            ast.col(0), ast.col(1), ast.col(2), ast.col(3), ast.col(4), $
            ast.unc(0), ast.unc(1), ast.unc(2), ast.unc(3), ast.unc(4), $
            ast.classId, $
            ast.cName, $
            ast.cID, $
            format='(I-6,",",A-15,",",A-3,10(",",F7.4),5(",",F7.4),",",I2,",",A1,",",I1)', $
            textout=demeoDir+'colors.csv', $
            comment='Num, Name, Class, rZ, rY, rJ, rH, rK, uZ, uY, uJ, uH, uK, cZ, cY, cJ, cH, cK, classId,complex,complexId', $
            /silent

  forprint, ast.num, ast.name, ast.class, $
            ast.ref(0), ast.ref(1), ast.ref(2), ast.ref(3), ast.ref(4), $
            ast.col(0), ast.col(1), ast.col(2), ast.col(3), ast.col(4), $
            ast.unc(0), ast.unc(1), ast.unc(2), ast.unc(3), ast.unc(4), $
            ast.classId, $
            format='(I-6,1x,A-15,1x,A-3,10(1x,F7.4),5(1x,F7.4),1x,I2,1x,A1,1x,I1)', $
            textout=demeoDir+'colors.dat', $
            comment='# Num Name Class rZ rY rJ rH rK uZ uY uJ uH uK cZ cY cJ cH cK classId complex complexId', $
            /silent


end