; NOTE:  This setup file shows some
; remnants of prior data reading-in and processing.
; At some point I simplified it so it just directly reads in
; the data ready to be plotted.


; for part (c)
;;; k=read_gdf('kat19603i-dec00u')
;;; w=where(k(0,*) gt 0,nw)
;;; temp=fltarr(2,nw) & temp(0,*)=k(5,w) & temp(1,*)=k(7,w)-k(5,w)
;;; av=avgbin(temp,/more)
av=read_gdf('data/av.gdf')

; for part (b):
;;; w=where(k(0,*) gt 0 and k(5,*) gt 2.4 and k(5,*) lt 2.6)
;;; plot_hist,k(7,w)-k(5,w),res2,/noplot
res2=read_gdf('data/res2.gdf')
;;; w=where(k(0,*) gt 0 and k(5,*) gt 3.0 and k(5,*) lt 3.5)
;;; plot_hist,k(7,w)-k(5,w),res3,/noplot
res3=read_gdf('data/res3.gdf')
t2=total(res2(1,*)) & t3=total(res3(1,*))
dr2=res2(0,1)-res2(0,0) & dr3=res3(0,1)-res3(0,0)

; for part (d):
gr=read_gdf('data/gr.trin19603i')

; for part (a):
;;; w=where(k(0,*) gt 0 and k(3,*) gt 0.1 and k(4,*) gt 0.1,nw)
;;; plot_hist,k(1,w),bins=0.01,res0,/noplot
res0=read_gdf('data/res0.gdf')
sine=sin(res0(0,*)+0.006) & wsin=where(sine gt 0)

;;; w=where(k(0,*) gt 0 and k(3,*) gt 0.02 and k(4,*) gt 0.02,nw)
;;; plot_hist,k(1,w),bins=0.01,res00,/noplot
res00=read_gdf('data/res00.gdf')
sine00=sin(res00(0,*)+0.006) & wsin00=where(sine00 gt 0)

;;; w=where(k(0,*) gt 0 and k(3,*) gt 0.237 and k(4,*) gt 0.237,nw)
;;; plot_hist,k(1,w),bins=0.01,res000,/noplot
res000=read_gdf('data/res000.gdf')
sine000=sin(res000(0,*)+0.006) & wsin000=where(sine000 gt 0)

