# used for offical lammps
#========================================================================================
variable skindistance equal 0.0001 #say 10% of max particles
variable interval_fix equal 1000
variable interval_dump equal 10000
#========================================================================================
units		si # (This means m, kg, s for base units.)
boundary	p p f  # (x and y are p; z is the rigid wall)
atom_style	bpm/sphere
dimension	3
newton		on off
echo            both
comm_modify     vel yes 
processors      * * * 
#========================================================================================
read_data       oneclump12.data
#========================================================================================
neighbor	${skindistance} bin 
group rubber type 1
group sand   type 2
neigh_modify exclude molecule/intra all  delay 0 every 1 one 200 page 2000
bond_style      bpm/rotational
bond_coeff      1 2e6 4e5 8e4 8e4 1.0e12 1.0e12 1.0e12 1.0e12 0.1 0.02 0.004 0.004 
special_bonds   lj 0 1 1 coul 1 1 1
#========================================================================================
pair_style granular 
pair_coeff * * hertz/material 2.0e7 0.1 0.5 tangential mindlin_rescale/force NULL 1.0 0.001
timestep	1e-3
run_style	verlet
#========================================================================================
compute     4 all property/atom id radius x y z vx vy vz fx fy fz 
compute     11 all bond/local dist force fx fy fz dx dy dz
compute     12 all property/local batom1 batom2 cutoff radius
#======================================================
fix 3 all wall/gran hooke/history 200000000.0 NULL 50.0 NULL 0.5 0 zplane -0.00501 NULL contacts
fix 4 all wall/gran hooke/history 200000000.0 NULL 50.0 NULL 0.5 0 zplane NULL 0.01261 wiggle z 0.1 0.1 contacts
#========================================================================================
dump        2 all custom ${interval_dump} r100d1_atom/dump*.atom id x y z diameter c_4[9] c_4[10] c_4[11]
dump        6 all local  ${interval_dump} r100d1_bond/dump*.bond c_12[1]  c_12[2] c_11[1] c_11[2] c_11[3] c_11[4] c_11[5] c_11[6] c_11[7] c_11[8] 
dump        7 all custom ${interval_dump} r100d1_contact/dump*.contact f_4[1] f_4[2] f_4[3] f_4[4] f_4[5] f_4[6] f_4[7] f_4[8]
#=========================================================================================
run		    999999999
#=========================================================================================