;$Id: pstress.pro,v 1.1 2021/01/24 05:50:49 brandenb Exp $
;
power,'_Str','hel_Str',k=k,spec1=stre1,spec2=stre2,i=n,tt=t,/noplot,/lks
nt=n_elements(t)
;
!p.multi=[0,1,3]
default,t1_stre,2.5
default,t2_stre,4.5
good=where(t ge t1_stre and t le t2_stre)
if good[0] eq -1 then begin
  it1_stre=nt-2
  it2_stre=nt-2
endif else begin
  it1_stre=min(good)
  it2_stre=max(good)
endelse
print,'AXEL: t1_stre, t2_stre=',t(it1_stre),t(it2_stre)
;
;plot_oi,k[1:*],stre2m[1:*]/stre1m[1:*],yr=[0,.4]
;
strer=stre2/stre1
pc_error_range_array,t[it1_stre:*],transpose(strer[*,it1_stre:*]),mean=strerm,error=strere
pc_error_range_array,t[it1_stre:*],transpose(stre1[*,it1_stre:*]),mean=stre1m,error=stre1e
pc_error_range_array,t[it1_stre:*],transpose(stre2[*,it1_stre:*]),mean=stre2m,error=stre2e
plot_oo,k[1:*],stre2m[1:*]
errplot,k[1:*],stre2m[1:*]-stre2e[1:*],stre2m[1:*]+stre2e[1:*]
;
plot_oo,k[1:*],strerm[1:*]
errplot,k[1:*],strerm[1:*]-strere[1:*],strerm[1:*]+strere[1:*]
;
plot_oi,k[1:*],strerm[1:*],yr=[-.1,1] 
errplot,k[1:*],strerm[1:*]-strere[1:*],strerm[1:*]+strere[1:*]
;
save,file='strem.sav',t,t1_stre,it1_stre,k, strerm,strere, stre1m,stre1e, stre2m,stre2e
spawn,'cvs add -kb strem.sav'
END