a question on "atoms"

LAMMPS LAMMPS Mailing List Mirror
1

I have used “atoms” and assigned varialble for it to denote the number atoms of the system as follows.
However, it showed the result as the value of the variable is zero.

I tested it using the most current lammps (windows,no-mpi, 1 July) execution file and lammps file of October 14, 2012 in linux system.

The output is “Natoms:0”.
The input file is as follows.

3d Lennard-Jones melt

units lj
atom_style atomic

lattice fcc 0.8442
region box block 0 10 0 10 0 10
create_box 1 box
create_atoms 1 box
mass 1 1.0

velocity all create 3.0 87287

pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0 2.5

variable N equal atoms
print “Natoms:${N}”

neighbor 0.3 bin
neigh_modify every 20 delay 0 check no

fix 1 all nve

#dump id all atom 50 dump.melt
#dump 1 all image 25 image.*.jpg type type &

axes yes 0.8 0.02 view 60 -30

#dump_modify 1 pad 3

thermo 50

run 250

Would it be concerned with my incorrect compling? Otherwise, will it be a bug?
Thanks for your answer in advance.

Cheers,

Young

2

this topic was discussed recently.
you version (windows and source) needs to be
updated, or you can access atoms only after
you had "thermo" output. it is not initialized
earlier. steve recently changed the code to
do the initialization so that you can use atoms
before and run or mimize commands. with
the older versions that you have, this is not
possible.

axel.

3

Thanks for your answer, axel.

Best regards,

Young