getLoadAmplitude

0001 function [mjd, ch1mag, ch2mag, ch3mag, ch4mag, ch5mag, ch6mag, ch7mag, ch8mag] = getLoadAmplitude(data,load,debug)
0002 d = data     
0003 % first apply alpha and load correction
0004 %--------------------------------------
0005 d = determineIndices(d);
0006 d = apparentAzEl(d);
0007 
0008 
0009 if(~isfield(d.antenna0.servo, 'equa')) 
0010   display('Calculating RA/DEC');
0011   long=-118.2822;
0012   lat=37.2339;
0013   
0014   az = d.antenna0.servo.apparent(:,1);
0015   el = d.antenna0.servo.apparent(:,2);
0016   jd=mjd2jd(d.antenna0.receiver.utc);
0017   [equa] = horiz_coo([pi/180*(az) pi/180*(el)],jd,[pi/180*(long) ...
0018         pi/180*(lat)],'e'); 
0019   d.antenna0.servo.equa=equa;
0020   clear equa;
0021   clear az;
0022   clear el;
0023 end
0024 
0025 s = 'FILTERED';
0026 d = assembleAlphaStreams(d,s); 
0027 d = applyAlpha(d, s);
0028 %-----------------------------------
0029 
0030 if (load==1)
0031 alphaType = 0; %FILTERED mode, see loadWrapper
0032 scaleType = 0;
0033 [d] = load_template_corr(d,alphaType,scaleType);
0034 end
0035 
0036 % Select survey data based on the feature flag
0037 d = framecut(d, bitsearch(d.array.frame.features, 0, 'any'));
0038 
0039 mjd = median(d.antenna0.receiver.utc);
0040 data = d.antenna0.receiver.data;
0041 % 8 data channels
0042 ch1 = data(:,1);
0043 ch2 = data(:,2);
0044 ch3 = data(:,3);
0045 ch4 = data(:,4);
0046 ch5 = data(:,5);
0047 ch6 = data(:,6);
0048 ch7 = data(:,7);
0049 ch8 = data(:,8);
0050 
0051 % Get rid of everything that has a signal in it (over 150% of the mode)
0052 m1 = mode(ch1);
0053 ind = abs(ch1 - m1) > 0.2;
0054 ch1(ind) = 0; 
0055 ch2(ind) = 0; 
0056 ch3(ind) = 0; 
0057 ch4(ind) = 0; 
0058 ch5(ind) = 0; 
0059 ch6(ind) = 0; 
0060 ch7(ind) = 0; 
0061 ch8(ind) = 0; 
0062 
0063 % FFT the data
0064 [f1 Y1] = psmej(ch1, 0)
0065 [f2 Y2] = psmej(ch2, 0)
0066 [f3 Y3] = psmej(ch3, 0)
0067 [f4 Y4] = psmej(ch4, 0)
0068 [f5 Y5] = psmej(ch5, 0)
0069 [f6 Y6] = psmej(ch6, 0)
0070 [f7 Y7] = psmej(ch7, 0)
0071 [f8 Y8] = psmej(ch8, 0)
0072 
0073 if (debug==1)
0074 figure(1)
0075 plot(ch1)
0076 figure(2)
0077 loglog(f1,Y1)
0078 keyboard
0079 end
0080 
0081 % Select data at the relevant frequencies (1.2Hz and harmonics, and 1.5Hz for background)
0082 disp('Range')
0083 range1 = find(f1 > 1.198 & f1 < 1.202);
0084 range2 = find(f1 > 2.398 & f1 < 2.402);
0085 range3 = find(f1 > 3.598 & f1 < 3.602);
0086 range4 = find(f1 > 1.498 & f1 < 1.502);
0087 range5 = find(f1 > 10.0 & f1 < 30.0);
0088 %range = [range1, range2, range3, range4];
0089 
0090 ch1mag = [max(Y1(range1)),max(Y1(range2)),max(Y1(range3)),mean(Y1(range4)),mean(Y1(range5))];
0091 ch2mag = [max(Y2(range1)),max(Y2(range2)),max(Y2(range3)),mean(Y2(range4)),mean(Y2(range5))];
0092 ch3mag = [max(Y3(range1)),max(Y3(range2)),max(Y3(range3)),mean(Y3(range4)),mean(Y3(range5))];
0093 ch4mag = [max(Y4(range1)),max(Y4(range2)),max(Y4(range3)),mean(Y4(range4)),mean(Y4(range5))];
0094 ch5mag = [max(Y5(range1)),max(Y5(range2)),max(Y5(range3)),mean(Y5(range4)),mean(Y5(range5))];
0095 ch6mag = [max(Y6(range1)),max(Y6(range2)),max(Y6(range3)),mean(Y6(range4)),mean(Y6(range5))];
0096 ch7mag = [max(Y7(range1)),max(Y7(range2)),max(Y7(range3)),mean(Y7(range4)),mean(Y7(range5))];
0097 ch8mag = [max(Y8(range1)),max(Y8(range2)),max(Y8(range3)),mean(Y8(range4)),mean(Y8(range5))];
0098 
0099 end