;$Id: psaff_comp.pro,v 1.4 2022/12/26 09:55:08 brandenb Exp $
if !d.name eq 'PS' then begin
  ysize=11
  thick=1.6
  device,xsize=18,ysize=ysize,yoffset=3
  !p.charthick=thick & !p.thick=thick & !x.thick=thick & !y.thick=thick
end
;
!x.margin=[6.8,.3]
!y.margin=[3.2,.3]
!p.charsize=1.2
;
siz=1.2
;
dir1='Hf0_t2em5_k180a_Hosk'
dir2='Hf0_t1em5_k180a_Hosk'
;
!x.title='!8t !7g!8k!6!s!d0!n!r!u2!n'
!y.title='!8I!6!dH!n(t)/!8I!6!dH!n(0)'
restore,'../'+dir1+'/saff.sav'
n=n_elements(r)
;
;  choice of point
;
i=180
i=205
i=305
k0=180.
eta=2e-5
s=1./ir[i,0]
st=eta*k0^2
xr=[5e-5,85.]
yr=[1e-1,2e0]
plot_oo,st*t,s*ir[i,*],xr=xr,yr=yr,li=2
xx=[.5,40.] & oplot,st*xx,1.26/xx^.11,li=3
xx=[.9,38.] & oplot,st*xx,0.96/xx^.16,li=3
xyouts, 5.,1.1,'!9A!8t!6!u-0.11!n'
xyouts,23.,0.48,'!9A!8t!6!u-0.16!n'
;
;  overplot second data set
;
restore,'../'+dir2+'/saff.sav'
n=n_elements(r)
eta=1e-5
s=1./ir[i,0]
st=eta*k0^2
oplot,st*t,s*ir[i,*]
;
;  inset
;
!p.charsize=1.1
!p.position=[.23,.27,.65,.65]
!x.title='!8k!6!d0!n !8R!6'
!y.title='!13I!6!dH!n(!8t!6,!8R!6)/!8I!6!dH!n(0)'
xr=[2e-3,1e0]
yr=[5e-2,2e0]
nt=n_elements(t)
plot_oo,/noerase,r[1:*],s*ir[1:*,0],xr=xr,yr=yr
for it=0,nt-1 do oplot,r[1:*],s*ir[1:*,it]
oplot,[1,1]*r[i],yr,li=3
arrow,/data,.03,.2,.2,.1
xyouts,.22,.09,'!8t!6'
;
cwd,run
print,'$mv idl.ps ~/tex/mhd/HallNonhel/fig/psaff_comp.eps'
!p.position=0
save,file='saff.sav',r,ir,t
spawn,'cvs add -kb saff.sav'
;
END