PASSKEy 使用 DVODPK 出现 ISTATE=-1、MXSTEP=500，程序继续运行，如何判断结果有效性

asdfghhjkl

各位老师好，我在使用 PASSKEy 进行 二维 DBD 仿真时，遇到化学动力学积分报错，想请教原因及处理方法。
1．模型与计算环境

• 二维笛卡尔坐标下的同轴反应器截面；
• 气隙中加入 12 个介质颗粒，颗粒半径 0.42 mm，相对介电常数为 4，最小几何间隙 50 μm；
• 不考虑催化表面反应；
• Intel Fortran / oneAPI 2023，OpenMP 8 核；
• 使用 LMEA，相关设置如下：
• Kinetic_implicit=1
• Kinetic_ILFA=0
• Kinetic_asynchronous=0
• Ebarlimit=100.d0
• Time_end=5.d-9
• 2．报错内容
• 在约 0.398 ns 时，多个网格单元出现：
  
• DVODPK- At current T (=R1), MXSTEP (=I1) steps
          taken on this call before reaching TOUT
  In above message, I1 = 500
  In above message, R1 = 0.3978834499471D-09
  ODE solver failed at cell 53205
  ISTATE = -1
  
• 报错后程序未退出，仍继续更新电场、粒子密度并输出结果，
  
• 停止作业前，日志累计有 164 条 ODE solver failed 记录。
  
• 3．目前排查情况
  编译时链接了：
  bin_vodpk.obin_dvode_openmp.obin_Kinetics.obin_Timestep.o通过 strings 确认，ODE solver failed 报错文本位于 bin_Kinetics.o。