fix bug with assigning molecule IDs in parallel

2024-02-09 11:04:07 -05:00
parent 524470f9fc
commit 1421d53b8c
1 changed files with 1 additions and 1 deletions
--- a/src/create_atoms.cpp
+++ b/src/create_atoms.cpp
@ -505,7 +505,7 @@ void CreateAtoms::command(int narg, char **arg)

    // molcreate = # of molecules I created

-    tagint molcreate = (atom->nlocal - nlocal_previous) / onemol->natoms;
+    tagint molcreate = (atom->nlocal - nlocal_previous) / onemol->natoms * onemol->nmolecules;

    // increment total bonds,angles,etc