update the values of functions related to temperature

Author: Xsmos

ares/static/ChemicalNetwork.py

@@ -203,6 +203,7 @@ def RateEquations(self, t, q, args):
         if not self.isothermal:
 
             for i, sp in enumerate(self.neutrals):
+                self.Beta, self.alpha, self.zeta, self.eta, self.psi, self.xi, self.omega = self.SourceIndependentCoefficients(q[-1], z).values()
                 elem = self.grid.parents_by_ion[sp]
 
                 heat += k_heat[i] * x[sp] * n[elem]          # photo-heating