read_data add problem

hi Steve
I can’t figure out if I am misusing the “read_data add” command or I found a bug.

Here is the min working example:

atom_style full
#atom_modify map hash
#atom_modify map array
boundary p p p
read_data test.data
read_data test.data add append shift 0 0 4

which outputs

……
read_data test.data
Reading data file …
orthogonal box = (-20 -20 -20) to (20 20 20)
1 by 1 by 1 MPI processor grid
reading atoms …
5 atoms
scanning bonds …
4 = max bonds/atom
scanning angles …
4 = max angles/atom
reading bonds …
4 bonds
reading angles …
4 angles
Finding 1-2 1-3 1-4 neighbors …
Special bond factors lj: 0 0 0
Special bond factors coul: 0 0 0
4 = max # of 1-2 neighbors
3 = max # of 1-3 neighbors
3 = max # of 1-4 neighbors
4 = max # of special neighbors
read_data test.data add append shift 0 0 4
Reading data file …
orthogonal box = (-20 -20 -20) to (20 20 24)
1 by 1 by 1 MPI processor grid
reading atoms …
5 atoms
scanning bonds …
ERROR on proc 0: Invalid atom ID in Bonds section of data file (…/atom.cpp:973)

this seems to trigger an error regarding a headroom variable on a map of sorts called “map_max_tag” which in this case is 10 whereas the atom ids are 11 and 15.

I can’t figure out why the new ids are 6-10 since there are only 5 atoms in the data file nor the logic in atom_map.

any help is appreciated.
thanks,
Reese
the little test data file is attached.

NiN4.data (869 Bytes)

ATT00001.htm (801 Bytes)

It’s probably a bug - another user has reported a similar

issue …

Steve