We can determine the angle between two subsequent slices by using arctan2 as defined in NumPy . We use this to measure the angular propagation of the MO nodal plane through the molecule. By taking the cumulated sum of all the angles in C, we can compare the helicality of MOs in different molecules quantitatively.