The "DynamoRegimePlots.py" file requires the results of our MCMC interior structure models for each moon (e.g., "MCMC2_Europa_06032024.csv"). However, these files are too large to upload directly to the AGU website and GitHub as well, even in a compressed format. If you want to reproduce Figure 3, then do the following:

1) Run "IcyMoonStructureMCMC.py" or "IoStructure.py" with the appropriate parameters as described in the Supplementary Information. Make sure to change the filename to save after each run (both the .h5 and .csv).

2) Move the .csv files to the same directory as "DynamoRegimePlots".

3) Run "DynamoRegimePlots.py". Ensure that the filenames are consistent with Step 1.