boundary p p p dimension 3 region total block -10 10 -10 10 -10 10 create_box 1 total create_atoms 1 single 0 0 0 fix prop all property/atom mol ghost yes variable cluster atom ((id-1)/10)+1 set id * mol v_cluster