Data and Figures used in "Fast Quasi‐Geostrophic Magneto‐Coriolis Modes in the Earth's Core"

Data, plotting routines and analysis routines to reproduce all results from the article Fast Quasi‐Geostrophic Magneto‐Coriolis Modes in the Earth's Core. The package uses the freely available code Mire.jl.

Prerequisites

Installed python3 with matplotlib ≥v2.1 and cartopy. A working Julia ≥v1.6.

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 FastMCM
allfigs() #or allfigs(true) to use LaTeX fonts

If loading FastMCM fails, due to a missing cartopy 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 data, run

using FastMCM
N=9
calc1dtm()
calcLHSRHS(N)
calcesol(N)
postprocess(N)

for a given N. Note, for large N (e.g. 35) this may take a lot of memory and a long time (~13h). 

If there are any issues, please don't hesitate to get in touch.