HELP: why this "read_restart" command does not work?

LAMMPS LAMMPS General Discussion
1

Hello,everyone!

I use the following methods to test if the ‘read_restart’ command does, but getting an absonant results.

First, (whole) it is the whole process. I grouped a region and setforce to 0 0 0. Then
①addforced on the region with the value 0.001，run 100, dump value of the force on the region is 0.001
②addforced on the region with the value 0.002, run 100, dump value of the force on the region is 0.003.
③only run 100 without any command, dump value kept the same
④addforce 0.001，dump value 0.004.
⑤addforce 0.002，dump value 0.006.

Second, (part1) i wrote the code whoes process only contains the first two steps, that is ①②, and used the command ‘write_restart’ command at last. During the part1 process, every result was the same as the whole process.

Third, (part2) ‘read_restart’ command was used to continue the process. However, after step ③, the result in this part2 was totally different from the resut in the whole part.

So, my question is why this “read_restart” command does not work? and how to correct use the restart command. (fix, pair_style, pair_coeff were all reused in the part2).

THANK YOU FOR YOUR HELP！

2

First off, the read_restart command does work. It is being tested regularly.
It is much more likely that there is an issue with your scripts.

In order to tell what goes wrong, however, it is necessary to have a (simplified) version of your input decks and the corresponding log files in order to try an reproduce it. Ideally this would be done on a small test system so the validation would be done quickly on a desktop.