當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

alvin510

銀蟲 (小有名氣)

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

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

我在用FORTRAN計(jì)算NLO過(guò)程中，編的一段程序遇到錯(cuò)誤，請(qǐng)高手幫忙看下，萬(wàn)分感謝。
其中的一段程序?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)上傳了上來(lái)。
再次感謝你的幫助！

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wè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

» 猜你喜歡

數(shù)一英一271專碩（085401）求調(diào)劑，可跨已經(jīng)有8人回復(fù)
295求調(diào)劑已經(jīng)有4人回復(fù)
0856求調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有10人回復(fù)
299求調(diào)劑已經(jīng)有6人回復(fù)
一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑已經(jīng)有6人回復(fù)
289求調(diào)劑已經(jīng)有12人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑，求導(dǎo)師收已經(jīng)有10人回復(fù)
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩已經(jīng)有5人回復(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過(guò)程中遇到如下錯(cuò)誤敢問(wèn)何解？不勝感激~ 已經(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)問(wè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é)果，大俠們看看哪里出了問(wèn)題？已經(jīng)有4人回復(fù)
求助Fortran編程已經(jīng)有8人回復(fù)
MAC ox 下怎么用fortran 編程？已經(jīng)有9人回復(fù)
求助fortran編程時(shí)怎么使用數(shù)據(jù)庫(kù) 已經(jīng)有4人回復(fù)
寫了一個(gè)fortran90的小程序，編譯通不過(guò)，請(qǐng)大俠幫忙已經(jīng)有59人回復(fù)
【求助】fortran編程問(wèn)題已經(jīng)有3人回復(fù)
【求助】運(yùn)行fortran程序出現(xiàn)錯(cuò)誤 : error LNK2019 怎么回事啊已經(jīng)有8人回復(fù)
【求助】同樣一個(gè)問(wèn)題用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)問(wèn)一下fortran的可視化編程已經(jīng)有7人回復(fù)
【分享】C/C++常用算法設(shè)計(jì)方法,fortran其他語(yǔ)言編程也可以參考下已經(jīng)有8人回復(fù)

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

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

jerkwin

專家顧問(wèn) (正式寫手)

專家經(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沒(méi)有定義啊

贊一下

回復(fù)此樓

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

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

查看全部 6 個(gè)回答

alvin510

銀蟲 (小有名氣)

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

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

贊一下

回復(fù)此樓

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

已閱回復(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沒(méi)有定義啊

謝謝你發(fā)現(xiàn)這個(gè)問(wèn)題。后來(lái)我自己又發(fā)現(xiàn)了幾個(gè)問(wèn)題一并做了更改，發(fā)現(xiàn)還是出現(xiàn)問(wè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的回帖

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

[考研] 求調(diào)劑 +7	爭(zhēng)取九點(diǎn)睡 2026-03-28	8/400	2026-03-28 21:07 by 爭(zhēng)取九點(diǎn)睡
[考研] 275求調(diào)劑 +10	Micky11223 2026-03-25	14/700	2026-03-28 15:48 by Micky11223
[考研] 085600 286分材料求調(diào)劑 +7	麻辣魷魚 2026-03-27	8/400	2026-03-28 12:17 by zllcz
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 材料求調(diào)劑一志愿哈工大324 +7	閆旭東 2026-03-28	9/450	2026-03-28 08:51 by Xu de nuo
[有機(jī)交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +9	vv迷 2026-03-22	9/450	2026-03-27 15:59 by 不吃魚的貓
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當(dāng)當(dāng)
[考研] 305求調(diào)劑 +5	哇盧卡庫(kù) 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機(jī) 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 一志愿南航 335分 \| 0856材料化工 \| GPA 4.07 \| 有科研經(jīng)歷 +6	cccchenso 2026-03-23	6/300	2026-03-25 22:25 by 544594351
[考研] 材料與化工304求B區(qū)調(diào)劑 +3	邱gl 2026-03-25	3/150	2026-03-25 19:03 by Ainin_
[考研] 材料專碩 335 分求調(diào)劑 +4	拒絕冷暴力 2026-03-25	4/200	2026-03-25 18:45 by haxia
[考研] 求b區(qū)院校調(diào)劑 +4	周56 2026-03-24	5/250	2026-03-25 17:12 by yishunmin
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 292求調(diào)劑 +4	鵝鵝鵝額額額額?/a> 2026-03-24	4/200	2026-03-24 16:41 by peike
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[基金申請(qǐng)] 請(qǐng)教下大家 2026年國(guó)家基金申請(qǐng)是雙盲審嗎？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch