function w_hat = qroj_hat(u_hat)
% QROJ_HAT projects onto negative wavenumbers
n = length(u_hat);
w_hat = [zeros(1,(n/2)) 0.5*u_hat(n/2+1) u_hat((n/2+2):n)]; 
end