;$Id: ppower_comp_decay.pro,v 1.2 2018/06/29 01:26:54 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=8,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
end
;
;  mv idl.ps ~/tex/roper/PencilGW/fig/ppower_comp_decay.ps
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,.2]
;
dir1='M144_k10d_exp'
dir2='M144_k10b_exp'
dir3='M144_k10c_exp'
;
power,'_GWs','_mag',k=k,spec1=spec1,spec2=spec1b,i=n,tt=t1,/noplot,datatopdir='../'+dir1+'/data'
power,'_GWs','_mag',k=k,spec1=spec2,spec2=spec2b,i=n,tt=t2,/noplot,datatopdir='../'+dir2+'/data'
power,'_GWs','_mag',k=k,spec1=spec3,spec2=spec3b,i=n,tt=t3,/noplot,datatopdir='../'+dir3+'/data'
;
!x.title='!8t!6/!8t!6!d*!n'
!y.title='!8E!6(!8k,t!6)'
!p.multi=0
!p.charsize=1.2
fo='(f6.1)'
siz=1.7
xx=1.2
th2=10
th=8
k0=10.
;
ikNy=71 ;(for 144)
yr=[2e-8,2e-1]
it=10 & itd2=0
!x.range=[1,50]
print,t1[it],t2[it+itd2],t3[it]
plot_io,t1,t1^2*spec1[ikNy,*],yr=yr
oplot,t2,t2^2*spec2[ikNy,*],col=122
oplot,t3,t3^2*spec3[ikNy,*],col=55
;
oplot,t1,t1^2*spec1[ikNy/2,*],li=1
oplot,t2,t2^2*spec2[ikNy/2,*],li=1,col=122
oplot,t3,t3^2*spec3[ikNy/2,*],li=1,col=55
;
oplot,t1,t1^2*spec1[ikNy/4,*],li=2
oplot,t2,t2^2*spec2[ikNy/4,*],li=2,col=122
oplot,t3,t3^2*spec3[ikNy/4,*],li=2,col=55
;
END