fig, ax = plt.subplots(1,1,figsize=(10, 7))
ax.plot(longTimes/(np.pi*1e7*1e6), hotSideTemp, label="Hot Side")
ax.plot(longTimes/(np.pi*1e7*1e6), coldSideTemp, label="Cold Side")
ax.set_xlabel("Time [Myr]", fontsize=27)
ax.set_ylabel("Teff [K]", fontsize=27)
ax.legend()