%
%   WPI_MICROS_2012.M    Matlab Script File with Homogeneized Densities and 
%                   Microscopic Cross Sections for the WPI LEU Fuel Assembly
%
%   This file contains 2-group microscopic data and assembly average nuclide 
%   densities for the WPI LEU fuel assembly.  The format of the variables 
%   is defined in the comments within the script file.  This file can be used
%   within other Matlab codes that need these data for various computations for 
%   the UMLRR with the WPI fuel assembly.
%
%   Data Ref: 1. micro xsec data from wpi2gvxs_v7n238.bcdlib (generated Aug. 2012)
%             2. assembly average densities from umlrrden.m and wpifdata.m
%
%   File prepared by J. R. White, UMass-Lowell (last update:  November 2016)
%
  
%
%   define the micro xsecs for 5 isotopes (from wpi2gvxs_v7n238.bcdlib)
%   units:  all cross sections have units of barns (where 1 b = 1e-24 cm^2)
%   format:  [ tr1       tr2;   fis1   fis2;   abs1      abs2; 
%             nufis1  nufis2;   tot1   tot2;   scat11  scat22;
%               0     scat12]
%   U235   WPI LEU fuel assy                    
  U235 = [1.28524E+01 3.79064E+02; 8.10212E+00 3.83641E+02; 1.24539E+01 4.49357E+02;
          1.98614E+01 9.34819E+02; 2.12437E+01 4.63669E+02; 8.79051E+00 1.43125E+01;
          0.00000E+00 5.48160E-03];
%   U238   full LEU fuel assy                    
  U238 = [9.37397E+00 1.05592E+01; 1.34232E-01 1.16138E-05; 3.14884E+00 1.86039E+00;
          3.74304E-01 2.89426E-05; 1.50334E+01 1.08353E+01; 1.18871E+01 8.97487E+00;
          0.00000E+00 3.79322E-03];
%   H1     full LEU fuel assy                    
    H1 = [2.72080E+00 2.51899E+01; 0.00000E+00 0.00000E+00; 3.46903E-03 2.45124E-01;
          0.00000E+00 0.00000E+00; 1.10472E+01 4.17276E+01; 1.04025E+01 4.14824E+01;
          0.00000E+00 6.40995E-01];
%   O-16   full LEU fuel assy                    
   O16 = [2.62927E+00 3.82666E+00; 0.00000E+00 0.00000E+00; 3.29889E-03 1.40274E-04;
          0.00000E+00 0.00000E+00; 3.36376E+00 4.00749E+00; 3.34533E+00 4.00735E+00;
          0.00000E+00 1.51274E-02];
%   Al27   full LEU fuel assy                    
  Al27 = [2.13458E+00 1.51099E+00; 0.00000E+00 0.00000E+00; 5.09056E-03 1.62161E-01;
          0.00000E+00 0.00000E+00; 2.69290E+00 1.56598E+00; 2.68437E+00 1.40382E+00;
          0.00000E+00 3.43263E-03];
%
%   now define the assembly averaged densities (data from umlrrden.m & leufdata.m)
%   units: N has units of atoms/b-cm  (multiply by 1e24 to convert to atoms/cm^3)
%   nuclide order: U235, U238, H1, O16, Al27
   nam = ['U235';'U238';'H   ';'O16 ';'Al27'];  nam = char(nam);
   den = [1.1866e-004 4.7606e-004 3.8788e-002 1.9394e-002 2.3371e-002];
%
%   end of data