old=mrdfits('56535/apVisit-r4-6897-56535-004.fits',1)
oldw=mrdfits('56535/apVisit-r4-6897-56535-004.fits',4)
olds=mrdfits('56535/apVisit-r4-6897-56535-004.fits',5)
oldt=mrdfits('56535/apVisit-r4-6897-56535-004.fits',7)
new=mrdfits('testlsf2/apVisit-r4-6897-56535-004.fits',1)
neww=mrdfits('testlsf2/apVisit-r4-6897-56535-004.fits',4)
news=mrdfits('testlsf2/apVisit-r4-6897-56535-004.fits',5)
newt=mrdfits('testlsf2/apVisit-r4-6897-56535-004.fits',7)
set_plot,'PS'
smcolor
!p.multi=[0,1,3]

device,file='bold-004.eps',/encap,/color,xsize=72,ysize=72,/in
for i=0,2 do begin
  if i eq 0 then xr=[15140,15820]
  if i eq 1 then xr=[15850,16440]
  if i eq 2 then xr=[16470,16960]
  plot,oldw[*,2-i],old[*,2-i],yr=[0,6500],xr=xr,xstyle=1,thick=2,charsize=4,xtickformat='(I5.5)'
  oplot,oldw[*,2-i],oldt[*,2-i]*6400,color=3
  oplot,oldw[*,2-i],olds[*,2-i],color=4
  oplot,oldw[*,2-i],old[*,2-i]
endfor
device,/close
ps2jpg,'bold-004.eps',/eps,/delete
device,file='btestlsf2-004.eps',/encap,/color,xsize=72,ysize=72,/in
for i=0,2 do begin
  if i eq 0 then xr=[15140,15820]
  if i eq 1 then xr=[15850,16440]
  if i eq 2 then xr=[16470,16960]
  plot,neww[*,2-i],new[*,2-i],yr=[0,6500],xr=xr,xstyle=1,thick=2,charsize=4,xtickformat='(I5.5)'
  oplot,neww[*,2-i],newt[*,2-i]*6400,color=3
  oplot,neww[*,2-i],news[*,2-i],color=4
  oplot,neww[*,2-i],new[*,2-i]
endfor
device,/close
ps2jpg,'btestlsf2-004.eps',/eps,/delete
end