@parameters
default,kin,0
pc_read_ts,obj=ts
;
!p.charsize=3
!p.multi=[0,1,3]
plot,ts.t,(ts.ggT2m+ts.ggX2m)/(ts.hhT2m+ts.hhX2m),ps=-1
;
!y.title='!8h!6!drms!n'
if kin then begin
  plot,ts.t,16*!pi*k0^2*sqrt(accum(ts.hhT2m+ts.hhX2m))/(.5*ts.urms^2),ps=-1
endif else begin
  plot,ts.t,16*!pi*k0^2*sqrt(accum(ts.hhT2m+ts.hhX2m))/(.5*ts.brms^2),ps=-1
endelse
print,'1) ',2*sqrt(6.)*!pi
;
!y.title='!13E!6!dGW!n'
if kin then begin
  plot,ts.t,accum(8*!pi*k0^2*(ts.ggT2m+ts.ggX2m)/(.5*ts.urms^2)^2,2),ps=-1
endif else begin
  plot,ts.t,accum(8*!pi*k0^2*(ts.ggT2m+ts.ggX2m)/(.5*ts.brms^2)^2),ps=-1
endelse
print,'2) ',!pi
;
!p.multi=0
;
END