Fortran子程序大数组疑问：-CB编译报错，声明wa(1)为何正常？

你遇到的核心问题是Fortran中64位超大数组与32位边界检查机制的兼容性冲突，下面一步步拆解原因和解决办法：

为什么wa(1_8:1_8*n1*n2)会触发边界检查错误？ #

当你用-CB选项启用数组边界检查时，ifort的部分版本在处理64位数组下标时，内部边界检查的存储逻辑可能默认使用32位整数：

• 你的数组大小是2^11 * 2^20 = 2^31，也就是2147483648，这个值超过了32位有符号整数的最大值2^31-1（即2147483647）。
• 当你在子程序中声明wa(1_8:1_8*n1*n2)时，64位的上界2147483648会被截断为32位有符号整数，结果是-2147483648（二进制补码的特性）。
• 后续你计算的下标ms是64位的正数（比如错误信息里的108003329），远大于这个被截断后的负数上限，自然触发了边界检查的严重错误。

你看到的“子程序中wa数组大小：0”也是这个问题导致的——编译器无法正确识别64位的数组边界，错误计算了数组大小。

为什么wa(1)能正常运行？ #

当你把数组声明为wa(1)时，它本质上是一个假定大小数组（assumed-size array）：

• Fortran标准中，假定大小数组的上界是未知的（编译器只知道第一个维度的下界是1，上界被隐含为*），此时-CB选项的边界检查会被弱化或者跳过——编译器不会严格验证你访问的下标是否超过实际数组的大小，因此不会触发错误。
• 但这种方式是“治标不治本”，它只是避开了边界检查，并没有真正解决64位数组的类型安全问题，后续如果数组操作出现越界，很难被检测到。

推荐的解决方案 #

1. 使用假定形状数组（最安全、标准的写法） #

把子程序中的数组声明改成假定形状数组，这是Fortran 90及以后的推荐写法，编译器会自动传递数组的完整描述信息（包括64位的上下界和大小）：

subroutine sub(n1,n2,wa)
    implicit none
    integer, intent(in) :: n1, n2
    real(4), intent(inout) :: wa(:)  ! 假定形状数组
    integer(8) :: i, j, ms
    print*, 'in sub, size of wa:', size(wa,kind=8)
    ms=0
    !$omp parallel default(shared) private(i,j,ms)
    !$omp do
    do i=1, n1
        do j=1, n2
            ms=(i-1)*n2+j
            wa(ms)=ms*1.d0
        enddo
    enddo
    !$omp end do nowait
    !$omp end parallel
    print*, 'size of wa:', size(wa,kind=8)
    return
end subroutine sub

这种写法下，-CB的边界检查会正确识别64位的数组边界，不会再出现截断错误，同时保留了类型安全。

2. 强制编译器使用64位整数 #

添加编译器选项-i8，强制ifort将默认整数类型设为64位，这样数组的边界值会用64位存储，避免截断：

ifort -openmp -CB -i8 test.f

这种方式可以保留你原来的显式下标声明，但需要注意代码中所有整数类型的兼容性（比如主程序中的n1、n2如果是32位整数，可能需要调整为64位）。

3. 显式使用假定大小数组（不推荐） #

如果你坚持用类似wa(1)的写法，更规范的方式是写成wa(*)，明确表示这是一个假定大小数组：

real(4), intent(inout) :: wa(*)

但如前所述，这种方式会失去边界检查的保护，仅适合临时调试场景。

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