[dfrayer@levi DATAdemo]$ pwd /home/scratch/dfrayer/DATAdemo [dfrayer@levi DATAdemo]$ sdfits './AGBT16A_054_02' 'sdfits_16a054_02_scans56_70' -scans=56:70 ./AGBT16A_054_02 Will attempt to fix bad lags in any ACS data (default) Gathering project info ... Scan 70 (VEGAS) integration 51, sampler A2_0, state [SIG, CAL ON ] written. Gathering project info ... "AGBT16A_054_02/ScanLog.fits" does not contain a valid spectrometer scan within the specified scan list. Gathering project info ... "AGBT16A_054_02/ScanLog.fits" does not contain a valid spectral processor scan within the specified scan list. Gathering project info ... "AGBT16A_054_02/ScanLog.fits" does not contain a valid DCR scan within the specified scan list. Gathering project info ... "AGBT16A_054_02/ScanLog.fits" does not contain a valid zpectrometer scan within the specified scan list. [dfrayer@levi DATAdemo]$ gbtidl IDL Version 7.1.1 (linux x86_64 m64). (c) 2009, ITT Visual Information Solutions Installation number: 15269-3. Licensed for use by: National Radio Astronomy Observatory Starting GBTIDL Display Device : X Visual Class : TrueColor Visual Depth : 24-Bit Color Table Size: 256 Number of Colors: 16777216 Decomposed Color: 0 -------------------------------------------------------------------- Welcome to GBTIDL v2.10.1 For news, documentation, enhancement requests, bug tracking, discussion, and contributions, visit: http://gbtidl.nrao.edu For help with a GBTIDL routine from the command line, use the procedure 'usage'. For example: usage,'show' ; gives the syntax of the procedure 'show' usage,'show',/verbose ; gives more information on 'show' -------------------------------------------------------------------- GBTIDL -> filein,'sdfits_example' GBTIDL -> summary Scan Source Vel Proc Seq RestF nIF nInt nFd Az El ------------------------------------------------------------------------------- 56 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 194.8 51.0 57 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 198.8 50.4 58 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 198.9 50.4 59 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 200.9 50.0 60 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 204.7 49.2 61 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 204.9 49.2 62 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 206.7 48.7 63 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 210.3 47.8 64 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 210.5 47.7 65 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 212.2 47.2 66 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 215.6 46.1 67 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 215.7 46.0 68 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 217.4 45.4 69 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 220.5 44.2 70 J085915.7+002329 3387.0 Track 1 1.420 1 51 1 220.7 44.1 ################################## dotri.pro script: pro dotri,sc1,sc2 ;;16A054 HI project ;;does ON-OFF-ON reduction for scans sc1-->sc2 ;;Session 1: sc1=16, sc2=39 sclear for ii=sc1,sc2,3 do begin print,'combining ON-OFF',ii,ii+1 getsigref,ii,ii+1,plnum=0 accum getsigref,ii,ii+1,plnum=1 accum print,'combining ON-OFF',ii+2,ii+1 getsigref,ii+2,ii+1,plnum=0 accum getsigref,ii+2,ii+1,plnum=1 accum endfor ave copy,0,10 return end ######################################3 GBTIDL -> dotri,56,70 % Compiled module: DOTRI. combining ON-OFF 56 57 Blanked spectra: ignored 1 integrations SigScan: 56 RefScan: 57 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.16 Blanked spectra: ignored 1 integrations SigScan: 56 RefScan: 57 (IF:0 FD:0 PL:1) units: Ta (K) Tsys: 16.73 combining ON-OFF 58 57 Blanked spectra: ignored 1 integrations SigScan: 58 RefScan: 57 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.16 Blanked spectra: ignored 1 integrations SigScan: 58 RefScan: 57 (IF:0 FD:0 PL:1) units: Ta (K) Tsys: 16.73 combining ON-OFF 59 60 Blanked spectra: ignored 1 integrations SigScan: 59 RefScan: 60 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.15 Blanked spectra: ignored 1 integrations SigScan: 59 RefScan: 60 (IF:0 FD:0 PL:1) units: Ta (K) Tsys: 16.74 combining ON-OFF 61 60 Blanked spectra: ignored 1 integrations SigScan: 61 RefScan: 60 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.15 Blanked spectra: ignored 1 integrations SigScan: 61 RefScan: 60 (IF:0 FD:0 PL:1) units: Ta (K) Tsys: 16.74 combining ON-OFF 62 63 Blanked spectra: ignored 1 integrations SigScan: 62 RefScan: 63 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.17 Blanked spectra: ignored 1 integrations SigScan: 62 RefScan: 63 (IF:0 FD:0 PL:1) units: Ta (K) Tsys: 16.78 combining ON-OFF 64 63 Blanked spectra: ignored 1 integrations SigScan: 64 RefScan: 63 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.17 Blanked spectra: ignored 1 integrations SigScan: 64 RefScan: 63 (IF:0 FD:0 PL:1) units: Ta (K) Tsys: 16.78 combining ON-OFF 65 66 Blanked spectra: ignored 1 integrations SigScan: 65 RefScan: 66 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.17 Blanked spectra: ignored 1 integrations SigScan: 65 RefScan: 66 (IF:0 FD:0 PL:1) units: Ta (K) Tsys: 16.84 combining ON-OFF 67 66 Blanked spectra: ignored 1 integrations SigScan: 67 RefScan: 66 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.17 Blanked spectra: ignored 1 integrations SigScan: 67 RefScan: 66 (IF:0 FD:0 PL:1) units: Ta (K) Tsys: 16.84 combining ON-OFF 68 69 Blanked spectra: ignored 1 integrations SigScan: 68 RefScan: 69 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.24 Blanked spectra: ignored 1 integrations SigScan: 68 RefScan: 69 (IF:0 FD:0 PL:1) units: Ta (K) Tsys: 16.93 combining ON-OFF 70 69 Blanked spectra: ignored 1 integrations SigScan: 70 RefScan: 69 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.24 Blanked spectra: ignored 1 integrations SigScan: 70 RefScan: 69 (IF:0 FD:0 PL:1) units: Ta (K) Tsys: 16.93 % ACCUMAVE: Average of : 20 spectra GBTIDL -> copy,0,10 GBTIDL -> gsmooth,20,/decimate GBTIDL -> setregion Use the left button to mark the region and the right button to exit GBTIDL -> nfit,5 GBTIDL -> bshape GBTIDL -> baseline GBTIDL -> getsigref,56,57 Blanked spectra: ignored 1 integrations SigScan: 56 RefScan: 57 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.16 GBTIDL -> gsmooth,20 GBTIDL -> getsigref,56,57 Blanked spectra: ignored 1 integrations SigScan: 56 RefScan: 57 (IF:0 FD:0 PL:0) units: Ta (K) Tsys: 16.16 GBTIDL -> gsmooth,20 Shows steps of getsigref: GBTIDL -> gettp,56 Blanked spectra: ignored 1 integrations Scan: 56 (IF:0 FD:0 PL:0) Tsys: 16.12 GBTIDL -> copy,0,2 GBTIDL -> gettp,57 osBlanked spectra: ignored 1 integrations Scan: 57 (IF:0 FD:0 PL:0) Tsys: 16.16 GBTIDL -> copy,0,3 GBTIDL -> show,2 GBTIDL -> oshow,3 GBTIDL -> subtract,2,3 GBTIDL -> freexy GBTIDL -> divide,0,3 GBTIDL -> scale,16.2 GBTIDL -> !g.s[0].units='Ta'