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

:param T: [K] temperature

:returns: [J/mol] The magnetic enthalpy 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