#Header
atom_style      granular
atom_modify     map array
hard_particles  yes
boundary        f f f
newton          off
communicate     single vel yes

units           si

neighbor        0.01 bin
neigh_modify    delay 0


#Simulation box
region      reg block -0.5 0.5 -0.5 0.5 0. 0.5 units box
create_box  2 reg


#Material properties
fix    m1  all  property/global youngsModulus peratomtype 1.e7 1.e7
fix    m2  all  property/global poissonsRatio peratomtype 0.4 0.4
fix    m3  all  property/global coefficientRestitution peratomtypepair 2 0.3 0.3 0.3 0.3
fix    m4  all  property/global coefficientFriction peratomtypepair 2 0.4 0.4 0.4 0.4
fix    m5  all  property/global coefficientDynamicFriction peratomtypepair 2 0.4 0.4 0.4 0.4


#Contact model
pair_style gran model hertz tangential history
pair_coeff	* *


#Time step
timestep    0.00001


#Gravity
fix    gravi  all  gravity 9.81 vector 0.0 0.0 -1.0


#Walls
fix   wall_x1  all  wall/gran model hertz tangential history primitive type 1 xplane -0.5
fix   wall_x2  all  wall/gran model hertz tangential history primitive type 1 xplane 0.5
fix   wall_y1  all  wall/gran model hertz tangential history primitive type 1 yplane -0.5
fix   wall_y2  all  wall/gran model hertz tangential history primitive type 1 yplane 0.5
fix   wall_z1  all  wall/gran model hertz tangential history primitive type 1 zplane 0.0
fix   wall_z2  all  wall/gran model hertz tangential history primitive type 1 zplane 0.5


#Spheres/Multispheres groups
group spheres initialize
group multispheres initialize


#Integrators
fix integr_sphere       spheres       nve/sphere
fix integr_multisphere  multispheres  multisphere


#Templates for distributions
fix    pts_sphere       spheres       particletemplate/sphere      1028283577  atom_type 1 density constant 2500 radius constant 0.02
fix    pts_multisphere  multispheres  particletemplate/multisphere 1192794983  atom_type 2 density constant 2500 nspheres 3 ntry 1000000  &
                                      spheres file cluster.multisphere scale 0.001 type 1


# Different distributions (activate ONLY one by commenting lines)
#fix    pdd  all  particledistribution/discrete 1207956223 1 pts_sphere 1.0
#fix    pdd  all  particledistribution/discrete 1207956223 1 pts_multisphere 1.0
fix    pdd  all  particledistribution/discrete 1207956223 2 pts_sphere 0.5 pts_multisphere 0.5

#Insertion
region ins_reg cylinder z 0.0 0.0 0.3 0.1 0.4 units box
fix    ins  all  insert/pack seed 2137351217 distributiontemplate  pdd  &
                 maxattempt 100  insert_every once  overlapcheck yes  all_in yes  &
                 vel constant 0. 0. -0.1  region ins_reg  particles_in_region 200


#Insert and run
run   1
dump  dmp  all  custom/vtk 1000 post/test_*.vtk  id type x y z vx vy vz fx fy fz omegax omegay omegaz radius
run   100000