# A. Pohl, 18 Feb 2025
# This is the code supporting the manuscript 'Spatial biases in oxygen-based Phanerozoic seawater temperature reconstructions'.

# CONTENT
# All files, code and data required to repeat the analysis and generate the figures.

# REQUIREMENTS
# In addition to a series of standard python packages, you will need a pygplates install (linked and used in step 1).
# Tested with Python 3.9.6 on MacBook Pro Apple M1 Max working under MacOS Sequoia 15.3

# HOW-TO:
# First you will need to download HadCM3 data. They are freely accessible and are thus not included here to limit total file size. All needed info and scripts to be found in directories: ./GCM_output/HadCM3_output_tfks)

# Guidelines to repeat main calculations:
# Start by running 'step0_generate_d18Osw_maps.py' to generate the content used in Fig. 1 and the NetCDF files needed in next step.
#     > python3 step0_generate_d18Osw_maps.py
# Then, run 'step1_plot_maps_and_trends.py' to generate the content used in Figs 2–7 and the .csv tables needed in next step.
#     > python3 step1_plot_maps_and_trends.py
# Finally, run 'step2_revisit_and_plot_StabIso_trend.py' to generate Fig. 8. 
#     > python3 step2_revisit_and_plot_StabIso_trend.py

# Sensitivity tests can be repeated by adjusting the user-params on top of the scripts before launching them.