"使用 Fortran 90 在输入文件中正确读取注释行"

Question

据我了解，Fortran 在从文件中读取数据时，会跳过以星号 (*) 开头的行，假设它们是注释。 好吧，我似乎在用我创建的一个非常简单的程序实现这种行为时遇到了问题。 这是我的简单 Fortran 程序：

  1       program test
  2 
  3       integer dat1
  4 
  5       open(unit=1,file="file.inp")
  6 
  7       read(1,*) dat1
  8 
  9 
 10       end program test

这是“file.inp”：

  1 *Hello
  2 1

我用

gfortran -g -o test test.f90

当我运行时，我收到错误：

At line 7 of file test.f90 (unit = 1, file = 'file.inp')
Fortran runtime error: Bad integer for item 1 in list input

当我运行删除了注释行的输入文件时，即：

1 1

代码运行良好。 因此，Fortran 正确解释该注释行似乎是一个问题。 这一定是我在这里想念的非常简单的东西，但我无法在谷歌上找到任何东西。

Answer 1

Fortran 不会自动跳过输入文件中的注释行。 您可以很容易地做到这一点，首先将行读入字符串，检查您的注释符号的第一个字符或在字符串中搜索该符号，然后如果该行不是注释，则对字符串进行“内部读取”以获取数值。

就像是：

use, intrinsic :: iso_fortran_env

character (len=200) :: line
integer :: dat1, RetCode

read_loop: do
   read (1, '(A)', isostat=RetCode)  line
    if ( RetCode == iostat_end)  exit ReadLoop
    if ( RetCode /= 0 ) then
      ... read error
      exit read_loop
    end if
    if ( index (line, "*") /= 0 )  cycle read_loop
    read (line, *) dat1
end do read_loop

Answer 2

Fortran 默认情况下不会忽略任何内容，除非您使用名称列表并且在这种情况下注释以感叹号开头。

"

Answer 3

我发现使用退格<\/a>语句比建议的解决方案更直观。 当在行首遇到注释字符“#”时，以下子例程将跳过该行。

subroutine skip_comments(fileUnit)
  integer, intent(in) :: fileUnit
  character(len=1) :: firstChar

  firstChar = '#'
  do while (firstChar .eq. '#')
    read(fileUnit, '(A)') firstChar
  enddo
  backspace(fileUnit)

end subroutine skip_comments