# Pour granular particles into chute container, then induce flow
echo both
atom_style	granular
atom_modify	map array
communicate	single vel yes

boundary	m m m
newton		off

units		si
processors	2 1 1

#read the restart file
read_restart 	../DEM/liggghts.restartCFDEM

#do not do this here, the simulation box is in the restart file!
#region		reg block -0.015 0.015 -0.015 0.015 -0.001 0.0554 units box
#create_box	1 reg

neighbor	0.0005 bin
neigh_modify	delay 0

#Material properties required for new pair styles

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

#pair style
pair_style 	gran/hertz/history #Hertzian without cohesion
pair_coeff	* *

#timestep, gravity
timestep	0.00001
fix		gravi all gravity 9.81 vector 0.0 0.0 -1.0

#walls
#walls
fix			zwalls1 all wall/gran/hertz/history   primitive type 1  zplane 0.0
fix			zwalls2 all wall/gran/hertz/history   primitive type 1  zplane 0.0553
fix         cylwalls all wall/gran/hertz/history primitive type 1 zcylinder 0.01385 0. 0. tangential_damping off

# change the particles density
set             group all density 2000

#cfd coupling
fix		cfd all couple/cfd couple_every 100 mpi
fix		cfd2 all couple/cfd/force

#apply nve integration to all particles that are inserted as single particles
fix		integr all nve/sphere

#center of mass
compute         centerOfMass all com

#compute total dragforce
#compute		dragtotal all reduce sum f_dragforce[1] f_dragforce[2] f_dragforce[3]

#screen output
compute		1 all erotate/sphere
thermo_style	custom step atoms ke c_1 vol c_centerOfMass[3] #c_dragtotal[1] c_dragtotal[2] c_dragtotal[3]
thermo		10
thermo_modify	lost ignore norm no
compute_modify	thermo_temp dynamic yes

#insert the first particles so that dump is not empty
dump 		myDump all stl 1 post/dump_*.stl
#run		1
dump		dmp all custom 5000 ../DEM/post/dump.liggghts_restart id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius 
undump		myDump

run		1