;ls ./2011/*.sav>2011fismlist.dat

bb = 295.
ee = 305.
ndays = ee-bb ;365
ntimes = ndays*1440L
nbins = 195
time = fltarr(ntimes)
fx_total = fltarr(ntimes)

close,1
dir = './2011/'
openr,1,'./2011fismlist.dat'
str = ' '

;need a container to hold one line of text
trash = 'garbage'
;let's skip the first few lines 
header = bb-1
for jj = 0,header-1 do readf,1,trash

for i = 0, ndays-1 do begin
    readf, 1, str
    fn = str ;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]

    fx_total(i*1440L:(i+1)*1440L-1) = total(fism_pred,2,/nan)  ;W/m2
    time(i*1440L:(i+1)*1440L-1)= i+bb+double(utc)/double(3600L*24.)   ;days(1.000~365.999)
endfor
close,1

set_plot,'ps'
device,file ='2011295_305_radiance.eps',/color
loadct, 39

plot, time, fx_total,xstyle=1,ystyle=1, xr=[bb,ee],yr=[0.059,0.063], $
  title = 'Solar irradiance'+strmid(string(DATE),5,4), $
  xtitle ='DOY', ytitle = 'W/m2 (0-195 nm)'

device,/close
set_plot,'x'

end