longTermTimeSeries

 script to plot data taken over a time series on the Southern RX 
 inputs - start_date, stop_date

0001 function longTermTimeSeries(timeStart,timeEnd)
0002 % script to plot data taken over a time series on the Southern RX
0003 % inputs - start_date, stop_date
0004 
0005 %%% CJC 27 April 2013 %%%%%
0006 %%% ACT 28 April 2013 - changed to have separate plots including lakeshore
0007 %%% temp sensors. Also uses LL averaged over frequency
0008 %%% ACt 29 April 2013 - takes start and stop date on the command line
0009 
0010 % takes the following inputs
0011 %timeStart='26-Apr-2013:17:00:49' ; % time start in C-BASS string
0012 %timeEnd = '27-Apr-2013:10:00:49' ; % time end in C-BASS string
0013 
0014 %timeStart='26-Apr-2013:22:00:49' ; % time start in C-BASS string%
0015 %timeEnd = '29-Apr-2013:05:00:49' ; % time end in C-BASS string
0016 
0017 
0018 scanLength = 0.5 ; % scan length in hours
0019 
0020 
0021 %%%%%%%%%%%%%%%
0022 dateStart = datenum(timeStart,'dd-mmm-yyyy:HH:MM:SS')
0023 dateEnd = datenum(timeEnd,'dd-mmm-yyyy:HH:MM:SS')
0024 
0025 currentTimeStart = dateStart;
0026 
0027 while(currentTimeStart <=dateEnd)
0028     
0029     currentTimeEnd = currentTimeStart+scanLength/24; %calculate end time
0030     
0031     currentTimeStartString = datestr(currentTimeStart,'dd-mmm-yyyy:HH:MM:SS') % convert to C-BASS string
0032     currentTimeEndString = datestr(currentTimeEnd,'dd-mmm-yyyy:HH:MM:SS')% convert to C-BASS string
0033     
0034     % read in the data
0035     d=read_arcSouth(currentTimeStartString,currentTimeEndString) 
0036     
0037     %convert to matlab times so we can format the axis
0038     time=d.antenna0.roach1.utc; %roach timestamps
0039     timeTemps=d.antenna0.thermal.utc; %temperature timestamps
0040     [year, month, day, hour, minute, second] = mjd2date(time);
0041     dateRoach=[year, month, day, hour, minute, second];
0042     [year, month, day, hour, minute, second] = mjd2date(timeTemps);
0043     dateTemps=[year, month, day, hour, minute, second];
0044     Roachtime=datenum(dateRoach);
0045     TempTime=datenum(dateTemps);
0046     
0047     
0048     % now plot the data making sure that we will just keep incrementing the
0049     % plot with time....
0050     
0051 %    subplot(411)
0052 chan = 46
0053 figure(1)
0054     plot(Roachtime,d.antenna0.roach1.LL(:,chan))
0055     hold all
0056     plot(Roachtime,d.antenna0.roach1.load2(:,chan))
0057     title('Intensity','Fontsize',16);
0058     xlabel('Time','Fontsize',14)
0059     ylabel('Backend Units','Fontsize',14)
0060     datetick('x','yyyy-mm-dd HH:MM:SS','keeplimits')
0061     h=legend('LCP','LoadL')
0062     grid on
0063     
0064 %    subplot(412)
0065 figure(20)
0066     plot(Roachtime,d.antenna0.roach1.Q(:,30))
0067     hold all
0068     plot(Roachtime,d.antenna0.roach1.U(:,30))
0069     title('Polarization','Fontsize',16);
0070     xlabel('Time','Fontsize',14)
0071     ylabel('Backend Units','Fontsize',14)
0072     datetick('x','yyyy-mm-dd HH:MM:SS','keeplimits')
0073     h=legend('Q','U')
0074     grid on
0075     
0076    % subplot(413)
0077 figure(30)   
0078  plot(TempTime,d.antenna0.thermal.ccTemperatureLoad)
0079     hold all
0080     title('Cold Load','Fontsize',16);
0081     xlabel('Time','Fontsize',14)
0082     ylabel('Temperature [K]','Fontsize',14)
0083     datetick('x','yyyy-mm-dd HH:MM:SS','keeplimits')
0084     grid on
0085     
0086 
0087     %subplot(414)
0088 figure(40)    
0089 plot(Roachtime(1:100:end),d.antenna0.thermal.lsTemperatureSensors)
0090     hold all
0091     title('Lakeshore','Fontsize',16);
0092     xlabel('Time','Fontsize',14)
0093     ylabel('Temperature [K]','Fontsize',14)
0094     datetick('x','yyyy-mm-dd HH:MM:SS','keeplimits')
0095     grid on
0096 
0097 
0098 
0099     %increment start time by scan length so that the next time will be at a
0100     %later time....
0101     currentTimeStart = currentTimeStart+scanLength/24; 
0102 
0103 end
0104 
0105 
0106 
0107 
0108 
0109