Hi Christina, Steve, et. al.
That (adding epsilon to the right-hand x boundary) did the trick, so presumably the unhappy guy was at x>=L_x to the numerical precision with which I wrote down the box size (using "ascii restarts"). I would have guessed that he would have been mapped to x=zero before the initial MPI broadcast, but I guess this doesn't happen.
This is confusing, since Domain::remap() looks like it should be doing exactly this... should take x=xmax -> x=xmin if there are PBCs in x. Domain::remap() looks like it's being called where I would expect:
ReadData::atoms()
Atom::unpack_data()
Domain::remap()
*puzzled*
Thanks,
Craig