Changed tign(B(j,1),B(j,2)) to min(tign(B(j,1),B(j,2)),ts*interval) in
the tign calculation for the case when the neighbor is being updated
after a few time_steps