;$Id: ppxyaver_xi_comp.pro,v 1.4 2020/08/01 19:03:16 brandenb Exp $
if !d.name eq 'PS' then begin
  ysize=11
  thick=3
  device,xsize=18,ysize=ysize,yoffset=3
  !p.charthick=thick & !p.thick=thick & !x.thick=thick & !y.thick=thick
end
;
dir1='zHf0_t2em5_k180a_rep' ;(Bz)
dir2='zHf1em3_t2em5_k180a'  ;(Cz)
dir2='yHf1em3_t2em5_k180b'  ;(Cy)
dir3='zHf1em2_t2em5_k180a'  ;(Dz)
;
siz=1.7
!x.margin=[8.2,0.7]
!y.margin=[3.2,0.2]
!p.charsize=siz
!p.multi=0
!x.title='!7g!6!d0!n!8k!6!s!d0!n!r!u2!n !8t!6'
!y.title='!7n!6'
;
t0=.1
r=-.43
k0=180.
eta0=2e-5
eta0eff=eta0*t0^r
s=eta0eff*k0^2
;
th1=2
th2=5
th3=7
;
col1=55
col2=122
if !d.name eq 'PS' then col3=0 else col3=255
;
xr=[2e-3,3.4e2]
yr=[4e-4,.18]
yr=[2e-3,.12]
plot_oo,xr,yr,/nodata
file='xyaver.sav'
;
restore,'../'+dir1+'/'+file
frac2=abmz/jbmz
iz1=findex(-6.,z-!pi)
iz2=findex(-3.,z-!pi)
iz3=findex(-1.,z-!pi)
;oplot,s*tgood,sqrt(abs(frac2(iz1,*))),li=1,col=col1,th=th1 & oplot,s*tgood,sqrt(frac2(iz1,*)>0.),col=col1,th=th1
;oplot,s*tgood,sqrt(abs(frac2(iz2,*))),li=1,col=col1,th=th2 & oplot,s*tgood,sqrt(frac2(iz2,*)>0.),col=col1,th=th2
oplot,s*tgood,sqrt(abs(frac2(iz3,*))),li=1,col=col1,th=th3 & oplot,s*tgood,sqrt(frac2(iz3,*)>0.),col=col1,th=th3
;
restore,'../'+dir2+'/'+file
frac2=abmz/jbmz
iz1=findex(-6.,z-!pi)
iz2=findex(-3.,z-!pi)
iz3=findex(-1.,z-!pi)
oplot,s*tgood,sqrt(abs(frac2(iz1,*))),li=1,col=col2,th=th1 & oplot,s*tgood,sqrt(frac2(iz1,*)>0.),col=col2,th=th1
oplot,s*tgood,sqrt(abs(frac2(iz2,*))),li=1,col=col2,th=th2 & oplot,s*tgood,sqrt(frac2(iz2,*)>0.),col=col2,th=th2
oplot,s*tgood,sqrt(abs(frac2(iz3,*))),li=1,col=col2,th=th3 & oplot,s*tgood,sqrt(frac2(iz3,*)>0.),col=col2,th=th3
;
restore,'../'+dir3+'/'+file
frac2=abmz/jbmz
iz1=findex(-6.,z-!pi)
iz2=findex(-3.,z-!pi)
iz3=findex(-1.,z-!pi)
oplot,s*tgood,sqrt(abs(frac2(iz1,*))),li=1,col=col3,th=th1 & oplot,s*tgood,sqrt(frac2(iz1,*)>0.),col=col3,th=th1
oplot,s*tgood,sqrt(abs(frac2(iz2,*))),li=1,col=col3,th=th2 & oplot,s*tgood,sqrt(frac2(iz2,*)>0.),col=col3,th=th2
oplot,s*tgood,sqrt(abs(frac2(iz3,*))),li=1,col=col3,th=th3 & oplot,s*tgood,sqrt(frac2(iz3,*)>0.),col=col3,th=th3
;
;  print,2e-5/.1^.43 = 5.38307e-05
;  print,2e-5/.1^.43*180.^2 = 1.74411
;  so eta0*k0^2*t0 = 0.174
;
xx0=.174
;
; for B (blue)
;xx=[.01,26] & oplot,xx,0.0070*(xx/xx0)^.6667,col=col1,li=3
;xx=[.20,40] & oplot,xx,0.0007*(xx/xx0)^.6667,col=col1,li=3
;xyouts,12,.004,siz=siz,'!6Run Bz',col=col1
;
; for C
;xx=[.01,3.4] & oplot,xx,0.043*(xx/xx0)^.6667,col=col2,li=3
;xx=[.03,5.0] & oplot,xx,0.006*(xx/xx0)^.6667,col=col2,li=3
;xyouts,.004,.008,siz=siz,'!6Cz',col=col2
;
; for D
;xx=[.01,2] & oplot,xx,0.300*(xx/xx0)^.6667,col=col3,li=3
;xx=[.08,4] & oplot,xx,0.060*(xx/xx0)^.6667,col=col3,li=3
;xyouts,.0033,.10,siz=siz,'!6Run Dz',col=col3
;
xx=[1.,200] & oplot,xx,0.0090*(xx/xx0)^.3333,col=col3,li=3
xx=[6.,200] & oplot,xx,0.0050*(xx/xx0)^.3333,col=col3,li=3
xx=[14.,220] & oplot,xx,0.0029*(xx/xx0)^.3333,col=col3,li=3
xyouts,1.6,.030,siz=siz,'!9A!8t!6!u1/3!n',col=col3
;
default,texdir,'HallMHD'
print,'$mv idl.ps ~/tex/mhd/'+texdir+'/fig/ppxyaver_xi_comp.ps'
END