loadct,0

plot,fltarr(n_elements(color))-20,xtitle='',ytitle='',title='',ps=3, $
  xr=[-1,2.5],yr=[-20,20],ytick_get=ytick,ystyle=4

loadct,6

colormm=minmax(color)
range=colormm[1]-colormm[0]
slope=(range/255.)
colorarr=findgen(255)/255.*200.

wlast=colormm[0]
for i=0,255 do begin

wplot=where(color lt (i*slope+colormm[0]) and color gt wlast)
wlast=(i*slope+colormm[0])

if wplot[0] ne -1 and n_elements(wplot) gt 1 then begin
oplot,color[wplot],aMag[wplot],color=colorarr[255-i],ps=8
endif

endfor

loadct,0

plot,fltarr(n_elements(color))-20,xtitle='Color (B-V)',ytitle='Absolute Visual Magnitude',title='HR-Diagram',ps=3, $
  xr=[-1,2.5],yr=[-20,20],ytickname=(-1)*strcompress(fix(ytick),/remo),/noerase