#### To make the error appears, the simulation must be computed with a multi-processors run (mpirun -n n with n>=2) ####

# Atom style (regardless the sub-styles defined in hybrid, this is the "hybrid + fix adapt" that make the problem occurs)
atom_style		hybrid granular
atom_modify		map array

# Fix boundaries
boundary		f f f
newton			off
communicate		single vel yes
processors		* * *
units			si
timestep		0.00005

# Simulation domain and atom types
region			reg block -0.05 0.05 -0.05 0.05 0. 0.05 units box
create_box		1 reg

# Bin size and check for a neighbor update every step
neighbor 		0.001 bin
neigh_modify	delay 0

# Add gravity
fix				grav	all gravity 9.81 vector 0.0 0.0 -1.0

# Pair style
pair_style		gran model hertz tangential history
pair_coeff		* *

# Material properties
fix 			m1		all property/global youngsModulus peratomtype 5.e6
fix 			m2		all property/global poissonsRatio peratomtype 0.3
fix 			m3		all property/global coefficientRestitution peratomtypepair 1 0.25
fix 			m4		all property/global coefficientFriction peratomtypepair 1 0.1

# Walls
fix				wall_x1	all wall/gran model hertz tangential history primitive type 1 xplane -0.05
fix				wall_x2	all wall/gran model hertz tangential history primitive type 1 xplane 0.05
fix				wall_y1	all wall/gran model hertz tangential history primitive type 1 yplane -0.05
fix				wall_y2	all wall/gran model hertz tangential history primitive type 1 yplane 0.05
fix				wall_z1	all wall/gran model hertz tangential history primitive type 1 zplane 0
fix				wall_z2	all wall/gran model hertz tangential history primitive type 1 zplane 0.05

# Particles template
fix				pts1	all  particletemplate/sphere  1939841531  atom_type 1  density constant 2500 radius constant 0.001

# Distribution of particles
fix				pdd1	all particledistribution/discrete 32452843  1 pts1 1.0

# Insertion region
region			ins_reg	block -0.05 0.05 -0.05 0.05 0. 0.05 units box

# Particle insertion
fix				ins		all insert/pack seed 32452867 distributiontemplate pdd1 &
				maxattempt 10 insert_every once overlapcheck yes &
				all_in yes vel constant 0.0 0.0 -0.0 region ins_reg &
				particles_in_region 1000

# Integrate using velocity verlet
fix				integr	all nve/sphere

# Dumping
run				1
#dump			dmpAtoms	all  custom/vtk  1000  post/particles_*.vtk  id  type  x y z  radius


##### Particles adapt command #####

compute			rad			all property/atom	radius
variable		dgrown		atom 1.001*2.*c_rad
# Remove the fix command and the issue disappears...
fix				grow		all adapt 50 atom diameter v_dgrown

run				10000