;$Id: ppxyaver_phase2.pro,v 1.7 2018/06/04 22:01:34 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=8,yoffset=3
  thick=4
  thick=1.6
  !p.charthick=thick & !p.thick=thick & !x.thick=thick & !y.thick=thick
endif else begin
  thick=1
end
;
;  reads the down-samples xyaver.dat file
;
;  mv idl.ps ~/tex/notes/MRIdynamo/fig/64x256c8_test.ps
;  mv idl.ps ~/tex/procs/catania07/fig/phase2_64x256c8_test.ps
;  mv idl.ps ~/tex/mhd/50yrMFT/fig/phase2_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=1.2
!x.title='!8t!6/!8T!6!drot!n'
!x.margin=[6.0,1.8]
!y.margin=[3.2,.5]
;
!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))
;
!p.multi=0
!x.title='!8B!dx!n!6/!8B!6!deq!n'
!y.title='!8B!dy!n!6/!8B!6!deq!n'
bxmz_tzm=total(bxmz_tz,2)/n_elements(zzz(goodz))/Beq
bymz_tzm=total(bymz_tz,2)/n_elements(zzz(goodz))/Beq
xr=[-1,1]*.04 & yr=[-1,1]*2
xr=[-1,1]*.025 & yr=[-1,1]*1.3
plot,bxmz_tzm,bymz_tzm,li=2,xr=xr,yr=yr
;
;for it=2,nt-1,5 do begin
;  oplot,bxmz_tzm(0:it),bymz_tzm(0:it)
;  print,it
;  wait,.01
;endfor
;
it2=nt-1
it1=fix(nt*.6)
oplot,bxmz_tzm(it1:it2),bymz_tzm(it1:it2)
;
it1=fix(nt*.7)
oplot,bxmz_tzm(it1:it2),bymz_tzm(it1:it2),thick=thick+2
;
it1=fix(nt*.8)
oplot,bxmz_tzm(it1:it2),bymz_tzm(it1:it2),thick=thick+4
;
it1=fix(nt*.9)
oplot,bxmz_tzm(it1:it2),bymz_tzm(it1:it2),thick=thick+8
;
ot=.2*tt*.02
;oplot,.030*cos(ot),1.74*cos(ot+.75*!pi),li=1
oplot,.020*cos(ot),1.2*cos(ot+.75*!pi),col=122
oplot,.015*cos(ot),1.1*cos(ot+.6*!pi),col=55
;
END