Wondering what could cause the "Did not assign all atoms correctly" error when using the "read_data" method (reading in ascii restart file).
The boundary conditions are "p s p". It's a 2d setup, so all z co-ords are set to zero. I'm checking now that all the y-values are within the initial range (but that shouldn't matter with "s" type boundaries, right?)
There are 435288 atoms in the ascii file. The number of atom lines in the data file and the "435288 atoms" line in the header agree. LAMMPS seems to successfully parse all the atoms, since it doesn't complain until the "Did not assign all atoms correctly" error which happens after the MPI broadcast of atoms to all procs. In the LAMMPS output, I get: