When trying to create bonds in "LIGGGHTS-WITH-BONDS" , I get the following error:
creating bond btw atoms 726 and 1741 (i has now 12 bonds) at step 1001
creating bond btw atoms 726 and 1743 (i has now 13 bonds) at step 1001
ERROR on proc 0: New bond exceeded bonds per atom in fix bond/create (../fix_bond_create_gran.cpp:485)
I tried changing the maximal amount of bonds in the atom style command so that the system does not build more bonds than the fix bond/create command allows. But this also does not work. The only workaround that worked for me was lowering the timesteps during which the bonds are created, but that means less bonds which is not very satisfying.
I attached the file that creates the system and saves a restart file and the file that loads it and creates the bonds.
