atom_style granular atom_modify map array boundary m p m newton off echo both communicate single vel yes units si region reg block -10 10 -10 10 -10 10 units box create_box 1 reg create_atoms 1 single 0 0 5 units box neighbor 1 bin neigh_modify delay 0 fix m1 all property/global youngsModulus peratomtype 5.e6 fix m2 all property/global poissonsRatio peratomtype 0.45 fix m3 all property/global coefficientRestitution & peratomtypepair 1 0.3 fix m4 all property/global coefficientFriction & peratomtypepair 1 0.5 fix m5 all property/global k_finnie & peratomtypepair 1 1.0 # New pair style pair_style gran/hertz/history 1 0 #Hertzian without cohesion pair_coeff * * timestep 0.01 # works # variable b equal 1.0 # fix flow all addforce NULL ${b} NULL # doesnt work # fix flow all addforce NULL v_b NULL # doesnt work # variable b atom type # fix flow all addforce NULL v_b[1] NULL # doesnt work variable b atom type fix flow all addforce NULL v_b NULL fix integr all nve dump dmp all custom 200 post/dump.tw & id type type x y z ix iy iz vx vy vz & fx fy fz omegax omegay omegaz radius run 1000