def Cp_mag(self, T):
tau = T / self.Tc_mag<EOL>if tau <= <NUM_LIT:1.0>:<EOL><INDENT>c = (self._B_mag*(<NUM_LIT:2>*tau**<NUM_LIT:3> + <NUM_LIT:2>*tau**<NUM_LIT:9>/<NUM_LIT:3> + <NUM_LIT:2>*tau**<NUM_LIT:15>/<NUM_LIT:5>))/self._D_mag<EOL><DEDENT>else:<EOL><INDENT>c = (<NUM_LIT:2>*tau**-<NUM_LIT:5> + <NUM_LIT:2>*tau**-<NUM_LIT:15>/<NUM_LIT:3> + <NUM_LIT:2>*tau**-<NUM_LIT>/<NUM_LIT:5>)/self._D_mag<EOL><DEDENT>result = R*math.log(self.beta0_mag + <NUM_LIT:1>)*c<EOL>return result<EOL>
Calculate the phase's magnetic contribution to heat capacity at the
specified temperature.

:param T: [K] temperature

:returns: [J/mol/K] The magnetic heat capacity of the compound phase.

Dinsdale, A. T. (1991). SGTE data for pure elements. Calphad, 15(4),
317–425. http://doi.org/10.1016/0364-5916(91)90030-N