;$Id: power_pol.pro,v 1.4 2020/08/17 17:07:35 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=15,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
!p.multi=[0,1,2]
@parameters
;
power,'_GWs','hel_GWs',k=k,spec1=grav1,spec2=grav2,i=n,tt=t,/noplot,/lks
;yr_power_pol=[3e-17,3e-9] & yr_pol=[0,1]
;yr_power_pol=[3e-17,3e-2] & yr_pol=[-1,1]
;
it1=1
it3=10
nt=n_elements(t)
it2=nt-1
col=grange(55,200,nt)
;
!x.title='!6'
!y.title='!13E!6!dGW!n(!8k!6)'
plot_oo,k[1:*],0*grav1[1:*,it2],yr=yr_power_pol
loadct,6
for it=it1,it2,it3 do oplot,k[1:*],grav1[1:*,it],col=fix(col(it))
loadct,5
;
!x.title='!8k!6'
!y.title='!13P!6!dGW!n(!8k!6)'
plot_oi,k[1:*],0*grav1[1:*,it2],yr=yr_pol
loadct,6
for it=it1,it2,it3 do oplot,k[1:*],grav2[1:*,it]/grav1[1:*,it],col=col(it)
loadct,5
;
cwd,run
print,"$mv idl.ps ~/tex/notes/GW_decay_compens/fig/pcomp_Ft_"+run+".ps"
END