;$Id: rspec_select_comp_AD.pro,v 1.4 2023/05/14 17:12:48 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=9,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
;
!x.margin=[8.0,0.3]
!y.margin=[3.1,0.3]
;
;  labels
;
si2=.8
siz=1.
xx1=180. & yy1=1.4e-6 & yy2=5e-12
;
;  list
;
dir1='Hf0_t2em5_k180a_Hosk_k2'
dir2='Hf0_t2em5_k180a_Hosk'
dir3='ADsig0_k180c'
;
power,/lk,'_mag','_saffman_mag',k=k,spec1=spec1a,spec2=spec2a,i=n,tt=ta,/noplot,datatopdir='../'+dir1+'/data'
power,/lk,'_mag','_saffman_mag',k=k,spec1=spec1b,spec2=spec2b,i=n,tt=tb,/noplot,datatopdir='../'+dir2+'/data'
power,/lk,'_mag','_saffman_mag',k=k,spec1=spec1c,spec2=spec2c,i=n,tt=tc,/noplot,datatopdir='../'+dir3+'/data'
;
!x.range=[1.,512.]
!y.range=[3e-11,6e-6]
!p.multi=[0,3,2]
;
help,ta
!x.title='!6'
!y.title='!6Sp(!8B!6)'
plot_oo,k,k,/nodata
for it=33,20,-3 do oplot,k,spec1a[*,it]
xx=[3.,10.] & oplot,xx,3e-8*xx^2
xyouts,3.0,1e-6,siz=si2,'!8k!6!u2!n'
xyouts,xx1,yy1,siz=siz,'!6(a)'
;
help,tb
!y.title='!6'
plot_oo,k,k,/nodata
for it=38,20,-4 do oplot,k,spec1b[*,it]
;xx=[3.,10.] & oplot,xx,3e-8*xx^2
xx=[2.,16.] & oplot,xx,2.3e-8*xx^1.5
xyouts,3.0,6e-7,siz=si2,'!8k!6!u3/2!n'
xyouts,xx1,yy1,siz=siz,'!6(b)'
;
help,tc
plot_oo,k,k,/nodata
for it=47,30,-3 do oplot,k,spec1c[*,it]
xx=[4.,22.] & oplot,xx,1e-8*xx^2
xyouts,5.0,7e-7,siz=si2,'!8k!6!u2!n'
xyouts,xx1,yy1,siz=siz,'!6(c)'
;
!x.title='!8k!6/!8k!6!d1!n'
!y.title='!6Sp(!8h!6)'
!y.range=[3e-17,3e-11]
help,ta
plot_oo,k,k,/nodata
for it=33,20,-3 do oplot,k,spec2a[*,it]
xx=[1.2,6.] & oplot,xx,2e-13*xx^2
xyouts,1.7,2e-12,siz=si2,'!8k!6!u2!n'
xyouts,xx1,yy2,siz=siz,'!6(d)'
;
help,tb
!y.title='!6'
plot_oo,k,k,/nodata
for it=38,20,-4 do oplot,k,spec2b[*,it]
xx=[1.2,6.] & oplot,xx,6e-14*xx^2
xyouts,2.0,1e-12,siz=si2,'!8k!6!u2!n'
xyouts,xx1,yy2,siz=siz,'!6(e)'
;
help,tc
plot_oo,k,k,/nodata
for it=47,30,-3 do oplot,k,spec2c[*,it]
xx=[1.2,16.] & oplot,xx,4e-14*xx^2
xyouts,3.5,2e-12,siz=si2,'!8k!6!u2!n'
xyouts,xx1,yy2,siz=siz,'!6(f)'
;
cwd,run
print,"$mv idl.ps ~/GitHub/Axel/Hosking-Shell/fig/rspec_select_comp_AD.eps"
END