###############################################################################
                            Project TBASEJRF021227
################################################################################
 
General Comments:

Goal is to study properties of L-band stability of continuum
source frequency response.

Used optical modules 1 and 3

o Spectrometer setup using 50 MHz BW, 9-level, 2 samplers, 1 quadrant,
  8192 lags per IF,  1N2-0A-50-9

 OR  CM  SF  AC sampler
 --  --  --  ----------
 1   1   1   0
 3   5   5   4

o LO setup: 

where IF1 =  3000 MHz
      LO2 = 13075 MHz
      LO3 = 10500 MHz
      IF3 =   425 MHz
      LO4 =   500 MHz
      IF4 =    75 MHz (AC BW=50 MHz)

2002_12_27_07:17:38A.fits scan: 1 1042+1203 2 spectra 2 phases 10 recs
   az: 120.2 el: 48.95 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_07:23:00A.fits scan: 2 1042+1203 2 spectra 2 phases 10 recs
   az: 120.2 el: 48.92 ra2000: 10.8 dec2000: 12.06 freq: 1405
show_avg_scan_diff(1, 2, spec_num=1, decimate=4, scale=1.0, linecolor=2, newplot=T, cf_offset=1405, usb=F)
show_avg_scan_diff(1, 2, spec_num=2, decimate=4, scale=1.0, linecolor=3, newplot=F, cf_offset=1405, usb=F)
pg.setyscale(0.3, 0.4)

2002_12_27_07:28:22A.fits scan: 3 1042+1203N 2 spectra 2 phases 10 recs
   az: 123.1 el: 50.78 ra2000: 10.71 dec2000: 12.13 freq: 1405
2002_12_27_07:33:45A.fits scan: 4 1042+1203N 2 spectra 2 phases 10 recs
   az: 123.1 el: 50.76 ra2000: 10.8 dec2000: 12.13 freq: 1405
2002_12_27_07:39:08A.fits scan: 5 1042+1203 2 spectra 2 phases 10 recs
   az: 126.4 el: 52.47 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_07:44:31A.fits scan: 6 1042+1203 2 spectra 2 phases 10 recs
   az: 126.3 el: 52.45 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_07:49:55A.fits scan: 7 1042+1203S 2 spectra 2 phases 10 recs
   az: 129.9 el: 54.08 ra2000: 10.71 dec2000: 12 freq: 1405
2002_12_27_07:55:19A.fits scan: 8 1042+1203S 2 spectra 2 phases 10 recs
   az: 129.8 el: 54.06 ra2000: 10.8 dec2000: 11.99 freq: 1405
2002_12_27_08:00:43A.fits scan: 9 1042+1203 2 spectra 2 phases 10 recs
   az: 133.5 el: 55.72 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_08:06:08A.fits scan: 10 1042+1203 2 spectra 2 phases 10 recs
   az: 133.4 el: 55.7 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_08:11:32A.fits scan: 11 1042+1203E 2 spectra 2 phases 10 recs
   az: 137.3 el: 57.17 ra2000: 10.72 dec2000: 12.06 freq: 1405
2002_12_27_08:16:57A.fits scan: 12 1042+1203E 2 spectra 2 phases 10 recs
   az: 137.3 el: 57.15 ra2000: 10.81 dec2000: 12.06 freq: 1405
2002_12_27_08:22:23A.fits scan: 13 1042+1203 2 spectra 2 phases 10 recs
   az: 141.7 el: 58.59 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_08:27:49A.fits scan: 14 1042+1203 2 spectra 2 phases 10 recs
   az: 141.7 el: 58.56 ra2000: 10.8 dec2000: 12.05 freq: 1405
2002_12_27_08:33:16A.fits scan: 15 1042+1203W 2 spectra 2 phases 10 recs
   az: 146.4 el: 59.87 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_08:38:42A.fits scan: 16 1042+1203W 2 spectra 2 phases 10 recs
   az: 146.4 el: 59.85 ra2000: 10.8 dec2000: 12.05 freq: 1405
2002_12_27_08:44:08A.fits scan: 17 1042+1203 2 spectra 2 phases 10 recs
   az: 151.2 el: 60.94 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_08:49:35A.fits scan: 18 1042+1203 2 spectra 2 phases 10 recs
   az: 151.2 el: 60.93 ra2000: 10.8 dec2000: 12.05 freq: 1405
2002_12_27_08:55:02A.fits scan: 19 1042+1203N 2 spectra 2 phases 10 recs
   az: 156.4 el: 61.95 ra2000: 10.71 dec2000: 12.13 freq: 1405
2002_12_27_09:00:30A.fits scan: 20 1042+1203N 2 spectra 2 phases 10 recs
   az: 156.3 el: 61.94 ra2000: 10.8 dec2000: 12.12 freq: 1405
2002_12_27_09:05:58A.fits scan: 21 1042+1203 2 spectra 2 phases 10 recs
   az: 161.9 el: 62.65 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_09:11:26A.fits scan: 22 1042+1203 2 spectra 2 phases 10 recs
   az: 161.9 el: 62.64 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_09:16:54A.fits scan: 23 1042+1203S 2 spectra 2 phases 10 recs
   az: 167.7 el: 63.14 ra2000: 10.71 dec2000: 11.99 freq: 1405
2002_12_27_09:22:23A.fits scan: 24 1042+1203S 2 spectra 2 phases 10 recs
   az: 167.7 el: 63.13 ra2000: 10.8 dec2000: 11.99 freq: 1405
2002_12_27_09:27:52A.fits scan: 25 1042+1203 2 spectra 2 phases 10 recs
   az: 173.6 el: 63.55 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_09:33:21A.fits scan: 26 1042+1203 2 spectra 2 phases 10 recs
   az: 173.6 el: 63.55 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_09:38:49A.fits scan: 27 1042+1203E 2 spectra 2 phases 10 recs
   az: 179.5 el: 63.68 ra2000: 10.72 dec2000: 12.06 freq: 1405
2002_12_27_09:44:18A.fits scan: 28 1042+1203E 2 spectra 2 phases 10 recs
   az: 179.5 el: 63.68 ra2000: 10.81 dec2000: 12.06 freq: 1405
2002_12_27_09:49:47A.fits scan: 29 1042+1203 2 spectra 2 phases 10 recs
   az: 185.7 el: 63.58 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_09:55:16A.fits scan: 30 1042+1203 2 spectra 2 phases 10 recs
   az: 185.6 el: 63.58 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_10:00:45A.fits scan: 31 1042+1203W 2 spectra 2 phases 10 recs
   az: 191.8 el: 63.24 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_10:06:13A.fits scan: 32 1042+1203W 2 spectra 2 phases 10 recs
   az: 191.7 el: 63.25 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_10:11:41A.fits scan: 33 1042+1203 2 spectra 2 phases 10 recs
   az: 197.4 el: 62.71 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_10:17:09A.fits scan: 34 1042+1203 2 spectra 2 phases 10 recs
   az: 197.4 el: 62.72 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_10:22:37A.fits scan: 35 1042+1203N 2 spectra 2 phases 10 recs
   az: 203 el: 62.04 ra2000: 10.71 dec2000: 12.12 freq: 1405
2002_12_27_10:28:05A.fits scan: 36 1042+1203N 2 spectra 2 phases 10 recs
   az: 202.9 el: 62.05 ra2000: 10.8 dec2000: 12.13 freq: 1405
2002_12_27_10:33:32A.fits scan: 37 1042+1203 2 spectra 2 phases 10 recs
   az: 208.2 el: 61.05 ra2000: 10.71 dec2000: 12.05 freq: 1405
2002_12_27_10:38:59A.fits scan: 38 1042+1203 2 spectra 2 phases 10 recs
   az: 208.1 el: 61.06 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_10:44:25A.fits scan: 39 1042+1203S 2 spectra 2 phases 10 recs
   az: 213.1 el: 59.9 ra2000: 10.71 dec2000: 11.99 freq: 1405
2002_12_27_10:49:51A.fits scan: 40 1042+1203S 2 spectra 2 phases 10 recs
   az: 213 el: 59.92 ra2000: 10.8 dec2000: 12 freq: 1405
2002_12_27_10:55:17A.fits scan: 41 1042+1203 2 spectra 2 phases 10 recs
   az: 217.7 el: 58.72 ra2000: 10.71 dec2000: 12.05 freq: 1405
2002_12_27_11:00:43A.fits scan: 42 1042+1203 2 spectra 2 phases 10 recs
   az: 217.7 el: 58.74 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_11:06:08A.fits scan: 43 1042+1203E 2 spectra 2 phases 10 recs
   az: 221.9 el: 57.39 ra2000: 10.72 dec2000: 12.05 freq: 1405
2002_12_27_11:11:33A.fits scan: 44 1042+1203E 2 spectra 2 phases 10 recs
   az: 221.9 el: 57.42 ra2000: 10.81 dec2000: 12.06 freq: 1405
2002_12_27_11:16:58A.fits scan: 45 1042+1203 2 spectra 2 phases 10 recs
   az: 226 el: 55.88 ra2000: 10.71 dec2000: 12.05 freq: 1405
2002_12_27_11:22:22A.fits scan: 46 1042+1203 2 spectra 2 phases 10 recs
   az: 226 el: 55.91 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_11:27:46A.fits scan: 47 1042+1203W 2 spectra 2 phases 10 recs
   az: 229.8 el: 54.27 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_11:33:10A.fits scan: 48 1042+1203W 2 spectra 2 phases 10 recs
   az: 229.8 el: 54.29 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_11:38:33A.fits scan: 49 1042+1203 2 spectra 2 phases 10 recs
   az: 233.2 el: 52.66 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_11:43:56A.fits scan: 50 1042+1203 2 spectra 2 phases 10 recs
   az: 233.2 el: 52.68 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_11:49:19A.fits scan: 51 1042+1203N 2 spectra 2 phases 10 recs
   az: 236.5 el: 50.98 ra2000: 10.71 dec2000: 12.12 freq: 1405
2002_12_27_11:54:42A.fits scan: 52 1042+1203N 2 spectra 2 phases 10 recs
   az: 236.5 el: 51.01 ra2000: 10.8 dec2000: 12.12 freq: 1405
2002_12_27_12:00:04A.fits scan: 53 1042+1203 2 spectra 2 phases 10 recs
   az: 239.4 el: 49.15 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_12:05:26A.fits scan: 54 1042+1203 2 spectra 2 phases 10 recs
   az: 239.4 el: 49.17 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_12:10:48A.fits scan: 55 1042+1203S 2 spectra 2 phases 10 recs
   az: 242.2 el: 47.26 ra2000: 10.71 dec2000: 11.99 freq: 1405
2002_12_27_12:16:10A.fits scan: 56 1042+1203S 2 spectra 2 phases 10 recs
   az: 242.2 el: 47.29 ra2000: 10.8 dec2000: 11.99 freq: 1405
2002_12_27_12:21:32A.fits scan: 57 1042+1203 2 spectra 2 phases 10 recs
   az: 244.9 el: 45.42 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_12:26:53A.fits scan: 58 1042+1203 2 spectra 2 phases 10 recs
   az: 244.9 el: 45.45 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_12:32:14A.fits scan: 59 1042+1203E 2 spectra 2 phases 10 recs
   az: 247.3 el: 43.55 ra2000: 10.72 dec2000: 12.06 freq: 1405
2002_12_27_12:37:35A.fits scan: 60 1042+1203E 2 spectra 2 phases 10 recs
   az: 247.3 el: 43.58 ra2000: 10.81 dec2000: 12.06 freq: 1405
2002_12_27_12:42:56A.fits scan: 61 1042+1203 2 spectra 2 phases 10 recs
   az: 249.8 el: 41.55 ra2000: 10.71 dec2000: 12.06 freq: 1405
2002_12_27_12:48:16A.fits scan: 62 1042+1203 2 spectra 2 phases 10 recs
   az: 249.7 el: 41.58 ra2000: 10.8 dec2000: 12.06 freq: 1405
2002_12_27_12:53:37A.fits scan: 63 1042+1203W 2 spectra 2 phases 10 recs
   az: 252.1 el: 39.52 ra2000: 10.71 dec2000: 12.06 freq: 1405


az[1] := 120.2; el[1] := 48.95;
az[2] := 120.2; el[2] := 48.92;
az[3] := 123.1; el[3] := 50.78;
az[4] := 123.1; el[4] := 50.76;
az[5] := 126.4; el[5] := 52.47;
az[6] := 126.3; el[6] := 52.45;
az[7] := 129.9; el[7] := 54.08;
az[8] := 129.8; el[8] := 54.06;
az[9] := 133.5; el[9] := 55.72;
az[10] := 133.4; el[10] := 55.70;
az[11] := 137.3; el[11] := 57.17;
az[12] := 137.3; el[12] := 57.15;
az[13] := 141.7; el[13] := 58.59;
az[14] := 141.7; el[14] := 58.56;
az[15] := 146.4; el[15] := 59.87;
az[16] := 146.4; el[16] := 59.85;
az[17] := 151.2; el[17] := 60.94;
az[18] := 151.2; el[18] := 60.93;
az[19] := 156.4; el[19] := 61.95;
az[20] := 156.3; el[20] := 61.94;
az[21] := 161.9; el[21] := 62.65;
az[22] := 161.9; el[22] := 62.64;
az[23] := 167.7; el[23] := 63.14;
az[24] := 167.7; el[24] := 63.13;
az[25] := 173.6; el[25] := 63.55;
az[26] := 173.6; el[26] := 63.55;
az[27] := 179.5; el[27] := 63.68;
az[28] := 179.5; el[28] := 63.68;
az[29] := 185.7; el[29] := 63.58;
az[30] := 185.6; el[30] := 63.58;
az[31] := 191.8; el[31] := 63.24;
az[32] := 191.7; el[32] := 63.25;
az[33] := 197.4; el[33] := 62.71;
az[34] := 197.4; el[34] := 62.72;
az[35] := 203.0; el[35] := 62.04;
az[36] := 202.9; el[36] := 62.05;
az[37] := 208.2; el[37] := 61.05;
az[38] := 208.1; el[38] := 61.06;
az[39] := 213.1; el[39] := 59.90;
az[40] := 213.0; el[40] := 59.92;
az[41] := 217.7; el[41] := 58.72;
az[42] := 217.7; el[42] := 58.74;
az[43] := 221.9; el[43] := 57.39;
az[44] := 221.9; el[44] := 57.42;
az[45] := 226.0; el[45] := 55.88;
az[46] := 226.0; el[46] := 55.91;
az[47] := 229.8; el[47] := 54.27;
az[48] := 229.8; el[48] := 54.29;
az[49] := 233.2; el[49] := 52.66;
az[50] := 233.2; el[50] := 52.68;
az[51] := 236.5; el[51] := 50.98;
az[52] := 236.5; el[52] := 51.01;
az[53] := 239.4; el[53] := 49.15;
az[54] := 239.4; el[54] := 49.17;
az[55] := 242.2; el[55] := 47.26;
az[56] := 242.2; el[56] := 47.29;
az[57] := 244.9; el[57] := 45.42;
az[58] := 244.9; el[58] := 45.45;
az[59] := 247.3; el[59] := 43.55;
az[60] := 247.3; el[60] := 43.58;
az[61] := 249.8; el[61] := 41.55;
az[62] := 249.7; el[62] := 41.58;
az[63] := 252.1; el[63] := 39.52;
center_on := seq(1,61,4);
num_pairs := length(center_on);
rdiffX := array(0.0, num_pairs, 2049);
sdiffX := array(0.0, num_pairs, 2049);
rdiffY := array(0.0, num_pairs, 2049);
sdiffY := array(0.0, num_pairs, 2049);
j := 0;
for (i in center_on) {
    j +:= 1;
    sn[j] := i;
    show_avg_scan_diff(i, i+1, spec_num=1, decimate=4, scale=1.0, linecolor=2, newplot=T, cf_offset=1405, usb=F)
    rdiffX[j, ] := rdiff;
    sdiffX[j, ] := sdiff;
    pg.setyscale(0.3, 0.4);
    show_avg_scan_diff(i, i+1, spec_num=2, decimate=4, scale=1.0, linecolor=2, newplot=T, cf_offset=1405, usb=F)
    rdiffY[j, ] := rdiff;
    sdiffY[j, ] := sdiff;
    pg.setyscale(0.3, 0.4);
}
jseq := seq(3,num_pairs,3);
pg.clear(); np := T; lc := 1;
lr1 := linear_response_az(100*.184/3.31, -83.1, az[sn[1]], el[sn[1]])
for (j in jseq) {
    lc +:= 1;
    lr :=  linear_response_az(100*.184/3.31, -83.1, az[sn[j]], el[sn[j]])
print lr1[1] / lr[1];
    pg.plotxy(freq, (lr1[1] / lr[1]) * rdiffX[j, ] / rdiffX[1, ], xtitle='Frequency in MHz', ytitle='Ratio of Spectral Power Density', title='(on - off) Ratio, Channel X, Beam Center, 1042+1203', linecolor=lc, newplot=np);
    if (np) pg.setyscale(0.95, 1.05);
    np := F;
}
# Xrcenter.ps
pg.clear(); np := T; lc := 1;
lr1 := linear_response_az(100*.184/3.31, -83.1, az[sn[1]], el[sn[1]])
for (j in jseq) {
    lc +:= 1;
    lr :=  linear_response_az(100*.184/3.31, -83.1, az[sn[j]], el[sn[j]])
    pg.plotxy(freq, (lr1[2] / lr[2]) * rdiffY[j, ] / rdiffY[1, ], xtitle='Frequency in MHz', ytitle='Ratio of Spectral Power Density', title='(on - off) Ratio, Channel Y, Beam Center, 1042+1203', linecolor=lc, newplot=np);
    if (np) pg.setyscale(0.95, 1.05);
    np := F;
}
# Yrcenter.ps
pg.clear(); np := T; lc := 1;
lr1 := linear_response_az(100*.184/3.31, -83.1, az[sn[1]], el[sn[1]])
for (j in jseq) {
    lc +:= 1;
    lr :=  linear_response_az(100*.184/3.31, -83.1, az[sn[j]], el[sn[j]])
    pg.plotxy(freq, (lr1[1] / lr[1]) * sdiffX[j, ] / sdiffX[1, ], xtitle='Frequency in MHz', ytitle='Ratio of Spectral Power Density', title='(on - off) / off Ratio, Channel X, Beam Center, 1042+1203', linecolor=lc, newplot=np);
    if (np) pg.setyscale(0.95, 1.05);
    np := F;
}
# Xscenter.ps
pg.clear(); np := T; lc := 1;
lr1 := linear_response_az(100*.184/3.31, -83.1, az[sn[1]], el[sn[1]])
for (j in jseq) {
    lc +:= 1;
    lr :=  linear_response_az(100*.184/3.31, -83.1, az[sn[j]], el[sn[j]])
    pg.plotxy(freq, (lr1[2] / lr[2]) * sdiffY[j, ] / sdiffY[1, ], xtitle='Frequency in MHz', ytitle='Ratio of Spectral Power Density', title='(on - off) / off Ratio, Channel Y, Beam Center, 1042+1203', linecolor=lc, newplot=np);
    if (np) pg.setyscale(0.95, 1.05);
    np := F;
}
# Yscenter.ps

rdiffX := array(0.0, 31, 2049);
rdiffY := array(0.0, 31, 2049);
sdiffX := array(0.0, 31, 2049);
sdiffY := array(0.0, 31, 2049);
for (i in 1:31) {
    sn := 1 + 2 * (i - 1);
    show_avg_scan_diff(sn, sn + 1, spec_num=1, decimate=4, scale=1.0,
		       linecolor=2, newplot=T, cf_offset=1405, usb=F)
    pg.setyscale(0, 0.4);
    rdiffX[i, ] := rdiff;
    sdiffX[i, ] := sdiff;
    show_avg_scan_diff(sn, sn + 1, spec_num=2, decimate=4, scale=1.0,
		       linecolor=2, newplot=T, cf_offset=1405, usb=F)
    pg.setyscale(0, 0.4);
    rdiffY[i, ] := rdiff;
    sdiffY[i, ] := sdiff;
}
write_value([rdiffX=rdiffX, sdiffX=sdiffX, rdiffY=rdiffY, sdiffY=sdiffY],
            'diffs.dat');

dd := read_value('diffs.dat');
lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[1], el[1]);
np := T;
for (i in seq(3, 31, 4)) {
    sn := 1 + 2 * (i - 1);
    lr3 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn])
    if (np) pg.clear();
    pg.plotxy(freq, (lr1[1] / lr3[1]) * dd.rdiffX[i, ] / dd.rdiffX[1, ],
              xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='Center scanN/scan1 (on - off) Ratio, Channel X, 1042+1203',
              linecolor=as_integer(i / 4) + 2, newplot=np);
    if (np) pg.setyscale(0.94, 1.02);
    np := F;
}
np := T;
for (i in seq(3, 31, 2)) {
    sn := 1 + 2 * (i - 1);
    lr3 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn])
    if (np) pg.clear();
    pg.plotxy(freq, (lr1[1] / lr3[1]) * dd.sdiffX[i, ] / dd.sdiffX[1, ],
              xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='Center scanN/scan1 (on - off) / off Ratio, Channel X, 1042+1203',
              linecolor=as_integer(i / 2) + 1, newplot=np);
    if (np) pg.setyscale(0.9, 1.1);
    np := F;
}
np := T;
for (i in seq(3, 31, 2)) {
    sn := 1 + 2 * (i - 1);
    lr3 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn])
    if (np) pg.clear();
    pg.plotxy(freq, (lr1[2] / lr3[2]) * dd.rdiffY[i, ] / dd.rdiffY[1, ],
              xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='Center scanN/scan1 (on - off) Ratio, Channel Y, 1042+1203',
              linecolor=as_integer(i / 2) + 1, newplot=np);
    if (np) pg.setyscale(0.96, 1.04);
    np := F;
}
np := T;
for (i in seq(3, 31, 2)) {
    sn := 1 + 2 * (i - 1);
    lr3 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn])
    if (np) pg.clear();
    pg.plotxy(freq, (lr1[2] / lr3[2]) * dd.sdiffY[i, ] / dd.sdiffY[1, ],
              xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='Center scanN/scan1 (on - off) / off Ratio, Channel Y, North, 1042+1203',
              linecolor=as_integer(i / 2) + 1, newplot=np);
    if (np) pg.setyscale(0.9, 1.1);
    np := F;
}

dd := read_value('diffs.dat');
diff_ratio := array(0.0, dd.rdiffX::shape[2]);
num_rats := 0
for (i in seq(1, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.rdiffX[i + 1, ] / dd.rdiffX[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
#              title='North/Center (on - off) Ratio, Channel X, 1042+1203',
              title='(N|S|E|W)/Center (on - off) Ratio, Channel X, 1042+1203',
              linecolor=2, newplot=T);
pg.setyscale(0.6, 0.7);
pg.plotxy(freq, beam_freq_response(freq, .64215), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.rdiffX::shape[2]);
num_rats := 0
for (i in seq(3, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.rdiffX[i + 1, ] / dd.rdiffX[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='South/Center (on - off) Ratio, Channel X, 1042+1203',
              linecolor=3, newplot=F);
pg.setyscale(0.5, 0.6);
pg.plotxy(freq, beam_freq_response(freq, .55897), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.rdiffX::shape[2]);
num_rats := 0
for (i in seq(5, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.rdiffX[i + 1, ] / dd.rdiffX[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='East/Center (on - off) Ratio, Channel X, 1042+1203',
              linecolor=5, newplot=F);
pg.setyscale(0.6, 0.7);
pg.plotxy(freq, beam_freq_response(freq, .631), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.rdiffX::shape[2]);
num_rats := 0
for (i in seq(7, 30, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.rdiffX[i + 1, ] / dd.rdiffX[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='West/Center (on - off) Ratio, Channel X, 1042+1203',
              linecolor=6, newplot=F);
pg.setyscale(0.55, 0.65);
pg.plotxy(freq, beam_freq_response(freq, .605), linecolor=4, newplot=F);
pg.setyscale(0.5, 0.7);
# Xroffset.ps
----------------------------------------

dd := read_value('diffs.dat');
diff_ratio := array(0.0, dd.rdiffY::shape[2]);
num_rats := 0
for (i in seq(1, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.rdiffY[i + 1, ] / dd.rdiffY[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
#              title='North/Center (on - off) Ratio, Channel Y, 1042+1203',
              title='(N|S|E|W)/Center (on - off) Ratio, Channel Y, 1042+1203',
              linecolor=2, newplot=T);
pg.setyscale(0.6, 0.7);
pg.plotxy(freq, beam_freq_response(freq, .663), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.rdiffY::shape[2]);
num_rats := 0
for (i in seq(3, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.rdiffY[i + 1, ] / dd.rdiffY[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='South/Center (on - off) Ratio, Channel Y, 1042+1203',
              linecolor=3, newplot=F);
pg.setyscale(0.5, 0.6);
pg.plotxy(freq, beam_freq_response(freq, .5435), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.rdiffY::shape[2]);
num_rats := 0
for (i in seq(5, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.rdiffY[i + 1, ] / dd.rdiffY[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='East/Center (on - off) Ratio, Channel Y, 1042+1203',
              linecolor=5, newplot=F);
pg.setyscale(0.6, 0.7);
pg.plotxy(freq, beam_freq_response(freq, .63), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.rdiffY::shape[2]);
num_rats := 0
for (i in seq(7, 30, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.rdiffY[i + 1, ] / dd.rdiffY[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='West/Center (on - off) Ratio, Channel Y, 1042+1203',
              linecolor=6, newplot=F);
pg.setyscale(0.55, 0.65);
pg.plotxy(freq, beam_freq_response(freq, .603), linecolor=4, newplot=F);
pg.setyscale(0.5, 0.7);
# Yroffset.ps
---------------------------------

dd := read_value('diffs.dat');
diff_ratio := array(0.0, dd.sdiffX::shape[2]);
num_rats := 0
for (i in seq(1, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.sdiffX[i + 1, ] / dd.sdiffX[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
#              title='North/Center (on - off) / off Ratio, Channel X, 1042+1203',
              title='(N|S|E|W)/Center (on - off) / off Ratio, Channel X, 1042+1203',
              linecolor=2, newplot=T);
pg.setyscale(0.6, 0.7);
pg.plotxy(freq, beam_freq_response(freq, .64215), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.sdiffX::shape[2]);
num_rats := 0
for (i in seq(3, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.sdiffX[i + 1, ] / dd.sdiffX[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='South/Center (on - off) / off Ratio, Channel X, 1042+1203',
              linecolor=3, newplot=F);
pg.setyscale(0.5, 0.6);
pg.plotxy(freq, beam_freq_response(freq, .55897), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.sdiffX::shape[2]);
num_rats := 0
for (i in seq(5, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.sdiffX[i + 1, ] / dd.sdiffX[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='East/Center (on - off) / off Ratio, Channel X, 1042+1203',
              linecolor=5, newplot=F);
pg.setyscale(0.6, 0.7);
pg.plotxy(freq, beam_freq_response(freq, .631), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.sdiffX::shape[2]);
num_rats := 0
for (i in seq(7, 30, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.sdiffX[i + 1, ] / dd.sdiffX[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='West/Center (on - off) / off Ratio, Channel X, 1042+1203',
              linecolor=6, newplot=F);
pg.setyscale(0.55, 0.65);
pg.plotxy(freq, beam_freq_response(freq, .605), linecolor=4, newplot=F);
pg.setyscale(0.5, 0.7);
# Xsoffset.ps
---------------------------------------

dd := read_value('diffs.dat');
diff_ratio := array(0.0, dd.sdiffY::shape[2]);
num_rats := 0
for (i in seq(1, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.sdiffY[i + 1, ] / dd.sdiffY[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
#              title='North/Center (on - off) / off Ratio, Channel Y, 1042+1203',
              title='(N|S|E|W)/Center (on - off) / off Ratio, Channel Y, 1042+1203',
              linecolor=2, newplot=T);
pg.setyscale(0.6, 0.7);
pg.plotxy(freq, beam_freq_response(freq, .663), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.sdiffY::shape[2]);
num_rats := 0
for (i in seq(3, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.sdiffY[i + 1, ] / dd.sdiffY[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='South/Center (on - off) / off Ratio, Channel Y, 1042+1203',
              linecolor=3, newplot=F);
pg.setyscale(0.5, 0.6);
pg.plotxy(freq, beam_freq_response(freq, .5435), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.sdiffY::shape[2]);
num_rats := 0
for (i in seq(5, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.sdiffY[i + 1, ] / dd.sdiffY[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='East/Center (on - off) / off Ratio, Channel Y, 1042+1203',
              linecolor=5, newplot=F);
pg.setyscale(0.6, 0.7);
pg.plotxy(freq, beam_freq_response(freq, .63), linecolor=4, newplot=F);

diff_ratio := array(0.0, dd.sdiffY::shape[2]);
num_rats := 0
for (i in seq(7, 30, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.sdiffY[i + 1, ] / dd.sdiffY[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
#pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='West/Center (on - off) / off Ratio, Channel Y, 1042+1203',
              linecolor=6, newplot=F);
pg.setyscale(0.55, 0.65);
pg.plotxy(freq, beam_freq_response(freq, .603), linecolor=4, newplot=F);
pg.setyscale(0.5, 0.7);
# Ysoffset.ps
---------------------------------

dd := read_value('diffs.dat');
diff_ratio := array(0.0, dd.rdiffX::shape[2]);
num_rats := 0
for (i in seq(3, 31, 8)) {
    sn := 1 + 2 * (i - 1);
    lr1 :=  linear_response_az(100*.184/3.31, -83.1, az[sn], el[sn]);
    lr2 :=  linear_response_az(100*.184/3.31, -83.1, az[sn + 1], el[sn + 1]);
    diff_ratio +:= (lr1[1] / lr2[1]) * dd.rdiffX[i + 1, ] / dd.rdiffX[i, ];
    num_rats +:= 1;
}
diff_ratio /:= num_rats;
pg.clear();
pg.plotxy(freq, diff_ratio, xtitle='Frequency in MHz',
              ytitle='Ratio of Spectral Power Density',
              title='South/Center (on - off) Ratio, Channel X, 1042+1203',
              linecolor=3, newplot=T)
pg.setyscale(0.5, 0.6);
pg.plotxy(freq, beam_freq_response(freq, .55897), linecolor=4, newplot=F);
spectrum_fft(diff_ratio, 50)
# XrFFTsouth.ps