Data and code used in "Satellite magnetic data reveal interannual waves in Earth's core"

Eigen mode solutions and code to obtain them for the results presented in Satellite magnetic data reveal interannual waves in Earth's core. The package uses the freely available code Mire.jl.

Prerequisites

Installed python3 with matplotlib ≥v2.1, cmocean and cartopy. A working Julia ≥v1.7.

Run

In the project folder run

julia --project=.

Then, from within the Julia REPL run

]instantiate

at first time, to install all dependencies.

After that, to compute all plots, run

using QGMCSat
allfigs()

They're automatically saved in the "figs" subfolder of the repository.

If loading QGMCSat fails, due to a missing cartopy or cmocean in the python version. Run (within Julia)

ENV["PYTHON"] = "python" #this should point to the python version that has cartopy installed
]build PyCall

To calculate all data, run

using QGMCSat
calculate_data()

This will take several hours/days depending on the machine (needs enough memory).

Individual data can be accessed directly through the .jld2 files from Julia. You can check out the individual figure functions to get an idea where which data is stored.

If there are any issues or questions, please don't hesitate to get in touch!