steppingLoads

0001 d = read_arcSouth('20-Apr-2013:21:41:00','20-Apr-2013:22:00:00');% on 20 April
0002 
0003 d=read_arcSouth('22-Apr-2013:13:54:43','22-Apr-2013:14:21:56')
0004 
0005 
0006 Q= mean(d.antenna0.roach1.Q(:,10:50),2);
0007 U= mean(d.antenna0.roach1.U(:,10:50),2);
0008 Il= mean(d.antenna0.roach1.LL(:,10:50),2);
0009 Ir= mean(d.antenna0.roach1.RR(:,10:50),2);
0010 IlLoad= mean(d.antenna0.roach1.load2(:,10:50),2);
0011 IrLoad= mean(d.antenna0.roach1.load1(:,10:50),2);
0012 time=24*3600*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1));
0013 
0014 time=d.antenna0.roach1.utc;
0015 time2=d.antenna0.thermal.utc;
0016 [year, month, day, hour, minute, second] = mjd2date(time);
0017 dateee=[year, month, day, hour, minute, second];
0018 [year, month, day, hour, minute, second] = mjd2date(time2);
0019 dateff=[year, month, day, hour, minute, second];
0020 ntime=datenum(dateee);
0021 time2=ntime;
0022 ntime=datenum(dateff);
0023 time3=ntime;
0024 subplot(211)
0025 plot(time2,Il);
0026 hold all
0027 plot(time2,Ir);
0028 plot(time2,IlLoad);
0029 plot(time2,IrLoad);
0030 legend('LCP','RCP','LCPLoad','RCPLoad')
0031 grid on
0032 xlabel('Time','Fontsize',14)
0033 ylabel('B.U','Fontsize',14)
0034 title('CBASS.S Step Cold Load 22 April','Fontsize',16);
0035 datetick('x','yyyy-mm-dd HH:MM:SS')
0036 subplot(212)
0037 plot(time3,d.antenna0.thermal.ccTemperatureLoad);
0038 title('Cold Load','Fontsize',16);
0039 xlabel('Time','Fontsize',14)
0040 ylabel('Temperature [K]','Fontsize',14)
0041 datetick('x','yyyy-mm-dd HH:MM:SS')
0042 grid on
0043 
0044 
0045 
0046 d=read_arcSouth('26-Apr-2013:13:19:49','26-Apr-2013:13:48:11')
0047 
0048 time=d.antenna0.roach1.utc; %roach timestamps
0049 timeTemps=d.antenna0.thermal.utc; %temperature timestamps
0050 [year, month, day, hour, minute, second] = mjd2date(time);
0051 dateRoach=[year, month, day, hour, minute, second];
0052 [year, month, day, hour, minute, second] = mjd2date(timeTemps);
0053 dateTemps=[year, month, day, hour, minute, second];
0054 Roachtime=datenum(dateRoach);
0055 TempTime=datenum(dateTemps);
0056 
0057 plot(d.antenna0.roach2.LL(t_12,:))
0058 hold all
0059 plot(d.antenna0.roach2.load2(t_12,:))
0060 
0061 subplot(211)
0062 plot(d.antenna0.roach2.LL(:,35));
0063 hold all
0064 plot(d.antenna0.roach2.load2(:,35));
0065 legend('LCP','LCPLoad')
0066 
0067 grid on
0068 subplot(212)
0069 temps=interp(d.antenna0.thermal.ccTemperatureLoad,20);
0070 plot(temps);
0071 grid on
0072 title('Cold Load Temperature','Fontsize',16);
0073 xlabel('Sample','Fontsize',14)
0074 ylabel('Temperature [K]','Fontsize',14)
0075 
0076 tempIndex(1,:)=[12 1.2e4]
0077 tempIndex(2,:)=[14 3.166e4]
0078 tempIndex(3,:)=[16 4.921e4]
0079 tempIndex(4,:)=[18 6.53e4]
0080 tempIndex(5,:)=[20 8.28e4]
0081 tempIndex(6,:)=[22 1.02e5]
0082 tempIndex(7,:)=[24 1.2e5]
0083 tempIndex(8,:)=[26 1.374e5]
0084 tempIndex(9,:)=[28 1.538e5]
0085 tempIndex(10,:)=[30 1.69e5]
0086 
0087 tempIndex(1,:)=[12 1.2e4]
0088 tempIndex(2,:)=[30 1.69e5]
0089 
0090 frequency = 5.5-[1:64]*500/64./1000
0091 figure
0092 plot(frequency,10*log10(d.antenna0.roach2.LL(tempIndex(1,2),:)),'r','Linewidth',3)
0093 hold all
0094 plot(frequency,10*log10(d.antenna0.roach2.load2(tempIndex(1,2),:)),'.-')
0095 for i=2:2
0096     plot(frequency,10*log10(d.antenna0.roach2.LL(tempIndex(i,2),:)),'g','Linewidth',3)
0097     plot(frequency,10*log10(d.antenna0.roach2.load2(tempIndex(i,2),:)),'.-')
0098 end
0099 title('Roach2','Fontsize',16);
0100 xlabel('Frequency [GHz]','Fontsize',14)
0101 ylabel('Power','Fontsize',14)
0102 legend('Tsky@12K','Tload@12K','Tsky@30K','Tload@30K')
0103 grid on
0104 
0105 frequency = [1:64]*500/64./1000+4.5
0106 figure
0107 plot(frequency,10*log10(d.antenna0.roach1.LL(tempIndex(1,2),:)),'r','Linewidth',3)
0108 hold all
0109 plot(frequency,10*log10(d.antenna0.roach1.load2(tempIndex(1,2),:)),'.-')
0110 for i=2:2
0111     plot(frequency,10*log10(d.antenna0.roach1.LL(tempIndex(i,2),:)),'g','Linewidth',3)
0112     plot(frequency,10*log10(d.antenna0.roach1.load2(tempIndex(i,2),:)),'.-')
0113 end
0114 title('Roach1','Fontsize',16);
0115 xlabel('Frequency [GHz]','Fontsize',14)
0116 ylabel('Power [dBm]','Fontsize',14)
0117 legend('Tsky@12K','Tload@12K','Tsky@30K','Tload@30K')
0118 grid on
0119 
0120 
0121 frequency = 5.5-[1:64]*500/64./1000
0122 figure
0123 plot(frequency,10*log10(d.antenna0.roach2.LL(tempIndex(1,2),:))-10*log10(d.antenna0.roach2.LL(tempIndex(i,2),:)),'r','Linewidth',3)
0124 hold all
0125 plot(frequency,10*log10(d.antenna0.roach2.load2(tempIndex(1,2),:))-10*log10(d.antenna0.roach2.load2(tempIndex(i,2),:)),'.-')
0126 
0127 title('Roach2- delta Power between 30K and 12K load','Fontsize',16);
0128 xlabel('Frequency [GHz]','Fontsize',14)
0129 ylabel('Power Difference [dB]','Fontsize',14)
0130 legend('Tsky','Tload')
0131 grid on
0132 
0133 frequency = [1:64]*500/64./1000+4.5
0134 figure
0135 plot(frequency,10*log10(d.antenna0.roach1.LL(tempIndex(1,2),:))-10*log10(d.antenna0.roach1.LL(tempIndex(i,2),:)),'r','Linewidth',3)
0136 hold all
0137 plot(frequency,10*log10(d.antenna0.roach1.load2(tempIndex(1,2),:))-10*log10(d.antenna0.roach1.load2(tempIndex(i,2),:)),'.-')
0138 
0139 title('Roach1- delta Power between 30K and 12K load','Fontsize',16);
0140 xlabel('Frequency [GHz]','Fontsize',14)
0141 ylabel('Power Difference [dB]','Fontsize',14)
0142 legend('Tsky','Tload')
0143 grid on
0144 
0145 
0146