The International Geomagnetic Reference Field (IGRF) is a standard mathematical description of the large-scale structure of the Earth's main magnetic field and its secular variation. It was created by fitting parameters of a mathematical model of the magnetic field to measured magnetic field data from surveys, observatories and satellites across the globe. The IGRF has been produced and updated under the direction of the International Association of Geomagnetism and Aeronomy (IAGA) since 1965. The International Geomagnetic Reference Field (IGRF) is a standard mathematical description of the large-scale structure of the Earth's main magnetic field and its secular variation. It was created by fitting parameters of a mathematical model of the magnetic field to measured magnetic field data from surveys, observatories and satellites across the globe. The IGRF has been produced and updated under the direction of the International Association of Geomagnetism and Aeronomy (IAGA) since 1965. The IGRF model covers a significant time span, and so is useful for interpreting historical data. (This is unlike the World Magnetic Model, which is intended for navigation in the next few years.) It is updated at 5-year intervals, reflecting the most accurate measurements available at that time. The current 12th edition of the IGRF model was released in 2015 and is valid from 1900 until 2020. For the interval from 1945 to 2010, it is 'definitive', meaning that future updates are unlikely to improve the model in any significant way. The IGRF models the geomagnetic field B → ( r , ϕ , θ , t ) {displaystyle {vec {B}}(r,phi , heta ,t)} as a gradient of a magnetic scalar potential V ( r , ϕ , θ , t ) {displaystyle V(r,phi , heta ,t)} B → ( r , ϕ , θ , t ) = − ∇ V ( r , ϕ , θ , t ) {displaystyle {vec {B}}(r,phi , heta ,t)=- abla V(r,phi , heta ,t)} The magnetic scalar potential model consists of the Gauss coefficients which define a spherical harmonic expansion of V {displaystyle V} V ( r , ϕ , θ , t ) = a ∑ ℓ = 1 L ∑ m = 0 ℓ ( a r ) ℓ + 1 ( g ℓ m ( t ) cos m ϕ + h ℓ m ( t ) sin m ϕ ) P ℓ m ( cos θ ) {displaystyle V(r,phi , heta ,t)=asum _{ell =1}^{L}sum _{m=0}^{ell }left({frac {a}{r}} ight)^{ell +1}left(g_{ell }^{m}(t)cos mphi +h_{ell }^{m}(t)sin mphi ight)P_{ell }^{m}left(cos heta ight)} where r {displaystyle r} is radial distance from the Earth's center, L {displaystyle L} is the maximum degree of the expansion, ϕ {displaystyle phi } is East longitude, θ {displaystyle heta } is colatitude (the polar angle), a {displaystyle a} is the Earth's radius, g ℓ m {displaystyle g_{ell }^{m}} and h ℓ m {displaystyle h_{ell }^{m}} are Gauss coefficients, and P ℓ m ( cos θ ) {displaystyle P_{ell }^{m}left(cos heta ight)} are the Schmidt normalized associated Legendre functions of degree l {displaystyle l} and order m {displaystyle m} . The Gauss coefficients are assumed to vary linearly over the time interval specified by the model.