A method for calculation of finite temperature thermodynamic properties for magnetic allotropes has been developed. The developed method is based on quasi-harmonic approximations using the Debye–Grüneisen model. Electronic contributions are included and magnetic disordering contributions are added with the empiric model often used in CALPHAD (CALculation of PHAse Diagrams) databases. The developed method also includes two new ways to estimate the Debye temperature (θD(V0)) which enables calculation of finite temperature thermodynamic properties also for allotropes that are dynamically unstable at 0 K. Also a new CALPHAD approach for estimation of Curie- (TC) or Néel- (TN) temperatures is presented. The developed method was used to calculate Gibbs energy as a function of temperature (G(T)) and isobaric heat capacity as a function of temperature (CP(T)) for the body-centred-cubic (bcc), the face-centred-cubic (fcc) and the hexagonal-closed-packed (hcp) allotropes of Fe, Co and Ni. The results were compared with the SGTE Unary database. By using the developed method the phase transitions in Fe and Co were predicted. The calculations of the metastable allotropes fcc-Fe, hcp-Fe and bcc-Ni in this work indicates that the descriptions of these allotropes in the SGTE Unary database need to be updated, especially regarding the magnetic parameters. Furthermore, contradictory to previously reported results in literature, the calculations in this work show that the magnetic ground state structure of fcc-Fe is the double layer antiferromagnetic structure (AFM-D) and not a spin-spiral. The final conclusion is that by using the developed method it is possible to calculate finite temperature thermodynamic properties, not only, for stable and metastable magnetic allotropes, but also for magnetic allotropes that are dynamically unstable at 0 K, which is very important for the development of thermodynamic databases.