plotAppliedCalibrationValues

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

  function plotAppliedNoiseDiodeTemps(start, stop, deltaT)

  This just calls the code in the pipeline that determines the noise diode
  temperatures.  
  
   start:  start time of your plot  [mjd]
   stop:   stop time of your plot   [mjd]
   deltaT:  interval on which you plot (default is 1 [day])

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

0001 function plotAppliedCalibrationValues(start, stop, deltaT)
0002 
0003 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0004 %
0005 %  function plotAppliedNoiseDiodeTemps(start, stop, deltaT)
0006 %
0007 %  This just calls the code in the pipeline that determines the noise diode
0008 %  temperatures.
0009 %
0010 %   start:  start time of your plot  [mjd]
0011 %   stop:   stop time of your plot   [mjd]
0012 %   deltaT:  interval on which you plot (default is 1 [day])
0013 %
0014 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0015 
0016 if(nargin==0)
0017   start = date2mjd(2014, 12, 0);
0018   stop  = date2mjd(2015, 01, 0);
0019   deltaT = 1;
0020 end
0021 
0022 if(nargin==2)
0023   deltaT = 1;
0024 end
0025 
0026 times = start:deltaT:stop;
0027 for m=1:length(times)
0028   d.array.frame.utc = times(m);
0029   d.antenna0.receiver.utc = times(m);
0030   d.par.autosave.flag = 1;
0031   [nv, nI(m,:), nIerr(m,:)] = getNoiseDiodeTemps(d);
0032   d2 = interpolateTau(d);
0033   nT(m,:) = d2.antenna0.receiver.tau;
0034 end
0035 
0036 figure()
0037 plot(times, nI);
0038 xlabel('Date [mjd]');
0039 ylabel('Noise diode temp applied [K]');
0040 legend('Diode 1', 'Diode 2');
0041 
0042 figure()
0043 plot(times, nT);
0044 xlabel('Date [mjd]');
0045 ylabel('Opacity Applied');
0046 
0047 return
0048