;$Id: psaff.pro,v 1.5 2022/11/24 06:07:48 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
;
n=512
siz=1.2
r=findgen(n)/n
pc_read_saffman,t,ir
!x.title='!8t !7g!8k!6!s!d0!n!r!u2!n'
!y.title='!8I!6!dH!n(t)/!8I!6!dH!n(0)'
;
;  choice of point
;
i=180
k0=180.
eta=2e-5
s=1./ir[i,0]
st=eta*k0^2
xr=[5e-5,2e2]
xr=[5e-5,90.]
yr=[1e-1,2e0]
plot_oo,st*t,s*ir[i,*],xr=xr,yr=yr
xx=[2.,50.] & oplot,st*xx,1.1/xx^.2
xyouts,5.,.8,'!9A!8t!6!u-0.2!n'
;
;  inset
;
!p.charsize=1.1
!p.position=[.23,.27,.65,.65]
!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
;
cwd,run
print,'$mv idl.ps ~/tex/mhd/HallNonhel/fig/psaff_'+run+'.eps'
!p.position=0
save,file='saff.sav',r,ir,t
spawn,'cvs add -kb saff.sav'
;
END