;
pro sflux, channel,ps=ps
;
;channel = 'G'
wavelength = 'all'
uarsday = 1050
units = 'photons/s/cm^2/A'

get_l3_hr, channel,wavelength,uarsday,l3hrdat,flux_units=units
;
; l3hrdat:
;   uday = int uday(m) array  of uars day numbers of the columns of flux data
;   wave = float array(n) of wavelengths of rows of flux data
;   flux(m,n) = fluxes
;   sf(m,n) = errors
;
nuflux = n_elements(l3hrdat.flux)
nuwave = n_elements(l3hrdat.wave)
print,'nuflux=',nuflux,' nuwave=',nuwave
;for i=0,nuflux-1 do begin
;  print,'i=',i,' wavelen=',l3hrdat.wave(i),' flux=',l3hrdat.flux(i)
;endfor
;
if keyword_set(ps) then begin
  set_plot,'PS'
  flnm = 'uars' + '.ch' + channel + '.ps'
  device,file=flnm
endif else begin
  set_plot,'X'
endelse
xtitle='WAVELENGTH (' + strcompress(string(l3hrdat.wave(nuflux-1))) + $
  ' to ' + strcompress(string(l3hrdat.wave(0))) + ')'
ytitle='UARS FLUX (' + units + ')'
toplab = 'UARS DAY ' + strcompress(string(uarsday)) + '  (CHANNEL ' + $
  channel + ')'
plot_io,l3hrdat.wave,l3hrdat.flux,xtitle=xtitle,ytitle=ytitle,title=toplab
; yrange=[1.e7,1.e12]	; appropriate for channel G
if keyword_set(ps) then begin
  device,/close
endif
return
end