From 5c9dbef47b5c8bb62ebc98e63e091805a122a990 Mon Sep 17 00:00:00 2001
From: ignis <ignis@postech.ac.kr>
Date: Sat, 28 Nov 2015 18:16:35 +0900
Subject: [PATCH] class template heThermo added init_qc() for qc_

---
 .../basic/heThermo/heThermo.C                     | 15 +++++++++++++++
 .../basic/heThermo/heThermo.H                     |  1 +
 2 files changed, 16 insertions(+)

diff --git a/src/thermophysicalModels/basic/heThermo/heThermo.C b/src/thermophysicalModels/basic/heThermo/heThermo.C
index 68fed6b9..5b7462e4 100644
--- a/src/thermophysicalModels/basic/heThermo/heThermo.C
+++ b/src/thermophysicalModels/basic/heThermo/heThermo.C
@@ -75,6 +75,21 @@ void Foam::heThermo<BasicThermo, MixtureType>::init()
     }
 
     this->heBoundaryCorrection(he_);
+
+    init_qc();
+}
+
+
+template<class BasicThermo, class MixtureType>
+void Foam::heThermo<BasicThermo, MixtureType>::init_qc()
+{
+    scalarField& qcCells = qc_.internalField();
+
+    forAll(qcCells, celli)
+    {
+        qcCells[celli] =
+            this->cellMixture(celli).QC();
+    }
 }
 
 
diff --git a/src/thermophysicalModels/basic/heThermo/heThermo.H b/src/thermophysicalModels/basic/heThermo/heThermo.H
index 68836137..2ae99b52 100644
--- a/src/thermophysicalModels/basic/heThermo/heThermo.H
+++ b/src/thermophysicalModels/basic/heThermo/heThermo.H
@@ -79,6 +79,7 @@ private:
 
         //- Initialize heThermo
         void init();
+        void init_qc();
 
 
 public: