longTermTimeSeriesFreq

 script to plot data taken over a time series on the Southern RX

0001 
0002 % script to plot data taken over a time series on the Southern RX
0003 
0004 %%% CJC 27 April 2013 %%%%%
0005 
0006 
0007 % takes the following inputs
0008 %timeStart='26-Apr-2013:17:00:49' ; % time start in C-BASS string
0009 %timeEnd = '27-Apr-2013:10:00:49' ; % time end in C-BASS string
0010 
0011 timeStart='27-Apr-2013:00:00:49' ; % time start in C-BASS string
0012 timeEnd = '27-Apr-2013:15:00:49' ; % time end in C-BASS string
0013 
0014 
0015 scanLength = 0.5 ; % scan length in hours
0016 
0017 
0018 %%%%%%%%%%%%%%%
0019 dateStart = datenum(timeStart,'dd-mmm-yyyy:HH:MM:SS')
0020 dateEnd = datenum(timeEnd,'dd-mmm-yyyy:HH:MM:SS')
0021 
0022 currentTimeStart = dateStart;
0023 
0024 while(currentTimeStart <=dateEnd)
0025     
0026     currentTimeEnd = currentTimeStart+scanLength/24; %calculate end time
0027     
0028     currentTimeStartString = datestr(currentTimeStart,'dd-mmm-yyyy:HH:MM:SS') % convert to C-BASS string
0029     currentTimeEndString = datestr(currentTimeEnd,'dd-mmm-yyyy:HH:MM:SS')% convert to C-BASS string
0030     
0031     % read in the data
0032     d=read_arcSouth(currentTimeStartString,currentTimeEndString) 
0033     
0034     %convert to matlab times so we can format the axis
0035     time=d.antenna0.roach1.utc; %roach timestamps
0036     timeTemps=d.antenna0.thermal.utc; %temperature timestamps
0037     [year, month, day, hour, minute, second] = mjd2date(time);
0038     dateRoach=[year, month, day, hour, minute, second];
0039     [year, month, day, hour, minute, second] = mjd2date(timeTemps);
0040     dateTemps=[year, month, day, hour, minute, second];
0041     Roachtime=datenum(dateRoach);
0042     TempTime=datenum(dateTemps);
0043     
0044     
0045     % now plot the data making sure that we will just keep incrementing the
0046     % plot with time....
0047     
0048     subplot(311)
0049     plot(Roachtime,d.antenna0.roach1.LLfreq(:,:))
0050     hold all
0051     plot(Roachtime,d.antenna0.roach1.load2freq(:,:))
0052     title('Intensity','Fontsize',16);
0053     xlabel('Time','Fontsize',14)
0054     ylabel('Backend Units','Fontsize',14)
0055     datetick('x','yyyy-mm-dd HH:MM:SS','keeplimits')
0056     h=legend('LCP','LoadL')
0057     grid on
0058     
0059     subplot(312)
0060     plot(Roachtime,d.antenna0.roach1.Q(:,30))
0061     hold all
0062     plot(Roachtime,d.antenna0.roach1.U(:,30))
0063     title('Polarization','Fontsize',16);
0064     xlabel('Time','Fontsize',14)
0065     ylabel('Backend Units','Fontsize',14)
0066     datetick('x','yyyy-mm-dd HH:MM:SS','keeplimits')
0067     h=legend('Q','U')
0068     grid on
0069     
0070     subplot(313)
0071     plot(TempTime,d.antenna0.thermal.ccTemperatureLoad)
0072     hold all
0073     title('Cold Load','Fontsize',16);
0074     xlabel('Time','Fontsize',14)
0075     ylabel('Temperature [K]','Fontsize',14)
0076     datetick('x','yyyy-mm-dd HH:MM:SS','keeplimits')
0077     grid on
0078     
0079     %increment start time by scan length so that the next time will be at a
0080     %later time....
0081     currentTimeStart = currentTimeStart+scanLength/24; 
0082 
0083 end
0084 
0085 
0086 
0087 
0088 
0089