逐行读取文本文件，但仅读取特定列

Question

我们如何在跳过其中的某些列的同时逐行读取特定文件？

例如，我有一个包含数据的文本文件，分为 5 列，但我只需要从中读取两列，它们可以是前两列或任何其他随机组合（我的意思是，需要一个可行的解决方案与列的任何组合，如第一和第三列）。

代码是这样的

        open(1, file=data_file)
        read (1,*) ! to skip first line, with metadata
        lmax = 0
        do while (.true.)
                ! read column 1 and 3 here, either write
                ! that to an array or just loop through each row
        end do
99      continue        
        close (1)

任何解释或示例都会有很大帮助。

Answer 1

High Performance Mark 的回答给出了简单的选择性列读取的要点：仍然读取列但将其转移到然后被忽略的变量。

为了扩展该答案，请考虑我们要读取五列行的第二列和第四列：

read(*,*) junk, x, junk, y

第一个值被转移到junk ，然后第二个被转移到x ，然后第三个（替换刚刚获得的那个）到junk ，最后第四个到y 。 第五个被忽略，因为我们已经用完了输入项并且传输语句终止（循环中的下一次读取将转到下一条记录）。

当然，当我们知道这是我们想要的那些列时，这很好。 当我们事先不知道时，让我们概括一下：

integer col1, col2   ! The columns we require, defined somehow (assume col1<col2)
<type>, dimension(nrows) :: x, y, junk(3)  ! For the number of rows
integer i

do i=1,nrows
  read(*,*) junk(:col1-1), x(i), junk(:col2-col1-1), y(i)
end do

在这里，我们将多个值（可能为零）转移到感兴趣的第一列之前，然后是感兴趣的值。 之后，更多要忽略的值（可能为零），然后是感兴趣的最终值。 跳过该行的其余部分。

这仍然是非常基本的，并且避免了许多潜在的需求复杂性。 在某种程度上，这是一种基本的方法，人们不妨考虑一下：

do i=1,nrows
  read(*,*) allofthem(:5)
  x(i) = allofthem(col1)
  y(i) = allofthem(col2)
end do

（该变量是逐行临时的）但多样性和选项都很好。

Answer 2

这很容易。 您只需从每一行读取 5 个变量，并忽略您不再使用的变量。 就像是

do i = 1, 100
    read(*,*) a(i), b, c(i), d, e
end do

这将在每次迭代时覆盖b 、 d和e中的值。

顺便说一句，你的行

99 continue

是多余的； 它不用作do循环的结束行，并且您不会从其他任何地方分支到它。 如果您从看不见的代码分支到它，您可以将标签99附加到下一行并删除continue语句。 通常，在现代 Fortran 中， continue是多余的； 特别是它在您的代码中似乎是多余的。