arrays - 正确使用 Fortran 流写入和读取整数数组

Question

这与我最近发表的有关 Fortran 流等的帖子有些相关:Converting data stored in Fortran 90 binaries to human readable format .

我试图将一个简单的整数数组写入一个文件，然后使用 Fortran 的 READ 函数读入我创建的二进制文件。我通过在我的 OPEN 指令中包含 ACCESS="STREAM" 来使用流。我有以下代码:

MODULE streamtest2subs

  IMPLICIT NONE

CONTAINS

SUBROUTINE writeUstream(myarray)
  IMPLICIT NONE
  INTEGER, INTENT(IN), DIMENSION(4,10) :: myarray
  INTEGER :: myvalue = 12345, mypos
  OPEN(UNIT=11, FILE="ustream.demo", STATUS="REPLACE", ACCESS="STREAM")
  WRITE(11) myarray
  CLOSE(UNIT=11)
END SUBROUTINE writeUstream

SUBROUTINE readUstream
  IMPLICIT NONE
  INTEGER :: test1, test2, test3
  INTEGER :: n
  OPEN(UNIT=42, FILE="ustream.demo", STATUS="OLD", ACCESS="STREAM")
  READ(42, POS=1) test1
  READ(42, POS=2) test2
  READ(42, POS=3) test3

  WRITE(*,*) "This is the output:"
  WRITE(*,*) test1
  WRITE(*,*) test2
  WRITE(*,*) test3
END SUBROUTINE readUstream

END MODULE streamtest2subs

PROGRAM streamtest2
  USE streamtest2subs
  IMPLICIT NONE
  INTEGER :: i, j, k
  INTEGER, DIMENSION(4,10) :: a

  WRITE(*,*) "This is my input array:"
  k=1
  DO i=1,4
    DO j=1,10
      a(i,j)=k
      WRITE(*, "(i3)", ADVANCE="NO") a(i,j)
      k=k+1
    END DO
    WRITE(*,*)
  END DO
  WRITE(*,*)

  CALL writeUstream(a)
  CALL readUstream
END PROGRAM streamtest2

但是，当我使用 gfortran 编译并运行它时，我得到以下输出:

 This is my input array:
  1  2  3  4  5  6  7  8  9 10
 11 12 13 14 15 16 17 18 19 20
 21 22 23 24 25 26 27 28 29 30
 31 32 33 34 35 36 37 38 39 40

 This is the output:
           1
   184549376
      720896

为什么输出这么复杂？ READ 是否正在将 ustream.demo 文件作为字符串而不是整数来读取？但是，当我将 test1、test2 和 test3 的类型更改为字符串时，我的输出只是一系列的三个空行。

我是否在 READ 中错误地使用了 POS 指令？我认为 POS 指定了输出中的字符编号(尽管我不确定数组的元素是否以任何方式分隔)；这是正确的吗？

非常感谢您的宝贵时间!

Answer 1

出于您在上一个问题中描述的目的，我认为您的编程比必要的更复杂。假设您想要一个没有额外记录结构的二进制文件，您想要未格式化的和流式的。您可以像编写文件一样阅读文件——您不需要使用 POS——除非目的是学习如何使用 POS。

我的测试示例基于一个 Fortran 程序，我在其中读取了一个由其他人编写的二进制文件，几乎可以肯定是使用 C 程序编写的。该文件由标题组成，后跟可变长度的数组。我打开文件:

open ( unit=75, file=FileName, status='old', access='stream', form='unformatted', action='read' )

我读了一个标题(一个具有许多子变量的用户定义类型的变量):

read (75) header

我将数组分配给读取到 header 的一个字段中的长度，然后读取数组:

allocate ( array (1:header % ArrayLen) )
read (75) array

然后我处理数组中的数据。然后我重复直到文件结尾(代码示例中未显示)。

非常简单...无需计算文件中的位置并使用 READ 的 POS 关键字。