在fortran中读取可变长度的实数值

Question

我有一个带有以下文本的文件：

-1.065211      246.0638 xlo xhi
-0.615       245.385 ylo yhi
-10 10 zlo zhi

我想将线上的数值读入3 * 2的实数矩阵中。
是否可以只读取一条记录的两个输入并转到下一行？ 就像是：
read(1,'(2F?.?/)') (matrix(i,1:2),i=1,3)
我在F中加了问号？ 因为我的数字长度是可变的。 换句话说，我只需要阅读自由格式的记录中的两项。 我了解您可以循环执行此操作。 我想知道一个优雅的一线客是否可行。

Answer 1

您可以用实现的目标

do i = 1, 3
   read(fd,*) matrix(i,:)
end do

假定fd已连接到包含您的数据的文件。

Answer 2

@Steve答案中的代码是可以完成您想要的最简单的代码。 如果您希望将其放在一行中 ，则只需

do i=1,3 ; read(fd,*) matrix(i,:) ; end do

如果您要寻找的是仅用一条语句 “仅读取一条记录的两个输入并转到下一行”的方法，就我所知，这是不可能的：由于您的实际值具有不同的格式，因此您需要使用列表导向的输入（格式中的星号）以读取具有相同语句的多行（所有其他格式的实数都需要输入正宽度，并且由于您的实数具有不同的宽度，因此会产生任何宽度规范）在运行时错误中）。

如果您只想使用一条语句，则需要明确告知程序读取每行末尾的两个字符串，即使它们已被丢弃（否则程序将尝试将字符串转换为实数，并且这将导致运行时错误）：

character(len=10) :: r
read(1,*) (a(i,:),r,r,i=1,3)

它不如@Steve的代码通用（要求每行恰好有2个字符串），并且实际上需要字符变量的附加声明。

最后，请注意，您不是按元素顺序读取数组：在Fortran中，数组在内存中的元素顺序为a(1,1) ， a(2,1) ， a(3,1) ， a(1,2)等。如果您计划操作的实际数组大得多，则索引的顺序可能会极大地影响性能，因为连续的内存访问可能会寻址内存的远处而不是连续的部分它。 在用Fortran编写程序时，就像使用其他任何语言（可能使用不同的顺序）一样，您需要考虑数组索引的顺序。