SPECFEM3D,poroelastic,stacey absorbing condition

I am trying to run specfem3D for an example for homogenous_poroelastic .This example originally has no absorption boundary, and I can successfully compile and run it. Now I want to use the stacey absorbing condition, it made some mistakes,

Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
Backtrace for this error:
#0 0x7f67ed4f2dfd in ???
#1 0x7f67ed4f2013 in ???
#2 0x7f67ed15020f in ???
#3 0x557c4dfb3c49 in compute_forces_poro_solid_part_
at src/specfem3D/compute_forces_poro_solid_part.f90:342
#4 0x557c4dfb022f in compute_forces_poroelastic_calling_
at src/specfem3D/compute_forces_poroelastic_calling_routine.f90:56
#5 0x557c4e01ebc6 in iterate_time_
at src/specfem3D/iterate_time.F90:253
#6 0x557c4df6fcd6 in xspecfem3d
at src/specfem3D/specfem3D.F90:394
#7 0x557c4df6fcd6 in main
at src/specfem3D/specfem3D.F90:32
Floating point exception (core dumped)

any suggestions?
Thanks and Kind Regards