LAMMPS (14 Aug 2013)
# LAMMPS benchmark of granular flow
# chute flow of 32000 atoms with frozen base at 26 degrees

variable	x index 1
variable	y index 1

units		lj
atom_style	sphere
boundary	p p fs
newton		off
communicate	single vel yes

read_data	data.chute
  orthogonal box = (0 0 0) to (40 20 37.2886)
  2 by 1 by 2 MPI processor grid
  32000 atoms
  32000 velocities

replicate	$x $y 1
replicate	2 $y 1
replicate	2 2 1
  orthogonal box = (0 0 0) to (80 40 37.2923)
  2 by 2 by 1 MPI processor grid
  128000 atoms

pair_style	gran/hooke/history 200000.0 NULL 50.0 NULL 0.5 0
pair_coeff	* *

neighbor	0.1 bin
neigh_modify	every 1 delay 0

timestep	0.0001

group		bottom type 2
3648 atoms in group bottom
group		active subtract all bottom
124352 atoms in group active
neigh_modify	exclude group bottom bottom

fix		1 all gravity 1.0 chute 26.0
fix		2 bottom freeze
fix		3 active nve/sphere

compute		1 all erotate/sphere
thermo_style	custom step atoms ke c_1 vol
thermo_modify	norm no
thermo		100

run		100
Memory usage per processor = 21.8161 Mbytes
Step Atoms KinEng 1 Volume 
       0   128000    3136556.5    6404.5051    119335.13 
     100   128000    3137168.3    6284.3873    119338.83 
Loop time of 1.01235 on 4 procs for 100 steps with 128000 atoms

Pair  time (%) = 0.508497 (50.2295)
Neigh time (%) = 0.0454835 (4.49288)
Comm  time (%) = 0.0667624 (6.59481)
Outpt time (%) = 0.00123399 (0.121894)
Other time (%) = 0.390371 (38.5609)

Nlocal:    32000 ave 32000 max 32000 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Nghost:    5463 ave 5463 max 5463 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Neighs:    115133 ave 115133 max 115133 min
Histogram: 4 0 0 0 0 0 0 0 0 0

Total # of neighbors = 460532
Ave neighs/atom = 3.59791
Neighbor list builds = 2
Dangerous builds = 0
