### Restart in.shock_2D file
read_restart shock.restart.150000


### Preliminaries ###
#units				si
#atom_style 	granular
#boundary 		p p f
newton 			off
communicate single vel yes
hard_particles yes


### DECLARE DOMAIN ###
## For 3D ##
region reg block -0.1 0.1 -0.01 0.01 0. 0.5 units box
#create_box 	2 reg
neighbor 		0.002 bin
neigh_modify delay 0


### MATERIAL PROPERTIES ###
fix 		m1 all property/global youngsModulus peratomtype 5.e6	5e99
fix 		m2 all property/global poissonsRatio peratomtype 0.35		0.35
fix 		m3 all property/global coefficientRestitution peratomtypepair 2 0.9 0.9 0.9 0.9
fix 		m4 all property/global coefficientFriction peratomtypepair 2 0.3 0.3 0.3 0.3
fix 		m5 all property/global characteristicVelocity scalar 2.
fix 		m6 all property/global cohesionEnergyDensity peratomtypepair 2 100000 0.1 0.1 0.1


### PAIR STYLE ###
pair_style gran model hertz tangential history cohesion sjkr #Hertzian with cohesion
pair_coeff 		* *


### PARTICLE DISTRIBUTIONS ###
## Soil
fix 				soil all particletemplate/sphere 15485863 atom_type 1 density constant 2500 radius constant 0.0015
fix 				soil_distribution all particledistribution/discrete 15485867 1 soil 1.0

## Impact
fix 				impact all particletemplate/sphere 32452843 atom_type 2 density constant 5000 radius constant 0.005
fix 				impact_distribution all particledistribution/discrete 32452867 1 impact 1.0

fix 1 				all nve/sphere
fix 2 				all gravity 5 vector 0.0 0.0 -1.0

### BOX WALLS ###

fix 				boxwalls_x1 all wall/gran model hertz tangential history primitive type 1 xplane -0.1
fix 				boxwalls_x2 all wall/gran model hertz tangential history primitive type 1 xplane +0.1
fix 				boxwalls_y1 all wall/gran model hertz tangential history primitive type 1 yplane -0.01
fix 				boxwalls_y2 all wall/gran model hertz tangential history primitive type 1 yplane +0.01
fix 				boxwalls_z1 all wall/gran model hertz tangential history primitive type 1 zplane 0.0


### PARTICLE INSERTION ###

## Region of insertion ##

group			soil1	type 1
group			impact1 type 2
group			layer1	type 1

region 			insertion_soil1 block -0.1 0.1 -0.01 0.01 0.4 0.5  units box
region			insertion_impact1 block -0.02 0.02 -0.01 0.01 0.4 0.5  units box
region			layer1 block -0.05 0.05 -0.005 0.005 0 0.1 units box
 
#fix ins_soil soil1 insert/pack seed 32452843 distributiontemplate soil_distribution vel constant 0. 0. 0. insert_every 5000 overlapcheck yes all_in yes particles_in_region 2000 region insertion_soil1


### THERMO SETTINGS ###

thermo_style custom step atoms ke vol
thermo 1000
thermo_modify lost ignore norm no


### COMPUTE ###
compute		soilStress	soil1	stress/atom
compute 		regionStress layer1 stress/atom


### EXECUTION and FURTHER SETTINGS ###
dump soil_data soil1 custom/vtk 500 post_Restart_2D/soil1_*.vtk id x y z vx vy vz fx fy fz c_soilStress[1] c_soilStress[2] c_soilStress[3] c_regionStress[1] c_regionStress[2] c_regionStress[3]
dump impact_data impact1 custom/vtk 500 post_Restart_2D/impact1_*.vtk id type x y z vx vy vz radius fx fy fz


#run 50000

# IMPACT

velocity all set 0. 0. 0.
timestep		0.000003
fix 				ins_impact impact1 insert/pack seed 49979693 distributiontemplate impact_distribution vel constant 0. 0. 0. insert_every 10000 overlapcheck yes all_in yes particles_in_region 1 region insertion_impact1
fix				4	soil1 gravity 1 vector 0.0 0.0 -1.0
fix				5	impact1 gravity 5 vector 0.0 0.0 -1.0
run 				300000

log log.liggghts_Jake