!p.multi=[0,1,1]

e=readfromfile('eigen.out')
e3=readfromfile('eigen-300.out')
e7=readfromfile('eigen-700.out')

print,e(1,0),e3(1,0)
print,e(1,1),e3(1,1)

th=3
sm=21

for i=0,1 do begin

  a=readfromfile('pc'+lzp(i+1,2)+'.out')
  a(1,*)=a(1,*)*e(1,i)
  a(1,*)=a(1,*)-avg(a(1,*))
  plot,a(0,*),smooth(a(1,*),sm),th=th

  b=readfromfile('pc'+lzp(i+1,2)+'.out-300')
  b(1,*)=b(1,*)*e3(1,i)
  b(1,*)=b(1,*)-avg(b(1,*))
  oplot,b(0,*),smooth(b(1,*),sm),col=2,th=th

  c=readfromfile('pc'+lzp(i+1,2)+'.out-700')
  corr=correlate(c(1,*),a(1,*))
  if (corr gt 0) then fac=1 else fac=-1
  print,corr,fac
  c(1,*)=fac*c(1,*)*e7(1,i)
  c(1,*)=c(1,*)-avg(c(1,*))
  oplot,c(0,*),smooth(c(1,*),sm),col=3,th=th

  oplot,1700*[1,1],makerange(a(1,*)),th=3,col=2
  oplot,1300*[1,1],makerange(a(1,*)),th=3,col=3
  oplot,1902*[1,1],makerange(a(1,*)),th=3

  abc_label,'PC: '+lzp(i,1)

  gettwogifs,out='~/np/wmc/e'+lzp(i+1,1)

endfor

end