This simple script performs a Boltzmann population from a .dat file containing two columns: the first column bears the IDs of each species, the second column bears the energies in Hartree (as from a gaussian calculation).
boltzmann.py file.dat 298.15 (you can put the temperature you want)