diff --git a/doc/sphinx/cython/zerodim.rst b/doc/sphinx/cython/zerodim.rst
index 301041160a3..6389f09ff02 100644
--- a/doc/sphinx/cython/zerodim.rst
+++ b/doc/sphinx/cython/zerodim.rst
@@ -42,6 +42,10 @@ Reactor
 ^^^^^^^
 .. autoclass:: Reactor
 
+MoleReactor
+^^^^^^^^^^^
+.. autoclass:: MoleReactor(contents=None, *, name=None, energy='on')
+
 IdealGasReactor
 ^^^^^^^^^^^^^^^
 .. autoclass:: IdealGasReactor(contents=None, *, name=None, energy='on')
@@ -54,6 +58,10 @@ ConstPressureReactor
 ^^^^^^^^^^^^^^^^^^^^
 .. autoclass:: ConstPressureReactor(contents=None, *, name=None, energy='on')
 
+ConstPressureMoleReactor
+^^^^^^^^^^^^^^^^^^^^^^^^
+.. autoclass:: ConstPressureMoleReactor(contents=None, *, name=None, energy='on')
+
 IdealGasConstPressureReactor
 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 .. autoclass:: IdealGasConstPressureReactor(contents=None, *, name=None, energy='on')
diff --git a/include/cantera/zeroD/MoleReactor.h b/include/cantera/zeroD/MoleReactor.h
index d42cf6b6986..4f33a495f13 100644
--- a/include/cantera/zeroD/MoleReactor.h
+++ b/include/cantera/zeroD/MoleReactor.h
@@ -38,11 +38,11 @@ class MoleReactor : public Reactor
 
 protected:
     //! Get moles of the system from mass fractions stored by thermo object
-    //! @param y[out] vector for system moles to be put into
+    //! @param y vector for moles to be put into
     virtual void getMoles(double* y);
 
     //! Set internal mass variable based on moles given
-    //! @param y[in] vector of moles of the system
+    //! @param y vector of moles of the system
     virtual void setMassFromMoles(double* y);
 
     virtual void evalSurfaces(double* LHS, double* RHS, double* sdot);
diff --git a/test/python/test_reactor.py b/test/python/test_reactor.py
index 8ddd133e2f0..89c08ba9ab0 100644
--- a/test/python/test_reactor.py
+++ b/test/python/test_reactor.py
@@ -220,10 +220,9 @@ def integrate(atol, rtol):
 
         n_baseline = integrate(1e-10, 1e-20)
         n_rtol = integrate(5e-7, 1e-20)
-        n_atol = integrate(1e-10, 1e-6)
-
-        self.assertTrue(n_baseline > n_rtol)
-        self.assertTrue(n_baseline > n_atol)
+        n_atol = integrate(1e-10, 1e-5)
+        self.assertTrue(n_baseline > n_rtol, f"{n_baseline} is less than {n_rtol}")
+        self.assertTrue(n_baseline > n_atol, f"{n_baseline} is less than {n_atol}")
 
     def test_advance_limits(self):
         P0 = 10 * ct.one_atm