Open boundaries don't work if `UpdateCallback` is not the last callback in the list

For some reason, we are getting 1e16 coordinates (which are used by inactive particles) into the neighborhood search update (which is only done for active particles):

┌ Info: 
│   point_coords =
│    2-element SVector{2, Float64} with indices SOneTo(2):
│     1.0e16
└     1.0e16
┌ Info: 
└   point in eachindex_y = true
ERROR: LoadError: particle coordinates are NaN or outside the domain bounds of the cell list
cell (3333333333333333505, 3333333333333333505) is out of bounds for cell grid of size [2:72, 2:37]

Start by reproducing the reported open-boundary case and trace the callback list into the neighborhood search update, focusing on how inactive-particle coordinates reach that update when UpdateCallback is not last. Done means the reported out-of-bounds cell-list error no longer occurs and open boundaries work regardless of callback ordering; the issue names no files or tests.