Fixed problem with repetitions>1 due to new data structures -- now
computes initial next_bond correctly