fixed bug computing heat capacities

Author: AlbertoCuadra

Solver/Functions/ComputeProperties.py

@@ -39,8 +39,6 @@ def __init__(self, app, SpeciesMatrix, p, T):
         self.ef = sum(SpeciesMatrix[:, 3])  # [kJ]
         self.ehT = sum(SpeciesMatrix[:, 4]) # [kJ]
         self.e = self.ef + self.ehT         # [kJ]
-        self.cP = sum(SpeciesMatrix[:, 5])  # [J/K]
-        self.cV = sum(SpeciesMatrix[:, 6])  # [J/K]
         self.S0 = sum(SpeciesMatrix[:, 7])  # [kJ/K]
         self.pv = sum(SpeciesMatrix[:, 8])  # [J/K]
         self.p = p                          # [bar]
@@ -50,6 +48,8 @@ def __init__(self, app, SpeciesMatrix, p, T):
         self.mi = sum(SpeciesMatrix[:, 10]) * 1e-3            # [kg]
         self.rho = self.mi / self.v * 1e3                     # [kg/m3]
         self.Yi = SpeciesMatrix[:, 10] / self.mi * 1e-3       # [-]
+        self.cP = sum(SpeciesMatrix[:, 5] * self.Yi)  # [J/K]
+        self.cV = sum(SpeciesMatrix[:, 6] * self.Yi)  # [J/K]
         self.W = 1/np.nansum(self.Yi / SpeciesMatrix[:, 11])  # []
         Ni = SpeciesMatrix[:, 0]                              # [mol]
         self.Xi = Ni / self.N                                 # [-]