A Fast Spherical Harmonics Transform for Global NWP and Climate Models

AbstractVery high-resolution spectral transform models are believed to become prohibitively expensive because of the relative increase in computational cost of the Legendre transforms compared to the gridpoint computations. This article describes the implementation of a practical fast spherical harmonics transform into the Integrated Forecast System (IFS) at ECMWF. Details of the accuracy of the computations, of the parallelization, and memory use are discussed. Results are presented that demonstrate the cost effectiveness and accuracy of the fast spherical harmonics transform, successfully mitigating the concern about the disproportionally growing computational cost. Using the new transforms, the first T7999 global weather forecast (equivalent to ≈2.5-km horizontal grid size) using a spectral transform model has been produced.

[1]  Pierre Bénard,et al.  Dynamical kernel of the Aladin–NH spectral limited‐area model: Revised formulation and sensitivity experiments , 2010 .

[2]  S. Gregory.  Redshifts and morphology of galaxies in the Coma cluster , 1975 .

[3]  C. Temperton.  Self-sorting mixed-radix fast Fourier transforms , 1983 .

[4]  S. Orszag.  Transform method for the calculation of vector-coupled sums: Application to the spectral form of the vorticity equation , 1970 .

[5]  Lars Isaksen,et al.  The IFS Model: A Parallel Production Weather Code , 1995, Parallel Comput..

[6]  Reiji Suda.  Fast Spherical Harmonic Transform Routine FLTSS Applied to the Shallow Water Test Set , 2005 .

[7]  Christiane Jablonowski,et al.  Block-structured adaptive meshes and reduced grids for atmospheric general circulation models , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[8]  S. L. Belousov,et al.  Tables of normalized associated Legendre polynomials , 1965 .

[9]  D. Williamson.  The Evolution of Dynamical Cores for Global Atmospheric Models(125th Anniversary Issue of the Meteorological Society of Japan) , 2007 .

[10]  Mats Hamrud,et al.  Revolutionizing Climate Modeling with Project Athena: A Multi-Institutional, International Collaboration , 2013 .

[11]  Colin J. Cotter,et al.  Mixed finite elements for numerical weather prediction , 2011, J. Comput. Phys..

[12]  D. S. Seljebotn.  WAVEMOTH–FAST SPHERICAL HARMONIC TRANSFORMS BY BUTTERFLY MATRIX COMPRESSION , 2011, 1110.4874.

[13]  Mark Tygert,et al.  Fast algorithms for spherical harmonic expansions, III , 2009, J. Comput. Phys..

[14]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[15]  Clive Temperton.  On Scalar and Vector Transform Methods for Global Spectral Models , 1991 .

[16]  John Thuburn,et al.  Horizontal grids for global weather and climate prediction models: a review , 2012 .

[17]  Paul N. Swarztrauber,et al.  On Computing the Points and Weights for Gauss-Legendre Quadrature , 2002, SIAM J. Sci. Comput..

[18]  Mark Tygert,et al.  Fast algorithms for spherical harmonic expansions, II , 2008, J. Comput. Phys..

[19]  Pierre Benard,et al.  Stability of Semi-Implicit and Iterative Centered-Implicit Time Discretizations for Various Equation Systems Used in NWP , 2003, physics/0304114.

[20]  A. Erdelyi,et al.  Fonctions Spheriques de Legendre et Fonctions Spheroidales. , 1959 .

[21]  Per-Gunnar Martinsson,et al.  On the Compression of Low Rank Matrices , 2005, SIAM J. Sci. Comput..

[22]  Reiji Suda,et al.  A fast spherical harmonics transform algorithm , 2002, Math. Comput..

[23]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[24]  Per-Gunnar Martinsson,et al.  An Accelerated Kernel-Independent Fast Multipole Method in One Dimension , 2007, SIAM J. Sci. Comput..

[25]  Philippe Courtier,et al.  Representation errors in various grid and spectral truncations for a symmetric feature on the sphere , 1996 .

[26]  Michael O'Neil,et al.  An algorithm for the rapid evaluation of special function transforms , 2010 .

[27]  Nils Wedi,et al.  High-Resolution Global Climate Simulations with the ECMWF Model in Project Athena: Experimental Design, Model Climate, and Seasonal Forecast Skill , 2012 .