# initialization and equlibration
include init_equil.in
#================================================================================
#The total simulation box is (-15,-5,-15) to (15,30,15) with periodic boundaries
#in all directions then I create a wall with 5 thickness (lj units) and a particle
#with diamete of 15 (lj units), I thermostat the fluid which is lj with 0.81 density
#================================================================================
reset_timestep		0
undump		2
unfix		1
unfix		2

region		wall block INF INF INF 0.0 INF INF units box
group		wall region wall
set		group wall type 2

region		particle sphere 0.0 12.5 0.0 7.50 side in units box
group		particle region particle
set		group particle type 3
velocity	particle set 0.0 0.0 0.0 units box 

group		flow subtract preset wall particle
group		mobile union wall flow

compute		pre2 flow temp/partial 1 1 1

fix		3 particle rigid single  
neigh_modify	exclude type 3 3

fix             1 mobile nve
fix             2 flow temp/berendsen 1.1 1.1 100
fix_modify      2 temp pre2

velocity	wall set 0.0 0.0 0.0 units box
fix		4 wall setforce 0.0 0.0 0.0

thermo_style    custom step temp c_pre2 etotal press vol
thermo		1

restart 1000000 3rd_restart.*
dump 	3 all atom 100000 rigid_particle.lammpstrj 

run	1000000
#==================================================================================================================================