hello everyone hope ur all well i am stuck on a problm i am trying to make a script to push or attract particles using an electrostatic field but i have a problme the same error keeps popping and i could not resolve the problm i know it is not a big problm for some of you but it is for me and i would apreciate any tips on how to improuve my script to reach my goal it is far from finnished thank you all
ERROR: Numeric index is out of bounds (../force.cpp:786)
#falling particles in electric field
atom_style hybrid charge granular
atom_modify map array
boundary f f f
newton off
communicate single vel yes
units si
region domain block -0.555 0.55 -0.15 0.15 -0.5 0.535 units box
create_box 1 domain
lattice sc 0.001
region factory block -0.14 0.10 -0.006 0.006 -0.002 0.205 units box
create_atoms 1 region factory
neighbor 0.002 bin
neigh_modify delay 0
#Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5
#New pair style
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *
timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
#distributions for insertion
fix pts all particletemplate/sphere 15485863 atom_type 1 density constant 1000 radius constant 0.004
fix pdd all particledistribution/discrete 15485867 1 pts 1.0
fix ins all insert/rate/region seed 32452843 distributiontemplate pdd &
nparticles 1000 massrate 1.66666667 insert_every 1000 &
overlapcheck yes vel constant 0. 0. -1. region factory ntry_mc 10000
# Geometry
fix wal11 all wall/gran model hertz tangential history primitive type 1 zplane 0.0
fix wall2 all wall/gran model hertz tangential history primitive type 1 xplane -0.15
fix wall3 all wall/gran model hertz tangential history primitive type 1 xplane 0.15
fix wall4 all wall/gran model hertz tangential history primitive type 1 yplane -0.01
fix wall5 all wall/gran model hertz tangential history primitive type 1 yplane 0.01
# Applied flow field
region electro block -0.15 -0.05 -0.01 0.01 0.1 0.2 units box
fix electro all efield 0 50000 0
set atom all charge 0.000001
#apply nve integration to all particles that are inserted as single particles
fix integr all nve/sphere
# Thermo settings
thermo_style custom step atoms ke cpu
thermo 10000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
# Run 1 step to check timestep and initialize
fix ctg all check/timestep/gran 1 0.01 0.01
run 1
unfix ctg
# Dump particle positions
dump dmp all custom/vtk 100 post/fall*.vtk id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
### Execution and further settings
# Run remainder
run 200000 upto