;Readin fism_60sec.sav data and calculate the integated solar
;irradiance(W/m2)
;UTC:second

restore,'./2003/FISM_60sec_2003301_v01_02.sav',/verbose
aa = total(fism_pred(*,0:194),2,/nan)
fp_mx = where(aa eq max(aa))
fp_mn = fp_mx-120 ;2 hours before
print, 'fp_mx=',fp_mx,'utc_mx=',utc(fp_mx)/3600.,'utc_mn=',utc(fp_mn)/3600.
fp_mn = 1
bb=fism_pred(fp_mx,*)-fism_pred(fp_mn,*)
bb = bb/fism_pred(fp_mn,*)

;stop

;loadct, 39
;set_plot,'ps'
;device,filename='2003301PI.eps',decomposed=0,/color
;device,/portrait
;device,filename='1989075PI.eps'
!p.multi =[0,1,3]
plot, utc/3600.,aa,xr=[0,24],yr=[0.04,0.12],xstyle=1,xtitle='Time(hour)',ytitle='W/m2 (0-190nm)',$
      title='Solar irradiance'+strmid(string(DATE),5,7),charsize = 1.8, pos=[0.1,0.7,0.9,0.9]
plots,utc(fp_mn)/3600.,aa(fp_mn),psym=5
plots,utc(fp_mx)/3600.,aa(fp_mx),psym=4
plot, fism_wv, fism_pred(fp_mn,*),/ylog,yr = [1.e-6, 1.e-2],ystyle=1,xtitle='Wavelength(nm)',ytitle='Solar Flux(W/m2/nm)',charsize = 1.8, pos=[0.1,0.4,0.9,0.6]
oplot, fism_wv, fism_pred(fp_mx,*),color = 250
xyouts,150,1.e-5,'Pre-Flare'
xyouts,150, 4.e-6,'Flare Peak',color = 250
plot, fism_wv, bb*100,/ylog,yr = [0.1,1.e4],ystyle=1,xtitle='Wavelength(nm)',ytitle='Percentage Increase(%)',charsize = 1.8, pos=[0.1,0.1,0.9,0.3]
;device,/close
;set_plot,'x'
!p.multi = 0.

;stop
;;save percentage increase to a data file
;nbins = n_elements(fism_wv)
;openw,1,'2002111PI.dat'
;;fmt='$(f9.2,5761e12.3)'
;for i=0,nbins-1 do begin
;   printf,1,fism_wv(i),bb(i)*100
;endfor
;close,1

;stop


;ls ./2003/*.sav>2003fismlist.dat
ndays = 1 ;365
ntimes = ndays*1440L
nbins = 195
fx_0_14 = fltarr(ntimes)
fx_14_25 = fltarr(ntimes)
fx_25_105 = fltarr(ntimes)
fx_105_120 = fltarr(ntimes)
fx_121 = fltarr(ntimes)
fx_122_195 = fltarr(ntimes)
time = fltarr(ntimes)
;
close,1
dir = './2003/'
openr,1,'./2003/2003fismlist.dat'
str = ' '

for i = 0, ndays-1 do begin
    readf, 1, str
    fn = dir+str
    print, fn
    fn = dir+'FISM_60sec_2003301_v01_02.sav'
    ;restore,'./2003/FISM_60sec_2003301_v01_02.sav',/verbose
    restore,fn,/verbose
;DATE, utc(1440) ;sec 
;fism_wv(195) ; 0.5:1:194.5 nm
;fism_pred(1440, 195) ;W/m2/nm, [time, bins]
;sradiance(ndays*1440),sum each of the columns of fism_pred.
;    xuv_flux(i*1440L:(i+1)*1440L-1) = total(fism_pred(*,0:24),2,/nan) ;W/m2
;    euv_flux(i*1440L:(i+1)*1440L-1) = total(fism_pred(*,25:120),2,/nan) ;W/m2
;    lyman_alpha(i*1440L:(i+1)*1440L-1) = fism_pred(*,121) ;121.5nm
;    fuv_flux(i*1440L:(i+1)*1440L-1) = total(fism_pred(*,122:194),2,/nan) ;W/m2
    fx_0_14(i*1440L:(i+1)*1440L-1) = total(fism_pred(*,0:13),2,/nan) ;W/m2
    fx_14_25(i*1440L:(i+1)*1440L-1) = total(fism_pred(*,14:24),2,/nan) ;W/m
    fx_25_105(i*1440L:(i+1)*1440L-1) = total(fism_pred(*,25:104),2,/nan) ;W/m2
    fx_105_120(i*1440L:(i+1)*1440L-1) = total(fism_pred(*,105:120),2,/nan) ;W/m2
    fx_121(i*1440L:(i+1)*1440L-1) = fism_pred(*,121) ;W/m2
    fx_122_195(i*1440L:(i+1)*1440L-1) = total(fism_pred(*,122:174),2,/nan) ;total(fism_pred(*,122:194),2,/nan) ;W/m2
    fx_total = total(fism_pred(*,0:174),2,/nan)
    time(i*1440L:(i+1)*1440L-1)= i+1+double(utc)/double(3600L*24.)   ;days(1.000~365.999)
endfor
close,1
;fx_total = fx_0_14+fx_14_25+fx_25_105+fx_105_120+fx_121+fx_122_195
stop

openw,12,'SolarEnergy2003301_new.dat',width = 160
;;fmt='$(f9.2,5761e12.3)'
printf, 12, 'time 0-14, 14-25, 25-105, 105-120, 121, 122-175 0-175 nm'
for i=0,n_elements(time)-1 do begin
   printf,12,time(i),fx_0_14(i), fx_14_25(i),fx_25_105(i), fx_105_120(i),fx_121(i),fx_122_195(i),fx_total(i)
endfor
close,12

stop


;set_plot,'ps'
;device,filename='sflux301.eps',decomposed=0,/color
;loadct,39
;plot, time,total,xr=[301, 302],yr=[0, 0.11],xtitle='time',ytitle='W/m2',title = 'Solar irradiance(0-195nm)'
;oplot, time, fx_0_14,color=50
;oplot, time, fx_14_25,color=100      
;oplot, time, fx_25_105,color=150
;oplot, time, fx_105_120, color = 200
;oplot, time, fx_121,color = 250, linestyle = 2
;oplot, time, fx_122_195, color = 0, linestyle = 2  
;xyouts, 301.8, 0.065,'total',color = 0 
;xyouts,301.8,0.001,'0-14nm',color=50  
;xyouts,301.8,0.00185,'14-25nm',color=100   
;xyouts,301.8,0.0035,'25-105nm',color=150  
;xyouts,301.8,0.00065,'105-120nm',color=200
;xyouts,301.8,0.009,'121nm',color=250
;xyouts,301.8,0.05,'122-195nm',color=0
;device,/close
;set_plot,'x'

stop

;set_plot,'ps'
;device,file ='2003_fism_solar_radiance.eps',/color
;loadct, 39
;plot, time, total,yr = [0.0, 0.1],xr=[0, 420],xstyle =1,title = 'Solar irradiance'+strmid(string(DATE),5,4),xtitle ='DOY', ytitle = 'W/m2 (0-195 nm)'
;oplot, time, xuv_flux, color = 50
;oplot, time, euv_flux, color = 250
;oplot, time, lyman_alpha, color = 150
;oplot, time, fuv_flux, color = 200
;;items = ['total','XUV','EUV','Lyman_'+Greek('alpha'),'FUV']
;;sym = [3,3,3,3,3]
;;legend,items,psym = sym, colors=[0,50,100,150,200]
;xyouts,365, 0.055,'total',color = 0
;xyouts,365, 0.001,'XUV',color = 50
;xyouts,365, 0.003,'EUV',color = 250
;xyouts,365, 0.008,'Lyman_'+Greek('alpha'),color = 150
;xyouts,365, 0.045,'FUV',color = 200
;device,/close
;set_plot,'x'

;set_plot,'ps'
;device,file ='2003_fism_percentage.eps',/color
;loadct, 39
;plot, time, fuv_flux/total_flux*100,yr = [0, 100],xr=[0, 420],xstyle =1,title = 'Percentage of total flux'+strmid(string(DATE),5,4),xtitle ='DOY', ytitle = 'percentage'
;oplot, time, xuv_flux/total_flux*100, color = 50
;oplot, time, euv_flux/total_flux*100, color = 250
;oplot, time, lyman_alpha/total_flux*100, color = 150
;xyouts,365, 2,'XUV',color = 50
;xyouts,365, 5,'EUV',color = 250
;xyouts,365, 13,'Lyman_'+Greek('alpha'),color = 150
;xyouts,365, 75,'FUV'
;device,/close
;set_plot,'x'

end