IGRF model Implementation (Calculating the Partial Derivatives of the Legendre Polynomials)
10 views (last 30 days)
Show older comments
I was implementing the IGRF-13 Geomagnetic model, and a parts of the equation (3 a, b, and c) requires calculating the Associated Legendre polynomials with Schmidt quasi-Normalization and its partial derivatives.
For the Associated Legendre polynomials, I used the function Legendre(n, x, 'sch') and it works perfectly. But now how would I calculate their partial derivatives?
Equations 19 a, b, and c provides some recursive formulae, but I am not getting the correct result possibly because I might have implemented them incorrectly. Even the sample MATLAB code provided in the paper is not producing the proper results.
I would really appreciate it if someone could provide a proper implementation for calculating the Partial derivatives of these Legendre Polynomials.
Thank you,
Abinay Brown
Reference: https://hanspeterschaub.info/Papers/UnderGradStudents/MagneticField.pdf
0 Comments
Answers (1)
Bjorn Gustavsson
on 7 Jun 2021
Have a look at the IGRF implementations available on the file exchange: international-geomagnetic-reference-field-igrf-model, international-reference-ionosphere-iri-model, earth-magnetic-field-with-igrf-13th-gen and igrf-magnetic-field. I've used one of the packages bby Drew Compston.
HTH
See Also
Categories
Find more on Oceanography and Hydrology in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!