Question

我编写了以下块来从外部数据文件中读取：

     open(unit=338,file='bounnodes.dat',form='formatted') 
      DO I=1,NQBOUN
         DO J=1,NUMBOUNNODES(I)
            read(338,2001) NODEBOUN(i,j)
            write(6,*) 'BOUNDARY NODES',  NODEBOUN(i,j)
         ENDDO
       ENDDO
     2001
     FORMAT(32I5)

据我所知，这应该从bounnodes.dat读取一个2 x 32阵列。但是，我在读取过程中收到错误end-of-file并打印出第一列。

我尝试使用相同的代码读取32 x 2数组，它读取第一列的32个元素，但为下一列输出0s。

你能解释一下发生了什么吗？我的格式错了吗？

Answer 1

Fortran中的每个读取语句都会进入下一条记录。这意味着普通文本文件中的新行。试试这个：

   DO I=1,NQBOUN
     DO J=1,NUMBOUNNODES(I)
        read(338,2001,advance='no') NODEBOUN(i,j)
        write(*,*) 'BOUNDARY NODES',  NODEBOUN(i,j)
     ENDDO
     read(338,*)
   ENDDO

其中NQBOUN是行数，NUMBOUNNODES（I）是行中的列数。（我总是有问题，32x2对2x32）

您可以使用隐含的操作

使其更短

   DO I=1,NQBOUN
        read(338,2001) ( NODEBOUN(i,j) , j=1,NUMBOUNNODES(I) )
        write(*,*) ( 'BOUNDARY NODES', NODEBOUN(i,j) , j=1,NUMBOUNNODES(I) )
   ENDDO

甚至

   DO I=1,NQBOUN
        read(338,2001) NODEBOUN(i,:)
        write(*,*) 'BOUNDARY NODES',  NODEBOUN(i,1:NUMBOUNNODES(I))
   ENDDO

所有这些都使用Fortran 90功能。

从fortran中的数据文件中读取列

1 个答案: