Add methods paper plots/code (Figures 5-7) (#567)

Author: Tom Wagg

examples/methods_paper_plots/fig_5_HR_diagram/README.md

@@ -0,0 +1,12 @@
+# How to recreate/edit the HR diagram from the COMPAS methods paper
+
+This folder contains everything you need to reproduce (or edit!) the HRD from the methods paper. You'll need to do the following:
+
+1. Run `python create_fig_5_grid.py` to create a grid of stars for the plot (you can change this to a custom range of masses or metallicities if you like)
+2. Run `python pythonSubmit.py` to run COMPAS for this grid
+3. Open `make_fig_5.ipynb` in Jupyter Lab/notebook and run the whole thing to create the figure (and learn how to change it)
+
+For reference, the changes from the *default* pythonSubmit.py are just:
+
+- use detailed output
+- use smaller time steps (to get smoother lines for the plot)

examples/methods_paper_plots/fig_5_HR_diagram/create_fig_5_grid.py

@@ -0,0 +1,15 @@
+""" This file creates a grid of stars for the HRD plot """
+
+# use solar and low metallicity
+metallicities = [0.0142, 0.001]
+
+# circular binaries with a separation so wide that they are effectively single
+eccentricity = 0.0
+a = 1e20
+
+with open("grid.txt", "w") as f:
+    # mass range: basically logspace from 0.5 to 50 but with cleaner numbers
+    masses = [0.5, 0.65, 0.80, 1.0, 1.3, 1.6, 2.0, 2.5, 3.2, 4.0, 5.2, 6.5, 8.2, 10.3, 13.0, 16.5, 21, 26, 33, 45, 65, 90, 115, 150]
+
+    grid_lines = ["--initial-mass-1 {} --initial-mass-2 {} --metallicity {} --eccentricity {} --semi-major-axis {} \n".format(masses[i], 0.1, metallicities[j], eccentricity, a) for j in range(len(metallicities)) for i in range(len(masses))]
+    f.writelines(grid_lines)

examples/methods_paper_plots/fig_5_HR_diagram/hrd.pdf

@@ -0,0 +1,12 @@
+# How to recreate/edit the maximum radius plot from the COMPAS methods paper
+
+This folder contains everything you need to reproduce (or edit!) the maximum radius plot from the methods paper. You'll need to do the following:
+
+1. Run `python create_fig_6_grid.py` to create a grid of stars for the plot (you can change this to a custom range of masses or metallicities if you like)
+2. Run `python pythonSubmit.py` to run COMPAS for this grid (note this is probably going to take about 10 minutes!)
+3. Open `make_fig_6.ipynb` in Jupyter Lab/notebook and run the whole thing to create the figure (and learn how to change it)
+
+For reference, the changes from the *default* pythonSubmit.py are just:
+
+- use detailed output
+- use smaller time steps (to get smoother lines for the plot)

examples/methods_paper_plots/fig_5_HR_diagram/make_fig_5.ipynb

@@ -0,0 +1,16 @@
+import numpy as np
+
+""" This file creates a grid of stars for the HRD plot """
+
+# use solar and low metallicity
+metallicities = [0.0142, 0.001]
+
+# circular binaries with a separation so wide that they are effectively single
+eccentricity = 0.0
+a = 1e20
+
+with open("grid.txt", "w") as f:
+    masses = np.round(np.logspace(np.log10(0.5), np.log10(150.0), 500), 3)
+
+    grid_lines = ["--initial-mass-1 {} --initial-mass-2 {} --metallicity {} --eccentricity {} --semi-major-axis {} \n".format(masses[i], 0.1, metallicities[j], eccentricity, a) for j in range(len(metallicities)) for i in range(len(masses))]
+    f.writelines(grid_lines)