Code for reproducing the analysis in the Nummelin et al (2020; Diagnosing the Scale and Space Dependent Horizontal Eddy Diffusivity at the Global Surface Ocean, JPO) results.
Workflow:
- Copy the git repository
- Get miniconda and setup the required python environment by running
$ conda env create -f environment.yml
- Go to data/raw/ and follow the instructions in the download.sh in order to download the data.
- Go to side_packages and use the fetch.sh MicroInverse and naturalHHD packages from GitHub.
- Once the data is downloaded and the required packages are in place, run the following
python bin/process_SST.py # Produce low-pass filtered SST python bin/process_eddy_atlas.py # Bin the Eddy trajectories to a grid python bin/run_inversion.py # Run the inversion python bin/run_HelmHoltzDecomposition.py # Run the Helmholtz-Hodge decomposition
- At this point all the major data processing is done.
The final step is to do some postprocessing and plot the figures in Nummelin et al. 2020.
python bin/plotting/run_plotting.py
Please note that this project requires a large amount of data and reproducing all the steps will take a significant amount of computing time and power. Please be in touch if you are interested in specific parts of the full project.
Project based on the cookiecutter science project template.