Fortran条件与积分代码错误排查及修正求助

先贴出对应编译报错场景的典型错误代码：

PROGRAM CALC_R
  IMPLICIT NONE
  INTEGER, PARAMETER :: T=2, M=3, Y=5
  REAL :: A(Y), B(Y), X(M), R(Y)
  INTEGER :: i, j

  ! 读取输入文件
  OPEN(UNIT=10, FILE='A.dat', STATUS='OLD')
  READ(10, *) A
  CLOSE(10)

  OPEN(UNIT=11, FILE='B.dat', STATUS='OLD')
  READ(11, *) B
  CLOSE(11)

  OPEN(UNIT=12, FILE='T.dat', STATUS='OLD')
  READ(12, *) T  ! 错误点1
  CLOSE(12)

  ! 初始化X数组
  X = [-2.0, -1.0, 0.0, 1.0, 2.0]  ! 错误点2

  ! 计算R值
  DO i = 1, Y
    DO j = 1, M  ! 错误点3
      R(i) = A(i) + B(i)*X(j) + T*M
    END DO
  END DO

  ! 输出R值
  DO i = 1, 6  ! 错误点4
    PRINT *, "R(", i, ") = ", R(i)
  END DO

END PROGRAM CALC_R

错误1：参数重定义（第9行） #

• 问题：T被声明为PARAMETER常量，意味着值不可修改，但后面又用READ语句给它赋值，编译直接触发“变量重定义/无法修改常量”错误。
• 修正：
  • 如果需要用T.dat的输入覆盖初始值：删掉PARAMETER，把T改成普通整型变量；
  • 如果T是固定参数：直接删掉读取T.dat的整个代码块。

修正后代码段二选一：

! 方案1：允许修改T的值
INTEGER :: T=2, M=3, Y=5
OPEN(UNIT=12, FILE='T.dat', STATUS='OLD')
READ(12, *) T
CLOSE(12)

! 方案2：固定T为参数，删除读取逻辑
INTEGER, PARAMETER :: T=2, M=3, Y=5

错误2：数组维度不兼容（第15行） #

• 问题：X被定义为长度M=3的数组，但赋值了5个元素（[-2,-1,0,1,2]），数组长度和元素数量不匹配，编译时报“维度不一致”错误。
• 修正：把M的值改成5，让X的长度和赋值元素数量对应：

INTEGER, PARAMETER :: T=2, M=5, Y=5  ! M改为5
REAL :: A(Y), B(Y), X(M), R(Y)
X = [-2.0, -1.0, 0.0, 1.0, 2.0]  ! 维度完全匹配

错误3：循环逻辑错误与越界风险（第19-21行） #

• 问题：内层循环j跑M次，但每次都把结果覆盖到同一个R(i)里，最后R(i)只保留最后一次循环的结果；如果需要存储所有(i,j)组合的R值，当前一维的R数组根本存不下，后续操作容易触发越界。
• 修正：把R改成二维数组，存储每个组合的计算结果：

REAL :: A(Y), B(Y), X(M), R(Y,M)  ! R定义为Y行M列的二维数组
! 计算R值
DO i = 1, Y
  DO j = 1, M
    R(i,j) = A(i) + B(i)*X(j) + T*M  ! 每个组合单独存值
  END DO
END DO

错误4：数组越界（第25行） #

• 问题：R的长度是Y=5，但循环i从1跑到6，超出数组最大下标，编译或运行时触发“数组越界”错误。
• 修正：把循环上限改成Y，二维数组则遍历所有维度：

! 二维数组输出所有R值
DO i = 1, Y
  DO j = 1, M
    PRINT *, "R(", i, ",", j, ") = ", R(i,j)
  END DO
END DO

! 一维数组输出版本
DO i = 1, Y
  PRINT *, "R(", i, ") = ", R(i)
END DO

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