use average k mode from binning in powerspectrum()

[jharno]

Ok, the pull request could not be done automatically. I have started manual merge, will need to test before applying to master branch.