Dear Lammps users,

I am testing some lattice system under displacement-control loading. The

steps are as follows:

Step 1. Run NVE assemble of the lattice for 1000 steps while holding the

displacements on the boundary to be 1.0.

Step 2. Increase the boundary displacements to be 2.0; Run NVE assemble

for another 1000 steps, where the initial configurations (i.e. partial

coordinates and velocities) are chosen to be the results from the last

integration in Step 1.

Step 3. Increase the boundary displacements to be 3.0; Run NVE assemble

for another 1000 steps, where the initial configurations (i.e. partial

coordinates and velocities) are chosen to be the results from the last

integration in Step 2.

Below is a part of my input script on applying 'variable loop' to run this

simulation.