converted References heading to rubric and only show for html

docs/source/methods/cmfd.rst

@@ -525,9 +525,9 @@ Once weight adjustment factors were calculated, each neutron's statistical
 weight in the source bank was modified according to its location and energy.
 Examples of CMFD simulations using OpenMC can be found in [HermanThesis]_.
 
-----------
-References
-----------
+.. only:: html
+
+   .. rubric:: References
 
 .. [BEAVRS] Nick Horelik, Bryan Herman. *Benchmark for Evaluation And Verification of Reactor
             Simulations*. Massachusetts Institute of Technology, http://crpg.mit.edu/pub/beavrs

docs/source/methods/cross_sections.rst

@@ -63,9 +63,9 @@ Other Methods
 A good survey of other energy grid techniques, including unionized energy grids,
 can be found in a paper by Leppanen_.
 
-----------
-References
-----------
+.. only:: html
+
+   .. rubric:: References
 
 .. [Brown] Forrest B. Brown, "New Hash-based Energy Lookup Algorithm for Monte
            Carlo codes," LA-UR-14-24530, Los Alamos National Laboratory (2014).

docs/source/methods/parallelization.rst

@@ -599,9 +599,9 @@ is actually independent of the number of nodes:
     E \left [ \Lambda_{j_{\text{max}}} \right ] = \sqrt{ \frac{N\sigma^2}{2\pi
     k^2}}.
 
-----------
-References
-----------
+.. only:: html
+
+   .. rubric:: References
 
 .. [Troubetzkoy] E. Troubetzkoy, H. Steinberg, and M. Kalos, "Monte Carlo
    Radiation Penetration Calculations on a Parallel Computer,"

docs/source/methods/physics.rst

@@ -1552,9 +1552,9 @@ default, the cutoff weight in OpenMC is :math:`w_c = 0.25` and the survival
 weight is :math:`w_s = 1.0`. These parameters vary from one Monte Carlo code to
 another.
 
-----------
-References
-----------
+.. only:: html
+
+   .. rubric:: References
 
 .. [Doyas] Richard J. Doyas and Sterrett T. Perkins, "Interpolation of Tabular
    Secondary Neutron and Photon Energy Distributions," *Nucl. Sci. Eng.*,

docs/source/methods/random_numbers.rst

@@ -62,9 +62,9 @@ Note that :eq:`lcg-skipahead` has the same general form as \eqref{eq:lcg}, so
 the idea is to determine the new multiplicative and additive constants in
 :math:`O(\log_2 N)` operations.
 
-----------
-References
-----------
+.. only:: html
+
+   .. rubric:: References
 
 .. [LEcuyer] P. L’Ecuyer, "Tables of Linear Congruential Generators of
    Different Sizes and Good Lattice Structures," *Math. Comput.*, **68**, 249

docs/source/methods/tallies.rst

@@ -475,9 +475,9 @@ normal distribution, we use an `unpublished rational approximation`_. After
 using the rational approximation, one iteration of Newton's method is applied to
 improve the estimate of the percentile.
 
-----------
-References
-----------
+.. only:: html
+
+   .. rubric:: References
 
 .. [George] E. E. Olusegun George and Meenakshi Sivaram, "A modification of the
    Fisher-Cornish approximation for the student t percentiles," Communication