Corrected fft.shift for Sf to only shift along the frequency axis, *NOT* ensemble (z) axis