fig, ax = plt.subplots(1,1,figsize=(10,7))
ax.invert_yaxis()
ax.invert_xaxis()
glob = pathlib.Path("../outputs.kttauTest/PopA+0.24/alpha-1.901/").rglob("*.iso")
for pID, path in enumerate(glob):
    iso, isoMeta = read_iso(path)
    print(isoMeta['M'])
    ageCut = iso['Age'].values < 12e9
    Log_T = iso['Log_T'][ageCut].values
    Log_g = iso['Log_g'][ageCut].values
    T = 10**Log_T
    ax.plot(T, Log_g)
ax.set_xlim(6100, 4000)