alvin510

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1784.3
散金: 423
紅花: 4
帖子: 288
在線: 225.4小時(shí)
蟲號(hào): 1245360
注冊(cè): 2011-03-26
專業(yè): 人工晶體

[求助] Fortran編程過程中遇到錯(cuò)誤，求大俠幫忙看下

我在用FORTRAN計(jì)算NLO過程中，編的一段程序遇到錯(cuò)誤，請(qǐng)高手幫忙看下，萬分感謝。
其中的一段程序?yàn)椋?br /> module def
implicit none
integer :: nk,nkpts,ndum,ns,nspins,nb,locat_str
character (len=80) :: dummy
real num_e
real fermi_energy
integer num_bands
real,allocatable :: weight(

real,allocatable :: kpoints(:,

end module

program main

!讀入 kpoints, weights 和 eigenvalues
use def
implicit none
integer,parameter :: fileid=10
character (len=80) :: filename="F:\FORTRAN\NLOC\CGA\unit_bands.txt"
logical alive

inquire (file=filename,exist=alive)
if(.not. alive) then
write(*,*) trim(filename), "doesn't exist!"
stop
end if

open(unit=10,file=filename)
read(unit=10,"(19XI6,/,26XI2,/,20XF7.1,/,22XI7,/,31XF13.8)"

nkpts,spins,num_e,num_bands,fermi_energy
allocate(kpoints(1,nkpts))
allocate(kpoints(2,nkpts))
allocate(kpoints(3,nkpts))
allocate(weight(nkpts))

do nk = 1, nkpts
read(fileid, '(a)') dummy
locat_str = index(dummy,'K-point')
read(dummy(locat_str+7

,*) ndum, kpoints(1,nk), kpoints(2,nk), kpoints(3,nk), weight(nk)
do ns = 1, nspins
read(fileid,*) dummy
do nb = 1, num_bands
read(fileid,*) eigenvalues(nb, nk, ns)
enddo
enddo
enddo
close(fileid)

stop
end program
其中調(diào)用了一個(gè)數(shù)據(jù)文件unit_bands.txt。這個(gè)文件我已經(jīng)上傳了上來。
再次感謝你的幫助！

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : unit_bands.txt

2013-07-18 18:42:01, 39.01 K

附件 2 : CGA-01.txt

2013-07-18 18:43:19, 1.17 K

» 猜你喜歡

312求調(diào)劑已經(jīng)有7人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有16人回復(fù)
280求調(diào)劑已經(jīng)有3人回復(fù)
266求調(diào)劑已經(jīng)有7人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有5人回復(fù)
環(huán)境工程調(diào)劑已經(jīng)有8人回復(fù)
302求調(diào)劑已經(jīng)有9人回復(fù)
297求調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有4人回復(fù)
焦慮已經(jīng)有10人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助：用高斯計(jì)算，總是出現(xiàn)2070錯(cuò)誤代碼，請(qǐng)大俠們幫幫忙已經(jīng)有5人回復(fù)
MATLAB, 代碼錯(cuò)誤，請(qǐng)大俠幫忙看看。已經(jīng)有5人回復(fù)
大俠們~本人編譯vasp.5.lib過程中遇到如下錯(cuò)誤敢問何解？不勝感激~ 已經(jīng)有7人回復(fù)
win8安裝Endnote失敗，顯示錯(cuò)誤1723，求高手大俠解決啊已經(jīng)有4人回復(fù)
VS20？？+intel visual fortran2011XE做并行計(jì)算的，能介紹一下經(jīng)驗(yàn)嗎？已經(jīng)有19人回復(fù)
急求fortran運(yùn)行錯(cuò)誤原因，在線等已經(jīng)有7人回復(fù)
求大俠幫忙解決一個(gè)B類錯(cuò)誤��！謝謝已經(jīng)有5人回復(fù)
Perl 生物信息學(xué)編程求助，運(yùn)行出現(xiàn)問題，煩請(qǐng)幫忙查找一下錯(cuò)誤，謝謝！已經(jīng)有3人回復(fù)
求高人幫忙看一下我這個(gè)Fortran程序的錯(cuò)誤出在哪里？已經(jīng)有12人回復(fù)
有會(huì)看fortran程序的，能幫忙看看這個(gè)程序嗎？已經(jīng)有10人回復(fù)
我編的Simpson積分法fortran程序給不出結(jié)果，大俠們看看哪里出了問題？已經(jīng)有4人回復(fù)
求助Fortran編程已經(jīng)有8人回復(fù)
MAC ox 下怎么用fortran 編程？已經(jīng)有9人回復(fù)
求助fortran編程時(shí)怎么使用數(shù)據(jù)庫已經(jīng)有4人回復(fù)
寫了一個(gè)fortran90的小程序，編譯通不過，請(qǐng)大俠幫忙已經(jīng)有59人回復(fù)
【求助】fortran編程問題已經(jīng)有3人回復(fù)
【求助】運(yùn)行fortran程序出現(xiàn)錯(cuò)誤 : error LNK2019 怎么回事啊已經(jīng)有8人回復(fù)
【求助】同樣一個(gè)問題用matlab和fortran編程計(jì)算，結(jié)果相差比較大，不知道為什么已經(jīng)有4人回復(fù)
【求助】用fortran怎么描述下面這樣的情況，我是新手，請(qǐng)大家?guī)兔?/a> 已經(jīng)有8人回復(fù)

【求助】ANSYS 中命令流編程（fortran）已經(jīng)有11人回復(fù)
【求助】請(qǐng)問一下fortran的可視化編程已經(jīng)有7人回復(fù)
【分享】C/C++常用算法設(shè)計(jì)方法,fortran其他語言編程也可以參考下已經(jīng)有8人回復(fù)

1樓 2013-07-18 18:43:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

alvin510

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1784.3
散金: 423
紅花: 4
帖子: 288
在線: 225.4小時(shí)
蟲號(hào): 1245360
注冊(cè): 2011-03-26
專業(yè): 人工晶體

大家還是看我附件中的代碼吧，網(wǎng)頁上的出了點(diǎn)錯(cuò)誤

贊一下

回復(fù)此樓

2樓2013-07-18 18:45:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jerkwin

專家顧問 (正式寫手)

專家經(jīng)驗(yàn): +14
應(yīng)助: 454 (碩士)
金幣: 20699.1
散金: 148
紅花: 81
帖子: 813
在線: 2648.3小時(shí)
蟲號(hào): 1023452
注冊(cè): 2010-05-19
專業(yè): 理論和計(jì)算化學(xué)
管轄: 分子模擬

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

eigenvalues沒有定義啊

贊一下

回復(fù)此樓

3樓2013-07-18 21:10:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

alvin510

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1784.3
散金: 423
紅花: 4
帖子: 288
在線: 225.4小時(shí)
蟲號(hào): 1245360
注冊(cè): 2011-03-26
專業(yè): 人工晶體

引用回帖:

3樓: Originally posted by jerkwin at 2013-07-18 21:10:42
eigenvalues沒有定義啊

謝謝你發(fā)現(xiàn)這個(gè)問題。后來我自己又發(fā)現(xiàn)了幾個(gè)問題一并做了更改，發(fā)現(xiàn)還是出現(xiàn)問題了。我的代碼如下：
module def
implicit none
integer :: nk,nkpts,ndum,ns,nspins,nb,locat_str
character (len=80) :: dummy
real num_e
real fermi_energy
integer num_bands
real,allocatable :: weight(

real,allocatable :: kpoints(:,

real,allocatable :: eigenvalues(:,:,

end module

program main

!讀入 kpoints, weights 和 eigenvalues
use def
implicit none
integer,parameter :: fileid=10
character (len=80) :: filename="F:\FORTRAN\NLOC\CGA\unit_bands.txt"
logical alive

inquire (file=filename,exist=alive)
if(.not. alive) then
write(*,*) trim(filename), "doesn't exist!"
stop
end if

open(unit=fileid,file=filename)
read(unit=fileid,"(19XI6,/,26XI2,/,20XF7.1,/,22XI7,/,31XF13.8)"

nkpts,nspins,num_e,num_bands,fermi_energy
!read(unit=fileid,"(19XI6)"

nkpts
!read(unit=fileid,"(26XI2)"

nspins
!read(unit=fileid,"(20XF7.1)"

num_e
!read(unit=fileid,"(22XI7)"

num_bands
!read(unit=fileid,"(31XF13.8)"

fermi_energy
allocate(kpoints(3,nkpts))
allocate(weight(nkpts))
allocate(eigenvalues(num_bands,nkpts,nspins))

do nk = 1, nkpts
read(fileid, '(a)') dummy
locat_str = index(dummy,'K-point')
read(dummy(locat_str+7

贊一下

回復(fù)此樓

4樓2013-07-18 22:47:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

alvin510

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1784.3
散金: 423
紅花: 4
帖子: 288
在線: 225.4小時(shí)
蟲號(hào): 1245360
注冊(cè): 2011-03-26
專業(yè): 人工晶體

提示出錯(cuò)的地方在
read(unit=fileid,"(19XI6,/,26XI2,/,20XF7.1,/,22XI7,/,31XF13.8)"

nkpts,nspins,num_e,num_bands,fermi_energy
這一行。
錯(cuò)誤 1 error #6357: This is not a valid specifier for a control list. ['(19XI6,/,26XI2,/,20XF7.1,/,22XI7,/,31XF13.] F:\FORTRAN\NLOC\CGA\CGA-01.f90 30

請(qǐng)幫我查查看原因，謝謝

贊一下

回復(fù)此樓

5樓2013-07-18 22:50:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jerkwin

專家顧問 (正式寫手)

專家經(jīng)驗(yàn): +14
應(yīng)助: 454 (碩士)
金幣: 20699.1
散金: 148
紅花: 81
帖子: 813
在線: 2648.3小時(shí)
蟲號(hào): 1023452
注冊(cè): 2010-05-19
專業(yè): 理論和計(jì)算化學(xué)
管轄: 分子模擬

【答案】應(yīng)助回帖

用你注釋掉的那種寫法，如
read(unit=fileid,×) nkpts

贊一下

回復(fù)此樓

6樓2013-07-18 22:54:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 312求調(diào)劑 +6	陌宸希 2026-03-16	7/350	2026-03-18 09:49 by zhukairuo
[考研] 化學(xué)工程321分求調(diào)劑 +13	大米飯！ 2026-03-15	16/800	2026-03-18 09:48 by zhukairuo
[考研] 工科材料085601 279求調(diào)劑 +5	困于星晨 2026-03-17	5/250	2026-03-18 09:39 by zhukairuo
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 070300化學(xué)319求調(diào)劑 +4	錦鯉0909 2026-03-17	4/200	2026-03-17 18:21 by 重科小霸王
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 274求調(diào)劑0856材料化工 +13	z2839474511 2026-03-11	14/700	2026-03-17 16:51 by share_joy
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4	小丁想進(jìn)步 2026-03-11	5/250	2026-03-17 16:05 by 外星文明
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 308 085701 四六級(jí)已過求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 材料工程調(diào)劑 +9	咪咪空空 2026-03-12	9/450	2026-03-13 22:05 by 星空星月
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 工科調(diào)劑 +4	Jiang191123！ 2026-03-11	4/200	2026-03-13 15:15 by Miko19