能否在Fortran FUNCTION中使用WRITE语句？遇程序卡顿问题求解释

Fortran标准完全允许在FUNCTION中使用WRITE语句，你的代码语法本身没有问题，程序出现的“卡顿”实际上是输出缓冲机制导致的视觉假象，而非代码错误。

问题根源 #

gfortran默认对终端输出（单元6对应标准输出stdout）采用行缓冲模式：只有当输出内容包含换行符，或者程序显式刷新缓冲区时，缓冲区里的内容才会被输出到终端。你的FORMAT语句中没有添加换行符，输出内容一直滞留在缓冲区中没有显示，让你误以为程序卡住了——实际上程序已经执行完函数逻辑，只是输出没被打印出来。

解决方案 #

有三种常见的解决方式：

1. 在FORMAT中添加换行符
   在格式串末尾加上/（Fortran的换行控制符），让输出自动触发缓冲刷新：

   INTEGER FUNCTION CURSED()
         I = 42
         WRITE(6, 2023) I
     2023 FORMAT('>>> In CURSED = ', I4, /)
         CURSED = I
         RETURN
         END
   

2. 手动刷新缓冲区
   在WRITE语句后调用CALL FLUSH(6)，强制将缓冲区内容输出到终端（gfortran支持该扩展命令，跨平台可用）：

   INTEGER FUNCTION CURSED()
         I = 42
         WRITE(6, 2023) I
         CALL FLUSH(6)
     2023 FORMAT('>>> In CURSED = ', I4)
         CURSED = I
         RETURN
         END
   

3. 编译时禁用缓冲
   使用gfortran -fflush-all-i/o编译选项，强制所有I/O操作立即刷新缓冲区，无需修改代码。

补充说明 #

Fortran标准从未禁止在函数中执行I/O操作，包括WRITE、READ等。你在MacOS和Windows平台遇到的一致现象，本质都是终端缓冲的共性行为，和平台本身关联不大，只要解决缓冲刷新问题就能正常输出。

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