0001
0002
0003
0004
0005
0006
0007
0008 [mean_power_ZP mean_power_P sim_beam sim_data sim_angles phi vals_sim FWHM_sim fit_angles simfile] = read_beam_simHart(3,100);
0009
0010
0011
0012 d=read_arcSouth('16-May-2013:10:15:12','16-May-2013:10:17:22')
0013
0014 az_scan = d.antenna0.servo.az;
0015 el_scan = d.antenna0.servo.el;
0016 az_src = az_scan.*cosd(el_scan);
0017
0018
0019
0020
0021 I(1,:) = d.antenna0.roach1.LLfreq-d.antenna0.roach1.LLfreq(3101);
0022 I(2,:) = d.antenna0.roach2.LLfreq-d.antenna0.roach2.LLfreq(3101);
0023 I(3,:) = d.antenna0.roach1.RRfreq-d.antenna0.roach1.RRfreq(3101);
0024 I(4,:) = d.antenna0.roach2.RRfreq-d.antenna0.roach2.RRfreq(3101);
0025 labels = {'LL1','LL2','RR1','RR2'}
0026
0027
0028 figure
0029 for chan=1:4
0030 subplot(2,2,chan)
0031 plot(az_src,I(chan,:)./max(I(chan,:)))
0032 hold all
0033 plot(sim_beam(:,1),sim_beam(:,2)./max(sim_beam(:,2)))
0034 ylim([0 1])
0035 xlim([-100 100])
0036 title(labels{chan})
0037 xlabel('Offset, deg')
0038 ylabel('Normalised amplitude')
0039 end
0040
0041
0042
0043
0044 LL = [d.antenna0.roach1.LL fliplr(d.antenna0.roach2.LL)];
0045 RR = [d.antenna0.roach1.RR fliplr(d.antenna0.roach2.RR)];
0046
0047 nd_L = 0.9
0048 nd_R = 0.63
0049
0050
0051 for chan=1:128;
0052 beam(:,chan) = LL(:,chan) - LL(3101,chan);
0053 beam(7088:7264,chan)=nan;
0054 beam_max = max(beam(:,chan));
0055 beam(:,chan) = beam(:,chan)/beam_max;
0056 baseline = nanmean(beam(2848:4777,chan));
0057 beam_nd = nanmean(beam(5741:6964,chan))- baseline;
0058 beam_sun= max(beam(:,chan))-baseline;
0059 nd = (beam_sun/beam_nd)*nd_L;
0060 beam_chansL(:,chan) = [beam_sun,beam_nd, nd];
0061
0062 end
0063
0064 for chan=1:128;
0065 beam(:,chan) = RR(:,chan) - RR(3101,chan);
0066 beam(7088:7264,chan)=nan;
0067 beam_max = max(beam(:,chan));
0068 beam(:,chan) = beam(:,chan)/beam_max;
0069 baseline = nanmean(beam(2848:4777,chan));
0070 beam_nd = nanmean(beam(5741:6964,chan))- baseline;
0071 beam_sun= max(beam(:,chan))-baseline;
0072 nd = (beam_sun/beam_nd)*nd_R;
0073 beam_chansR(:,chan) = [beam_sun,beam_nd, nd];
0074
0075 end
0076 plot(beam_chansL(3,:))
0077 hold all
0078 plot(beam_chansR(3,:))
0079 legend('LL','RR')
0080 xlabel('Channel')
0081 ylabel('Inferred Sun Temp, K (ND=1K)')
0082 title('Roach Channel number vs inferred Sun Temp (nd_L = 0.9K, nd_R=0.63K)')
0083 LL1_mean = nanmean(beam_chansL(3,1:64))
0084 LL2_mean = nanmean(beam_chansL(3,65:128))
0085 RR1_mean = nanmean(beam_chansR(3,1:64))
0086 RR1_mean = nanmean(beam_chansR(3,65:128))