;  
;  $Id: pgrowth.pro,v 1.2 2013/10/18 09:16:15 brandenb Exp $
;  for plotting brms and local growth rate
;  a summary of the results are in roberts.dat in this directory
;
t1=30. & t2=1e9
!p.multi=[0,1,2]
pc_read_ts,o=ts
tt=ts.t
brms=ts.brms
tmax=max(tt)
lam=deriv(tt,alog(brms))
default,t1,(1<(tmax/2.))
default,t2,(4<(tmax/2.))
@parameters
;
;  brms
;
plot_io,tt,brms,xr=[0.,2*t2]
;
;  growth rate
;
good=where(tt gt t1 and tt lt t2)
if good(0) ne -1 then begin
  lamm=mean(lam(good))
  print,'lam=',lamm
endif
;
yr=[-1,1]*6.*max(lamm)
plot,tt,lam,yr=yr,xtit='!8t!6',ytit='',xr=[0.,2*t2]
oplot,tt,tt*0
;
good=where(tt gt t1 and tt lt t2)
if good(0) ne -1 then begin
  oplot,tt(good),lam(good),col=122
  lamm=mean(lam(good))
  oplot,tt(good),tt(good)-tt(good)+lamm,li=3,col=188
  print,'lam=',lamm
endif else begin
  print,'time series still too short for getting growth rate'
endelse
;
END