Practical Considerations for Computing Dimensional Spectra from Gridded Data

Practical Considerations for Computing Dimensional Spectra from Gridded Data AbstractSpectra are often computed from gridded data to determine the horizontal-scale dependence of quantities such as kinetic energy, vertical velocity, or perturbation potential temperature. This paper discusses several important considerations for the practical computation of such spectra. To ensure that the sum of the spectral energy densities in wavenumber space matches the sum of the energies in the physical domain (the discrete Parseval relation), the constant coefficient multiplying the spectral energy density must properly account for the way the discrete Fourier transform pair is normalized. The normalization factor appropriate of many older FORTRAN-based fast Fourier transforms (FFTs) differs from that in Matlab and Python’s numpy.fft, and as a consequence, the correct scaling factor for the kinetic energy (KE) spectral density differs between one-dimensional FFTs computed using these two approaches by a factor equal to the square of the number of physical grid points. A common algorithm used to compute two-dimensional spectra as a function of the total-wavenumber magnitude sums the contributions from all pairs of x- and y-component wavenumbers whose vector magnitude lies with a series of bins. This approach introduces systematic short-wavelength noise, which can be largely eliminated though a simple multiplicative correction. One- and two-dimensional spectra will differ by a constant if computed for flows in which the KE spectral density decreases as a function of the wavenumber to some negative power. This constant is evaluated and the extension of theoretical results to numerically computed FFTs is examined. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Monthly Weather Review American Meteorological Society

Practical Considerations for Computing Dimensional Spectra from Gridded Data

Loading next page...
 
/lp/ams/practical-considerations-for-computing-dimensional-spectra-from-hTFWL0gWbS
Publisher
American Meteorological Society
Copyright
Copyright © American Meteorological Society
ISSN
1520-0493
D.O.I.
10.1175/MWR-D-17-0056.1
Publisher site
See Article on Publisher Site

Abstract

AbstractSpectra are often computed from gridded data to determine the horizontal-scale dependence of quantities such as kinetic energy, vertical velocity, or perturbation potential temperature. This paper discusses several important considerations for the practical computation of such spectra. To ensure that the sum of the spectral energy densities in wavenumber space matches the sum of the energies in the physical domain (the discrete Parseval relation), the constant coefficient multiplying the spectral energy density must properly account for the way the discrete Fourier transform pair is normalized. The normalization factor appropriate of many older FORTRAN-based fast Fourier transforms (FFTs) differs from that in Matlab and Python’s numpy.fft, and as a consequence, the correct scaling factor for the kinetic energy (KE) spectral density differs between one-dimensional FFTs computed using these two approaches by a factor equal to the square of the number of physical grid points. A common algorithm used to compute two-dimensional spectra as a function of the total-wavenumber magnitude sums the contributions from all pairs of x- and y-component wavenumbers whose vector magnitude lies with a series of bins. This approach introduces systematic short-wavelength noise, which can be largely eliminated though a simple multiplicative correction. One- and two-dimensional spectra will differ by a constant if computed for flows in which the KE spectral density decreases as a function of the wavenumber to some negative power. This constant is evaluated and the extension of theoretical results to numerically computed FFTs is examined.

Journal

Monthly Weather ReviewAmerican Meteorological Society

Published: Sep 6, 2017

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 12 million articles from more than
10,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Unlimited reading

Read as many articles as you need. Full articles with original layout, charts and figures. Read online, from anywhere.

Stay up to date

Keep up with your field with Personalized Recommendations and Follow Journals to get automatic updates.

Organize your research

It’s easy to organize your research with our built-in tools.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

Monthly Plan

  • Read unlimited articles
  • Personalized recommendations
  • No expiration
  • Print 20 pages per month
  • 20% off on PDF purchases
  • Organize your research
  • Get updates on your journals and topic searches

$49/month

Start Free Trial

14-day Free Trial

Best Deal — 39% off

Annual Plan

  • All the features of the Professional Plan, but for 39% off!
  • Billed annually
  • No expiration
  • For the normal price of 10 articles elsewhere, you get one full year of unlimited access to articles.

$588

$360/year

billed annually
Start Free Trial

14-day Free Trial