为什么直接访问 I/O 在英特尔可视路由器中工作不正确

Question

这个程序

  program test
  real a(10)
  open(1, file='f1',access='direct', recl=20)
  do i=1, 10
      a(i) = i-1
  end do
  write(1, rec=1)(a(i),i=1,5)
  write(1, rec=2)(a(i),i=6,10)
  close(1)
  open(1, file='f1',access='direct',recl=8)
  read(1, rec=4)(a(i),i =5,9,4)
  print*,a
  end

在视觉Fortran中工作不正确（不正确）：

  0.0000000E+00   1.000000       2.000000       3.000000       9.000000
   5.000000       6.000000       7.000000      0.0000000E+00   9.000000

WATCOM的结果（正确）：

  0.0000000   1.0000000       2.0000000       3.0000000       6.0000000
  5.0000000   6.0000000       7.0000000       7.0000000       9.0000000

为什么？

Answer 1

在编写代码时，您不需要知道记录以任何单位计算的长度到底有多长。设置recl=20或recl=5并根据编译器的确切行为和确切的real大小是未来问题的秘诀。

您应该向编译器询问记录的大小，无论它使用什么单位：

integer :: recl5, recl2
inquire(iolength=recl5) (a(i),i=1,5)
open(1, file='f1',access='direct', recl=recl5)

...

inquire(iolength=recl2) (a(i),i =5,9,4)
open(1, file='f1',access='direct',recl=recl2)

注意：在 Fortran 90 中，使用子数组表示法可以缩短它。

注意2：对于某些编译器，您的程序将永远无法工作，因为编译器可以在文件中使用记录结束标记，并且您必须始终使用相同的recl打开它。但这并不常见。

Answer 2

Visual Fortran 中 RECL 的默认单位长度是一个单词（4 字节）。如果您使用选项"使用字节作为未格式化文件的 RECL= 单位"（/assume：byterecl）进行编译，您将获得所需的内容。