;$Id: ppzaver_pdft.pro,v 1.1 2013/07/18 15:47:29 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=24,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.7
!x.margin=[8.8,.5]
!y.margin=[3.2,.5]
;
restore,'pzaver_pdf.sav'
s=size(StokesQ_xx)
nt=s[2]
print,nt
;nt=9
!x.range=0
!y.range=0
cwd,run
;
kurtQ=fltarr(nt)
kurtU=fltarr(nt)
!p.multi=[0,1,2]
!p.title='!6'+run
!x.title='!6'
!y.title='!6kurt(!8Q!6)'
it=0 & dx=StokesQ_xx(1,it)-StokesQ_xx(0,it)
for it=0,nt-1 do begin
  dx=StokesQ_xx(1,it)-StokesQ_xx(0,it)
  kurtQ(it)=dx*total(StokesQ_xx(*,it)^4*StokesQ_yy(*,it))/ $
           (dx*total(StokesQ_xx(*,it)^2*StokesQ_yy(*,it)))^2
endfor
;
!p.title='!6'
!x.title='!8t!6'
!y.title='!6kurt(!8Q!6)'
it=0 & dx=StokesU_xx(1,it)-StokesU_xx(0,it)
for it=0,nt-1 do begin
  dx=StokesU_xx(1,it)-StokesU_xx(0,it)
  kurtU(it)=dx*total(StokesU_xx(*,it)^4*StokesU_yy(*,it))/ $
           (dx*total(StokesU_xx(*,it)^2*StokesU_yy(*,it)))^2
endfor
;
plot,kurtQ,ps=-1
plot,kurtU,ps=-1
;
END