;plot percent increase for 61 solar flares
files = file_search('./*PI*dat')
nfiles = n_elements(files)
print, 'number of files = ',nfiles
PI = fltarr(195,nfiles)

;!p.multi = [0,1,2]
loadct, 39
set_plot,'ps'
device,/landscape,filename='PercentIncrease.eps',/color
DEVICE, DECOMPOSED = 0 
test = read_dat(files(0))
wv = test.wv
percent = test.percent
PI(*,0)= percent
plot, wv, percent,/ylog,yr = [0.1,1.e4],ystyle=1,xtitle='Wavelength(nm)',ytitle='Percentage Increase(%)',title='61 X-class Flares',pos = [0.1,0.3,.9,0.9],charsize = 1.5
for i = 1,nfiles-1 do begin
test = read_dat(files(i))
percent = test.percent
PI(*,i) = percent
oplot, wv, percent
endfor
average = total(PI,2)/nfiles
oplot,wv,average,color = 150, thick= 3

test2003301 = read_dat('./2003301PI.dat') ; Average of 34 events
percent2003301 = test2003301.percent
oplot,wv,percent2003301,color = 250, thick =3
legend,['Average','2003301'],/fill,psym=[-3,-3],colors=[150,250],/right,charsize = 1.5
device,/close
set_plot,'x'

stop
openw,1,'Average_61PI.dat'
;fmt='$(f9.2,5761e12.3)'
for i=0,194 do begin
   printf,1,wv[i],average[i]
endfor
close,1
stop

test2 = read_dat('./Average_34PI.dat') ; Average of 34 events
ave_percent = test2.percent
oplot,wv,ave_percent,color = 250, thick =2
end