當(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編程過程中遇到錯(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

» 猜你喜歡

329求調(diào)劑已經(jīng)有8人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有8人回復(fù)
349求調(diào)劑已經(jīng)有5人回復(fù)
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有7人回復(fù)
291求調(diào)劑已經(jīng)有4人回復(fù)
0856材料化工調(diào)劑總分330 已經(jīng)有14人回復(fù)
352分-085602-一志愿985 已經(jīng)有5人回復(fù)
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩已經(jīng)有8人回復(fù)
求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分已經(jīng)有4人回復(fù)
279求調(diào)劑已經(jīng)有4人回復(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è): 人工晶體

引用回帖:

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的回帖

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

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è): 人工晶體

提示出錯(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ā)表

[考研] 297求調(diào)劑 +7	田洪有 2026-03-26	7/350	2026-03-29 09:19 by lbsjt
[考研] 321求調(diào)劑 +7	璞玉~~ 2026-03-25	8/400	2026-03-29 06:41 by 544594351
[碩博家園] 求調(diào)劑有機(jī)化學(xué)考研356分 +3	Nadiums 2026-03-25	4/200	2026-03-29 00:50 by 544594351
[考研] 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。 +7	nnnnnnn5 2026-03-25	7/350	2026-03-29 00:32 by 我是小康
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3	llss0711 2026-03-28	5/250	2026-03-29 00:23 by MRMFYTL
[考研] 295材料工程專碩求調(diào)劑 +7	1428151015 2026-03-27	7/350	2026-03-28 19:58 by S240
[考研] 求調(diào)劑 +3	QiMing7 2026-03-25	4/200	2026-03-28 14:30 by QiMing7
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9	丹青奶蓋 2026-03-26	10/500	2026-03-28 07:45 by barnett0632
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 求調(diào)劑 +6	林之夕 2026-03-24	6/300	2026-03-27 08:38 by hypershenger
[考研] 材料調(diào)劑 +8	匹克i 2026-03-23	8/400	2026-03-27 08:11 by hypershenger
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 327求調(diào)劑 +7	prayer13 2026-03-23	7/350	2026-03-26 20:48 by 不吃魚的貓
[考研] 325求調(diào)劑 +3	Aoyijiang 2026-03-23	3/150	2026-03-26 20:46 by 不吃魚的貓
[考研] 一志愿南航 335分 \| 0856材料化工 \| GPA 4.07 \| 有科研經(jīng)歷 +6	cccchenso 2026-03-23	6/300	2026-03-25 22:25 by 544594351
[考研] 一志愿北化315 求調(diào)劑 +3	akrrain 2026-03-24	3/150	2026-03-24 19:35 by 了了了了。。
[考研] 化工專碩求調(diào)劑 +3	question挽風(fēng) 2026-03-24	3/150	2026-03-24 18:48 by jhhcooi
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋懷