帮我看看这个文件读入错误

straggler

在附件中的Fortran程序中，第51-60行有一个文件读入语句：

2.       OPEN(2,FILE='f',FORM='UNFORMATTED')
3.       READ (2) (F(J),J=1,NEQ)
4.     CLOSE (2)
5.     open(22,FILE='f1.dat',FORM='FORMATTED')
6.     write (22,*) (F(J),J=1,NEQ)
7.       CLOSE (22)
8. c    write(*,*)'f=',f
9.       DO J=1,NEQ
10.       U(J)=F(J)
11.     enddo

复制代码

在DEBUG时，我发现数组U(J)的元素的值全为NAN，我在ENDDO后面设置了一个断点来跟踪是
否文件读入有问题，跟踪后发现，F(J)的值也为NAN，这时我想可能是OPEN语句要打开的对
象文件'f'有问题，为了检测该数据文件是否有问题，我在READ语句后面加入了一个写语句
，让它输出一个有格式文件f1.dat。结果发现，f1.dat文件中却写入成功，也就是说，REA
D语句已经成功从文件f中读入了数据，那为什么F(J)的值为NAN？请高手指点一二，不胜感
激。附件中是程序文件solv.for及数据文件f。

freepu

程序简化到如下程度，没有发现问题

2.       implicit real*8 (a-h,o-z)
3.       integer, parameter :: NEQ=10
4.       real(4) F(NEQ), U(NEQ)
7. C.......OPEN F FILE
8.       OPEN(2,FILE='f',FORM='UNFORMATTED')
9.       READ (2) (F(I),I=1,NEQ)
10.       CLOSE (2)
11.       open(22,FILE='f1.dat',FORM='FORMATTED')
12.       write (22,*) (F(I),I=1,NEQ)
13.       CLOSE (22)
14. c  write(*,*)'f=',f
15.       DO I=1,NEQ
16.       U(I)=F(I)
17.      enddo
18. c  
19.      end

复制代码

2.   -1.236702     -0.1950102      0.0000000E+00  0.0000000E+00  -1.236700   
3. -0.2106352      -1.236700     -0.2106352      0.0000000E+00  0.0000000E+00

复制代码

straggler

谢谢FREEPU版主，我再看看是不是别的什么原因。

freepu

straggler wrote:
谢谢FREEPU版主，我再看看是不是别的什么原因。

不用客气，偶可不是什么版主啊;)
可是从傲雪来？