如何在Fortran中检测空行并按空行分割数据块计算均值？

Fortran默认的格式化输入会自动跳过空白字符（包括空行），但可以通过读取整行字符串的方式捕获空行，进而实现按空行分隔处理数据块。以下是具体实现思路和代码示例：

核心思路 #

使用READ(unit, '(A)') line语句读取每一行，空行会被读取为全空格的字符串，通过trim(line) == ''即可判断空行。同时维护状态变量跟踪是否处于数据块中，累加块内数值，遇到空行或文件结束时计算均值。

代码示例 #

program block_mean
    implicit none
    character(len=256) :: line  ! 足够长的字符串存储每行内容
    integer :: iostat, num, sum_val, count_val
    logical :: in_block

    ! 初始化状态
    in_block = .false.
    sum_val = 0
    count_val = 0

    ! 打开目标文件（替换为你的文件名）
    open(unit=10, file='data.txt', status='old', action='read', iostat=iostat)
    if (iostat /= 0) then
        print *, '文件打开失败'
        stop
    end if

    ! 逐行遍历文件
    do
        ! 读取整行，捕获空行
        read(10, '(A)', iostat=iostat) line
        ! 文件结束时退出循环
        if (iostat /= 0) exit

        ! 判断当前行是否为空行
        if (trim(line) == '') then
            ! 若处于数据块中，计算并输出均值
            if (in_block) then
                print *, '块均值：', real(sum_val)/count_val
                sum_val = 0
                count_val = 0
                in_block = .false.
            end if
        else
            ! 非空行，转换为整数并累加
            read(line, *) num
            sum_val = sum_val + num
            count_val = count_val + 1
            in_block = .true.
        end if
    end do

    ! 处理文件末尾的最后一个数据块
    if (in_block) then
        print *, '块均值：', real(sum_val)/count_val
    end if

    close(10)
end program block_mean

关键细节说明 #

• 读取整行：'(A)'格式会读取当前行的所有字符（包括空格），空行被读取为全空格字符串，trim()后变为空串，以此识别空行。
• 状态跟踪：in_block变量标记是否正在处理一个数据块，避免空行连续出现时重复计算。
• 文件结束处理：循环结束后需检查是否还有未处理的最后一个块（文件末尾无空行的情况）。
• 错误处理：示例中包含了文件打开的错误判断，若需要更严谨的数值转换错误处理，可以在read(line, *) num时添加iostat参数捕获转换失败的情况。

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