`read_abinit_WFK` does not store k-points

[brainandforce]

I've noticed that when read_abinit_WFK is called, the k-points are simply missing from the data, all of them being the zero vector with weight 1. There are some issues with the way PlanewaveWavefunction is implemented, and the fix will involve some refactoring of that code.