Fast and accurate evaluation of geomagnetic field elements at arbitrary points in space
Autor: | Pencho Petrushev, Kamen G. Ivanov |
---|---|
Rok vydání: | 2020 |
Předmět: |
Physics
Reference ellipsoid 010504 meteorology & atmospheric sciences Fortran business.industry Spherical harmonics 010103 numerical & computational mathematics 01 natural sciences Ellipsoid Computational science Geophysics Earth's magnetic field Software Geochemistry and Petrology Harmonics 0101 mathematics MATLAB business computer 0105 earth and related environmental sciences computer.programming_language |
Zdroj: | Geophysical Journal International. 224:181-190 |
ISSN: | 1365-246X 0956-540X |
DOI: | 10.1093/gji/ggaa439 |
Popis: | SUMMARY An algorithm and software are developed for fast and accurate evaluation of the elements of the geomagnetic field represented in high-degree (>720) solid spherical harmonics at many scattered points in the space above the surface of the Earth. The algorithm is based on representation of the geomagnetic field elements in solid ellipsoidal harmonics and application of tensor product needlets. Open source FORTRAN and MATLAB realizations of this algorithm that rely on data from the Enhanced Magnetic Models 2015, 2017 (EMM2015, EMM2017) have been developed and extensively tested. The capabilities of the software are demonstrated on the example of the north, east and down components of the geomagnetic field as well as the derived horizontal intensity, total intensity, inclination and declination. For the range from −417 m under the Earth reference ellipsoid up to 1000 km above it the FORTRAN and MATLAB versions of the software run 465 and 189 times faster than the respective FORTRAN and MATLAB versions of the software using the standard spherical harmonic series method, while the accuracy is less than 1 nT and the memory (RAM) usage is 9 GB. |
Databáze: | OpenAIRE |
Externí odkaz: |