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