import cclib data = cclib.io.ccopen('output.log').parse() energy = data.scfenergies[-1] total_charge = sum(data.atomcharges) for atom1, atom2, atom3, angle in data.geometry: cclib.plot.energy_level_diagram(data.moenergies, filename='energy_levels.png')