lukesSouthColdLoadTest

% Read Data

0001 %% Read Data
0002 d = read_arcSouth('25-Feb-2014:23:55:20','26-Feb-2014:01:43:25')
0003 
0004 %% Get params we want
0005 LL = [ d.antenna0.roach1.LL fliplr(d.antenna0.roach2.LL) ];
0006 RR = [ d.antenna0.roach1.RR fliplr(d.antenna0.roach2.RR) ];
0007 load1 = [ d.antenna0.roach1.load1 fliplr(d.antenna0.roach2.load1) ];
0008 load2 = [ d.antenna0.roach1.load2 fliplr(d.antenna0.roach2.load2) ];
0009 loadT = d.antenna0.thermal.ccTemperatureLoadFast;
0010 
0011 nchan = 128;
0012 chan = 1;
0013 %% When do the loads and the LL/RR cross?
0014 while chan<127
0015     difference = LL(:,chan)-load2(:,chan);
0016     difference(find(difference<300 & difference>-300)) = 0;
0017     ind = find( not(difference) );
0018     temp(chan) = median( loadT(ind) );
0019     chan = chan + 1;
0020 end
0021 
0022 chan = 120
0023     difference = LL(:,chan)-load2(:,chan);
0024     difference(find(difference<100 & difference>-100)) = 0;
0025     ind = find( not(difference) );
0026     temp(chan) = median( loadT(ind) );
0027     chan = chan + 1;