diff --git a/bmad/doc/cover-page.tex b/bmad/doc/cover-page.tex
index 276100828a..cb2d88b02f 100644
--- a/bmad/doc/cover-page.tex
+++ b/bmad/doc/cover-page.tex
@@ -3,7 +3,7 @@
 
 \begin{flushright}
 \large
-  Revision: July 7, 2024 \\
+  Revision: July 10, 2024 \\
 \end{flushright}
 
 \pdfbookmark[0]{Preamble}{Preamble} 
diff --git a/bmad/doc/elements.tex b/bmad/doc/elements.tex
index 8864443d22..54da59077f 100644
--- a/bmad/doc/elements.tex
+++ b/bmad/doc/elements.tex
@@ -1266,7 +1266,7 @@ \section{Crystal}
 \index{follow_diffracted_beam}\index{thickness}
 Attributes specific to a \vn{crystal} element are:
 \begin{example}
-  b_param            = <Real>       ! b parameter
+  b_param            = <Real>       ! b parameter for photons with the reference energy.
   crystal_type       = <String>     ! Crystal material (\sref{s:cryst.list}) and reflection plane.
   psi_angle          = <Real>       ! Rotation of H-vector about the surface normal.
   thickness          = <Real>       ! Thickness of crystal for Laue diffraction.
diff --git a/bmad/doc/synchrotron-radiation.tex b/bmad/doc/synchrotron-radiation.tex
index 933ec7ca60..0c4b74fe1b 100644
--- a/bmad/doc/synchrotron-radiation.tex
+++ b/bmad/doc/synchrotron-radiation.tex
@@ -247,9 +247,9 @@ \section{Synchrotron Radiation Integrals}
 \label{s:synch.ints}
 \index{synchrotron radiation!integrals}
 
-the synchrotron radiation integrals can be used to compute emittances, the energy spread,
+The synchrotron radiation integrals can be used to compute emittances, the energy spread,
 etc. However, using the 6D damped and stochastic transport matrices (\sref{s:map.rad}) has
-a number of advantages
+a number of advantages:
 \begin{itemize}
 %
 \item
diff --git a/regression_tests/nonlin_test/small_ring.bmad b/regression_tests/nonlin_test/small_ring.bmad
index 8bb24d6fd6..2f3e3ab54e 100644
--- a/regression_tests/nonlin_test/small_ring.bmad
+++ b/regression_tests/nonlin_test/small_ring.bmad
@@ -1,6 +1,7 @@
 parameter[p0c] = 1e9
 bmad_com[radiation_damping_on] = f
 bmad_com[radiation_fluctuations_on] = f
+bmad_com[spin_tracking_on] = T
 
 q1: quad, l = 0.5, k1 = 0.6, tilt = 0.05
 q2: quad, l = 0.5, k1 = -0.5