Question

我遇到Fortran 90代码问题。代码使用read语句读取大量数据。错误消息指向此段代码的最后一行：

        open(1,file=trim(filename),action='read',form='unformatted',status='old')
        read(1) np    ! np is a 6 element array of integers*4
        npart = np(species)

        if (sum(np(1:species-1))>0) then
            read(1) (xempty,yempty,zempty,i=1,sum(np(1:species-1)))
        end if
        allocate(x(npart),y(npart),z(npart),mass(npart))

        write(*,*) "number of particles: ", npart
        read(1) (x(i),y(i),z(i),i=1,npart

运行代码时出现错误消息。它是：

Fortran runtime error: I/O past end of record on unformatted file

我认为这个问题与Fortran 90中的记录长度限制有关，但我不知道如何解决这个问题，以及如何在gfortran中查找或修改记录限制。如果我使用较小的文件名，则不会显示错误消息。换句话说，当npart较小时。

我正在使用以下gfortran编译器：

Using built-in specs.
COLLECT_GCC=gfortran
COLLECT_LTO_WRAPPER=/usr/lib/gcc/x86_64-linux-gnu/5/lto-wrapper
Target: x86_64-linux-gnu
Configured with: ../src/configure -v --with-pkgversion='Ubuntu 5.4.0-
6ubuntu1~16.04.4' --with-bugurl=file:///usr/share/doc/gcc-
5/README.Bugs --enable-languages=c,ada,c++,java,go,d,fortran,objc,obj-
c++ --prefix=/usr --program-suffix=-5 --enable-shared --enable-linker-
build-id --libexecdir=/usr/lib --without-included-gettext --enable-
threads=posix --libdir=/usr/lib --enable-nls --with-sysroot=/ --
enable-clocale=gnu --enable-libstdcxx-debug --enable-libstdcxx-
time=yes --with-default-libstdcxx-abi=new --enable-gnu-unique-object -
-disable-vtable-verify --enable-libmpx --enable-plugin --with-system-
zlib --disable-browser-plugin --enable-java-awt=gtk --enable-gtk-cairo 
--with-java-home=/usr/lib/jvm/java-1.5.0-gcj-5-amd64/jre --enable-
java-home --with-jvm-root-dir=/usr/lib/jvm/java-1.5.0-gcj-5-amd64 --
with-jvm-jar-dir=/usr/lib/jvm-exports/java-1.5.0-gcj-5-amd64 --with-
arch-directory=amd64 --with-ecj-jar=/usr/share/java/eclipse-ecj.jar --
enable-objc-gc --enable-multiarch --disable-werror --with-arch-32=i686 
--with-abi=m64 --with-multilib-list=m32,m64,mx32 --enable-multilib --
with-tune=generic --enable-checking=release --build=x86_64-linux-gnu -
-host=x86_64-linux-gnu --target=x86_64-linux-gnu
Thread model: posix
gcc version 5.4.0 20160609 (Ubuntu 5.4.0-6ubuntu1~16.04.4)

输入文件是二进制格式，它是使用C代码创建的。我不可能知道用于创建文件的编译器。代码的最后一行：read（1）（x（i），y（i），z（i），i = 1，npart）读取实际（4字节）格式的粒子的位置。此外，当我使用粒子数较少的输入文件时，换句话说，当npart较小时，不会出现错误。

Fortran运行时错误：I / O在未格式化文件上的记录结束

0 个答案: