[lammps-users] variable command

[Joaquin Rodriguez Nieva <joaquinrodnie@hotmail.com>]

Hello,

I am having problems with the "variable" command (lammps-7Sep09, serial). I am trying to group the atoms that are inside a complicated geometry (in the example I chose a cylinder of radius 2 to make it simpler) but I cannot use the command "variable x equal x[id]" to extract the coordinates of atom id. I am using the following loop and the error sign is under it (the error also appears without the loop):

###########LOOP#############
variable i loop $N
label loop

variable x1 equal x[$i]
variable y1 equal y[$i]
variable dd equal (${x1}**2+${y1}**2)

if ${dd} <= 2.0 then "group Cyl id $i"

next i
jump in.spt loop
############################

ERROR: Indexed per-atom vector in variable formula without atom map

Thank you and best regards,

Joaquín R.N.

---

Hotmail: Trusted email with powerful SPAM protection. Sign up now.