That loop is creating a load of empty dictionaries, then your second loop is overwriting those dictionaries with new ones that have your desired data already contained within.

Combining the two modules (nestedWilsonLib and spin_resolved_analysis), one can compute the spin-resolved Wilson loop and the nested spin-resolved Wilson loop spectrum for a given PythTB tight-binding ...