rerun read_dump issue

Dear Lammps Users,

I am trying to perform a rerun for a simple LJ system.
I'm using a native atom dump file, which I wrote using the following command:

dump  myDump3 all atom ${out_freq} traj.file #dump data

However, when I launch the rerun I get the following error:

ERROR: Illegal read_dump command (../read_dump.cpp:606)

I suppose the read_dump is referred to the rerun command, since in my run I don't have any read_dump command.

Here's the input file:

#LJ bench run
echo both
read_restart    restart.start

pair_style lj/cut 2.5
pair_coeff * * 1 1 
pair_modify shift yes

neighbor        3.23 bin    #binning neighbor-list style, skin length (2.0 default)
neigh_modify    every 10 delay 0 check yes #modify neigh lists every N step

timestep 	0.001

rerun traj.file first 500000 dump xs ys zs 

And this is how the dump trajectory called by rerun looks like:

ITEM: TIMESTEP
500000
ITEM: NUMBER OF ATOMS
920
ITEM: BOX BOUNDS pp pp pp
0 10
0 10
0 10
ITEM: ATOMS id type xs ys zs
685 1 0.0784089 0.121576 0.0251055
361 1 0.395694 0.397481 0.135961
531 1 0.169525 0.394286 0.0922789
327 1 0.0620085 0.233832 0.0355341

I'm trying to follow the manual, but I'm new to Lammps.
Does anyone know how should I correct my input?

Thank you and best regards

Claudio

Dear Lammps Users,
I am trying to perform a rerun for a simple LJ system.
I'm using a native atom dump file, which I wrote using the following command:

dump myDump3 all atom ${out_freq} traj.file dump data

However, when I launch the rerun I get the following error:

ERROR: Illegal read_dump command (../read_dump.cpp:606)

I suppose the read_dump is referred to the rerun command, since in my run I don't have any read_dump command.

​exactly. ...and you should look at the read_dump documentation.
it explains that for reading x, y, or z​ coordinates, you should use x, y,
or z, respectively, regardless of whether they are written out plain,
scaled or unwrapped. LAMMPS will reconstruct the best data from it
transparently.

so change your rerun command to:


rerun traj.file first 500000 dump x y z

and it should work.

axel.

Ok, thank you for the reply.

Best

Claudio