Japan Geoscience Union Meeting 2016

Presentation information


Symbol S (Solid Earth Sciences) » S-GD Geodesy

[S-GD23] Geodesy General Contributions / Global Geodetic Observing System

Mon. May 23, 2016 5:15 PM - 6:30 PM Poster Hall (International Exhibition Hall HALL6)

Convener:*Koji Matsuo(Geospatial Information Authority of Japan), Takahito Kazama(Graduate School of Science, Kyoto University), Ryoji Kawabata(Geospatial Information Authority of Japan)

5:15 PM - 6:30 PM

[SGD23-P08] Comparison of computational methods of associated Legendre functions

*Takeshi Enomoto1 (1.Disaster Prevention Research Institute, Kyoto University)

Keywords:Sphere, Numerical method, Spectral method

Spherical harmonics composed of trigonometric and associated Legendre functions are used in geophysics and in other disciplines in science that deal with phenomena on the spherical surface. In meteorology, spherical harmonics are used to expand the prognostic variables of the atmosphere to compute the dynamical process in an atmospheric general circulation model or to analyze the energy spectrum. Recent increase in computing power allows us to use a large truncation wave number to achieve a horizontally high resolution. However, the values of associated Legendre functions of high order and degree cannot be computed accurately with the traditionally-used three-point recurrence in double precision. Alternatively, underflows can be avoided with the four-point recurrence in double precision or the three-point recurrence in extended floating arithmetic. Comparison of the two methods shows that the former and the latter have advantages in accuracy and speed, respectively. In addition, a method is shown to improve accuracy of the Fourier expansion of Legendre polynomials used along with the four-point recurrence.