Fortran95循环批量读取.dat文件编译报错，求代码修正方案

核心错误原因 #

编译错误是因为你在do i=1,10循环体内，使用read(filename, '(I0)') i试图修改循环变量i——Fortran语法明确禁止在循环体内重新定义循环控制变量，这会触发编译器报错。另外，这个read语句本身逻辑错误：你需要把整数i转换成字符串来生成文件名，应该用write而非read。

修正后的完整代码 #

program read_dat_files
    implicit none
    integer :: i, j, ios, unit_num
    character(len=10) :: filename
    ! 根据实际数据类型调整，这里假设是实数数组
    real, dimension(10, 10) :: aa  

    do i = 1, 10
        ! 把整数i格式化为字符串，直接生成目标文件名
        write(filename, '(I0,".dat")') i  

        ! 使用独立的文件单元号，避免和循环变量冲突
        unit_num = 10 + i  
        open(unit=unit_num, file=filename, status='old', iostat=ios)
      
        if (ios /= 0) then
            print *, 'Error opening file ', trim(filename)
        else
            do j = 1, 10
                read(unit_num, *) aa(i, j)
            end do
            close(unit_num)
        end if
    end do

    ! 可选：验证读取结果
    print *, 'Read data successfully:'
    do i = 1, 10
        print *, aa(i, :)
    end do
end program read_dat_files

关键修改点 #

• 替换文件名生成逻辑：用write(filename, '(I0,".dat")') i直接生成1.dat、2.dat这类文件名，替代错误的read语句，同时避免修改循环变量i。
• 显式声明数组aa：必须提前声明数组的类型和维度，否则编译器无法识别该变量。
• 使用独立文件单元号：不用循环变量i作为文件单元号，改用unit_num = 10 + i（选择10以上的安全单元号，避免和标准输入输出单元冲突），提升代码安全性。
• 添加implicit none：强制显式声明所有变量，避免隐式类型带来的潜在错误（Fortran95推荐写法）。

其他注意事项 #

• 如果你的.dat文件中数据类型不是实数，需要把aa的类型改成对应类型（比如integer）。
• 如果文件数量不是固定10个，可以通过iostat判断文件是否存在，循环直到找不到文件为止，增强代码灵活性。

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