Fortran旧代码更新：数组声明从A(*)改A(:)引发段错误求助

这是Fortran新旧语法过渡时非常容易踩的坑，我来给你把底层逻辑和解决办法讲清楚：

核心区别：假定大小数组 vs 假定形状数组 #

你修改前后的两种数组声明，本质是完全不同的特性：

• double precision A(*)是假定大小数组（assumed-size array），属于Fortran77的旧特性。它只告诉编译器「这是一个一维的双精度数组」，但不会传递数组的实际大小/边界信息给子程序。调用时主程序只传数组的起始地址，子程序里可以通过索引访问元素（但编译器不会帮你检查越界，全靠自己保证）。
• double precision A(:)是假定形状数组（assumed-shape array），Fortran90引入的现代特性。它要求编译器能获取数组的完整形状（包括各维度的上下界、总大小），这就必须通过**显式接口（explicit interface）**来传递这些额外的形状描述信息。

你的段错误根源 #

你修改了子程序的数组声明，但没有给主程序提供子程序的显式接口——这就导致调用时，主程序只传了数组的起始地址，而子程序却期望获取形状描述信息，结果拿到的是垃圾数据。当子程序尝试用这些错误的形状信息访问数组时，就会触发内存越界，直接导致段错误。

解决办法（按推荐程度排序） #

1. 用模块包裹子程序（最规范的现代写法） #

把你的子程序放到一个模块里，主程序通过use语句引用模块。这样编译器会自动生成并检查显式接口，完全不需要手动处理：

module array_ops
  implicit none
  contains
    subroutine my_subroutine(A)
      ! 现代假定形状数组声明
      double precision, intent(inout) :: A(:)
      ! 子程序逻辑示例：打印数组大小
      print *, "Array size in subroutine: ", size(A)
    end subroutine my_subroutine
end module array_ops

program main
  use array_ops
  implicit none
  integer, parameter :: arr_size = 10
  double precision :: A(arr_size)
  
  ! 初始化数组
  A = 1.0d0
  ! 调用子程序，接口自动生效
  call my_subroutine(A)
end program main

2. 手动编写显式接口块 #

如果不想用模块，就在主程序（或主程序能访问的作用域）里手动写出子程序的接口块，告诉编译器参数的类型和形状：

program main
  implicit none
  integer, parameter :: arr_size = 10
  double precision :: A(arr_size)
  
  ! 手动声明子程序的显式接口
  interface
    subroutine my_subroutine(A)
      double precision, intent(inout) :: A(:)
    end subroutine my_subroutine
  end interface
  
  A = 1.0d0
  call my_subroutine(A)
end program main

! 子程序可以单独放在外部
subroutine my_subroutine(A)
  double precision, intent(inout) :: A(:)
  print *, "Array size in subroutine: ", size(A)
end subroutine my_subroutine

额外提醒 #

使用假定形状数组的好处远不止语法更现代：编译器会帮你检查调用时的数组维度是否匹配，还能在子程序里直接用size(A)、ubound(A)等函数获取数组信息，不需要额外传递大小参数，代码更简洁安全。

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