;$Id: pspec_irro.pro,v 1.2 2026/07/27 01:13:14 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=11,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
  col0=0
endif else begin
  col0=255
end
;
siz=1.6
xout=18.
!p.charsize=1.7
!x.margin=[8.3,.5]
!y.margin=[3.2,.5]
!p.multi=0
;
@parameters
power,'o','_kin',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot,/lks
power,'_lr','_ro',k=k,spec1=spec3,spec2=spec4,i=n,tt=t,/noplot,/lks
pc_read_param,obj=param,/param2
pc_read_ts,obj=ts
nu=param.nu
;
;-----------------------------------------------------------------------------
k0=2.5
!x.title='!8k!6/!8k!d!7m!n!6'
!y.title='!6Sp(!8u!6!dirro!n) !8k!6!d0!n/!8c!6!s!ds!n!r!u2!n'
!x.range=[2.8e-2,3.3e1]
default,w,.01
yr_spec2=[1e-16,2e-1]
plot_oo,k[1:*],abs(spec2[1:*,0]-.5*spec2[1:*,0]/k[1:*]^2),yr=yr_spec2,/nodata
;
for j=0,n_elements(iit)-1 do begin
  it=iit[j]
  if ct[j] eq 6 then loadct,6
  if ct[j] eq 5 then loadct,5
  rhom=ts.rhom(findex(t[j],ts.t))
  epsK=ts.epsK(findex(t[j],ts.t))
  knu=(epsK/(rhom*nu^3))^.25
  ;oplot,k/knu,abs(spec2[*,it]-.5*spec2[1:*,0]/k[1:*]^2)*k[1:*]^1.6667/epsK^.66667,li=li[j],col=col[j]
  oplot,k/knu,abs(spec2[*,it]-.5*spec2[1:*,0]/k[1:*]^2)/k0,li=li[j],col=col[j]
  print,t[it]
  wait,w
endfor
siz=1.7
;xx=[10.,50.]/16. & oplot,xx,3e-4/xx^(11./3.),li=3
;xyouts,siz=siz,.26,3e-3,'!9A!8k!6!u-11/3!n'
xx=[10.,50.]/16. & oplot,xx,3e-6/xx^(11./3.),li=3
xyouts,siz=siz,.20,1e-6,'!9A!8k!6!u-11/3!n'
;
xx=.05 & dx=xx & yy=3e-9 & siz=1.6 & fac=16.
legend,xx,dx,yy,1,siz=siz,'!8t!6 = 0, 10!u-4!n, 0.1'
legend,xx,dx,yy/fac,0,siz=siz,'!8t!6 = 0.6'
legend,xx,dx,yy/fac^2,0,siz=siz,'!8t!6 = 0.9',col=55 & loadct,6
legend,xx,dx,yy/fac^3,0,siz=siz,'!8t!6 = 1.3',col=144 & loadct,5
legend,xx,dx,yy/fac^4,0,siz=siz,'!8t!6 = 1.6',col=155
legend,xx,dx,yy/fac^5,0,siz=siz,'!8t!6 = 2.5',col=122
;
cwd,run
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/pspec_irro_"+run+".eps"
END