Convert optdone from bool to list

These changes allows Gaussian calculations with multiple geometry optimizations (e.g. Opt=ModRedundant and possibly IRCs) and ORCA relaxed surface scans to be completely parsed.

The optdone list stores the index of atomcoords, scfenergies, and geovalues when an individual optimization has converged.