;$Id: pspec_select_hor.pro,v 1.1 2024/03/10 18:15:10 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
end
;
;  mv idl.ps ~/tex/roper/GW/fig/pspecm_scl2.ps
;  mv idl.ps ~/tex/prop/computing/pdc/18/fig/pspecm_scl3n.ps
;
@parameters
!p.charsize=.8
!x.margin=[8.2,0.5]
!y.margin=[3.1,0.4]
!p.multi=[0,2,1]
circ_sym,1.3,1
siz=.9
si2=.9
;
default,regime,2
print,"$sed.csh data/param2.nml"
pc_read_param,obj=param
pc_read_param,obj=param2,/param2
;
default,iread,0
if iread eq 0 then begin
  power,'_mag','_Str',k=k,spec1=spec1,spec2=stre1,i=n,tt=t,/noplot,/lks
  power,'_GWh','_GWs',k=k,spec1=grah1,spec2=grav1,i=n,tt=t,/noplot,/lks
  iread=1
endif
s=1d0/k[1] ;(scaling factor)
nk=n_elements(k)
sGW=s/6.
;
ni=n_elements(iii)
li=indgen(ni) mod 5
col=indgen(ni)/5
;
col(where(col eq 1))=55
col(where(col eq 2))=122
;
;  panel 1
;
!x.title='!8k!6'
!y.title='!8E!6!dM!n(!8k,t!6)'
xr=minmax(k[1:*])
xr2=[1e2,1.4e4]
yr=[2e-27,1]
yr=[1e-27,1e-3]
plot_oo,xr,yr,/nodata
for j=0,n_elements(iii)-1 do begin
  i=iii[j]
  oplot,k[1:*],s*spec1(1:*,i),li=li[j],col=col[j]
  print,'i,li,col,t=',i,li[j],col[j],t(i)-1.
endfor
if regime eq 2 then begin
  xx=[120.,800.] & oplot,xx,2e-15*(xx/100.)^4
  xyouts,900.,2e-11,'!9A!8k!6!u4!n',siz=siz
  i=99 & oplot,k[50]*[1,1],s*spec1(50,i)*[1,1],col=122,ps=8
endif else begin
  xx=[120.,1400.] & oplot,xx,5e-14*(xx/100.)^4
  xyouts,1500.,8e-11,'!9A!8k!6!u4!n',siz=siz
  if i eq 1920 then oplot,k[1]*[1,1]*1.1,s*spec1(1,i)*[1,1],col=122,ps=8
endelse
mu=param.mu5_const
lam=param2.lambda5
eta=param2.eta
lameta2=lam*eta^2
oplot,xr,xr*0+16.*mu*eta^2,li=1,thick=6
oplot,xr2,xr2*0+1.0*mu/lam,li=3
oplot,[1,1]*mu/2.,yr,li=2
oplot,[1,1]*mu,yr,li=1
xyouts,15000.,4e-6,'!6(a)',siz=siz
;
;  panel 3
;
!y.title='!8E!6!dGW!n(!8k,t!6)'
i=iii[0]
yr=[1e-37,1e-6]
yr=[1d-45,1e-12] ;(regime I)
yr=[1d-45,1e-10] ;(regime II)
plot_oo,k[1:*],k(1:*),yr=yr,/nodata
for j=0,n_elements(iii)-1 do begin
  i=iii[j]
  oplot,k[1:*],sgw*grav1(1:*,i),li=li[j],col=col[j]
endfor
if regime eq 2 then begin
  xx=[120.,4400.] & oplot,xx,1e-30*(xx/100.)^2
  xyouts,4800.,9e-29,'!9A!8k!6!u2!n',siz=si2
  xx=180. & arrow,xx,1e-18,xx,1e-12,/data
endif else begin
  xx=[120.,7400.] & oplot,xx,1e-23/(xx/100.)^2
  xyouts,1500.,2e-29,'!9A!8k!6!u-2!n',siz=si2
endelse
oplot,[1,1]*mu,yr,li=1
xyouts,15000.,4e-14,'!6(b)',siz=siz
;
cwd,run
print,'$mv idl.ps ~/Overleaf/Axel/CosmoBField/fig/pspec_select_hor_'+run+'.eps'
!p.multi=0
END