Global spherical harmonic paleomagnetic field model LSMOD.1 describes the magnetic field evolution from 50 to 30 ka BP based on published paleomagnetic sediment records and volcanic data. The time interval includes the Laschamp (~41 ka BP) and Mono Lake (~34 ka BP) excursions. The model is given with Fortran source code to obtain spherical harmonic magnetic field coefficients for individual epochs and to obtain time series of magnetic declination, inclination and field intensity from 49.95 to 30 ka BP for any location on Earth. For details see M. Brown, M. Korte, R. Holme, I. Wardinski and S. Gunnarson (2018): Earth's magnetic field is probably not reversing. PNAS, 115, 5111-5116.

File overviewLSMOD.1 -- ASCII file containing the time-dependent model by a list of spline basis knot points and spherical harmonic coefficients for these knot points.LSfield.f -- Fortran source code to obtain time series predictions of declination, inclination and intensity from the model file.LScoefs.f -- Fortran source code to obtain the spherical harmonic coefficients for an individual age from the time-dependent model file. The data are licenced under the Creative Commons Attribution 4.0 International Licence (CC BY 4.0) and the Fortran Codes under the Apache License, Version 2.0. The Fortran source code should work with any standard Fortran 77 or higher compiler. Each of the two program files can be compiled separately, all required subroutines are included in the files. The model file, LSMOD.1 or LSMOD.2, is read in by the executable program and has to be in the same directory. The programs work with interactive input, which will be requested when running the program.