;$Id: ppkft.pro,v 1.11 2019/01/28 10:32:03 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=10,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/tina/decdynamo/fig/ppkft.ps
;  This routine is supposed to be run from the Run A directory,
;  which is: KH1152tnuk4b_sig1b_M01
;
!p.charsize=1.5
!x.margin=[8.4,.5]
!y.margin=[3.2,.5]
;
dir1='KH1152tnuk4b_sig1b_M01'
dir2='KH1152tnuk4b_sig1b_M01_rep'
dir3='KH1152tnuk4c_sig1b_M01'
dir4='KH1152tnuk4d_sig1b_M01'
dir5='KH1152tnuk4e_sig1b_M01'  ;(Run H)
file='pkf.sav'
;
k0=60.
nexp=1.05
urms0=1.750E-01
tau=1./(urms0*k0)
s=2./(urms0^2/k0)
s2=2./(urms0^2/k0^nexp)
!x.title='!8t!6/!7s!6!d0!n'
!y.title='!6<!8B!6!u2!n>!7n!6!s!dM!n!r!u!7b!6!n / (!7q!6!d0!n!8u!6!s!u2!n!r!drms0!n !7n!6!d0!n)'
;
xr=[.1,56000]
restore,'../'+dir2+'/'+file
t1=t & EM1=EM & kI1=kI
restore,'../'+dir1+'/'+file
plot_oo,[t1,t[0],.5*(t[0]+t[1]),t[1:*]]/tau,s*[EM1,EM[0],.52*(EM[0]+EM[1]),EM[1:*]]/[kI1,kI[0],.5*(kI[0]+kI[1]),kI[1:*]],xr=xr,yr=[3e-4,1.4e-2]
oplot,[t1,t[0],.5*(t[0]+t[1]),t[1:*]]/tau,s2*[EM1,EM[0],.52*(EM[0]+EM[1]),EM[1:*]]/[kI1,kI[0],.5*(kI[0]+kI[1]),kI[1:*]]^nexp,li=1
print,t
;
restore,'../'+dir3+'/'+file
oplot,t/tau,s*EM/kI,col=122,li=1,thick=8
;oplot,t/tau,s2*EM/kI^nexp,col=122,li=1,thick=1
;
restore,'../'+dir4+'/'+file
oplot,t/tau,s*EM/kI,col=55,li=2,thick=5
;oplot,t/tau,s2*EM/kI^nexp,col=55,li=2,thick=1
;
loadct,6
restore,'../'+dir5+'/'+file
oplot,t/tau,s*EM/kI,col=122,li=3,thick=5
xyouts,500.,.00034,siz=siz,'Run H',col=122
loadct,5
;
xyouts,.60,.0089,siz=siz,'Run A'
xyouts,4000,.0008,siz=siz,'Run A!c!7b!6=0.05'
xyouts,40.,.006,siz=siz,'F',col=122
xyouts,4.,.0036,siz=siz,'Run G',col=55
;
print,'$mv idl.ps ~/tex/tina/decdynamo/fig/ppkft.ps'
nu=1e-5*(t>1.)^(-.43)
nu=2e-5*(t>1.)^(-.43)
print,"plot_oi,t/tau,s*EK^.5/kI/nu,yr=[7e4,1e5]"
;
;  inset
;
!p.charsize=.9
!p.charthick=2.4
!p.position=[.30,.3,.5,.5]
;
xr=[.0,10]
yr=[1e-3,1.4e-2]
restore,'../'+dir2+'/'+file
t1=t & EM1=EM & kI1=kI
restore,'../'+dir1+'/'+file
plot_io,[t1,t[0],.5*(t[0]+t[1]),t[1:*]]/tau,s*[EM1,EM[0],.52*(EM[0]+EM[1]),EM[1:*]]/[kI1,kI[0],.5*(kI[0]+kI[1]),kI[1:*]],xr=xr,yr=yr,/noerase
oplot,[t1,t[0],.5*(t[0]+t[1]),t[1:*]]/tau,s2*[EM1,EM[0],.52*(EM[0]+EM[1]),EM[1:*]]/[kI1,kI[0],.5*(kI[0]+kI[1]),kI[1:*]]^nexp,li=1
print,t
;
restore,'../'+dir3+'/'+file
oplot,t/tau,s*EM/kI,col=122,li=1,thick=8
;oplot,t/tau,s2*EM/kI^nexp,col=122,li=1,thick=1
;
restore,'../'+dir4+'/'+file
oplot,t/tau,s*EM/kI,col=55,li=2,thick=5
;oplot,t/tau,s2*EM/kI^nexp,col=55,li=2,thick=1
;
;  slope
;
xx=[.5,1.2]
oplot,xx,.0030*exp(.45*xx)
;
print,'$mv idl.ps ~/tex/tina/decdynamo/fig/ppkft.ps'
nu=1e-5*(t>1.)^(-.43)
nu=2e-5*(t>1.)^(-.43)
print,"plot_oi,t/tau,s*EK^.5/kI/nu,yr=[7e4,1e5]"
;
!p.position=0
END