Segmentation fault

Please check the following script. It is running fine with atom_style sph but with hybrid sph peri it gives Segmentation fault error.

Many thanks

echo screen
processors 2 2 2
atom_style hybrid sph peri
atom_modify map array

newton off
boundary f f f
units si
variable h equal 0.0002
variable c equal 14.80 # soundspeed for Tait’s EOS
variable dt equal 0.0000001*{h}/{c} # CFL criterion
region box block 0 0.04 0 0.04 0 0.02 units box
create_box 2 box
region drop sphere 0.02 0.02 0.0065 0.00135
region water block 0.015 0.025 0.015 0.025 0.0046 0.005
lattice bcc 0.00015
create_atoms 1 region drop
lattice sc 0.00015
create_atoms 2 region water
#lattice sc 1.0
set region drop type 1
set region water type 2
group drop type 1
group water type 2
mass 1 1.75972385E-9 #1.0767E-07
mass 2 1.76714587E-9 #1.09721E-07
set group drop sph/rho 995.8
set group water sph/rho 1000
pair_style hybrid/overlay sph/rhosum 1 sph/taitwater peri/pmb
pair_coeff * * sph/taitwater 995.8 {c} 1 {h}
pair_coeff 1 1 sph/rhosum {h}
pair_coeff 2 2 peri/pmb 6.2696893e24 0.00045001 0.00015 0.0
fix gfix water gravity 9.81 vector 0 0 -1
fix gfix drop gravity 9.81 vector 0 0 -1
velocity drop ramp vz -2.67 -2.67 z 0.0185 0.0115
velocity water set 0.0 0.0 0.0 sum yes units box
compute rho_peratom all sph/rho/atom
compute e_peratom all sph/e/atom
compute esph all reduce sum c_e_peratom
compute ke all ke
variable etot equal c_esph+c_ke+f_gfix
fix dtfix all dt/reset 1 NULL {dt} {h} units box
fix integrate_drop_fix drop sph
fix integrate_bc_fix water sph/stationary
fix F1 water nve
compute C1 water damage/atom
thermo_style custom step time atoms c_ke c_esph v_etot f_gfix press f_dtfix
thermo 10
thermo_modify lost ignore norm no
dump D1 all custom 100 …/perisph/post/dump.drop id type xs ys zs vx vy vz fx fy fz
variable skin equal 0.3*{h}
neighbor {skin} bin
neigh_modify every 1 delay 0 check no
run 5000

When reporting bugs or problems, please always report your LAMMPS version, your compilers, how you configured/installed LAMMPS and what platform you are running on.

If you compiled/installed LAMMPS from source code, please apply the following change, recompile, and try again.

diff --git a/src/PERI/atom_vec_peri.cpp b/src/PERI/atom_vec_peri.cpp
index f0cbd6eb04..e393bc28b6 100644
--- a/src/PERI/atom_vec_peri.cpp
+++ b/src/PERI/atom_vec_peri.cpp
@@ -88,12 +88,13 @@ void AtomVecPeri::grow_pointers()
 
 void AtomVecPeri::create_atom_post(int ilocal)
 {
+  const auto xinit = atom->x;
   vfrac[ilocal] = 1.0;
   rmass[ilocal] = 1.0;
   s0[ilocal] = DBL_MAX;
-  x0[ilocal][0] = x[ilocal][0];
-  x0[ilocal][1] = x[ilocal][1];
-  x0[ilocal][2] = x[ilocal][2];
+  x0[ilocal][0] = xinit[ilocal][0];
+  x0[ilocal][1] = xinit[ilocal][1];
+  x0[ilocal][2] = xinit[ilocal][2];
 }
 
 /* ----------------------------------------------------------------------
@@ -103,10 +104,11 @@ void AtomVecPeri::create_atom_post(int ilocal)
 
 void AtomVecPeri::data_atom_post(int ilocal)
 {
+  const auto xinit = atom->x;
   s0[ilocal] = DBL_MAX;
-  x0[ilocal][0] = x[ilocal][0];
-  x0[ilocal][1] = x[ilocal][1];
-  x0[ilocal][2] = x[ilocal][2];
+  x0[ilocal][0] = xinit[ilocal][0];
+  x0[ilocal][1] = xinit[ilocal][1];
+  x0[ilocal][2] = xinit[ilocal][2];
 
   if (rmass[ilocal] <= 0.0) error->one(FLERR, "Invalid mass in Atoms section of data file");
 }

Problem solved. Many thanks @akohlmey.