diff --git a/ares/static/ChemicalNetwork.py b/ares/static/ChemicalNetwork.py
index 81e0eb706..1081670dd 100755
--- a/ares/static/ChemicalNetwork.py
+++ b/ares/static/ChemicalNetwork.py
@@ -201,7 +201,7 @@ def RateEquations(self, t, q, args):
         heat = 0.0
         cool = 0.0
         if not self.isothermal:
-
+            self.Beta, self.alpha, self.zeta, self.eta, self.psi, self.xi, self.omega = self.SourceIndependentCoefficients(q[-1], z).values()
             for i, sp in enumerate(self.neutrals):
                 elem = self.grid.parents_by_ion[sp]