;$Id: ppxyaver_comp.pro,v 1.4 2020/08/01 19:03:16 brandenb Exp $
if !d.name eq 'PS' then begin
  ysize=18
  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)
dir3='zHf1em2_t2em5_k180a'  ;(Dz)
dir4='yHf1em3_t2em5_k180b'  ;(Cy)
;
siz=1.7
!x.margin=[8.2,0.7]
!y.margin=[3.2,0.2]
!p.charsize=siz
!p.multi=[0,1,2]
!x.title='!6'
!y.title='!7u!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,3e2]
yr=[4e-4,1.8]
plot_oo,xr,yr,/nodata
file='xyaver.sav'
;
restore,'../'+dir1+'/'+file
frac2=abmz*jbmz/b2mz^2
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/b2mz^2
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/b2mz^2
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
;
restore,'../'+dir4+'/'+file
frac2=abmz*jbmz/b2mz^2
iz1=findex(-6.,z-!pi)
iz2=findex(-3.,z-!pi)
iz3=findex(-1.,z-!pi)
loadct,6
oplot,s*tgood,sqrt(abs(frac2(iz1,*))),li=1,col=col3,th=th1 & oplot,s*tgood,sqrt(frac2(iz1,*)>0.),col=col2,th=th1
oplot,s*tgood,sqrt(abs(frac2(iz2,*))),li=1,col=col3,th=th2 & oplot,s*tgood,sqrt(frac2(iz2,*)>0.),col=col2,th=th2
oplot,s*tgood,sqrt(abs(frac2(iz3,*))),li=1,col=col3,th=th3 & oplot,s*tgood,sqrt(frac2(iz3,*)>0.),col=col2,th=th3
loadct,5
;
;  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.0020*(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 Cz
xx=[.01,3.4] & oplot,xx,0.043*(xx/xx0)^.6667,col=col2,li=3
;xx=[.03,5.0] & oplot,xx,0.012*(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 CZ
loadct,6
xx=[.004,8.8] & oplot,xx,0.013*(xx/xx0)^.6667,col=col2,li=3
xx=[.03,5.0] & oplot,xx,0.0033*(xx/xx0)^.6667,col=col2,li=3
xx=[3.,100.] & oplot,xx,0.0003*(xx/xx0)^.6667,col=col2,li=3
xyouts,.004,.008,siz=siz,'!6Cz',col=col2
loadct,5
;
; for D
xx=[.01,2] & oplot,xx,0.300*(xx/xx0)^.6667,col=col3,li=3
;xx=[.08,4] & oplot,xx,0.120*(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
xyouts,.003,.75,siz=siz,'!6(a)'
;
;-----------------------------------------------------------------------------
!x.title='!7g!6!d0!n!8k!6!s!d1!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
;
restore,'../'+dir4+'/'+file
frac2=abmz/jbmz
iz1=findex(-6.,z-!pi)
iz2=findex(-3.,z-!pi)
iz3=findex(-1.,z-!pi)
loadct,6
oplot,s*tgood,sqrt(abs(frac2(iz1,*))),li=1,col=col3,th=th1 & oplot,s*tgood,sqrt(frac2(iz1,*)>0.),col=col2,th=th1
oplot,s*tgood,sqrt(abs(frac2(iz2,*))),li=1,col=col3,th=th2 & oplot,s*tgood,sqrt(frac2(iz2,*)>0.),col=col2,th=th2
oplot,s*tgood,sqrt(abs(frac2(iz3,*))),li=1,col=col3,th=th3 & oplot,s*tgood,sqrt(frac2(iz3,*)>0.),col=col2,th=th3
loadct,5
;
;  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
loadct,6
xx=[.8,60] & oplot,xx,0.0046*(xx/xx0)^.3333,col=col2,li=3,th=5
xx=[.7,10] & oplot,xx,0.0022*(xx/xx0)^.3333,col=col2,li=3,th=5
xx=[8.,80] & oplot,xx,0.0012*(xx/xx0)^.3333,col=col2,li=3,th=5
loadct,5
xyouts,1.6,.030,siz=siz,'!9A!8t!6!u1/3!n',col=col3
xyouts,.003,.075,siz=siz,'!6(b)'
;
default,texdir,'HallMHD'
print,'$mv idl.ps ~/tex/mhd/'+texdir+'/fig/ppxyaver_comp.ps'
!p.multi=0
END