variable skindistance equal 0.0001 #say 10% of max particles
variable pstart equal 0 # use 0 
variable pend equal 1000000
variable pdamp equal 1000
#========================================================================================
units		si # (This means m, kg, s for base units.)
boundary	p p p 
atom_style	bpm/sphere
dimension	3
newton		on off
echo            both
comm_modify     vel yes 
processors      * * * 
#========================================================================================
read_data       rubber_sand100small.data
#read_restart   restart_filers100d68.106400000
#========================================================================================
neighbor	${skindistance} bin 
group rubber type 1
group sand   type 2
neigh_modify exclude molecule/intra all  delay 0 every 1 one 20000000 page 2000000000
bond_style      bpm/rotational
bond_coeff      1 1.2e7 1.2e7  3e6 3e6 1.0e12 1.0e12 1.0e12 1.0e12 0.1 0.1 0.005 0.005
special_bonds   lj 0 1 1 coul 1 1 1
#========================================================================================
pair_style granular 
pair_coeff * * hertz/material 1.2e7 0.1 0.5 tangential mindlin_rescale/force NULL 1.0 0.001
timestep	1e-8
run_style	verlet
#========================================================================================
compute     7 all pressure NULL pair 
#========================================================================================
fix         9 all press/berendsen x ${pstart} ${pend} ${pdamp} y ${pstart} ${pend} ${pdamp} z ${pstart} ${pend} ${pdamp}
fix	        3 all nve/bpm/sphere
#========================================================================================
restart     100000 restart_filers100d68
#=========================================================================================
run		    999999999
#=========================================================================================