Fortran 77读取外部文件时如何跳过#开头的注释行？

你的代码存在的核心问题 #

• 硬编码跳过行完全失效：你手动指定的跳过行号是固定值，但文件中注释行是每约4000行出现一次，位置可能不固定，这种写法既不灵活，也容易因文件内容变化直接出错。
• 读取逻辑顺序错误：你先跳了几行后直接循环读固定N行，但中间遇到的注释行仍会被read(7,*) a,b,c,d,e读取，程序尝试把#开头的字符串解析成实数，直接触发"Bad real number"错误。
• 错误的写操作：write(7, *) event是往输入数据文件里写内容，会直接破坏你的原始数据，属于完全错误的操作。
• 固定行数计数不合理：你用固定的N=40010作为总行数，但因为存在注释行，实际有效数据行数少于这个值，最终计算的平均值会不准确。

通俗易懂的解决方案 #

核心思路是逐行读取，先判断是否为注释行，再处理数据，不需要提前知道注释行位置，也不用固定行数：

program mean_analysis
    implicit none
    
    integer :: count_lines  ! 计数有效数据行数
    double precision :: a, b, c, d, e
    double precision :: suma, sumb, sumc, sumd, sume
    character(100) :: line  ! 存储读取的每一行内容
    logical :: is_comment   ! 判断是否为注释行
    integer :: io_status    ! 文件读取状态码

    ! 初始化总和和计数器
    suma = 0.0d0
    sumb = 0.0d0
    sumc = 0.0d0
    sumd = 0.0d0
    sume = 0.0d0
    count_lines = 0

    ! 打开输入文件，指定status="old"确保文件存在
    open(unit=7, file="zpc_initial_momenta.dat", status="old", iostat=io_status)
    if (io_status /= 0) then
        print *, "Error opening input file!"
        stop  ! 打开失败直接终止程序
    end if

    ! 循环读取直到文件结束
    do
        ! 读取一行到字符变量，同时获取读取状态
        read(7, '(A)', iostat=io_status) line
        ! 检查是否读到文件末尾，是就退出循环
        if (io_status /= 0) exit

        ! 判断是否为注释行：第一个字符是#就跳过
        is_comment = (line(1:1) == '#')
        if (is_comment) cycle  ! 跳过当前循环，处理下一行

        ! 不是注释行，解析成5个数值
        read(line, *) a, b, c, d, e

        ! 累加总和，计数有效行数
        suma = suma + a
        sumb = sumb + b
        sumc = sumc + c
        sumd = sumd + d
        sume = sume + e
        count_lines = count_lines + 1
    end do

    close(7)

    ! 检查是否有有效数据，避免除以0的错误
    if (count_lines == 0) then
        print *, "No valid data found!"
        stop
    end if

    ! 写入结果到输出文件，用格式化输出让结果更整齐
    open(unit=8, file="outputanalysis.dat", status="replace")
    write(8, '(5F12.6)') suma/count_lines, sumb/count_lines, &
                          sumc/count_lines, sumd/count_lines, sume/count_lines
    close(8)

    ! 打印提示信息，方便确认运行结果
    print *, "Analysis done! Processed ", count_lines, " data lines."
    print *, "Results saved to outputanalysis.dat"

    end program

关键细节说明 #

• 用iostat检查文件状态：每次读取都判断是否读到文件末尾，或文件打开失败，避免程序莫名崩溃。
• 逐行判断注释行：读取每一行到字符变量后，检查第一个字符是否为#，是的话直接跳过，不解析数值。
• 动态计数有效行数：用count_lines统计实际处理的数据行数，最后用这个数计算平均值，避免固定N带来的误差。
• 移除错误写操作：删掉了往输入文件写内容的代码，保护你的原始数据。
• 格式化输出：用(5F12.6)格式输出平均值，让结果排版更整齐易读。

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