@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