The problem with all the ones I discussed in the last commit (as well as with the A1 array) was that the k array was being overwritten in femsim.f90. I changed this and it now correctly reads in everything (I hope) (not just k). We get reasonable output now! I don't know if it's actually correct though.