Rcvr Ch  X  Y
ODM	 2  3
CM	 1  5
SG	 1  5
Sampler	 0  1
IFCF	 3000

Scan
  1	test scan 1 minute

  2	800 MHz BW, 4300 MHz CF, Balanced before scan
2002_12_30_22:08:46A.fits scan: 2 Rx4_6_baselines 2 spectra 2 phases 10 recs

  3	NGC7027, xonoff, on, CF=4300
  4	NGC7027, xonoff, off, CF=4300
2002_12_30_22:21:15A.fits scan: 3 NGC7027 2 spectra 2 phases 10 recs
   az: -70.24 el: 61.19 ra2000: 21.12 dec2000: 42.24 freq: 4300
2002_12_30_22:26:34A.fits scan: 4 NGC7027 2 spectra 2 phases 10 recs
   az: -70.25 el: 61.23 ra2000: 21.21 dec2000: 42.23 freq: 4300
show_avg_scan_diff(3, 4, spec_num=1, decimate=1, linecolor=2, newplot=T)

  5	false start
2002_12_30_22:32:49A.fits scan: 5 NGC7027 2 spectra 2 phases 1 recs
   az: -70.01 el: 60.08 ra2000: 21.21 dec2000: 42.24 freq: 4700

  6	NGC7027, xonoff, on, CF=4700
  7	NGC7027, xonoff, off, CF=4700
2002_12_30_22:33:47A.fits scan: 6 NGC7027 2 spectra 2 phases 10 recs
   az: -69.73 el: 58.89 ra2000: 21.12 dec2000: 42.24 freq: 4700
2002_12_30_22:39:06A.fits scan: 7 NGC7027 2 spectra 2 phases 10 recs
   az: -69.76 el: 58.92 ra2000: 21.21 dec2000: 42.23 freq: 4700

  8	NGC7027, xonoff, on, CF=5100
  9	NGC7027, xonoff, off, CF=5100
2002_12_30_22:45:41A.fits scan: 8 NGC7027 2 spectra 2 phases 10 recs
   az: -69.18 el: 56.7 ra2000: 21.12 dec2000: 42.24 freq: 5100
2002_12_30_22:51:00A.fits scan: 9 NGC7027 2 spectra 2 phases 10 recs
   az: -69.2 el: 56.73 ra2000: 21.21 dec2000: 42.23 freq: 5100
show_record_dsmror(9,  phase=1, rec_num=[1,2], spec_num=1, decimate=1, linecolor=2, newplot=T )
pg.setyscale(-0.04, 0.02)
for (i in 3:10) show_record_dsmror(9,  phase=1, rec_num=[1,i], spec_num=1, decimate=1, linecolor=1+i/2, newplot=F)

 10	NGC7027, xonoff, on, CF=5500
 11	NGC7027, xonoff, off, CF=5500
2002_12_30_22:57:23A.fits scan: 10 NGC7027 2 spectra 2 phases 10 recs
   az: -68.56 el: 54.56 ra2000: 21.12 dec2000: 42.24 freq: 5500
2002_12_30_23:02:42A.fits scan: 11 NGC7027 2 spectra 2 phases 10 recs
   az: -68.59 el: 54.6 ra2000: 21.21 dec2000: 42.23 freq: 5500
show_avg_scan_diff(3, 4, spec_num=1, decimate=1, linecolor=2, newplot=T)
pg.setyscale(0, 0.6)
pg.setxscale(3900, 5900)
show_avg_scan_diff(6, 7, spec_num=1, decimate=1, linecolor=2, newplot=F)
show_avg_scan_diff(8, 9, spec_num=1, decimate=1, linecolor=2, newplot=F)
show_avg_scan_diff(10, 11, spec_num=1, decimate=1, linecolor=2, newplot=F)
# compositeX.ps
# channel 2 unuseable due to spectrometer defect

Change to 200 MHz BW, CF=5000
 12	test scan
2002_12_30_23:15:17A.fits scan: 12 NGC7027 2 spectra 2 phases 1 recs

 13	cband021230.g
# sc.setParameter('projectId', 'TBASEJRF021230')
# cband_200mhz();
# rx4_6.setParameter('yrIfFilterSwitch', 2);
# rx4_6.setParameter('xlIfFilterSwitch', 2);
# rx4_6.setParameter('tuning_frequency', 5000.0);
# rx4_6.prepare();
focusonoff('HighElevEighth', '23:00:00', '+30:00:00', 302,
		        X=[0.0,0.0], Y=[0.3,0.0], Z=[0.0,0.0],
		        coord_mode='J2000', overhead=30.0)
2002_12_30_23:27:34A.fits scan: 13 HighElevEighth 2 spectra 2 phases 10 recs
   az: 256.6 el: 66.89 ra2000: 23 dec2000: 30 freq: 5000
2002_12_30_23:32:56A.fits scan: 14 HighElevEighth 2 spectra 2 phases 10 recs
   az: 256.6 el: 66.91 ra2000: 23.09 dec2000: 30 freq: 5000

focusonoff('HighElevQuarter', '23:00:00', '+30:00:00', 302,
		        X=[0.0,0.0], Y=[0.6,0.0], Z=[0.0,0.0],
		        coord_mode='J2000', overhead=30.0)
2002_12_30_23:38:17A.fits scan: 15 HighElevQuarter 2 spectra 2 phases 10 recs
   az: 259.3 el: 64.83 ra2000: 23 dec2000: 30 freq: 5000
2002_12_30_23:43:38A.fits scan: 16 HighElevQuarter 2 spectra 2 phases 10 recs
   az: 259.2 el: 64.86 ra2000: 23.09 dec2000: 30 freq: 5000
show_avg_scan_diff(13, 14, spec_num=1, decimate=16, offset=0.002, linecolor=2, newplot=T, usb=F)
x8Hsdiff := sdiff
pg.setyscale(-0.01, 0.01) 
show_avg_scan_diff(13, 14, spec_num=2, decimate=16, offset=0.004, linecolor=3, newplot=F, usb=F)
y8Hsdiff := sdiff
show_avg_scan_diff(15, 16, spec_num=1, decimate=16, offset=-0.006, linecolor=2, newplot=F, usb=F)
x4Hsdiff := sdiff
show_avg_scan_diff(15, 16, spec_num=2, decimate=16, offset=-0.006, linecolor=3, newplot=F, usb=F)
y4Hsdiff := sdiff
# high8th_qrt.ps

focusonoff('LowElevEighth', '23:30:00', '-42:00:00', 302,
		        X=[0.0,0.0], Y=[0.3,0.0], Z=[0.0,0.0],
		        coord_mode='J2000', overhead=30.0)
2002_12_30_23:52:08A.fits scan: 17 LowElevEighth 2 spectra 2 phases 10 recs
   az: 198.3 el: 6.625 ra2000: 23.5 dec2000: -42 freq: 5000
2002_12_30_23:57:27A.fits scan: 18 LowElevEighth 2 spectra 2 phases 10 recs
   az: 198.3 el: 6.641 ra2000: 23.59 dec2000: -41.99 freq: 5000

focusonoff('LowElevQuarter', '23:30:00', '-42:00:00', 302,
		        X=[0.0,0.0], Y=[0.6,0.0], Z=[0.0,0.0],
		        coord_mode='J2000', overhead=30.0)
2002_12_31_00:02:46A.fits scan: 19 LowElevQuarter 2 spectra 2 phases 10 recs
   az: 200.2 el: 5.939 ra2000: 23.5 dec2000: -42.01 freq: 5000
2002_12_31_00:08:05A.fits scan: 20 LowElevQuarter 2 spectra 2 phases 10 recs
   az: 200.2 el: 5.964 ra2000: 23.59 dec2000: -41.99 freq: 5000

show_avg_scan_diff(17, 18, spec_num=1, decimate=16, offset=-0.002, linecolor=2, newplot=T, usb=F)
x8Lsdiff := sdiff
pg.setyscale(-0.01, 0.01) 
show_avg_scan_diff(17, 18, spec_num=2, decimate=16, offset=0.002, linecolor=3, newplot=F, usb=F)
y8Lsdiff := sdiff
show_avg_scan_diff(19, 20, spec_num=1, decimate=16, offset=-0.010, linecolor=2, newplot=F, usb=F)
x4Lsdiff := sdiff
show_avg_scan_diff(19, 20, spec_num=2, decimate=16, offset=-0.006, linecolor=3, newplot=F, usb=F)
y4Lsdiff := sdiff

focusonoff('MedElevEighth', '23:50:00', '-37:00:00', 302,
		        X=[0.0,0.0], Y=[0.3,0.0], Z=[0.0,0.0],
		        coord_mode='J2000', overhead=30.0)
2002_12_31_00:13:42A.fits scan: 21 MedElevEighth 2 spectra 2 phases 10 recs
   az: 200.3 el: 11.18 ra2000: 23.83 dec2000: -37 freq: 5000
2002_12_31_00:19:02A.fits scan: 22 MedElevEighth 2 spectra 2 phases 10 recs
   az: 200.3 el: 11.2 ra2000: 23.93 dec2000: -36.99 freq: 5000

focusonoff('MedElevQuarter', '23:50:00', '-37:00:00', 302,
		        X=[0.0,0.0], Y=[0.6,0.0], Z=[0.0,0.0],
		        coord_mode='J2000', overhead=30.0)
2002_12_31_00:24:22A.fits scan: 23 MedElevQuarter 2 spectra 2 phases 10 recs
   az: 202.4 el: 10.42 ra2000: 23.83 dec2000: -37.01 freq: 5000
2002_12_31_00:29:42A.fits scan: 24 MedElevQuarter 2 spectra 2 phases 10 recs
   az: 202.3 el: 10.44 ra2000: 23.93 dec2000: -36.99 freq: 5000

show_avg_scan_diff(21, 22, spec_num=1, decimate=16, offset=-0.001, linecolor=2, newplot=T, usb=F)
x8Msdiff := sdiff
pg.setyscale(-0.01, 0.01) 
show_avg_scan_diff(21, 22, spec_num=2, decimate=16, offset=0.003, linecolor=3, newplot=F, usb=F)
y8Msdiff := sdiff
show_avg_scan_diff(23, 24, spec_num=1, decimate=16, offset=-0.006, linecolor=2, newplot=F, usb=F)
x4Msdiff := sdiff
show_avg_scan_diff(23, 24, spec_num=2, decimate=16, offset=-0.003, linecolor=3, newplot=F, usb=F)
y4Msdiff := sdiff
# med8th_qrt.ps

spectrum_fft((x4Hsdiff + x4Msdiff + x4Lsdiff) / 3, 200);
pg.setyscale(-0.01, 0.25) 
pg.setxscale(-1, 140) 
spectrum_fft((y4Hsdiff + y4Msdiff + y4Lsdiff) / 3, 200, linecolor=3, newplot=F);
# avgspec4ab.ps

spectrum_fft((x8Hsdiff + x8Msdiff + x8Lsdiff) / 3, 200);
pg.setyscale(-0.01, 0.25) 
pg.setxscale(-1, 140) 
spectrum_fft((y8Hsdiff + y8Msdiff + y8Lsdiff) / 3, 200, linecolor=3, newplot=F);
# avgspec8ab.ps

spectrum_fft(x4Lsdiff, 200);
pg.setyscale(-0.01, 0.5) 
pg.setxscale(-1, 140) 
spectrum_fft(x4Msdiff, 200, linecolor=4, newplot=F);
spectrum_fft(x4Hsdiff, 200, linecolor=3, newplot=F);
# allspec4a.ps
pg.setxscale(-1, 25)
# allspec4ahx.ps

spectrum_fft(y4Lsdiff, 200);
pg.setyscale(-0.01, 0.5) 
pg.setxscale(-1, 140) 
spectrum_fft(y4Msdiff, 200, linecolor=4, newplot=F);
spectrum_fft(y4Hsdiff, 200, linecolor=3, newplot=F);
# allspec4b.ps
pg.setxscale(-1, 25)
# allspec4bhx.ps

spectrum_fft(x8Lsdiff, 200);
pg.setyscale(-0.01, 0.5) 
pg.setxscale(-1, 140) 
spectrum_fft(x8Msdiff, 200, linecolor=4, newplot=F);
spectrum_fft(x8Hsdiff, 200, linecolor=3, newplot=F);
# allspec8a.ps
pg.setxscale(-1, 25)
# allspec8ahx.ps

spectrum_fft(y8Lsdiff, 200);
pg.setyscale(-0.01, 0.5) 
pg.setxscale(-1, 140) 
spectrum_fft(y8Msdiff, 200, linecolor=4, newplot=F);
spectrum_fft(y8Hsdiff, 200, linecolor=3, newplot=F);
# allspec8b.ps
pg.setxscale(-1, 25)
# allspec8bhx.ps

spectrum_fft(y8Lsdiff, 200);
pg.clear(); pg.plotxy(ind(spec_fft), spec_fft)
pg.setxscale(80, 180)
pg.setyscale(-0.01, 0.2)
maskoff := [80:110,160:180]
maskon := [111:159]
off := sum(spec_fft[maskoff]) / length(maskoff)
print sum(spec_fft[maskon] - off)
# 0.40818545

spectrum_fft(y8Msdiff, 200);
pg.clear(); pg.plotxy(ind(spec_fft), spec_fft)
pg.setxscale(80, 180)
pg.setyscale(-0.01, 0.2)
maskoff := [80:110,160:180]
maskon := [111:159]
off := sum(spec_fft[maskoff]) / length(maskoff)
print sum(spec_fft[maskon] - off)
# 0.344785286

spectrum_fft(y8Hsdiff, 200);
pg.clear(); pg.plotxy(ind(spec_fft), spec_fft)
pg.setxscale(80, 180)
pg.setyscale(-0.01, 0.2)
maskoff := [80:110,160:180]
maskon := [111:159]
off := sum(spec_fft[maskoff]) / length(maskoff)
print sum(spec_fft[maskon] - off)
# 1.03848205

# total power change: X H: 1.00, M: 1.36, L: 1.73
# total power change: Y H: 1.00, M: 1.32, L: 1.65

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_00:37:52A.fits scan: 25 NGC7027 2 spectra 2 phases 10 recs
   az: 298.7 el: 36.68 ra2000: 21.12 dec2000: 42.23 freq: 5000
2002_12_31_00:43:11A.fits scan: 26 NGC7027 2 spectra 2 phases 10 recs
   az: 298.6 el: 36.71 ra2000: 21.21 dec2000: 42.23 freq: 5000
show_avg_scan_diff(25, 26, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.35, 0.45)
x25sdiff := sdiff;
x25rdiff := rdiff;
show_avg_scan_diff(25, 26, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F, usb=F)
y25sdiff := sdiff;
y25rdiff := rdiff;
# Tsrc_o_Tsys7027.ps

xyz := xyz_save := ant.getValue('xyzSubrSegments');
xyz[1].length.seconds := 302;
xyz[1].Y.pos +:= 0.3;
ant.setParameter('xyzSubrSegments', xyz);
xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_00:48:30A.fits scan: 27 NGC7027 2 spectra 2 phases 10 recs
   az: 299.6 el: 34.86 ra2000: 21.12 dec2000: 42.25 freq: 5000
2002_12_31_00:53:49A.fits scan: 28 NGC7027 2 spectra 2 phases 10 recs
   az: 299.6 el: 34.89 ra2000: 21.21 dec2000: 42.22 freq: 5000
show_avg_scan_diff(27, 28, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.35, 0.45)
x27sdiff := sdiff;
x27rdiff := rdiff;
show_avg_scan_diff(27, 28, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y27sdiff := sdiff;
y27rdiff := rdiff;

xyz[1].Y.pos +:= 0.3;
ant.setParameter('xyzSubrSegments', xyz);
xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_00:59:08A.fits scan: 29 NGC7027 2 spectra 2 phases 10 recs
   az: 300.6 el: 33.06 ra2000: 21.12 dec2000: 42.25 freq: 5000
2002_12_31_01:04:27A.fits scan: 30 NGC7027 2 spectra 2 phases 10 recs
   az: 300.5 el: 33.08 ra2000: 21.21 dec2000: 42.22 freq: 5000
show_avg_scan_diff(29, 30, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.35, 0.45)
x29sdiff := sdiff;
x29rdiff := rdiff;
show_avg_scan_diff(29, 30, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y29sdiff := sdiff;
y29rdiff := rdiff;

ant.setParameter('xyzSubrSegments', xyz_save);
xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_01:09:46A.fits scan: 31 NGC7027 2 spectra 2 phases 10 recs
   az: 301.5 el: 31.27 ra2000: 21.12 dec2000: 42.25 freq: 5000
2002_12_31_01:15:05A.fits scan: 32 NGC7027 2 spectra 2 phases 10 recs
   az: 301.5 el: 31.3 ra2000: 21.21 dec2000: 42.22 freq: 5000
show_avg_scan_diff(31, 32, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.35, 0.45)
x31sdiff := sdiff;
x31rdiff := rdiff;
show_avg_scan_diff(31, 32, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y31sdiff := sdiff;
y31rdiff := rdiff;

xonoff('3C48', '01:37:41.30', '+33:09:35.4', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_01:22:57A.fits scan: 33 3C48 2 spectra 2 phases 10 recs
   az: 252.3 el: 76.25 ra2000: 1.628 dec2000: 33.16 freq: 5000
2002_12_31_01:28:22A.fits scan: 34 3C48 2 spectra 2 phases 10 recs
   az: 252.3 el: 76.27 ra2000: 1.72 dec2000: 33.16 freq: 5000
show_avg_scan_diff(33, 34, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.35, 0.45)
x33sdiff := sdiff;
x33rdiff := rdiff;
show_avg_scan_diff(33, 34, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y33sdiff := sdiff;
y33rdiff := rdiff;

pg.clear(); pg.plotxy(freq, x33sdiff / x31sdiff, linecolor=2, newplot=T)
pg.setyscale(0.80, 1.0)
pg.plotxy(freq, y33sdiff / y31sdiff, linecolor=3, newplot=F)
pg.plotxy(freq, x33rdiff / x31rdiff, linecolor=2, newplot=F)
pg.plotxy(freq, y33rdiff / y31rdiff, linecolor=3, newplot=F)
pg.plotxy(freq, x33rdiff / x35rdiff - 0.05, linecolor=2, newplot=F)
pg.plotxy(freq, y33rdiff / y35rdiff - 0.06, linecolor=3, newplot=F)
# 3C48rat.ps
spectrum_fft(x33sdiff / x31sdiff - [1:2049] / 5e4, 200)
pg.setxscale(-2, 120)
spectrum_fft(y33sdiff / y31sdiff - [1:2049] / 5e4, 200, linecolor=3, newplot=F)
pg.setyscale(-0.05, 1.7) 
# fft3C48srat.ps
spectrum_fft(x33rdiff / x31rdiff - [1:2049] / 5e4, 200)
pg.setxscale(-2, 120)
spectrum_fft(y33rdiff / y31rdiff - [1:2049] / 5e4, 200, linecolor=3, newplot=F)
pg.setyscale(-0.05, 1.7) 
# fft3C48rrat.ps

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_01:36:22A.fits scan: 35 NGC7027 2 spectra 2 phases 10 recs
   az: 304.1 el: 26.88 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_01:41:40A.fits scan: 36 NGC7027 2 spectra 2 phases 10 recs
   az: 304.1 el: 26.91 ra2000: 21.21 dec2000: 42.22 freq: 5000
show_avg_scan_diff(35, 36, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.35, 0.45)
x35sdiff := sdiff;
x35rdiff := rdiff;
show_avg_scan_diff(35, 36, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y35sdiff := sdiff;
y35rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_01:46:59A.fits scan: 37 NGC7027 2 spectra 2 phases 10 recs
   az: 305.2 el: 25.17 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_01:52:17A.fits scan: 38 NGC7027 2 spectra 2 phases 10 recs
   az: 305.2 el: 25.19 ra2000: 21.21 dec2000: 42.22 freq: 5000
show_avg_scan_diff(37, 38, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.30, 0.45)
x37sdiff := sdiff;
x37rdiff := rdiff;
show_avg_scan_diff(37, 38, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y37sdiff := sdiff;
y37rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_01:57:36A.fits scan: 39 NGC7027 2 spectra 2 phases 10 recs
   az: 306.3 el: 23.47 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_02:02:54A.fits scan: 40 NGC7027 2 spectra 2 phases 10 recs
   az: 306.3 el: 23.5 ra2000: 21.21 dec2000: 42.22 freq: 5000
show_avg_scan_diff(39, 40, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.30, 0.45)
x39sdiff := sdiff;
x39rdiff := rdiff;
show_avg_scan_diff(39, 40, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y39sdiff := sdiff;
y39rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_02:08:13A.fits scan: 41 NGC7027 2 spectra 2 phases 10 recs
   az: 307.4 el: 21.81 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_02:13:31A.fits scan: 42 NGC7027 2 spectra 2 phases 10 recs
   az: 307.4 el: 21.84 ra2000: 21.21 dec2000: 42.22 freq: 5000
show_avg_scan_diff(41, 42, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.30, 0.45)
x41sdiff := sdiff;
x41rdiff := rdiff;
show_avg_scan_diff(41, 42, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y41sdiff := sdiff;
y41rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_02:18:50A.fits scan: 43 NGC7027 2 spectra 2 phases 10 recs
   az: 308.6 el: 20.17 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_02:24:08A.fits scan: 44 NGC7027 2 spectra 2 phases 10 recs
   az: 308.6 el: 20.19 ra2000: 21.21 dec2000: 42.22 freq: 5000
show_avg_scan_diff(43, 44, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.30, 0.45)
x43sdiff := sdiff;
x43rdiff := rdiff;
show_avg_scan_diff(43, 44, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y43sdiff := sdiff;
y43rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_02:29:27A.fits scan: 45 NGC7027 2 spectra 2 phases 10 recs
   az: 309.8 el: 18.55 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_02:34:45A.fits scan: 46 NGC7027 2 spectra 2 phases 10 recs
   az: 309.8 el: 18.58 ra2000: 21.21 dec2000: 42.22 freq: 5000
show_avg_scan_diff(45, 46, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.30, 0.45)
x45sdiff := sdiff;
x45rdiff := rdiff;
show_avg_scan_diff(45, 46, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y45sdiff := sdiff;
y45rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_02:40:04A.fits scan: 47 NGC7027 2 spectra 2 phases 10 recs
   az: 311 el: 16.97 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_02:45:22A.fits scan: 48 NGC7027 2 spectra 2 phases 10 recs
   az: 311 el: 16.99 ra2000: 21.21 dec2000: 42.22 freq: 5000
show_avg_scan_diff(47, 48, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.30, 0.45)
x47sdiff := sdiff;
x47rdiff := rdiff;
show_avg_scan_diff(47, 48, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y47sdiff := sdiff;
y47rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_02:50:40A.fits scan: 49 NGC7027 2 spectra 2 phases 10 recs
   az: 312.3 el: 15.42 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_02:55:58A.fits scan: 50 NGC7027 2 spectra 2 phases 10 recs
   az: 312.2 el: 15.44 ra2000: 21.21 dec2000: 42.23 freq: 5000
show_avg_scan_diff(49, 50, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.30, 0.45)
x49sdiff := sdiff;
x49rdiff := rdiff;
show_avg_scan_diff(49, 50, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y49sdiff := sdiff;
y49rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_03:01:16A.fits scan: 51 NGC7027 2 spectra 2 phases 10 recs
   az: 313.6 el: 13.89 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_03:06:33A.fits scan: 52 NGC7027 2 spectra 2 phases 10 recs
   az: 313.5 el: 13.92 ra2000: 21.21 dec2000: 42.23 freq: 5000
show_avg_scan_diff(51, 52, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.30, 0.45)
x51sdiff := sdiff;
x51rdiff := rdiff;
show_avg_scan_diff(51, 52, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y51sdiff := sdiff;
y51rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_03:11:51A.fits scan: 53 NGC7027 2 spectra 2 phases 10 recs
   az: 314.9 el: 12.41 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_03:17:08A.fits scan: 54 NGC7027 2 spectra 2 phases 10 recs
   az: 314.8 el: 12.44 ra2000: 21.21 dec2000: 42.23 freq: 5000
show_avg_scan_diff(53, 54, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.30, 0.45)
x53sdiff := sdiff;
x53rdiff := rdiff;
show_avg_scan_diff(53, 54, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y53sdiff := sdiff;
y53rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_03:22:26A.fits scan: 55 NGC7027 2 spectra 2 phases 10 recs
   az: 316.2 el: 10.96 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_03:27:43A.fits scan: 56 NGC7027 2 spectra 2 phases 10 recs
   az: 316.2 el: 10.99 ra2000: 21.21 dec2000: 42.23 freq: 5000
show_avg_scan_diff(55, 56, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.30, 0.45)
x55sdiff := sdiff;
x55rdiff := rdiff;
show_avg_scan_diff(55, 56, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y55sdiff := sdiff;
y55rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_03:33:01A.fits scan: 57 NGC7027 2 spectra 2 phases 10 recs
   az: 317.6 el: 9.551 ra2000: 21.12 dec2000: 42.23 freq: 5000
2002_12_31_03:38:19A.fits scan: 58 NGC7027 2 spectra 2 phases 10 recs
   az: 317.6 el: 9.574 ra2000: 21.21 dec2000: 42.23 freq: 5000
show_avg_scan_diff(57, 58, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.20, 0.45)
x57sdiff := sdiff;
x57rdiff := rdiff;
show_avg_scan_diff(57, 58, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y57sdiff := sdiff;
y57rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_03:43:37A.fits scan: 59 NGC7027 2 spectra 2 phases 10 recs
   az: 319 el: 8.18 ra2000: 21.12 dec2000: 42.24 freq: 5000
2002_12_31_03:48:54A.fits scan: 60 NGC7027 2 spectra 2 phases 10 recs
   az: 319 el: 8.204 ra2000: 21.21 dec2000: 42.23 freq: 5000
show_avg_scan_diff(59, 60, spec_num=1, decimate=16, offset=0, linecolor=2, newplot=T, usb=F)
pg.setyscale(0.20, 0.45)
x59sdiff := sdiff;
x59rdiff := rdiff;
show_avg_scan_diff(59, 60, spec_num=2, decimate=16, offset=0, linecolor=3, newplot=F, usb=F)
y59sdiff := sdiff;
y59rdiff := rdiff;

xonoff('NGC7027', '21:07:01.60', '+42:14:10.0', 302, coord_mode='J2000',
		        overhead=30.0)
2002_12_31_03:54:12A.fits scan: 61 NGC7027 2 spectra 2 phases 10 recs
   az: 320.5 el: 6.854 ra2000: 21.12 dec2000: 42.24 freq: 5000

pg.clear(); pg.plotxy(freq, x27sdiff / x25sdiff, linecolor=2, newplot=T)
pg.setyscale(0.81, 1.02)
pg.plotxy(freq, x29sdiff / x25sdiff + 0.03, linecolor=2, newplot=F)
pg.plotxy(freq, x31sdiff / x25sdiff + 0.00, linecolor=2, newplot=F)
pg.plotxy(freq, x35sdiff / x25sdiff + 0.03, linecolor=2, newplot=F)
pg.plotxy(freq, x37sdiff / x25sdiff + 0.033, linecolor=2, newplot=F)
pg.plotxy(freq, x39sdiff / x25sdiff + 0.045, linecolor=2, newplot=F)
pg.plotxy(freq, x41sdiff / x25sdiff + 0.042, linecolor=2, newplot=F)
pg.plotxy(freq, x43sdiff / x25sdiff + 0.025, linecolor=2, newplot=F)
pg.plotxy(freq, x45sdiff / x25sdiff + 0.053, linecolor=2, newplot=F)
pg.plotxy(freq, x47sdiff / x25sdiff + 0.067, linecolor=2, newplot=F)
pg.plotxy(freq, x49sdiff / x25sdiff + 0.093, linecolor=2, newplot=F)
pg.plotxy(freq, x51sdiff / x25sdiff + 0.105, linecolor=2, newplot=F)
pg.plotxy(freq, x53sdiff / x25sdiff + 0.131, linecolor=2, newplot=F)
pg.plotxy(freq, x55sdiff / x25sdiff + 0.085, linecolor=2, newplot=F)
pg.plotxy(freq, x57sdiff / x25sdiff + 0.207, linecolor=2, newplot=F)
pg.plotxy(freq, x59sdiff / x25sdiff + 0.22, linecolor=2, newplot=F)
pg.plotxy(freq, x27sdiff / x25sdiff - 0.165, linecolor=3, newplot=F)
# xsdiffrat.ps

pg.clear(); pg.plotxy(freq, y27sdiff / y25sdiff, linecolor=2, newplot=T)
pg.setyscale(0.81, 1.02)
pg.plotxy(freq, y29sdiff / y25sdiff + 0.022, linecolor=2, newplot=F)
pg.plotxy(freq, y31sdiff / y25sdiff - 0.01, linecolor=2, newplot=F)
pg.plotxy(freq, y35sdiff / y25sdiff + 0.02, linecolor=2, newplot=F)
pg.plotxy(freq, y37sdiff / y25sdiff + 0.027, linecolor=2, newplot=F)
pg.plotxy(freq, y39sdiff / y25sdiff + 0.032, linecolor=2, newplot=F)
pg.plotxy(freq, y41sdiff / y25sdiff + 0.028, linecolor=2, newplot=F)
pg.plotxy(freq, y43sdiff / y25sdiff + 0.006, linecolor=2, newplot=F)
pg.plotxy(freq, y45sdiff / y25sdiff + 0.038, linecolor=2, newplot=F)
pg.plotxy(freq, y47sdiff / y25sdiff + 0.047, linecolor=2, newplot=F)
pg.plotxy(freq, y49sdiff / y25sdiff + 0.069, linecolor=2, newplot=F)
pg.plotxy(freq, y51sdiff / y25sdiff + 0.08, linecolor=2, newplot=F)
pg.plotxy(freq, y53sdiff / y25sdiff + 0.103, linecolor=2, newplot=F)
pg.plotxy(freq, y55sdiff / y25sdiff + 0.06, linecolor=2, newplot=F)
pg.plotxy(freq, y57sdiff / y25sdiff + 0.185, linecolor=2, newplot=F)
pg.plotxy(freq, y59sdiff / y25sdiff + 0.193, linecolor=2, newplot=F)
pg.plotxy(freq, y27sdiff / y25sdiff - 0.155, linecolor=3, newplot=F)
# ysdiffrat.ps

elev := [34.86, 33.06, 31.27, 26.88, 25.17, 23.47, 21.81,
	 20.17, 18.55, 16.97, 15.42, 13.89, 12.41, 10.96, 9.55, 8.18];
avgxs := avgys := avgxr := avgyr := array(0.0, length(elev))
mask := (length(freq) / 10):(9 *  length(freq) / 10)
avgxs[1] := sum((x27sdiff / x25sdiff)[mask]) / length(mask);
avgxs[2] := sum((x29sdiff / x25sdiff)[mask]) / length(mask);
avgxs[3] := sum((x31sdiff / x25sdiff)[mask]) / length(mask);
avgxs[4] := sum((x35sdiff / x25sdiff)[mask]) / length(mask);
avgxs[5] := sum((x37sdiff / x25sdiff)[mask]) / length(mask);
avgxs[6] := sum((x39sdiff / x25sdiff)[mask]) / length(mask);
avgxs[7] := sum((x41sdiff / x25sdiff)[mask]) / length(mask);
avgxs[8] := sum((x43sdiff / x25sdiff)[mask]) / length(mask);
avgxs[9] := sum((x45sdiff / x25sdiff)[mask]) / length(mask);
avgxs[10] := sum((x47sdiff / x25sdiff)[mask]) / length(mask);
avgxs[11] := sum((x49sdiff / x25sdiff)[mask]) / length(mask);
avgxs[12] := sum((x51sdiff / x25sdiff)[mask]) / length(mask);
avgxs[13] := sum((x53sdiff / x25sdiff)[mask]) / length(mask);
avgxs[14] := sum((x55sdiff / x25sdiff)[mask]) / length(mask);
avgxs[15] := sum((x57sdiff / x25sdiff)[mask]) / length(mask);
avgxs[16] := sum((x59sdiff / x25sdiff)[mask]) / length(mask);
pg.clear(); pg.plotxy((1/sin(pi*elev/180)) / (1/sin(pi*36.68/180)), avgxs,
                      xtitle='Secant(z) Ratio',
                      ytitle='Source Power Ratio', linecolor=2, newplot=T)
avgys[1] := sum((y27sdiff / y25sdiff)[mask]) / length(mask);
avgys[2] := sum((y29sdiff / y25sdiff)[mask]) / length(mask);
avgys[3] := sum((y31sdiff / y25sdiff)[mask]) / length(mask);
avgys[4] := sum((y35sdiff / y25sdiff)[mask]) / length(mask);
avgys[5] := sum((y37sdiff / y25sdiff)[mask]) / length(mask);
avgys[6] := sum((y39sdiff / y25sdiff)[mask]) / length(mask);
avgys[7] := sum((y41sdiff / y25sdiff)[mask]) / length(mask);
avgys[8] := sum((y43sdiff / y25sdiff)[mask]) / length(mask);
avgys[9] := sum((y45sdiff / y25sdiff)[mask]) / length(mask);
avgys[10] := sum((y47sdiff / y25sdiff)[mask]) / length(mask);
avgys[11] := sum((y49sdiff / y25sdiff)[mask]) / length(mask);
avgys[12] := sum((y51sdiff / y25sdiff)[mask]) / length(mask);
avgys[13] := sum((y53sdiff / y25sdiff)[mask]) / length(mask);
avgys[14] := sum((y55sdiff / y25sdiff)[mask]) / length(mask);
avgys[15] := sum((y57sdiff / y25sdiff)[mask]) / length(mask);
avgys[16] := sum((y59sdiff / y25sdiff)[mask]) / length(mask);
pg.plotxy((1/sin(pi*elev/180)) / (1/sin(pi*36.68/180)), avgys,
          xtitle='Secant(z) Ratio',
          ytitle='Source Power Ratio', linecolor=3, newplot=F)
# avgsdiffrat.ps

pg.clear(); pg.plotxy(freq, x27rdiff / x25rdiff, linecolor=2, newplot=T)
pg.setyscale(0.81, 1.02)
pg.plotxy(freq, x29rdiff / x25rdiff + 0.02, linecolor=2, newplot=F)
pg.plotxy(freq, x31rdiff / x25rdiff - 0.016, linecolor=2, newplot=F)
pg.plotxy(freq, x35rdiff / x25rdiff - 0.015, linecolor=2, newplot=F)
pg.plotxy(freq, x37rdiff / x25rdiff - 0.025, linecolor=2, newplot=F)
pg.plotxy(freq, x39rdiff / x25rdiff - 0.023, linecolor=2, newplot=F)
pg.plotxy(freq, x41rdiff / x25rdiff - 0.043, linecolor=2, newplot=F)
pg.plotxy(freq, x43rdiff / x25rdiff - 0.075, linecolor=2, newplot=F)
pg.plotxy(freq, x45rdiff / x25rdiff - 0.07, linecolor=2, newplot=F)
pg.plotxy(freq, x47rdiff / x25rdiff - 0.072, linecolor=2, newplot=F)
pg.plotxy(freq, x49rdiff / x25rdiff - 0.062, linecolor=2, newplot=F)
pg.plotxy(freq, x51rdiff / x25rdiff - 0.0687, linecolor=2, newplot=F)
pg.plotxy(freq, x53rdiff / x25rdiff - 0.068, linecolor=2, newplot=F)
pg.plotxy(freq, x55rdiff / x25rdiff - 0.16, linecolor=2, newplot=F)
pg.plotxy(freq, x57rdiff / x25rdiff - 0.04, linecolor=2, newplot=F)
pg.plotxy(freq, x59rdiff / x25rdiff - 0.083, linecolor=2, newplot=F)
pg.plotxy(freq, x27rdiff / x25rdiff - 0.165, linecolor=3, newplot=F)
# xrdiffrat.ps

pg.clear(); pg.plotxy(freq, y27rdiff / y25rdiff, linecolor=2, newplot=T)
pg.setyscale(0.81, 1.02)
pg.plotxy(freq, y29rdiff / y25rdiff + 0.018, linecolor=2, newplot=F)
pg.plotxy(freq, y31rdiff / y25rdiff - 0.018, linecolor=2, newplot=F)
pg.plotxy(freq, y35rdiff / y25rdiff - 0.010, linecolor=2, newplot=F)
pg.plotxy(freq, y37rdiff / y25rdiff - 0.014, linecolor=2, newplot=F)
pg.plotxy(freq, y39rdiff / y25rdiff - 0.021, linecolor=2, newplot=F)
pg.plotxy(freq, y41rdiff / y25rdiff - 0.041, linecolor=2, newplot=F)
pg.plotxy(freq, y43rdiff / y25rdiff - 0.077, linecolor=2, newplot=F)
pg.plotxy(freq, y45rdiff / y25rdiff - 0.066, linecolor=2, newplot=F)
pg.plotxy(freq, y47rdiff / y25rdiff - 0.074, linecolor=2, newplot=F)
pg.plotxy(freq, y49rdiff / y25rdiff - 0.072, linecolor=2, newplot=F)
pg.plotxy(freq, y51rdiff / y25rdiff - 0.079, linecolor=2, newplot=F)
pg.plotxy(freq, y53rdiff / y25rdiff - 0.082, linecolor=2, newplot=F)
pg.plotxy(freq, y55rdiff / y25rdiff - 0.17, linecolor=2, newplot=F)
pg.plotxy(freq, y57rdiff / y25rdiff - 0.055, linecolor=2, newplot=F)
pg.plotxy(freq, y59rdiff / y25rdiff - 0.108, linecolor=2, newplot=F)
pg.plotxy(freq, y27rdiff / y25rdiff - 0.17, linecolor=3, newplot=F)
# yrdiffrat.ps

mask := (length(freq) / 10):(9 *  length(freq) / 10)
avgxr[1] := sum((x27rdiff / x25rdiff)[mask]) / length(mask);
avgxr[2] := sum((x29rdiff / x25rdiff)[mask]) / length(mask);
avgxr[3] := sum((x31rdiff / x25rdiff)[mask]) / length(mask);
avgxr[4] := sum((x35rdiff / x25rdiff)[mask]) / length(mask);
avgxr[5] := sum((x37rdiff / x25rdiff)[mask]) / length(mask);
avgxr[6] := sum((x39rdiff / x25rdiff)[mask]) / length(mask);
avgxr[7] := sum((x41rdiff / x25rdiff)[mask]) / length(mask);
avgxr[8] := sum((x43rdiff / x25rdiff)[mask]) / length(mask);
avgxr[9] := sum((x45rdiff / x25rdiff)[mask]) / length(mask);
avgxr[10] := sum((x47rdiff / x25rdiff)[mask]) / length(mask);
avgxr[11] := sum((x49rdiff / x25rdiff)[mask]) / length(mask);
avgxr[12] := sum((x51rdiff / x25rdiff)[mask]) / length(mask);
avgxr[13] := sum((x53rdiff / x25rdiff)[mask]) / length(mask);
avgxr[14] := sum((x55rdiff / x25rdiff)[mask]) / length(mask);
avgxr[15] := sum((x57rdiff / x25rdiff)[mask]) / length(mask);
avgxr[16] := sum((x59rdiff / x25rdiff)[mask]) / length(mask);
pg.clear(); pg.plotxy((1/sin(pi*elev/180)) / (1/sin(pi*36.68/180)), avgxr,
                      xtitle='Secant(z) Ratio',
                      ytitle='Source Power Ratio', linecolor=2, newplot=T)
pg.setyscale(0.85, 1.05)
avgyr[1] := sum((y27rdiff / y25rdiff)[mask]) / length(mask);
avgyr[2] := sum((y29rdiff / y25rdiff)[mask]) / length(mask);
avgyr[3] := sum((y31rdiff / y25rdiff)[mask]) / length(mask);
avgyr[4] := sum((y35rdiff / y25rdiff)[mask]) / length(mask);
avgyr[5] := sum((y37rdiff / y25rdiff)[mask]) / length(mask);
avgyr[6] := sum((y39rdiff / y25rdiff)[mask]) / length(mask);
avgyr[7] := sum((y41rdiff / y25rdiff)[mask]) / length(mask);
avgyr[8] := sum((y43rdiff / y25rdiff)[mask]) / length(mask);
avgyr[9] := sum((y45rdiff / y25rdiff)[mask]) / length(mask);
avgyr[10] := sum((y47rdiff / y25rdiff)[mask]) / length(mask);
avgyr[11] := sum((y49rdiff / y25rdiff)[mask]) / length(mask);
avgyr[12] := sum((y51rdiff / y25rdiff)[mask]) / length(mask);
avgyr[13] := sum((y53rdiff / y25rdiff)[mask]) / length(mask);
avgyr[14] := sum((y55rdiff / y25rdiff)[mask]) / length(mask);
avgyr[15] := sum((y57rdiff / y25rdiff)[mask]) / length(mask);
avgyr[16] := sum((y59rdiff / y25rdiff)[mask]) / length(mask);
pg.plotxy((1/sin(pi*elev/180)) / (1/sin(pi*36.68/180)), avgyr,
          xtitle='Secant(z) Ratio',
          ytitle='Source Power Ratio', linecolor=3, newplot=F)
# avgrdiffrat.ps