The Role of Vibrational Anharmonicity in the Computational Study of Thermal Spin Crossover