I have optimized a set of Fe complexes and I want to compare computed bond lengths and angles to theoretical ones. I have started doing this manually but I was just wondering if there is a way to do this automatically as this is a boring task. I would appreciate it if anyone has a software or python script that can do this automatically.