<div><div><font face="Courier" size="4"><br></font></div><div><font face="Courier" size="4">Dear all,</font></div><div><font face="Courier" size="4"><br></font></div><div><font face="Courier" size="4">These days, I used the 'fft2df' function to calculate the spectra of a two-dimensional physical space [f(ny,nx)]. Then I got a (ny * nx/2+1) shaped array [F(ny,nx/2+1)]. However, how to get the wave numbers corresponding to the spectral-spatial array F? Is zonal wave number k = i/nx, meridional wave number l = j/ny, and total wave number K = sqrt(k^2+l^2),(j=0,1,...,ny-1; i=0,1,...,nx/2)?&nbsp;</font></div><div><font face="Courier" size="4">And I also wonder why return array from 'fft2df' is a (ny * nx/2+1) shape. Why not (ny/2+1 * nx/2+1) shape? Because from (ny/2+2) onwards, all meridional wave numbers will be larger than Nyquist frequency (= 2*sampling rate) and therefore have to be cut off.&nbsp;</font></div><div><font face="Courier" size="4">Could someone can answer my doub
 ts above? Thank you very much!</font></div><div><font face="Courier" size="4"><br></font></div><div><font face="Courier" size="4"><br></font></div><div><font face="Courier" size="4">Yours sincerely,</font></div><div><font face="Courier" size="4">Yong-Jie Huang</font></div></div><br><br><br>