Home > cbassSouthFunctions > ACTFunctions > HartRAO_Beam > plot_sun_scan.m

plot_sun_scan

PURPOSE ^

%

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

%
Plot roach data for sun scans -  averaged over all frequencies

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 
0002 %%
0003 %Plot roach data for sun scans -  averaged over all frequencies
0004 
0005 %% Note Roach 2 frequencies are inverted i.e frequency couns down (5.5-->5GHz)from channel 65:128
0006 % Read in simulated beam
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 % Read in azimuth scan data
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 %Get Averaged Intensity (LL,RR on each of raoch 1 and 2) - subtract a zero
0019 %level - not perfect but will do for now to see things
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 % Plot the data
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 % To plot the beam channel by channel to check the amplitude vs noise
0043 % diode amplitude
0044 LL = [d.antenna0.roach1.LL fliplr(d.antenna0.roach2.LL)];%so actual frequency counts up on 2nd roach
0045 RR = [d.antenna0.roach1.RR fliplr(d.antenna0.roach2.RR)];%so actual frequency counts up on 2nd roach
0046 
0047 nd_L = 0.9 % Kelvin
0048 nd_R = 0.63 % Kelvin
0049 
0050 %Left polarization
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 %Right polarization
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))%so actual frequency counts up
0085 RR1_mean = nanmean(beam_chansR(3,1:64))
0086 RR1_mean = nanmean(beam_chansR(3,65:128))%so actual frequency counts up

Generated on Sun 14-Jun-2015 17:12:45 by m2html © 2005