# Pour granular particles into chute container, then induce flow
units si
#read_restart restart.519000.bed
#read_dump dump.pourCohe.519000 519000 x y
atom_style	 sphere
boundary	p fm p
newton          off
communicate     single vel yes
#comm_modify vel yes
dimension 2
region          reg block -1 1 -1 1 -0.1 0.1 units box
create_box	 2 reg
#mass            2 0.01
#
neighbor        0.2 bin
neigh_modify    delay 0 page 80000000 one 800000
#

#
pair_style	gran/hooke/history 20000.0 NULL 2500.0 1000 0.5 1
pair_coeff	* *
#
timestep        0.0001
#
fix             1 all nve/sphere
fix             2 all gravity 9.8 vector 0 -1 0
fix             ylower all wall/gran 5000.0 NULL 10000.0 NULL 0.5 1 yplane -1.0 -0.5
#               ID group-ID wall/gran Kn Kt gamma_n gamma_t xmu dampflag wallstyle args keyword values
region         slab1 block -1.0 1.0 -0.5 0.11 -0.10 0.10 units box
region         slab2 block -1.0 1.0 0.55 0.60 -0.1 0.1 units box
region          bedding block -1.0 1.0 -0.5 -0.4 -0.1 0.1 units box
##
group shen type 1
group salt type 2
group aab region slab1
group bed region bedding
#
#
compute         1 all erotate/sphere
thermo_style    custom step atoms ke c_1 vol temp
thermo          1000
thermo_modify   lost ignore norm no
compute_modify  thermo_temp dynamic yes
#
#velocity bed set 0.0 NULL NULL sum no units box
#
fix 5 all enforce2d
fix damp aab viscous 0.0025
fix stop bed freeze
#
shell mkdir dump
shell cd dump
dump     	        mydmp all custom 100 dump.pourCohe.* id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
#
restart 1000 restart.*.bed
#run 10000
unfix damp
# 
fix             ins21 all pour 6000 2 10719 region slab2 diam range 0.0008 0.0008 dens 2800 2800 vel 5.26 5.26 -3.04
restart 1000 restart.*.bed
run 400000