def ignition_delay(states, species):
    """
    This function computes the ignition delay from the occurrence of the
    peak in species' concentration.
    """
    i_ign = states(species).Y.argmax()
    return states.t[i_ign]