fig, ax = plt.subplots(1,1, figsize=(10, 7))
ax.plot(10**isoEAtAge['log_Teff'], isoEAtAge['log_L'])
ax.plot(10**isoAAtAge['log_Teff'], isoAAtAge['log_L'])
ax.set_ylim(2.0,2.15)
ax.set_xlim(4600,4700)
ax.invert_xaxis()