default,tgood,.2
;
pc_read_ts,obj=ts
tt=ts.t
ppm=ts.ppm
TTm=ts.TTm
rhom=ts.rhom
plot,tt,ppm,ps=-1
;
tgood=tgood < max(tt)
good=where(tt ge tgood)
ppm0=mean(ppm(good))
print,'ppm0=',ppm0
oplot,tt(good),0*tt(good)+ppm0,col=122,li=2
;
print,'tgood=',tgood
END