;;;;
;;;
;;;  need load ssw packages
pro read_CBSmm_flux, fname, datastr, headm, headext
   ;;;  datastr   data structure
   ;;;  Fcal is center of frequency,in MHz with 100MHz band;
   
   mrd_head, fname, head
   headm = fitshead2struct( head )
   datastr = mrdfits( fname, 1, head1 )
   headext = fitshead2struct( head1 )
   
end



;;;;==============================================================


pro test_plot_CBSmm_flux

fn = file_search( 'D:\sswdb\CSOradio\CBSmm\*.fits')
read_CBSmm_flux, fn[0], datastr, headm, headext
dateobs = headext.Date_obs
freq = datastr.Fcal / 1000.0

trange = [ datastr.time00[0],datastr.time00[-1] ]
frange = [1000, 4000]
titless = 'CBSmm Flux vs Time (' + strmid( dateobs, 0, 10 ) + ' )'  
window, 0 
loadct,39
;;; plot 35.25GHz flux
plot,datastr.time00, datastr.flux00, xstyle=5,ystyle=1, xrange=trange, yrange = frange, $
      ytitle='Flux (SFU)', color=0, background= 255 ,title=titless 
;;; plot 35.75GHz flux
oplot,datastr.time01, datastr.flux01, color= 25 
;;; plot 36.25GHz flux
oplot,datastr.time02, datastr.flux02, color= 50
;;; plot 36.75GHz flux
oplot,datastr.time03, datastr.flux03, color= 75
;;; plot 37.25GHz flux
oplot,datastr.time04, datastr.flux04, color= 100
;;; plot 37.75GHz flux
oplot,datastr.time05, datastr.flux05, color= 125
;;; plot 38.25GHz flux
oplot,datastr.time06, datastr.flux06, color= 150
;;; plot 38.75GHz flux
oplot,datastr.time07, datastr.flux07, color= 175
;;; plot 39.25GHz flux
oplot,datastr.time08, datastr.flux08, color= 200
;;; plot 39.75GHz flux
oplot,datastr.time09, datastr.flux09, color= 225


;;plot time axis
timeaxis,color=0, title='Time (UT)'



stop
end