PRO JAN20spectPLOT, PS=ps, AGU=agu, ENERGY_EDGES=energy_edges
;
; Reads, plots and ratios spectra for Jan 20, 2005 event as provided by GS, 11-16-05
;
; 17-Nov-05		Initial version
; 26-Nov-05	gh	Add /AGU keyword to produce plot for AGU poster
;				Add ENERGY_EDGES keyword to print average over energy values
;
IF KEYWORD_SET(ps) THEN set_plot, 'ps'
datin 		= READ_ASCII('JAN20spectrumGS.txt', DATA_START=6)
kev 		= datin.field1[0,*]
totspect 	= datin.field1[1,*]
bremspect 	= datin.field1[2,*]
nucfraction	= 1 - bremspect/totspect
IF KEYWORD_SET(agu) THEN BEGIN
	!P.MULTI=[0,1,2]
	PLOT, kev, nucfraction, /XLOG, YRANGE=[0,1], TITLE='Nuclear Fraction of Solar Flux', $
		XTITLE='Energy (keV)', YTITLE='Fraction', CHARSIZE=1.5, CHARTHICK=3, THICK=3, XTHICK=3, YTHICK=3
		XYOUTS, 120, 0.90, '20 Jan 2005'
		ENDIF ELSE BEGIN
	!P.MULTI=[0,2,2]
	PLOT,  kev, totspect, /XLOG, /YLOG,      XTITLE='Energy(keV)',TITLE='20 Jan 2005 06:44:40-06:50:00'
	OPLOT, kev, bremspect
	PLOT,  kev, (totspect-bremspect)>1.e-5,/XLOG,/YLOG, XTITLE='Energy(keV)',YTITLE='nuclear', TITLE='20 Jan 2005 06:44:40-06:50:00'
	PLOT,  kev, nucfraction, /XLOG, XTITLE='Energy(keV)', YTITLE='Nuclear Fraction of Solar Flux', $
			TITLE='20 Jan 2005 06:44:40-06:50:00'
	PLOT, kev, 1 - bremspect/totspect, XTITLE='Energy(keV)', YTITLE='Nuclear Fraction of Solar Flux', $
			TITLE='20 Jan 2005 06:44:40-06:50:00'
	!P.MULTI=[0,1,1]
ENDELSE
IF KEYWORD_SET(ps) THEN DEVICE,/CLOSE
nrg = N_ELEMENTS(energy_edges)
IF nrg EQ 0 THEN RETURN
FOR n = 0, nrg-2 DO BEGIN
	ok = WHERE(kev GE energy_edges(n) AND kev LE energy_edges(n+1), nok)
	IF nok EQ 0 THEN CONTINUE
	avfract = AVERAGE(nucfraction[ok])
	fmax = MAX(nucfraction[ok])
	fmin = MIN(nucfraction[ok])
	PRINT, energy_edges(n), energy_edges(n+1), fmin, fmax, avfract
ENDFOR
RETURN
END