read_netcdf,'2012_0122-0124_fism_flare_37.nc',d,attr
day=2012022
index=where(d.date eq day)
wave=0.5*(d.wave1+d.wave2)
color8
plot_io,wave,d.sp_flux[*,index],/xstyle,xrange=[-50,2000],color=5

LOOP: repeat begin
  print, 'Enter 1 for next day'
  read,dummy
  if dummy eq 0 then stop
  day=day+1L
  str=string(day)
  index=where(d.date eq day,count)
  plot_io,wave,d.sp_flux[*,index[0]],/xstyle,xrange=[-50,2000],/nodata
  for i=0,count-1 do oplot,wave,d.sp_flux[*,index[i]],color=5
  xyouts,1200,1.e8,str
  print,d.sp_flux[*,index[0]]
endrep until day eq 2012025
end