;finds plots area vs. time from text file
;checkvar, ps, 0
initialize, ps

;reads ascii file and returns it as a text array
c = rd_ascii('flux_output_jul23_diff_det_6-12_0030.txt')

;removes top line which was:
;box  date  stime  etime  flux  area  centr  peak  st.dev  peak.val  Elow  Ehi  detectors
c = c[1:*]

;makes 1 line string two rows
b = str2cols(c)

;area in arcsec^2
areas = float(reform(b[5,*]))

;changes arcseconds to cm
areas = areas * 72500000.^2

;detectors = reform(b[15,*])
detectors = b[15,*]
;detectors = ['9','8-9','7-9','6-9','5-9','4-9','3-9','2-9','1-9']
;print, detectors

plot, detectors, areas, xtitle = 'FWHM resolution in arc sec', ytitle = 'area in cm!u2!n', $
title='log plot Area vs. detector  6-12 keV 00:30:00-00:31:00'

if ps gt 0 then device,/close
set_plot,'win'

end