# for t in range(0, len(star_formation_times) - 1):
    #     print(star_formation_times[t], node["time"], star_formation_times[t+1], node.ds.all_data()[("gas", "H2_fraction")])

    #     if star_formation_times[t] < node["time"] and star_formation_times[t+1] > node["time"]:
    #         time_before = node["time"] - star_formation_times[t]
    #         time_after = star_formation_times[t+1] - node["time"]

    #         if 0.33 > min(time_before, time_after):
    #             if node.ds.all_data().max(("gas", "H2_fraction")) <= 1E-4:
    #                 return False

    #             # @TODO - Check mass ratio vs cosmic mean

    #         break