为什么在 Fortran 中使用命令 PRINT 会覆盖输入文件？

我正在编写代码并使用 Fortran 中的输入和输出功能。代码看起来像这样（仅用于简化):

PROGRAM TEST

  REAL, DIMENSION(1000):: A
  REAL:: B
  INTEGER::T

 !Defining input and output  
  OPEN(UNIT=1, FILE='input.dat', STATUS='OLD')
  OPEN(UNIT=2, FILE='output.dat', STATUS='NEW')  

 !Reading from file "input.dat"  
  READ(1,*) (A(I),I=1,1000)

 !Just for initial condition
  B=0.0  

  DO T=1, 10
    PRINT *, 'Step =', T 
        DO I=1, 1000     
           B=B+1.0     
           A(I)=A(I)/B  
        END DO
  END DO

 !Writing results into file "output.dat"
   DO I=1, 1000
      WRITE (2,100) I, A(I)
   END DO 
   100 FORMAT (' ',T3, I12, T17, F14.4)   

END PROGRAM TEST

我使用的是 Gfortran 5.3，结果并不像我预期的那样。我期望获得变量的结果T程序运行时屏幕上（或 Ubuntu 操作系统中的终端）上显示的变量I and A(I)被写入文件output.dat。我的变量没有问题I and A(I)，因为它们已成功写入文件output.dat。问题出在变量上T，它没有出现在终端上，但被写入文件中input.dat。好吧，即使是文件中的前一个文件input.dat没有被覆盖。有人能给我建议吗？

仅供参考，我还尝试过其他编译器（使用 Windows 操作系统），例如：

1. Microsoft Fortran Powerstation（非常旧的）：但它的工作方式正如我预期的那样。
2. MinGW-w64（Windows 的 GCC 版本）：但无法正常工作。

这可能是因为对于平台/编译器/编译器版本/编译器选项的特定组合，单元 1 是控制台的预连接单元。

您的 OPEN 语句将该单元定向到您的输入文件。因此，隐式寻址该单元的 PRINT 语句然后将其输出定向到同一文件。

使用不同的单元编号 - 选择大于 10 的值通常对于编译器预连接的单元是安全的。为了进一步安全，您可以使用INQUIRE(UNIT=unit_number, EXIST=some_logical_variable)语句来检查特定单元是否在 OPEN 语句之前连接到文件 - 如果是，则选择不同的单元号。理想情况下，如果您正在编写 Fortran 2008，则可以使用 NEWUNIT 说明符。

（不要将单位编号的值硬编码到输入/输出语句中 - 它们应始终由变量或命名常量表示，以便可以轻松地在一处设置/更改该值。）