#author: Jayson S. Bazar
#Intruder impacting a granular media

atom_style       sphere
atom_modify      map array
dimension        2
boundary         p p p
newton           off
comm_modify      vel yes
units            si


# create geometry

lattice          hex 0.7
region           box block 0 0.3  0 0.3 -0.001 0.001 units box
create_box       2 box

neighbor         0.003 bin
neigh_modify     delay 0 every 1 check no

#atom regions

region           media block 0.1 INF INF INF INF INF

create_atoms     1 single 0 0.15 0
create_atoms     2 region media 

#set

set              type 1 diameter 0.05 
set              type 2  

# potentials 

pair_style       gran/hertz/history 20000.0 NULL 50.0 NULL 0.5 1
pair_coeff	 * *

#groups

group            1 type 1
group            2 region media
group            2 type 2




fix		 1 all nve/sphere
fix              2 1 drag 0.15 0.25 0.0 1.0  0.0
fix              3 viscous 0.001
fix              4 all enforce2d
#fix             4 all lb/fluid 1 1 0.5 0.0009
#fix             5 all lb/viscous



compute          1 all erotate/sphere

thermo		 1000

dump		 1 all atom 400 dump.mod_intruder


dump		 2 all image 400 image.*.jpg type type &
		 zoom 1.6 adiam 1.5

dump_modify	 2 pad 5

dump		 3 all movie 100 movie.mpg type type &
		 zoom 1.6 adiam 1.5
dump_modify	 3 pad 5

run		 6000