;$Id: ppxyaver_phase.pro,v 1.3 2009/09/17 05:10:30 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=18,yoffset=3
  !p.charthick=4 & !p.thick=4 & !x.thick=4 & !y.thick=4
end
;
;  reads the down-samples xyaver.dat file
;
;  mv idl.ps ~/tex/notes/MRIdynamo/fig/64x512a.ps
;  mv idl.ps ~/tex/notes/MRIdynamo/fig/64x256c7.ps
;  mv idl.ps ~/tex/notes/MRIdynamo/fig/64x256c8.ps
;  mv idl.ps ~/tex/notes/MRIdynamo/fig/64x256c8_test.ps
;  mv idl.ps ~/tex/procs/catania07/fig/phase_64x256c8_test.ps
;
;restore,'xyaver.sav'
;
nt=n_elements(tt)
nz=n_elements(zzz)
;
Trot=2.*!pi/Omega
Hdisc=cs0/Omega*sqrt(2.)
Hdisc=cs0/Omega
!p.charsize=2.4
!x.title='!8t!6/!8T!6!drot!n'
;
!p.multi=[0,1,3]
z2d=zzz#replicate(1.,nt)
;
;zgood=1.9*Hdisc
;good=where(z2d gt zgood)
;goodz=where(zzz gt zgood)
;
zgood1=-.5*Hdisc & zgood2=+.5*Hdisc
zgood1=-.2*Hdisc & zgood2=+.2*Hdisc
good =where(z2d gt zgood1 and z2d lt zgood2)
goodz=where(zzz gt zgood1 and zzz lt zgood2)
;
fact=30.
Beq=.06
bxmz_tz=transpose(reform(bxmz(good),n_elements(zzz(goodz)),nt))
bymz_tz=transpose(reform(bymz(good),n_elements(zzz(goodz)),nt))
;
!y.title='!8z!6/!8H!6'
contour,bxmz_tz,tt/Trot,zzz(goodz)/Hdisc,nlev=30,/fil
contour,bymz_tz,tt/Trot,zzz(goodz)/Hdisc,nlev=30,/fil
;
!y.title='!8B!di!n!6/!8B!6!deq!n'
plot,tt/Trot,total(bymz_tz,2)/n_elements(zzz(goodz))/Beq
oplot,tt/Trot,total(bxmz_tz,2)/n_elements(zzz(goodz))/Beq*fact,li=2
oplot,tt/Trot,tt*0.
;
END