数据读取问题

Kieran

大家好，

我想请教大家一个问题。我在读取一个文件（INPUT.dat）的数据，其中文件的内容如下。
Hamiltonian file with magnetization along z axis
ac_111_0_90_wannier90_hr.dat
Output file with magnetization along z axis
ac_111_0_90_wannier90.wout
Hamiltonian file with magnetization along different direction
ac_111_0_85_wannier90_hr.dat
ac_111_0_95_wannier90_hr.dat
ac_111_90_85_wannier90_hr.dat
ac_111_90_95_wannier90_hr.dat
Angle difference in unit of degree
5.0
k point mesh
50 50
Velocity direction of Hamitlonian (0 - along x directin; 1 - along y direction)
0
Femi energy value in unit of eV
-0.5415
Energy window for Fermi surface in unit of eV
0.01
Broadening value in unit of eV
0.05

这个文件中的第11行是一个浮点数5.0。我读取该文件的部分程序如下。

[Fortran] 纯文本查看 复制代码

CHARACTER (LEN=50)    :: na(6)
DOUBLE PRECISION       :: an
INTEGER                       :: km(2)
INTEGER                       :: vd
DOUBLE PRECISION       :: fermi
DOUBLE PRECISION       :: wf
DOUBLE PRECISION       :: bv

OPEN (UNIT=3, FILE='INPUT.dat', STATUS='OLD')

READ (UNIT=3, FMT=*)
READ (UNIT=3, FMT='(a)') na(1)
READ (UNIT=3, FMT=*)
READ (UNIT=3, FMT='(a)') na(2)
READ (UNIT=3, FMT=*)
READ (UNIT=3, FMT='(a)') na(3)
READ (UNIT=3, FMT='(a)') na(4)
READ (UNIT=3, FMT='(a)') na(5)
READ (UNIT=3, FMT='(a)') na(6)
READ (UNIT=3, FMT=*)
READ (UNIT=3, FMT='(a)') an
READ (UNIT=3, FMT=*)
READ (UNIT=3, FMT=*) km
READ (UNIT=3, FMT=*)
READ (UNIT=3, FMT=*) vd
READ (UNIT=3, FMT=*)
READ (UNIT=3, FMT=*) fermi
READ (UNIT=3, FMT=*)
READ (UNIT=3, FMT='(a)') wf
READ (UNIT=3, FMT=*)
READ (UNIT=3, FMT=*) bv
write (unit=*,fmt=*) na(6), an

但是，当我查看屏幕上输出的变量时，看到的却是
ac_111_90_95_wannier90_hr.dat                       6.013568264329796E-154
也就是说变量an的值是6.013568264329796E-154，可是我的文件里对应于an变量的数据是5.0。


请问为什么会是这样呢？能麻烦大家给些建议，怎么解决这个问题吗？谢谢啦。

Kieran

fcode 发表于 2020-4-30 08:06
我亲爱的朋友，你的代码太麻烦了。

改成这样，简单又易读。

非常感谢你的建议。

下次问问题前，我得多仔细想想。

fcode

我亲爱的朋友，你的代码太麻烦了。

改成这样，简单又易读。

[Fortran] 纯文本查看 复制代码

program www_fcode_cn
  use ISO_C_Binding
  implicit none
  integer :: i 
  real :: an
  Open(12,File="INPUT.dat")
  Do i = 1 , 10
    Read(12,*)
  End Do
  read(12,*) an
  write(*,*) an
  Close(12)
end program www_fcode_cn

Kieran

对不起，我发现我的代码问题了。

READ (UNIT=3, FMT='(a)') an，这个语句的读取格式不对。应该是READ (UNIT=3, FMT=*) an

实在抱歉，问了这么乌龙的问题。我想删除这个帖子，但不知道怎么删除。所以就只能在这里回复下自己。耽误大家宝贵时间了。实在抱歉。