gfortran编译已初始化数组报may be used uninitialized警告问题

问题背景 #

编译未改动的Fortran代码时，gfortran抛出-Wmaybe-uninitialized警告。代码逻辑为：局部可分配数组完成内存分配后先整体赋值为0，再将传入的rarray数组值拷贝到该数组的对应切片位置。
相关代码片段：

SUBROUTINE advection(rbox0,rfrac,rarray,na)
  IMPLICIT NONE

  INTEGER,  INTENT(IN)   :: na
  REAL, INTENT(IN)   :: rfrac 
  REAL, INTENT(IN)   :: rbox0 
  REAL, INTENT(INOUT) :: rarray(0:na)

  ! local arrays
  REAL, ALLOCATABLE :: zarray1(:),zarray2(:)


  ALLOCATE(zarray1(-4:2*na))
  zarray1(:)=0.0
  zarray1(0:na)=rarray ! zarray stores the original array
  rarray(:)=0.0
  rarray(0)=(1.0-rbox0)*zarray1(0)

  ! etc etc 

编译环境为Mac平台Homebrew安装的GNU Fortran (Homebrew GCC 11.3.0_1) 11.3.0版本，抛出的警告信息如下：

vectri.f90:1101:38:

 1101 |       rarray(0)=(1.0-rbox0)*zarray1(0)
      |                                      ^
Warning: 'MEM <real(kind=4)[0:]> [(real(kind=4)[0:] *)_77][4]' may be used uninitialized [-Wmaybe-uninitialized]

已尝试的调整：

• 最初使用无切片的全数组赋值写法zarray1=0.0和rarray=0.0，警告出现后改为带(:)的切片全量赋值写法，警告未消除。
• 这类未初始化警告通常对应真实的内存访问问题，不能直接忽略，但无法定位编译器触发警告的根因。

根因说明 #

该警告是GCC 11.x版本gfortran前端的静态分析误报，代码本身不存在未初始化内存访问：

• gfortran 11版本的静态分析器在追踪可分配数组跨切片赋值的内存状态时存在缺陷，没有正确识别zarray1(0:na)=rarray语句已经完成了zarray1(0)位置的赋值，错误判定该偏移位置的内存仍为未初始化状态。
• 无论是全数组赋值还是带(:)的切片赋值，都没有触发旧版本分析器的内存状态更新逻辑，因此调整写法无法消除警告。

处理方式 #

1. 先确认代码实际无运行时问题：编译时添加-fcheck=mem -g选项，运行编译后的程序。如果真的存在未初始化内存访问，Fortran运行时会直接抛出明确的内存错误；如果程序正常运行无报错，即可确认是编译器误报。
2. 消除警告的可选方案：
   • 升级GCC版本到12及以上，该版本修复了大量Fortran前端的未初始化检测误报，该场景下不会再抛出错误警告。
   • 若暂时无法升级编译器，可在确认逻辑正确的前提下，对zarray1的目标访问段增加一次显式赋值，或使用编译诊断pragma屏蔽对应代码段的-Wmaybe-uninitialized警告。
   • 分配数组时直接通过数组构造器完成全量初始化，例如zarray1 = [ (0.0, i = -4, 2*na) ]，部分场景下可以让旧版本分析器正确识别初始化状态。

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