;$Id: pzaver_pdf.pro,v 1.2 2014/03/02 21:14:38 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
;
npdf=100
StokesQ_xx=fltarr(npdf,nt) & StokesQ_yy=fltarr(npdf,nt)
StokesU_xx=fltarr(npdf,nt) & StokesU_yy=fltarr(npdf,nt)
;
for it=0,nt-1 do begin
  print,it
  pdf,o.StokesQmxy(*,*,it),xx,yy,ndim=npdf
  StokesQ_xx(*,it)=xx & StokesQ_yy(*,it)=yy
  ;
  pdf,o.StokesUmxy(*,*,it),xx,yy,ndim=npdf
  StokesU_xx(*,it)=xx & StokesU_yy(*,it)=yy
endfor
;
save,file='pzaver_pdf.sav',StokesQ_xx,StokesQ_yy,StokesU_xx,StokesU_yy
END