能否向Fortran的INTENT(INOUT)子例程传递维度未定的数组？

你遇到的错误核心原因是：[k,i1,i2]是数组构造器表达式，它生成的是一个临时匿名数组，属于"右值"（只能被读取，无法被赋值修改）。而子例程中intent(inout)的形参要求实参必须是可修改的变量实体（比如命名数组、可分配数组等），临时数组无法满足这个要求，因此编译器报错。

方案1：使用命名数组中转（兼容性最好） #

先把需要修改的变量打包到一个命名数组里，传递给子例程后再将修改后的值拆回原变量，示例代码：

主程序 #

program main
  integer :: k, i1, i2
  integer, dimension(3) :: int_array  ! 声明与变量数量匹配的数组

  ! 将变量值写入数组
  int_array = [k, i1, i2]
  ! 调用子例程修改数组
  call load_checkpoint(int_array)
  ! 将修改后的值写回原变量
  k = int_array(1)
  i1 = int_array(2)
  i2 = int_array(3)
end program main

子例程（无需修改） #

subroutine load_checkpoint(ints)
  integer, intent(inout) :: ints(:)
  ! 读取外部文件并修改变量的逻辑示例
  open(unit=10, file='checkpoint.txt', status='old')
  read(10, *) ints  ! 假设文件中整数个数与数组长度一致
  close(10)
end subroutine load_checkpoint

方案2：动态分配数组（适配未知数量的整数） #

如果文件中整数个数n不确定，可以使用可分配数组，让子例程先读取n再分配空间，这种方式通用性更强：

主程序 #

program main
  integer, allocatable :: int_array(:)
  
  call load_checkpoint(int_array)
  
  ! 使用修改后的数组，比如输出验证
  if (allocated(int_array)) then
    print *, "读取到的整数：", int_array
  end if
end program main

子例程（调整为intent(out)并处理分配） #

subroutine load_checkpoint(ints)
  integer, intent(out), allocatable :: ints(:)
  integer :: n, i
  
  open(unit=10, file='checkpoint.txt', status='old')
  read(10, *) n  ! 先读取文件中的整数个数
  allocate(ints(n))
  read(10, *) (ints(i), i=1,n)  ! 读取所有整数
  close(10)
end subroutine load_checkpoint

注意事项 #

• 确保外部文件的格式与读取逻辑匹配：比如方案1中文件需包含与数组长度一致的整数；方案2中第一行需是整数个数，后续行是对应数量的整数。
• 若使用方案2，需确保编译器支持Fortran 2003及以上标准（主流编译器如GCC、Intel Fortran均支持）。

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