; Load the data.
filein,"data/TSCAL_220105_W.raw.vegas"
; Look at its contents.
summary
; We have four Nod scans at different frequencies,
; covering the range of frequencies available to the
; W-band RX.
; The W-band RX does not have a noide diode (currently),
; and at this point in time GBTIDL does not know
; how to calibrate data without a noise diode.
; So, we need to load staff contributed data reduction
; routines. Hopefully the new software will have this
; capability out of the box (i.e., no need to load 
; additional scripts to work with this kind of data).
; The first script determines the counts-to-Kelvin
; conversion factor from observations of a hot and cold 
; loads.
.compile /home/astro-util/projects/4mm/PRO/calseq_4ch.pro
.compile /home/astro-util/projects/4mm/PRO/calseq.pro
; This one handles the Nods using the derived conversion 
; factors.
; For more details read GBT memo #302.
.compile /home/astro-util/projects/4mm/PRO/wnod.pro
.compile /home/astro-util/projects/4mm/PRO/wonoff_gain.pro
; Derive conversion factors and system temperatures.
calseq_4ch,23
; Put the gains in a vector. These are the values
; reported by calseq_4ch,23 (lovely user experience).
gains=[9.7830632e-06,2.7649641e-06,6.2953872e-06,1.0926164e-06]
; Now calibrate the Nods.
wnod,24,gains=gains
; This will calibrate and average all the on-source data.
; Copy the results to a different data container.
copy,0,10
; Repeat for the rest of the frequencies/scans.
calseq_4ch,26
gains=[1.0146364e-05,3.0414764e-06,8.7260057e-06,4.9196119e-06]
wnod,27,gains=gains
copy,0,11
calseq_4ch,29
gains=[8.7696096e-06,2.2162095e-07,1.2493783e-06,1.8103832e-07]
wnod,30,gains=gains
copy,0,12
calseq_4ch,32
gains=[3.2101523e-05,2.3262233e-06,9.6356663e-06,3.4275235e-07]
wnod,33,gains=gains
copy,0,13
; Now plot them all together.
freeze
copy,10,0
show
copy,11,0
oshow
copy,12,0
oshow
copy,13,0
oshow
; We can see that the temperature scale looks reasonable
; (all callibrated data has similar antenna temperatures),
; and that there are gaps in the frequency coverage.
; Save the plot for future reference.
write_ps,"outputs/TSCAL_220105_W_noded_spectra.ps"