t=5.
tmin=1e-3
pc_read_ts,o=ts
good=where(ts.t gt t1)
!p.multi=[0,2,2]
;
E=.5*mean(ts.urms^2)
tgood=ts.t(good)
tcomp=1./(tgood > tmin)
tcomp=tcomp*0.+1.
;
plot,ts.t(good),ts.hrms(good)/EM*tcomp
oplot,ts.t(good),accum(ts.hrms(good)/EM*tcomp)
;
plot,ts.t(good),ts.hrms(good)/EM
oplot,ts.t(good),accum(ts.hrms(good)/EM)
;
tcomp2=tcomp^2
plot,ts.t(good),tcomp2*ts.gg2m(good)/(32.*!pi*EM^2)
oplot,ts.t(good),accum(tcomp2*ts.gg2m(good)/(32.*!pi*EM^2))
;
plot,ts.t(good),ts.gg2m(good)/(32.*!pi*EM^2)
oplot,ts.t(good),accum(ts.gg2m(good)/(32.*!pi*EM^2))
;
print
print,'hrms,gg2m=',.5*max(ts.hrms(good)/EM*tcomp),.5*max(ts.gg2m(good)*tcomp^2)/(32.*!pi*EM^2)
;
!p.multi=0
END