c语言调用fortran数组,正确使用Fortran流来编写和读取整数数组

这与我最近发布的关于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指定了输出中的字符编号(尽管我不确定数组中的元素是否以任何方式分隔)。它是否正确？

非常感谢您的宝贵时间！