def get_nat_orb_occ(lines):
    """
    Find the natural orbital occupations
    """
    nat_occ = []
    start = False
    for line in lines:
        if start:
            if line.strip():
                nat_occ.append(abs(float(line.split('=')[-1].strip())))
            else:
                break
        elif line == 'Natural Orbital Occupation Numbers:\n':
            start = True

    return nat_occ