;$Id: rspec_select_hoskM.pro,v 1.3 2026/10/27 15:01:38 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=19,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
  col0=0
endif else begin
  col0=255
endelse
;
;  for Run S (but file name was bad); to be run in:
;  1024a_mu01_k002o_rr
;
siz=1.7
siz2=1.4
!p.position=0
!p.charsize=1.7
!p.multi=[0,1,2]
!x.margin=[9.5,2.0]
!y.margin=[3.2,0.3]
;
@parameters
power,'_mag','_saffman_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot,/lks
k0=60.
w=.2
w=.02
xr=[k[1],max(k)]
th=[0.,2.,2.,4.,8.]
k1=k & k1[0]=1. & k1=1./k1 & k1[0]=0.
xxi=intarr(7)
circ_sym,1.3,1
default,csize,.2
default,csize2,.5
xlab=940.
close,1
openw,1,'xi_list.txt'
;
;---------------------------------------------
;  panel 1
;
!x.title='!6'
!y.title='!8E!6!dM!n(!8k!6,!8t!6)'
!y.title='!6Sp(!8B!6) !8k!6!d0!n/<!8B!6!s!d0!n!r!u2!n>'
col=[col0,col0,col0,55,122,155,122]
li=[2,0,1,0,0,0,0]
fact=1./total(spec1(*,0))
print,'fact=',fact
yr=[2e-10,1e-3]*fact
yr=[1e-7,2e-2]
plot_oo,xr,yr,/nodata
for j=0,n_elements(iit)-1 do begin
  if j eq 4 then loadct,6
  if j eq 5 then loadct,5
  i=iit[j]
  ;
  ;if j eq 4 then loadct,6
  ;if j eq 5 then loadct,5
  ;oplot,k,spec1(*,i),col=col[j],li=li[j]
  ;
  oplot,k,fact*spec1(*,i),col=col[j],li=li[j]
  ;xi=total(k1*spec1(*,i))/total(spec1(*,i))
  ;xxi[j]=xi
  ;print,j,i,t(i),xi
  ;printf,1,j,i,t(i),xi
  ;ik=findex(1./xi,k)
  ;
  ;circ_sym,1.3,1
  ;if j eq 4 then loadct,6
  ;if j eq 5 then loadct,5
  ;if j ge 3 then oplot,[1,1]*k[ik],[1,1]*reform(spec1(ik,i)),col=col[j],ps=8
  ;if j eq 4 then loadct,5
  ;wait,w
end
close,1
;xx=[.052,.5] & oplot,xx,fact*4.3e-5*xx^0.,li=3
xyouts,xlab,5.0e-3,siz=siz,'!6(a)'
arrow,/data,4.3,3e-7,1.8,1e-4
xx=[3.,66.] & oplot,xx,2e-4*xx,li=3
xyouts,5.,3e-3,'!9A!8k!6!u1!n',siz=siz
xx=[14.,80.] & oplot,xx,1.4e-2/xx^(7./3.),li=3
xyouts,14.,4e-7,'!9A!8k!6!u-7/3!n',siz=siz
print
;
;---------------------------------------------
;  panel 2
;
lam=2e4
!x.title='!8k!6/!8k!6!d1!n'
!y.title='!6Sp(!8h!6)'
!y.title='!6Sp(!8h!6) !8k!6!d0!n/<!8h!6!s!d0!n!r!u2!n>'
!y.title='!6(2!7p!6!u2!n/!8k!6!s!d1!n!r!u2!n) !6Sp(!8h!6) !8k!6!s!d0!n!r!u4!n/!13E!6!u2!n'
fact_old=1./total(spec2(*,0))
fact=2.*!pi^2*k0^4/total(spec1(*,0))^2
yr_hel=[3e-14,2e-7]*fact
s=fact/fact_old
yr_hel=[5e-5,.1]*s
plot_oo,xr,yr_hel,/nodata
print,'fact=',fact
for j=0,n_elements(iit)-1 do begin
  if j eq 4 then loadct,6
  if j eq 5 then loadct,5
  i=iit[j]
  oplot,k,fact*spec2(*,i),col=col[j],li=li[j]
  print,i,t(i)
  wait,w
end
xx=[1.7,30.] & oplot,xx,s*6e-5*xx,li=3
xyouts,5.,s*1.4e-4,'!9A!8k!6!u1!n',siz=siz
xyouts,xlab,s*4e-2,siz=siz,'!6(b)'
;
;!p.charsize=1.0
;!p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
;!x.title='!6'
;xx=[.052,.5] & oplot,xx,s*8.3e-8*xx^0.,li=3
;
;!p.position=[.35,.66,.54,.78]
;!y.title='!8E!6!dM!n(1,!8t!6)'
;plot_oo,t,spec1(1,*),/noerase,xr=[2e-1,3e2],yr=[1e-8,1e-6]
;xx=[10.,240.] & oplot,xx,5e-9*xx^.667,li=3
;xyouts,19.,1.7e-8,'!9A!8t!6!u2/3!n',siz=siz*.7
;
!p.position=0
cwd,run
print,"$mv idl.ps ~/Overleaf/Vinay/Helicity-Gauge/figs/rspec_select_hoskM.eps"
;spawn,'cvs add Sp_htot_k012.txt'
END