;$Id: psat.pro,v 1.3 2007/02/12 10:43:46 brandenb Exp $
;
;  plots urms and brms and prints the mean values near the end
;
@info
default,t1,max(tt)/2.
;
good=where(tt gt t1)
urmsm=mean(urms(good))
brmsm=mean(brms(good))
;
!p.multi=[0,1,2]
;
plot,tt,urms
oplot,tt(good),0.*tt(good)+urmsm,col=122
;
plot,tt,brms
oplot,tt(good),0.*tt(good)+brmsm,col=122
;
print,'urmsm,brmsm=',urmsm,brmsm
!p.multi=0
END