这问题有很多种可能，是比较麻烦的一个运行时错误了，需要很高的耐心去解决。

可能是未分配的数组，直接传递给了子程序，而子程序直接开始使用。

也可能是其他原因。

你可以认真检查 assemb_norm.for 的 41 行，这一行使用了哪些数组。
以及 nbody.f90 的 263 行，这一行调用了哪个函数，传入了哪些参数？