def read_SCF(g09_file):
    """
    Opens a g09 log file and returns the final SCF energy

    Parameters
    ----------
    g09_file: Path to g09 log file
        File path

    Returns
    ----------
    SCF: float
        Final SCF energy

    """
    energies=[]
    with open(g09_file) as f:
        for line in f:
            # Ensures line is not blank
            if line.strip():
                if " SCF Done:" in line:
                    energies.append(float(line.split()[4]))
        f.close()
    SCF=energies[-1]
    return SCF