e_fermi is wrong for k-points and fixed occupations

The code clearly selects the value at Gamma, which is wrong in general.