def Cp(self, phase, T):
if phase not in self._phases:<EOL><INDENT>raise Exception("<STR_LIT>" %<EOL>(phase, self.formula))<EOL><DEDENT>return self._phases[phase].Cp(T)<EOL>
Calculate the heat capacity of a phase of the compound at a specified
temperature.

:param phase: A phase of the compound, e.g. 'S', 'L', 'G'.
:param T: [K] temperature

:returns: [J/mol/K] Heat capacity.