for i in range(number_of_bands // 2):
                        energies[i, int(kpoint_counter)] = float(words[5 + 2 * i])
                        occupancy[i, int(kpoint_counter)] = float(words[4 + 2 * i])