; Load the file.
filein,"data/TGBT21A_501_11.raw.vegas.fits"
; Display a summary of its contents.
summary
; Calibrate the position switched data in scans 152 and 153.
; This will only calibrate the data for the YY polarization,
; and the first spectral window (with the redshifted HI line).
getps,152
; Smooth the data using a Gaussian convolution kernel with a
; standard deviation of 5 channels. Keep only every 5th sample
; (otherwise the channels will be highly correlated, changing 
; the noise properties -- this should be the default).
;gsmooth,5,/decimate
; Change the x axis to velocity.
; By default this will apply a Doppler correction based on the 
; configuration used by the user at the time of observing.
; The data will be converted to the velocity frame specified 
; in the observing configuration as well.
setxunit,"km/s"
; After consulting with their project friend, the user learns
; that the default value of TCAL in the SDFITS is out of date,
; so their data should be scaled by a factor of 1.2 to remove
; a systematic offset.
;scale,1.2
; Save the result as a text file for further analysis.
write_ascii,"outputs/NGC2415_HI.ascii"
fileout,"outputs/NGC2415_HI.fits"
keep
; Now process the remaining spectral windows. No signal is 
; expected, but we'll look at them just in case.
; First, lets freeze the plotter (auto update off).
;freeze
; Change the x axis units back to GHz to avoid crowding.
;setxunit,"GHz"
; Loop over spectral windows doing the same as for the HI line.
;for k=1,4,1 do begin & getps,152,ifnum=k & gsmooth,5,/decimate & scale,1.2 & oshow & endfor
; We do not save these spectra.
; Zoom in to the OH line (this can also be done interactively 
; using the plotter).
;setxy,1.64,1.66,-0.5,1
; Nothing obvious in there.