filein,"data/TGBT21A_501_11.raw.vegas.fits"
fileout,"outputs/TGBT21A_501_11_ifnum_0.fits"
;; Now get the data for spectral window 0 scan 152 integration 1 to 3.
;for i=0,2 do begin $
;    chunk = getchunk(scan=152,ifnum=0,int=i,count=nchunk) &
;    ; Save it to the output file.
;    putchunk,chunk & $
;    ;data_free, chunk & $
;    ; Same for scan 153.
;    ;chunk = getchunk(scan=153,ifnum=0,int=i,count=nchunk) &
;    ;putchunk,chunk & $
;end
; Bloody putchunk does not seem to work properly in a for loop!
chunk = getchunk(scan=152,ifnum=0,int=0,count=nchunk)
putchunk,chunk
chunk = getchunk(scan=152,ifnum=0,int=1,count=nchunk)
putchunk,chunk
chunk = getchunk(scan=152,ifnum=0,int=2,count=nchunk)
putchunk,chunk
chunk = getchunk(scan=153,ifnum=0,int=0,count=nchunk)
putchunk,chunk
chunk = getchunk(scan=153,ifnum=0,int=1,count=nchunk)
putchunk,chunk
chunk = getchunk(scan=153,ifnum=0,int=2,count=nchunk)
putchunk,chunk


exit