LAMMPS (5 Oct 2016)
# Pour 2d granular particles into container

dimension	2
atom_style	sphere
boundary	f fm p
newton		off
comm_modify	vel yes

region		reg block 0 100 0 50 -0.5 0.5 units box
create_box	1 reg
Created orthogonal box = (0 0 -0.5) to (100 50 0.5)
  1 by 1 by 1 MPI processor grid

neighbor	0.2 bin
neigh_modify	delay 0

# IMPORTANT NOTE: these values of stiffness (4000) and timestep (0.001)
# are used in this example file to produce a quick simulation and movie.
# More appropriate values for realistic simulations would be
# k = 4.0e5 and dt = 0.0001, as in bench/in.chute (if it were Hertzian).

pair_style      gran/hertz/history 4000.0 NULL 100.0 NULL 0.5 0
pair_coeff	* *

timestep	0.001

fix             1 all nve/sphere
fix		2 all gravity 1.0 spherical 0.0 -180.0
fix		xwalls all wall/gran hertz/history 4000.0 NULL 100.0 NULL 0 1 		xplane 0 100
fix		ywalls all wall/gran hertz/history 4000.0 NULL 100.0 NULL 0 1 		yplane 0 NULL

region          slab block 1.0 99.0 30 34.5 -0.5 0.5 units box
fix             ins all pour 1000 1 4767548 vol 0.4 10 	        diam range 0.5 1.0 region slab
Particle insertion: 224 every 3000 steps, 1000 by step 12001

fix             3 all enforce2d

compute		1 all erotate/sphere
thermo_style	custom step atoms ke c_1 vol
thermo		1000
thermo_modify	lost ignore norm no
compute_modify	thermo_temp dynamic yes

#dump		id all atom 250 dump.pour

#dump		2 all image 500 image.*.jpg type type #		zoom 1.6 adiam 1.5
#dump_modify	2 pad 5

#dump		3 all movie 500 movie.mpg type type #		zoom 1.6 adiam 1.5
#dump_modify	3 pad 5

run		25000
Neighbor list info ...
  2 neighbor list requests
  update every 1 steps, delay 0 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 1.2
  ghost atom cutoff = 1.2
  binsize = 0.6 -> bins = 167 84 2
Memory usage per processor = 2.75186 Mbytes
Step Atoms KinEng c_1 Volume 
       0        0           -0            0         5000 
    1000      224    201.77464            0         5000 
    2000      224    372.42868            0         5000 
    3000      224    599.60221            0         5000 
    4000      448    1084.3752            0         5000 
    5000      448    1592.4543            0         5000 
    6000      448    1763.3153    37.274939         5000 
    7000      672     1805.206    89.331853         5000 
    8000      672    1778.0015    111.58381         5000 
    9000      672    1592.6805    97.550311         5000 
   10000      896     1631.962    76.905078         5000 
   11000      896    1590.1527    68.644552         5000 
   12000      896    1564.4707    60.119562         5000 
   13000     1000     1490.072    56.982972         5000 
   14000     1000    1266.4669    57.579371         5000 
   15000     1000    1034.1055    46.079374         5000 
   16000     1000    748.70283    43.933154         5000 
   17000     1000    462.65166    37.522016         5000 
   18000     1000    327.40481    29.450285         5000 
   19000     1000    223.24319    31.981144         5000 
   20000     1000    54.792042    23.449101         5000 
   21000     1000    25.805484    16.352528         5000 
   22000     1000    20.617765    11.756965         5000 
   23000     1000    16.115382    7.9435503         5000 
   24000     1000    11.971003    5.9411126         5000 
   25000     1000    10.404117    4.3811155         5000 
Loop time of 3.00132 on 1 procs for 25000 steps with 1000 atoms

Performance: 719682.852 tau/day, 8329.663 timesteps/s
99.9% CPU use with 1 MPI tasks x no OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 1.2554     | 1.2554     | 1.2554     |   0.0 | 41.83
Neigh   | 0.4288     | 0.4288     | 0.4288     |   0.0 | 14.29
Comm    | 0.0078034  | 0.0078034  | 0.0078034  |   0.0 |  0.26
Output  | 0.0003016  | 0.0003016  | 0.0003016  |   0.0 |  0.01
Modify  | 1.177      | 1.177      | 1.177      |   0.0 | 39.21
Other   |            | 0.132      |            |       |  4.40

Nlocal:    1000 ave 1000 max 1000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    2374 ave 2374 max 2374 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 2374
Ave neighs/atom = 2.374
Neighbor list builds = 2089
Dangerous builds = 0
Total wall time: 0:00:03
