Change_box

Dear Lammps Users,

I developed an input script to deform the simulation box (2D) and then perform stress calculations on a GRANULAR packing. I am using the September 2013 lammps version. But I am getting an error.

ERROR :

change_box all triclinic
ERROR: Cannot change_box after reading restart file with per-atom info (…/change_box.cpp:52)

Dear Lammps Users,

I developed an input script to deform the simulation box (2D) and then
perform stress calculations on a GRANULAR packing. I am using the September
2013 lammps version. But I am getting an error.

ERROR :

change_box all triclinic
ERROR: Cannot change_box after reading restart file with per-atom info
(../change_box.cpp:52)
-------------------------------------------------------------------------------------------------------------------
Input script was :

read_restart ./restart

pair_style gran/hooke/history 1.0 0.0 1.0 0.0 0.0 1
pair_coeff * *

neighbor 0.2 bin
neigh_modify every 1 delay 0 check no

velocity all set 0.0 0.0 0.0 units box

reset_timestep 0

change_box all triclinic
fix lebc all deform 1 xy erate ${mystrainrate} remap x
------------------------------------------------------------------------------------------------------------------------
I found some suggestions in the mail list to fix this error but i couldn't
understand very much.

i.e.,
Cannot change_box after reading restart file with per-atom info This is
because the restart file info cannot be migrated with the atoms. You can get
around this by performing a 0-timestep run which will assign the restart
file info to actual atoms. I am not clear myself what does it means by
0-timestep run means ?

it means to use the command:

run 0

that will execute the setup phase of a run, but not actually run any
timesteps that advance the positions.

axel.