pc_read_ts,obj=ts
Q=ts.brms/ts.urms
nt=n_elements(ts.t)
print,Q[0],Q[nt-1]
END