%
%   TAVE_EXPT_022707.M    Experimental data for the average pool temperature in the 
%             UMLRR for a specific pool heatup experiment performed on Feb. 27, 2007
%
%   Note:  Experimental data are stored in the Te array at 1 minute intervals.  This 
%          contains 300 minutes of data while the secondary coolant loop is turned 
%          off.  During this time period the reactor was at about 93% of full power
%          according to the primary power indicators.  The initial pool temperature
%          was 67.8 F and the air temperature in containment was nearly constant at 
%          about 70.3 F during the full experiment.  However, all these measurements
%          are noisy, the thermocouple data are not really accurate to 0.1 degree,
%          and there are often some minor losses due to convection heat transfer from
%          within the primary loop.  Thus, a curve fit to help minimize some of this
%          uncertainty is used to define the efffective power and initial pool
%          temperature more accurately (this is done interactively with Matlab's
%          Basic Fitting Tool and the results are indicated at the end of this file).
%
%   File generated by J. R. White, UMass-Lowell (Feb. 27, 2007)
%

%  
%   raw temperture data from UMLRR Online averaged over 1 minute intervals
      Te = [6.78E+01 6.78E+01 6.82E+01 6.81E+01 6.82E+01 6.83E+01 6.84E+01 ...
            6.84E+01 6.85E+01 6.87E+01 6.88E+01 6.89E+01 6.89E+01 6.91E+01 ...
            6.92E+01 6.92E+01 6.93E+01 6.93E+01 6.95E+01 6.95E+01 6.96E+01 ...
            6.96E+01 6.96E+01 6.98E+01 6.99E+01 7.00E+01 7.00E+01 7.01E+01 ...
            7.03E+01 7.04E+01 7.03E+01 7.04E+01 7.04E+01 7.05E+01 7.04E+01 ...
            7.05E+01 7.08E+01 7.09E+01 7.10E+01 7.10E+01 7.10E+01 7.11E+01 ...
            7.12E+01 7.13E+01 7.13E+01 7.16E+01 7.16E+01 7.18E+01 7.18E+01 ...
            7.19E+01 7.20E+01 7.20E+01 7.21E+01 7.22E+01 7.23E+01 7.24E+01 ...
            7.24E+01 7.24E+01 7.26E+01 7.27E+01 7.27E+01 7.28E+01 7.28E+01 ...
            7.30E+01 7.31E+01 7.31E+01 7.31E+01 7.32E+01 7.32E+01 7.33E+01 ...
            7.36E+01 7.36E+01 7.37E+01 7.37E+01 7.38E+01 7.39E+01 7.39E+01 ...
            7.39E+01 7.41E+01 7.42E+01 7.42E+01 7.43E+01 7.44E+01 7.44E+01 ...
            7.45E+01 7.47E+01 7.47E+01 7.48E+01 7.49E+01 7.49E+01 7.51E+01 ...
            7.51E+01 7.52E+01 7.53E+01 7.53E+01 7.54E+01 7.54E+01 7.55E+01 ...
            7.56E+01 7.56E+01 7.57E+01 7.57E+01 7.58E+01 7.59E+01 7.61E+01 ...
            7.61E+01 7.62E+01 7.62E+01 7.64E+01 7.65E+01 7.66E+01 7.66E+01 ...
            7.67E+01 7.67E+01 7.68E+01 7.70E+01 7.70E+01 7.71E+01 7.71E+01 ...
            7.72E+01 7.73E+01 7.74E+01 7.75E+01 7.76E+01 7.76E+01 7.77E+01 ...
            7.78E+01 7.78E+01 7.79E+01 7.80E+01 7.80E+01 7.81E+01 7.81E+01 ...
            7.82E+01 7.83E+01 7.83E+01 7.84E+01 7.86E+01 7.86E+01 7.85E+01 ...
            7.86E+01 7.88E+01 7.88E+01 7.88E+01 7.89E+01 7.90E+01 7.91E+01 ...
            7.92E+01 7.93E+01 7.94E+01 7.94E+01 7.96E+01 7.97E+01 7.97E+01 ...
            7.98E+01 7.99E+01 8.00E+01 8.00E+01 8.00E+01 8.01E+01 8.01E+01 ...
            8.01E+01 8.02E+01 8.04E+01 8.05E+01 8.06E+01 8.06E+01 8.07E+01 ...
            8.08E+01 8.09E+01 8.10E+01 8.11E+01 8.12E+01 8.12E+01 8.13E+01 ...
            8.14E+01 8.15E+01 8.16E+01 8.17E+01 8.17E+01 8.17E+01 8.18E+01 ...
            8.19E+01 8.20E+01 8.20E+01 8.21E+01 8.22E+01 8.22E+01 8.22E+01 ...
            8.24E+01 8.26E+01 8.26E+01 8.26E+01 8.27E+01 8.28E+01 8.28E+01 ...
            8.30E+01 8.31E+01 8.31E+01 8.31E+01 8.32E+01 8.32E+01 8.33E+01 ...
            8.34E+01 8.36E+01 8.36E+01 8.37E+01 8.37E+01 8.38E+01 8.38E+01 ...
            8.39E+01 8.39E+01 8.40E+01 8.40E+01 8.41E+01 8.42E+01 8.43E+01 ...
            8.44E+01 8.45E+01 8.45E+01 8.46E+01 8.47E+01 8.47E+01 8.47E+01 ...
            8.49E+01 8.49E+01 8.50E+01 8.51E+01 8.52E+01 8.52E+01 8.54E+01 ...
            8.54E+01 8.54E+01 8.56E+01 8.55E+01 8.57E+01 8.57E+01 8.57E+01 ...
            8.59E+01 8.60E+01 8.60E+01 8.61E+01 8.62E+01 8.63E+01 8.63E+01 ...
            8.64E+01 8.64E+01 8.65E+01 8.66E+01 8.66E+01 8.68E+01 8.69E+01 ...
            8.69E+01 8.70E+01 8.71E+01 8.72E+01 8.72E+01 8.73E+01 8.74E+01 ...
            8.74E+01 8.75E+01 8.75E+01 8.76E+01 8.77E+01 8.77E+01 8.79E+01 ...
            8.79E+01 8.79E+01 8.80E+01 8.81E+01 8.82E+01 8.83E+01 8.83E+01 ...
            8.84E+01 8.84E+01 8.85E+01 8.86E+01 8.87E+01 8.87E+01 8.88E+01 ...
            8.89E+01 8.90E+01 8.90E+01 8.91E+01 8.92E+01 8.92E+01 8.92E+01 ...
            8.94E+01 8.95E+01 8.95E+01 8.95E+01 8.95E+01 8.96E+01 8.97E+01 ...
            8.98E+01 8.98E+01 9.00E+01 9.00E+01 9.01E+01 9.02E+01 9.03E+01];
      Nt = length(Te);           % length of temperature array
      te = 0:Nt-1;               % time vector - data taken at 1 minute intervals
      Tair = 70.3;               % air temperature (F) (this was nearly constant)
      edate = 'Feb. 27, 2007';   % date for experiment
%
%   Note: The effective core heatup rate and initial pool temperature were obtained
%   from a linear fit to the first 60 minutes of data. Assuming negligible losses
%   from the pool surface during this period, the pool temperature should increase
%   linearly according to T(t) = To + (P/mc)*t, where P/mc is the pool heatup rate
%   (with units of F/min).  From a linear fit and known approximate values for
%   the mass, m, and heat capacity, c, of the pool water, we can compute the
%   net power level of the core.  This value may be somewhat different from the 
%   value noted above, and it tries to take account of the small losses in the
%   primary loop that may occur and some of the uncertainty in the measured data -- 
%   assuming, of course, negligible losses from the pool surface during this time.
      To = 67.9;                 % estimated initial pool temp from curve fit (F)
      hurate = 0.0799;           % estimated pool heatup rate from curve fit (F/min)
%
%   end of data file
