
red344 = redindices(344)
red492 = redindices(492)
redall = [red344,red492]

splot,redall.ew_mg2,redall.modelcolor[0],ps=4
soplot,redall[i].ew_mg2,redall[i].modelcolor[0],ps=1,color='red',syms=2

i344=where(redall.plate EQ 344 AND redall.redsample)
i492=where(redall.plate EQ 492 AND redall.redsample)
i344=where(redall.plate EQ 344)
i492=where(redall.plate EQ 492)

!p.multi = [0,1,2]
xrange = [0,3]
yrange = [-12,-2]
model_kcorr = redall.modelcolor[0] - 2*redall.z
model_dered = model_kcorr  - redall.reddening[0] + redall.reddening[1]
djs_plot,model_kcorr[i344],redall[i344].ew_mg2,ps=1,color='red', $
 xr=xrange,yr=yrange
djs_oplot,model_kcorr[i492],redall[i492].ew_mg2,ps=4
djs_plot,model_dered[i344],redall[i344].ew_mg2,ps=1,color='red', $
 xr=xrange,yr=yrange
djs_oplot,model_dered[i492],redall[i492].ew_mg2,ps=4

;------------------------------------------------------------------------------
plate = [414,415,416,428,464,453,488,489,491,492]
redall = redindices(plate[0])
for i=1,n_elements(plate)-1 do redall = [redall, redindices(plate[i])]
save, file='redall.ss'

restore, 'redall.ss'
i1 = where(redall.reddening[0] LT 0.25 AND redall.redsample)
i2 = where(redall.reddening[0] GT 0.25 AND redall.redsample)
djs_plot,redall[i1].ew_mg2,redall[i1].modelcolor[0],ps=4
djs_oplot,redall[i2].ew_mg2,redall[i2].modelcolor[0],ps=4,color='red'

model_kcorr = redall.modelcolor[0] - 2*redall.z
model_dered = model_kcorr  - redall.reddening[0] + redall.reddening[1]
!p.multi = [0,1,2]
yrange=[1.0,2.5]
djs_plot,redall[i1].ew_mg2,model_kcorr[i1],ps=1,yr=yrange,/ystyle
djs_oplot,redall[i2].ew_mg2,model_kcorr[i1],ps=1,color='red'
djs_plot,redall[i1].ew_mg2,model_dered[i1],ps=1,yr=yrange,/ystyle
djs_oplot,redall[i2].ew_mg2,model_dered[i1],ps=1,color='red'

; Find residuals to Mg-color
acoeff = linfit(redall[i1].ew_mg2,model_dered[i1])
cresid = model_kcorr - (acoeff[0] + acoeff[1] * redall.ew_mg2)
ugpredict = redall.reddening[0] - redall.reddening[1]
j = where(redall.redsample)
!p.multi=0
djs_plot,ugpredict[j],cresid[j],ps=1, $
 xtitle='SFD (u-g)',ytitle='Residual (u-g)',charsize=2
djs_oplot, [0,1], [0,1]

;----------------------
; Try a dispersion-color relation instead of Mg-color...
i=where(redall.reddening[0] LT 0.25 AND redall.redsample)
acoeff=linfit(redall[i].vdisp,redall[i].modelcolor[0])
splot,alog10(redall[i].vdisp),redall[i].modelcolor[0],ps=1
splot,alog10(redall[i].vdisp),redall[i].modelcolor[0] $
- acoeff[0] - acoeff[1]*alog10(redall[i].vdisp),ps=1
cresid = redall.modelcolor[0] - (acoeff[0] + acoeff[1] * alog10(redall.vdisp))
ugpredict = redall.reddening[0] - redall.reddening[1]
j = where(redall.redsample)
splot,ugpredict[j],cresid[j],ps=4

