;$Id: ppxyaver.pro,v 1.12 2020/07/28 20:28:35 brandenb Exp $
if !d.name eq 'PS' then begin
  ysize=24
  default,thick,4
  device,xsize=18,ysize=ysize,yoffset=3
  !p.charthick=thick & !p.thick=thick & !x.thick=thick & !y.thick=thick
end
;
siz=1.7
!x.margin=[9.0,0.7]
!y.margin=[3.2,0.2]
!p.multi=[0,1,3]
!p.charsize=3.2
;
file='xyaver.sav'
restore,file
eta0=2e-5 & k0=180.
eps0=eta0^3*k0^2
B20=eta0^2
j1=4
j3=5
w=0.
xout=-6.12
for j=j1,n_elements(tgood)-1,j3 do print,tgood[j]
;
;-----------------------------------------------------------------------------
;  panel 1
;
!x.title='!6'
!y.title='!7e!6/[!7e!6]'
yr=[1e2,1e10]
plot_io,z-!pi,yr,/nodata
for j=j1,n_elements(tgood)-1,j3 do begin
  oplot,z-!pi,epsMmz[*,j]/eps0
  print,tgood[j]
  wait,w
endfor
xyouts,xout,6e2,siz=siz,'!6(a)'
xyouts,-3.8,400,siz=siz,'!8t!6'
arrow,/data,-3.0,1e9,-4.,1e3
;
;-----------------------------------------------------------------------------
;  panel 2
;
!y.title='!8B!6!drms!n/[!8B!6!drms!n]'
yr=[2e2,2e4]
plot_io,z-!pi,yr,/nodata
for j=j1,n_elements(tgood)-1,j3 do begin
  oplot,z-!pi,sqrt(b2mz[*,j]/B20)
  wait,w
endfor
xyouts,xout,3e2,siz=siz,'!6(b)'
xyouts,-4.0,3e2,siz=siz,'!8t!6'
arrow,/data,-3.5,1e4,-4.2,3e2
;
;-----------------------------------------------------------------------------
;  panel 3
;
frac2=abmz*jbmz/b2mz^2
iz1=findex(-6.,z-!pi)
iz2=findex(-3.,z-!pi)
iz3=findex(-1.,z-!pi)
!x.title='!8k!6!d1!n!8z!6'
!y.title='!7u!6'
plot_io,z-!pi,z,yr=[2e-2,2.0],/nodata
for j=j1,n_elements(tgood)-1,j3 do begin
  oplot,z-!pi,sqrt(abmz[*,j]*jbmz[*,j])/b2mz[*,j]
  wait,w
endfor
xyouts,xout*1.01,1.26,siz=siz*.9,'!6(c)'
xyouts,-5.0,1.2,siz=siz,'!8t!6'
arrow,/data,-3.5,.03,-5.2,1.3
;
cwd,run
default,texdir,'HallMHD'
print,'$mv idl.ps ~/tex/mhd/'+texdir+'/fig/ppxyaver.ps'
print,'$mv idl.ps ~/tex/mhd/'+texdir+'/fig/ppxyaver_'+run+'.ps'
print,'.r ppxyaver_parametric'
!p.multi=0
END