Kieran 发表于 2020-4-30 02:32:19

数据读取问题

大家好，

我想请教大家一个问题。我在读取一个文件（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。我读取该文件的部分程序如下。
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 发表于 2020-4-30 03:23:47

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

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

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

fcode 发表于 2020-4-30 08:06:53

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

改成这样，简单又易读。

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 发表于 2020-4-30 13:52:34

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

改成这样，简单又易读。


非常感谢你的建议。

下次问问题前，我得多仔细想想。
页: [1]
查看完整版本: 数据读取问题