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

这与我最近发表的关于Fortran stream之类的文章有些关联：

我试图将一个简单的整数数组写入一个文件，然后使用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的类型更改为string时，我的输出只是一系列三个空行

我是否错误地使用了

READ

中的

POS

指令？我认为

POS

指定了输出中的字符编号（尽管我不确定数组的元素是否以任何方式分隔）；这是正确的吗

---

非常感谢您抽出时间

问题在于您的读取语句，

POS

是从文件开始的字节偏移量，与写入文件的原始

a

数组的数组索引没有直接关系

因此，在您的代码中，如果

test1

是一个4字节的整数，Fortran read将通过读取文件的1到4字节来构造它。类似地，

test2

将从字节2到5构造，而值实际上存储在字节5到8中

当我将读取行更改为

READ(42, POS=1) test1
READ(42, POS=5) test2
READ(42, POS=9) test3

输出变为

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
      11
      21

我想这更像是你想要的。请注意，如果最初将值作为2字节或8字节整数写入文件，则在读取文件时，必须考虑

POS

偏移量中的这些大小

另外，值得学习用于快速查询二进制文件的命令行工具

od

。例如，我做的第一件事就是用

$ od -t d4 ustream.demo 
0000000                 1              11              21              31
0000020                 2              12              22              32
0000040                 3              13              23              33
0000060                 4              14              24              34
0000100                 5              15              25              35
0000120                 6              16              26              36
0000140                 7              17              27              37
0000160                 8              18              28              38
0000200                 9              19              29              39
0000220                10              20              30              40
0000240

---

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

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

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

我读取了一个标题（一个包含许多子变量的用户定义类型的变量）：

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

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

然后我处理数组中的数据。 然后重复，直到文件结束（代码示例中未显示）

很简单。。。无需计算文件中的位置，使用读取的POS关键字

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