* The implicit assumption that the gamma point is the first k point was eliminated.
* The opening of the kpts file to determine the number of k points in the full BZ
  was eliminated (It is not yet tested if the new approach always works).
* Some variables have been removed, some have been renamed.
* Some debugging code has been removed.
* An error message was improved.