Calculate quasi-harmonic free energies from Gaussian output files with temperature and other corrections