【已解决】Fortran77老代码中的do循环，在Fortran95中能否使用？

前言

老师给我了一个Fortran77的老代码，这段代码在他的电脑上可以顺利运行，在我这里就会出现各种各样的奇怪问题，考虑到老代码新编译器没法好好运行，所以我想把它改写成Fortran95的，但是里面有一些奇怪的语法不知道怎么改。

---

一、奇怪的运行错误

在导师电脑上可以成功运行，在我这里就会变成NaN

奇怪的是，我想看看哪里开始出错的，加了一句‘write(6,*) d(20,20,20)’在一个函数之前，想看看这个值的变化，这个值和加这一句的位置是随手选的，奇怪的是，不再出现NaN了，虽然结果是错的，但是能够顺利运行下去了，而且换个位置，或者改成‘write(5)’，输出去别处，NaN就又出现了。

二、这是什么写法？

代码如下：

      ixjx=ix*jx
      ixjxkx=ix*jx*kx

      do