;$Id: pzaver.pro,v 1.2 2013/04/02 13:06:09 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=12,ysize=21,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/koen/RFP/fig/pfrac.ps
;
!p.charsize=1.0
!x.margin=[8.8,.5]
!y.margin=[3.2,.5]
;
;!p.multi=[0,2,1]
!p.multi=[0,1,2]
;!p.multi=0
pc_read_zaver,o=o
t=o.t
nt=n_elements(t)
it2=nt-1
;surface,o.StokesQmxy[*,*,it2]/o.StokesImxy[*,*,it2]
x=findgen(512)/512.*2*!pi & x=x-mean(x)
y=findgen(512)/512.*2*!pi & y=y-mean(y)
lev=grange(-1,1,30)*1.1*max(abs(o.bzmxy[*,*,it2]))
contour,o.bzmxy[*,*,it2],x,y,lev=lev,/fil
velovect1,o.bxmxy[*,*,it2],o.bymxy[*,*,it2],x,y,/over
;contour,o.StokesImxy[*,*,it2],x,y,nlev=30,/fil
PI=o.StokesQmxy^2+o.StokesUmxy^2
RM=-(o.StokesQ1mxy*o.StokesUmxy-o.StokesU1mxy*o.StokesQmxy)/PI
contour,RM[*,*,it2],x,y,nlev=30,/fil
chi=.5*atan(o.StokesUmxy,o.StokesQmxy)
Bxrec=PI*cos(chi)
Byrec=PI*sin(chi)
velovect1,Bxrec[*,*,it2],Byrec[*,*,it2],x,y,/over,/nohead
;vel_a,o.StokesQmxy[*,*,it2],o.StokesUmxy[*,*,it2],x,y,nvec=1500,len=.002,/noerase
;vel_a,o.bxmxy[*,*,it2],o.bymxy[*,*,it2],x,y,nvec=500,len=.1,tit=''
;vel,o.StokesUmxy[*,*,it2],o.StokesQmxy[*,*,it2],nvec=500,len=.1,tit=''
;contour,o.StokesImxy[*,*,it2],tit='I',nlev=30,/fil
;contour,o.StokesQmxy[*,*,it2],tit='Q',nlev=30,/fil
;contour,o.StokesUmxy[*,*,it2],tit='U',nlev=30,/fil
END