function m = mom(u, k)
% MOM calculates the momentum
  m = real(ifft(abs(fft(u.')).^2)./k);
end