% function to calculate the inverse Fourier Transform of the frequency series Sw
function invfourtran = ift(Sw,dw)

  invfourtran = fftshift(fft(ifftshift(Sw)))*dw/sqrt(2*pi);

return