;$Id: rspec_select_hoskM_hallk2.pro,v 1.1 2023/04/30 09:36:21 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=6,yoffset=3
  !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6
  col0=0
endif else begin
  col0=255
endelse
;
;  to be run from:
;  Hf0_t2em5_k180a_Hosk
;  Hf0_t2em5_k180a_Hosk_k2
;
siz=1.0
!p.position=0
!p.charsize=1.0
!p.multi=[0,2,1]
!x.margin=[8.0,1.8]
!y.margin=[3.2,0.3]
;
@parameters
power,'_mag','_saffman_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot,/lks
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
close,1
openw,1,'xi_list.txt'
;
;---------------------------------------------
;  panel 1
;
!x.title='!6'
!x.title='!8k!6/!8k!6!d0!n'
!y.title='!6Sp(!8B!6)'
!y.title='!8E!6!dM!n(!8k!6)  and  !8kH!6!dM!n(!8k!6)/2'
yr=[3e-14,1e-4]
plot_oo,xr,yr,/nodata
col=[col0,col0,col0,55,122,155,122]
li=[2,0,1,0,0,0,0]
for j=0,n_elements(iit)-1 do begin
  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]
  ;
    circ_sym,csize,1
    oplot,k,spec1(*,i)
  ;
  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)
  ;oplot,[1,1]/xi,yr,col=col[j]
  ;
  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=[.06,.4] & oplot,xx,1.7e-6*xx^2,li=3
xyouts,.13,1.0e-10,siz=siz,'!9A!8k!6!u4!n'
xx=[.08,.7] & oplot,xx,1.3e-5*xx^0.3333,li=3
xyouts,.12,1.0e-5,siz=siz,'!9A!8k!6!u1/3!n'
;xx=[.5,1.7] & oplot,xx,4e-9/xx^2.5,li=3
;xyouts,.58,8.0e-10,siz=siz*.8,'!9A!8k!6!u-5/2!n'
xyouts,5.5,2.0e-5,siz=siz,'!6(a)'
arrow,/data,.04,1e-11,.04,5e-8
print
;
;---------------------------------------------
;  panel 2
;
lam=2e4
yr_hel=[3e-18,8e-11]
!y.title='!6Sp(!8h!6!dtot!n)'
plot_oo,xr,yr_hel,/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]
  oplot,k,spec2(*,i),col=col[j],li=li[j]
  print,i,t(i)
  wait,w
end
xx=[.064,.1] & oplot,xx,3.8e-6*xx^4,li=3
xx=[.10,.4] & oplot,xx,1.2e-9*xx^2,li=3
xyouts,.08,2.0e-11,siz=siz,'!9A!8k!6!u4!n'
xyouts,.26,4.0e-11,siz=siz,'!9A!8k!6!u2!n'
;xx=[.64,1.8] & oplot,xx,6e-12/xx^4,li=3
;xyouts,.94,3.5e-13,siz=siz*.8,'!9A!8k!6!u-4!n'
xyouts,5.3,7.2e-12,siz=siz,'!6(b)'
;
cwd,run
print,"$mv idl.ps ~/GitHub/Axel/Hosking-Shell/fig/rspec_select_hoskM_AD.eps"
print,"$mv idl.ps ~/GitHub/Axel/Hosking-Shell/fig/rspec_select_hoskM_hall_k2.eps"
print,"$mv idl.ps ~/GitHub/Axel/Hosking-Shell/fig/rspec_select_hoskM_hall_k4.eps"
END