Fortran三维数组t初始化异常：输出全为0和NaN的问题求助

我之前调试Fortran老代码的时候也碰到过类似的情况，尤其是处理USGS这类遗留开源脚本时，很容易踩数组初始化的坑。咱们拆解下这个问题的核心：三维数组t打印出0和NaN，但一维数组x正常，本质和Fortran的数组初始化规则、内存使用逻辑直接相关，具体原因和解决办法如下：

可能的原因 #

• 未显式初始化的内存残留：
  Fortran对局部数组（比如子程序内声明的数组）默认不会自动初始化任何值——数组会直接占用内存里的“脏数据”。一维数组x(nx0)只有101个元素，刚好占用了之前内存里相对干净的区域（比如其他程序执行后留下的0）；而三维数组t有101401=4040个元素，尺寸更大，覆盖到了之前程序运算残留的NaN区域（比如曾经出现过除以0、无效浮点运算产生的NaN，未被清理）。

• 数组存储布局的影响：
  Fortran是列优先存储，三维数组t(nx0, nd0, nindex)的存储顺序是先填满nx0维度，再填nd0，最后nindex。如果内存某块区域刚好留存了NaN，三维数组的大尺寸更大概率会覆盖到这块区域，而一维数组可能刚好避开了。

• 作用域与编译器行为差异：
  如果x是模块级的全局数组，Fortran会默认初始化为0；但如果t是子程序的局部数组，就没有这个默认初始化逻辑——这也是两类数组表现不同的常见原因。

解决办法 #

针对这个问题，最直接有效的方式是显式初始化数组，有几种适配不同Fortran版本的写法：

1. 声明时直接初始化（适配Fortran90及以上）：

   real t(nx0, nd0, nindex) = 0.0
   

2. 代码内整体赋值（适配Fortran90及以上）：

   t = 0.0  ! Fortran支持数组整体赋值，无需循环
   

3. 循环赋值（适配老版本Fortran，比如Fortran77）：

   do k = 1, nindex
     do j = 1, nd0
       do i = 1, nx0
         t(i,j,k) = 0.0
       end do
     end do
   end do
   

另外，如果NaN的出现是因为之前的代码存在非法浮点运算（比如除以0、对数输入负数等），可以在编译时开启浮点异常检测（比如GCC的-ffpe-trap=invalid,zero,overflow选项），定位到产生NaN的代码，从根源解决问题。

内容的提问来源于stack exchange，提问作者zigzs