longTermPowerSpectraSpecific

 this is a script to plot out data over a long period of time from the
 Southern System
CJC 10 January 2013

0001 % this is a script to plot out data over a long period of time from the
0002 % Southern System
0003 %CJC 10 January 2013
0004 
0005 
0006 %%%%%%%
0007 ylimTemp=[11.995 12.005];
0008 ylimPow = [1e-5 1e2]
0009 ylimTimePow=[-200 200]
0010 scheduleString='12K Load Sky Stare'
0011 d=read_arcSouth('21-Apr-2013:13:15:30','21-Apr-2013:13:20:30');% eccosorb in place -
0012 d = read_arcSouth('29-Dec-2012:23:40:00','29-Dec-2012:23:55:00')
0013    
0014     [dataFreqRoach,dataTimeRoach,dataTimeTempLs,dataTimeTempCryo]=generatePowerSpectra(d,0);
0015     string=['29-Dec-2012:23:40:00->29-Dec-2012:23:55:00']
0016     figure
0017     subplot(414)
0018     loglog(dataFreqRoach(:,1),dataFreqRoach(:,4))
0019     hold all
0020     loglog(dataFreqRoach(:,1),dataFreqRoach(:,5))
0021     grid on
0022     title('Polarisation')
0023     subplot(413)
0024     plot(dataTimeRoach(:,1),dataTimeRoach(:,2))
0025     hold all
0026     plot(dataTimeRoach(:,1),dataTimeRoach(:,3))
0027     %legend('RCP','LCP')
0028     xlabel('Time [Hr]')
0029     ylabel('BU')
0030     titlestring=['Time Series: ',scheduleString]
0031     title(titlestring)
0032     ylim(ylimTimePow)
0033     grid on
0034     subplot(412)
0035     plot(dataTimeTempCryo(:,1),dataTimeTempCryo(:,2))
0036     xlabel('Time [Hr]')
0037     ylabel('Temperature [K]')
0038     title('Load/Sky Temperature')
0039     grid on
0040     ylim(ylimTemp)
0041     subplot(411)
0042     loglog(dataFreqRoach(:,1),dataFreqRoach(:,2))
0043     hold all
0044     loglog(dataFreqRoach(:,1),dataFreqRoach(:,3))
0045     grid on
0046     legend('RCP','LCP','Location','SouthWest')
0047     xlabel('Time [Hr]')
0048     ylabel('BU')
0049     titlestring=['Noise Power Spec: ',string,' : ',scheduleString]
0050     title(titlestring)
0051     ylim(ylimPow)
0052     filename=['29-Dec-2012:23:40:00',scheduleString];
0053     print('-depsc',filename)
0054     close all
0055