Internal function
time_bp_to_index(time_bp, time_steps)
time_bp
vector of times BP
time_steps
time steps for which reconstructions are available
the indeces for the relevant time steps