gfortran-10中Fortran数组deallocate后shape仍返回原维度问题咨询

核心结论 #

你观测到的是gfortran编译器对未定义行为的具体实现表现，符合Fortran语言标准约定，不属于编译器异常。 #

具体原因 #

• Fortran标准明确规定：仅当allocatable数组处于已分配状态时，shape()、size()、lbound()、ubound()等维度属性查询操作才合法。数组被deallocate释放后，调用上述查询函数的行为属于未定义行为，编译器可返回任意值、也可触发运行时报错，不同编译器的实现逻辑不受约束。
• gfortran 10的底层实现中，allocatable数组的描述符（存储数组维度、各维长度、数据指针等元信息的结构）在执行deallocate时，只会将「已分配标记位」设为.false.、释放数据段内存，不会主动清空描述符中存储的维度长度字段，因此调用shape(array)时会读到之前分配的旧值。
• 若切换其他编译器（如Intel Fortran），相同操作大概率会直接触发运行时崩溃或返回乱码值，都属于未定义行为的正常表现。

最佳实践 #

所有对allocatable数组的维度属性查询，都必须前置allocated(array)判断，只有确认返回.true.时再执行后续查询和操作，避免依赖未定义行为带来的兼容性问题。

! 正确写法示例
if (allocated(array)) then
    print *, shape(array)
    ! 其他数组操作
else
    print *, "数组未分配，无法查询维度"
end if

内容的提问来源于stack exchange，提问作者mgo.fis