diff --git a/src/thermophysicalModels/reactionThermo/mixtures/basicSpecieMixture/basicSpecieMixture.C b/src/thermophysicalModels/reactionThermo/mixtures/basicSpecieMixture/basicSpecieMixture.C
index 0ff307384..64af9f76c 100644
--- a/src/thermophysicalModels/reactionThermo/mixtures/basicSpecieMixture/basicSpecieMixture.C
+++ b/src/thermophysicalModels/reactionThermo/mixtures/basicSpecieMixture/basicSpecieMixture.C
@@ -96,15 +96,19 @@ Foam::tmp<Foam::volScalarField> Foam::basicSpecieMixture::Qc() const
                 Y[0].mesh()
             ),
             Y[0].mesh(),
-            dimensionedScalar("zero", dimless, 0)
+            dimensionedScalar("zero", dimCurrent*dimTime/dimMass, 0)
         )
     );
 
     volScalarField& Qc = tQc.ref();
 
+    dimensionedScalar Qci("Qci", dimCurrent*dimTime/dimMass, 0);
+
     forAll(Y, i)
     {
-        Qc += Y[i]*this->Qc(i);
+        Qci.value() = this->Qc(i);
+        Qc += Y[i]*Qci;
+        // Qc += Y[i]*this->Qc(i);
     }
 
     return tQc;
@@ -126,15 +130,19 @@ Foam::tmp<Foam::volScalarField> Foam::basicSpecieMixture::Qc2() const
                 Y[0].mesh()
             ),
             Y[0].mesh(),
-            dimensionedScalar("zero", dimless, 0)
+            dimensionedScalar("zero", sqr(dimCurrent*dimTime)/dimMass, 0)
         )
     );
 
     volScalarField& Qc = tQc.ref();
 
+    dimensionedScalar Qc2i("Qc2i", sqr(dimCurrent*dimTime)/dimMass, 0);
+
     forAll(Y, i)
     {
-        Qc += Y[i]*this->Qc2(i);
+        Qc2i.value() = this->Qc2(i);
+        Qc += Y[i]*Qc2i;
+        //Qc += Y[i]*this->Qc2(i);
     }
 
     return tQc;