;$Id: pgrah_comp_hor.pro,v 1.5 2026/01/13 09:20:18 brandenb Exp $
if !d.name eq 'PS' then begin
  ;device,xsize=18,ysize=12,yoffset=3
  device,xsize=18,ysize=6,yoffset=3
  !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6
end
;
siz=.9
!p.charsize=1.0
!x.margin=[6.8,.3]
!y.margin=[3.1,.2]
;
!p.multi=[0,2,1]
!x.range=[90.,27000.]
;
dirK1='K512sig0_k6_ramp1a'
;dirK2='K512sig01_k6_ramp1a'
dirK2='K512sig01_k6_ramp1c'
dirK3='K512sig03_k6_ramp1a'
dirK4='K512sig05_k6_ramp1a'
dirK5='K512sig1_k6_ramp1a'
;
dirM1='M512sig0_k6_ramp1a'
;dirM2='M512sig01_k6_ramp1b'
dirM2='M512sig01_k6_ramp1c'
dirM3='M512sig03_k6_ramp1a'
dirM4='M512sig05_k6_ramp1a'
dirM5='M512sig1_k6_ramp1a'
;
col2=55
col3=122
col4=155
col5=122
;
xout=110.
yout=.90
;
;-----------------------------------------------------------------------------
!x.title='!8k!6'
!y.title='!13P!8!dh!n!6(!8k!6)'
;
restore,"../"+dirK1+"/grahm.sav"
plot_oi,k[1:*],grahrm[1:*],yr=[-.0,1] 
grahre[1:*]=grahre[1:*]+2.5/k[1:*]^.5
errplot,k[1:*],grahrm[1:*]-grahre[1:*],grahrm[1:*]+grahre[1:*]
oplot,!x.range,[0.,0.]
print,'K1 ',mean(grahrm[11:13])
;
restore,"../"+dirK2+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col2
grahre[1:*]=grahre[1:*]+2.5/k[1:*]^.5
errplot,k[1:*],grahrm[1:*]-grahre[1:*],grahrm[1:*]+grahre[1:*],col=col2,li=1
print,'K2 ',mean(grahrm[11:13])
;
loadct,6
restore,"../"+dirK3+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col3
grahre[1:*]=grahre[1:*]+2.5/k[1:*]^.5
errplot,k[1:*],grahrm[1:*]-grahre[1:*],grahrm[1:*]+grahre[1:*],col=col3,li=1
loadct,5
print,'K3 ',mean(grahrm[11:13])
;
restore,"../"+dirK4+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col4
grahre[1:*]=grahre[1:*]+2.5/k[1:*]^.5
errplot,k[1:*],grahrm[1:*]-grahre[1:*],grahrm[1:*]+grahre[1:*],col=col4,li=1
print,'K4 ',mean(grahrm[11:13])
;
restore,"../"+dirK5+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col5
grahre[1:*]=grahre[1:*]+2.5/k[1:*]^.5
errplot,k[1:*],grahrm[1:*]-grahre[1:*],grahrm[1:*]+grahre[1:*],col=col5,li=1
print,'K5 ',mean(grahrm[11:13])
;
xyouts,xout,yout,siz=siz,'!6(c)'
;-----------------------------------------------------------------------------
!y.title='!6'
;
restore,"../"+dirM1+"/grahm.sav"
plot_oi,k[1:*],-grahrm[1:*],yr=[-.0,1],/nodata
grahre[1:*]=grahre[1:*]+2.0/k[1:*]^.5
errplot,k[1:*],-grahrm[1:*]-grahre[1:*],-grahrm[1:*]+grahre[1:*]
oplot,!x.range,[0.,0.]
print,'M1 ',mean(grahrm[11:13])
;
restore,"../"+dirM2+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col2
grahre[1:*]=grahre[1:*]+2.2/k[1:*]^.5
errplot,k[1:*],grahrm[1:*]-grahre[1:*],grahrm[1:*]+grahre[1:*],col=col2,li=1
print,'M2 ',mean(grahrm[11:13])
;
loadct,6
restore,"../"+dirM3+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col3
grahre[1:*]=grahre[1:*]+2.2/k[1:*]^.5
errplot,k[1:*],grahrm[1:*]-grahre[1:*],grahrm[1:*]+grahre[1:*],col=col3,li=1
loadct,5
print,'M3 ',mean(grahrm[11:13])
;
restore,"../"+dirM4+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col4
grahre[1:*]=grahre[1:*]+2.2/k[1:*]^.5
errplot,k[1:*],grahrm[1:*]-grahre[1:*],grahrm[1:*]+grahre[1:*],col=col4,li=1
print,'M4 ',mean(grahrm[11:13])
;
restore,"../"+dirM5+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col5
grahre[1:*]=grahre[1:*]+2.2/k[1:*]^.5
errplot,k[1:*],grahrm[1:*]-grahre[1:*],grahrm[1:*]+grahre[1:*],col=col5,li=1
print,'M5 ',mean(grahrm[11:13])
;
xyouts,xout,.06,siz=siz,'!6(d)'
;-----------------------------------------------------------------------------
!p.charsize=.7
!x.range=[300.,27000.]
!p.charthick=1.3 & !p.thick=1.3 & !x.thick=1.3 & !y.thick=1.3
x1=.37 & x2=.87 & dx=.11 & y1=.68 & y2=.95
;
!p.position=[x1,y1,x1+dx,y2]
plot_oo,k[1:*],grahrm[1:*],yr=[.05,1],/noerase,/nodata
;
restore,"../"+dirK2+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col2
;
loadct,6
restore,"../"+dirK3+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col3
loadct,5
;
restore,"../"+dirK5+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col5
xx=[3e3,1.8e4] & oplot,xx,5e1/xx^.55
;
;-----------------------------------------------------------------------------
!p.position=[x2,y1,x2+dx,y2]
plot_oo,k[1:*],grahrm[1:*],yr=[.05,1],/noerase,/nodata
xx=[3e3,1.8e4] & oplot,xx,7e1/xx^.55
;xyouts,2.8e3,.4,siz=siz,'!9A!8k!6!u-0.55!n'
;
restore,"../"+dirM2+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col2
;
loadct,6
restore,"../"+dirM3+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col3
loadct,5
;
restore,"../"+dirM5+"/grahm.sav"
oplot,k[1:*],grahrm[1:*],col=col5
;
!p.position=0
print,"$mv idl.ps ~/tex/teach/COSMOMAG26/2_relic_gravitational_waves/fig/pgrah_comp_hor.ps"
print,"$mv idl.ps ~/tex/sayan/SVT/fig/pgrah_comp_hor.ps"
END