## Example 3:	drop rigid hexapods falling under gravity
#  New tricks: add rigid bodies
#  Commands that appears in previous examples will not be repeatedly explained.
#  More command details in the documentation https://lammps.sandia.gov/doc/Commands.html

atom_style	sphere
atom_modify	map array
boundary	f f fm
newton		off
comm_modify	vel yes cutoff 2.5

region		reg block -8 8 -8 8 0 15 units box
create_box	1 reg

neighbor	0.5 bin
neigh_modify	delay 0

pair_style	gran/hertz/history 200.0 200.0 5.0 5.0 0.5 1
pair_coeff	* *

timestep	0.01

fix		prop all property/atom mol ghost yes
molecule	hexapod hexapod.template
fix		hexarigid all rigid/small molecule mol hexapod
neigh_modify	exclude molecule/intra all

fix		1 all gravity 1.0 spherical 0.0 -180.0

region		reg2 block -8 8 -8 8 0 15 units box open 6
fix		2 all wall/gran/region hertz/history 800.0 200.0 5.0 5.0 0.5 1 region reg2

region		pourinregion block -5 5 -5 5 5 12 units box
fix		ins all pour 10 0 123456789 id next vol 0.8 10 &
		region pourinregion mol hexapod rigid hexarigid

compute		keT all ke/rigid hexarigid
compute		erT all erotate/rigid hexarigid

thermo_style	custom step atoms c_keT c_erT
thermo_modify	lost warn norm no

thermo		100

variable	colors string "red green blue yellow white"
variable	mol2 atom mol%5
dump		2 all image 100 drop_*.jpg v_mol2 type zoom 1.5
dump_modify	2 pad 5 amap 0 5 sa 1 5 ${colors}

run		2000