This package provides Python code to calculate magnetic fields of Andrew Jackson's GUFM1 historical magnetic field model. This code was largely based off of the Fortran code eval_field.f found on his website, but is completely re-written in native Python.

To use, simply clone the git directory, then either run the code from the directory directly, or run python setup.py install to allow for importing the method in any python project. A few different routines are currently available, allowing users to calculate radial field strength at any spatial location, rms field strength over the whole core, and other metrics of core power spectrum.

Users should be able to be calculate parameters for any year from 1590 to 1990 from the main gufm1_data.txt file, but you can also use data files from individual years downloaded directly from his website if desired.

For examples, see the [GUFM1 Demonstration](examples/GUFM1 Demonstration.ipynb) ipython notebook.

Details of the GUFM1 model can be found on the website, or in the paper

Jackson, A., Jonkers, A. R., & Walker, M. R. (2000). Four centuries of geomagnetic secular variation from historical records. Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, 358(1768), 957-990.