plot,[0,1],/nodata,xr=[1000,2000],yr=[-.03,.03],/ys

cols=[!p.color,indgen(5)+1]

w=15

for i=3,0,-1 do begin

  a=readfromfile('pc'+lzp(i+1,2)+'.out')
  ys=a(0,*)
  a1=smooth1(a(1,*),w=w,ys=ys)
  oplot,ys,a1,col=cols(i),th=3

endfor

gettwogifs,out='~/np/wmc/a'

end