One issue causing the jagged spectra at high wavenumbers is under sampling there. For example consider the picture:

[![enter image description here][1]][1]

You don't want to sample from the red zones as they will become increasingly under-sampled as you move past a radius of size $|k_{x}|$.


  [1]: https://i.sstatic.net/zlCnN.png