Using an ocean numerical model, we studied how the polar transition zone is impacted by surface buoyancy fluxes and the equation of state for seawater.
The reproducibility is made through the snakemake workflow (Python). Output of the model runs are available on zenodo for only rerunning analyses. However, the runs themselves are also reproducible with the help of snakemake (slightly harder as the model runs need to be conducted on HPC).
I tried hard to make it reproducible, so hopefully this paper can serve as an example on how reproducibility can be achieved.
I think that being reproducible with only few commands to type in a terminal is quite an achievment. At least in my field, where I usually see code published along with paper, but with almost no documentation on how to rerun it.
One of the hardest point for reproducibility is getting the proper environment with all installed packages. This can be the tricky part, and this may be where the analyses become non-reproducible, unfortunately. However, I think that with a little effort, it is possible to fix the potential environment issues.
Resources
Associated event
