;$Id: pkt2_one.pro,v 1.10 2020/06/20 19:58:52 brandenb Exp $
if !d.name eq 'PS' then begin
  ysize=18
  default,thick,4
  device,xsize=18,ysize=ysize,yoffset=3
  !p.charthick=thick & !p.thick=thick & !x.thick=thick & !y.thick=thick
end
;
cwd,run
siz=1.7
xstyle=0
@parameters
default,iread,0
default,ihydro,0
default,thick1,3
default,thick2,3
default,fourthird,1.
  iread=0
if iread eq 0 then begin
  power,'_mag','hel_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot
  power,'_Lor','hel_Lor',k=k,spec1=spec1h,spec2=spec2h,i=n,tt=t,/noplot
  power,'_2Lor','hel_2Lor',k=k,spec1=spec1c,spec2=spec2c,i=n,tt=t,/noplot
  iread=1
endif
nt=n_elements(t)
;
default,siz,1.7
!p.multi=[0,1,2]
half='!s!u 1!n!r!s-!r!d 2!n'
!y.title='!8E!6!dM!n!6(!8k,t!6)'
!x.title='!6'
!y.title='!6'
!x.margin=[5.0,.3]
!y.margin=[3.2,.3]
!p.charsize=siz
tilde='!9!s!aA!n!r!6'
;
default,w,.1
default,k0,60.
default,vA0,1.
s=k0/(vA0^2*k0)
xr=[.98,520]/k0
tauA=1./(k0*vA0)
print,'tauA=',tauA
pc_read_param,/param2,o=param
;eta=param.eta
E=t
D=t
II=t
IIK=t
;
default,xstyle,1
if xstyle then begin
  plot_oo,xr,yr1_one,/nodata,xtickf='logticks_exp',ytickf='logticks_exp'
endif else begin
  plot_oo,xr,yr1_one,/nodata
endelse
;
;  compute 1/k
;
k1=k & k1(0)=1. & k1=1./k1 & k1(0)=0.
EEM=total(spec1(*,it)) & xi=total(spec1(*,it)*k1)/EEM
;
oplot,k/k0,abs(s*spec2(*,it)*k/2.),thick=thick1,col=55,li=1
oplot,k/k0,s*spec2(*,it)*k/2.,thick=thick1,col=55
oplot,k/k0,s*spec1(*,it),thick=thick2,col=122
oplot,[1,1]/(k0*xi),yr1_one,li=3
;xx=[.0065,.04] & oplot,xx,.80*xx^5
;xx=[.14,1.2] & oplot,xx,1e-10/xx^2.333
;print,1e-5*xx^5
;xyouts,.012,3e-9,'!9A!8k!6!u5!n'
;xyouts,.7,7e-10,'!9A!8k!6!u-7/3!n'
default,yout_labela,1e-3
default,yout_labela,1e-3
xyouts,.0067,yout_labela,'(a)'
if run eq 'Hf0_t2em5_k180c' then begin
  xyouts,.50,1e-6,'!8E!6(!8k!6)',col=122
  xyouts,.038,1e-10,'!8kH!6(!8k!6)/2',col=55
endif else begin
  xyouts,.30,5e-6,'!8E!6(!8k!6)',col=122
  xyouts,.12,1e-7,'!8kH!6(!8k!6)/2',col=55
endelse
;
;  plot II
;
!x.title='!8k!6/!8k!6!d0!n'
default,xstyle,1
if xstyle then begin
  plot_oo,xr,yr2_one,/nodata,xtickf='logticks_exp',ytickf='logticks_exp'
endif else begin
  plot_oo,xr,yr2_one,/nodata
endelse
oplot,[1,1]/(k0*xi),yr2_one,li=3
;
oplot,k/k0,abs(s*spec1h(*,it)),li=1,thick=thick2
oplot,k/k0,   -s*spec1h(*,it), li=0,thick=thick2
loadct,6
s3=s*k/2.
oplot,k/k0,abs(s3*spec2h(*,it)),li=1,thick=thick2,col=122
oplot,k/k0,   -s3*spec2h(*,it), li=0,thick=thick2,col=122
loadct,5
;
;  plot III
;
;s2=s*eta
;oplot,k/k0,abs(s2*spec2c(*,it)*k/2.),thick=thick1,col=55,li=1
;oplot,k/k0,s2*spec2c(*,it)*k/2.,thick=thick1,col=55
;oplot,k/k0,s2*spec1c(*,it),thick=thick2,col=122
;oplot,k/k0,s2*spec1(*,it)*k^2,thick=thick2*2,col=122,li=2
;xx=[.0065,.04] & oplot,xx,4e4*xx^7
;xx=[.14,1.2] & oplot,xx,4.2e-6/xx^.3333
;xyouts,.013,1e-7,'!9A!8k!6!u7!n'
;xyouts,1.3,2e-6,'!9A!8k!6!u-1/3!n'
default,yout_labelb,1e-4
xyouts,.0067,yout_labelb,'(b)'
if run eq 'Hf0_t2em5_k180c' then begin
  xyouts,.60,4e-8,'!8T!dE!n!6(!8k!6)'
  loadct,6
  xyouts,.035,1e-12,'!8kT!dH!n!6(!8k!6)/2',col=122
  loadct,5
endif else begin
  xyouts,.14,3e-6,'!8T!dE!n!6(!8k!6)'
  loadct,6
  xyouts,.07,8e-9,'!8kT!dH!n!6(!8k!6)/2',col=122
  loadct,5
endelse
;
print,'it,t=',it,t[it]
;
cwd,run
print,'$mv idl.ps ~/tex/mhd/HallMHD/fig/pkt2_one_'+run+'.ps'
!p.multi=0
END