 |
 |
 |
| 24小時熱門版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2455 | 回復: 5 | |||
alvin510銀蟲 (小有名氣)
|
[求助]
Fortran編程過程中遇到錯誤，求大俠幫忙看下
|
|
我在用FORTRAN計算NLO過程中，編的一段程序遇到錯誤，請高手幫忙看下，萬分感謝。 其中的一段程序為： 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_energyallocate(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 其中調用了一個數據文件unit_bands.txt。這個文件我已經上傳了上來。 再次感謝你的幫助！ |
回復此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內容負責。
本內容由用戶自主發(fā)布，如果其內容涉及到知識產權問題，其責任在于用戶本人，如對版權有異議，請聯系郵箱：xiaomuchong@tal.com - 附件 1 : unit_bands.txt
- 附件 2 : CGA-01.txt
2013-07-18 18:42:01, 39.01 K
2013-07-18 18:43:19, 1.17 K
» 猜你喜歡
26博士申請
已經有3人回復
-
268求調劑
已經有3人回復
-
302求調劑
已經有5人回復
-
311求調劑
已經有10人回復
-
被我言中：新模板不強調格式了，假專家開始管格式了
已經有4人回復
-
303求調劑
已經有5人回復
-
考研求調劑
已經有4人回復
-
085601專碩，總分342求調劑，地區(qū)不限
已經有4人回復
-
085601材料工程專碩求調劑
已經有5人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調劑
已經有6人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求助：用高斯計算，總是出現2070錯誤代碼，請大俠們幫幫忙
已經有5人回復
- MATLAB, 代碼錯誤，請大俠幫忙看看。
已經有5人回復
- 大俠們~本人編譯vasp.5.lib過程中遇到如下錯誤 敢問何解？不勝感激~
已經有7人回復
- win8安裝Endnote失敗，顯示錯誤1723，求高手大俠解決啊
已經有4人回復
- VS20？？+intel visual fortran2011XE做并行計算的，能介紹一下經驗嗎？
已經有19人回復
- 急求fortran運行錯誤原因，在線等
已經有7人回復
- 求大俠幫忙解決一個B類錯誤！！謝謝
已經有5人回復
- Perl 生物信息學編程求助，運行出現問題，煩請幫忙查找一下錯誤，謝謝！
已經有3人回復
- 求高人幫忙看一下我這個Fortran程序的錯誤出在哪里？
已經有12人回復
- 有會看fortran程序的，能幫忙看看這個程序嗎？
已經有10人回復
- 我編的Simpson積分法fortran程序給不出結果，大俠們看看哪里出了問題？
已經有4人回復
- 求助Fortran編程
已經有8人回復
- MAC ox 下怎么用fortran 編程？
已經有9人回復
- 求助fortran編程時怎么使用數據庫
已經有4人回復
- 寫了一個fortran90的小程序，編譯通不過，請大俠幫忙
已經有59人回復
- 【求助】fortran編程問題
已經有3人回復
- 【求助】運行fortran程序出現錯誤 : error LNK2019 怎么回事啊
已經有8人回復
- 【求助】同樣一個問題用matlab和fortran編程計算，結果相差比較大，不知道為什么
已經有4人回復
- 【求助】用fortran怎么描述下面這樣的情況，我是新手，請大家?guī)兔?/a>
已經有8人回復
- 【求助】ANSYS 中命令流編程（fortran）
已經有11人回復
- 【求助】請問一下fortran的可視化編程
已經有7人回復
- 【分享】C/C++常用算法設計方法,fortran其他語言編程也可以參考下
已經有8人回復
alvin510
銀蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 1784.3
- 散金: 423
- 紅花: 4
- 帖子: 288
- 在線: 225.4小時
- 蟲號: 1245360
- 注冊: 2011-03-26
- 專業(yè): 人工晶體
2樓2013-07-18 18:45:13
jerkwin
專家顧問 (正式寫手)
-
專家經驗: +14 - 應助: 454 (碩士)
- 金幣: 20699.1
- 散金: 148
- 紅花: 81
- 帖子: 813
- 在線: 2648.3小時
- 蟲號: 1023452
- 注冊: 2010-05-19
- 專業(yè): 理論和計算化學
- 管轄: 分子模擬
3樓2013-07-18 21:10:42
alvin510
銀蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 1784.3
- 散金: 423
- 紅花: 4
- 帖子: 288
- 在線: 225.4小時
- 蟲號: 1245360
- 注冊: 2011-03-26
- 專業(yè): 人工晶體
|
謝謝你發(fā)現這個問題。后來我自己又發(fā)現了幾個問題一并做了更改，發(fā)現還是出現問題了。我的代碼如下： 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_energyallocate(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 ,*) 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 |
4樓2013-07-18 22:47:56
alvin510
銀蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 1784.3
- 散金: 423
- 紅花: 4
- 帖子: 288
- 在線: 225.4小時
- 蟲號: 1245360
- 注冊: 2011-03-26
- 專業(yè): 人工晶體
5樓2013-07-18 22:50:43
jerkwin
專家顧問 (正式寫手)
-
專家經驗: +14 - 應助: 454 (碩士)
- 金幣: 20699.1
- 散金: 148
- 紅花: 81
- 帖子: 813
- 在線: 2648.3小時
- 蟲號: 1023452
- 注冊: 2010-05-19
- 專業(yè): 理論和計算化學
- 管轄: 分子模擬
6樓2013-07-18 22:54:07
相關版塊跳轉
我要訂閱樓主 alvin510 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調劑 +5 | 呼呼呼。。。。 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 085601專碩，總分342求調劑，地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|
|
[考研] 301求調劑 +5 | yy要上岸呀 2026-03-17 | 5/250 |
|
|
|
[考研] 312求調劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
[考研] 本人考085602 化學工程 專碩 +16 | 不知道叫什么！ 2026-03-15 | 18/900 |
|
|
|
[考研] 求調劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學。調劑到任何專業(yè)都可以 +4 | 小丁想進步 2026-03-11 | 5/250 |
|
|
|
[考研] 0854可跨調劑，一作一項核心論文五項專利，省、國級證書40+數一英一287 +3 | 小李0854 2026-03-16 | 3/150 |
|
|
|
[考研] 285化工學碩求調劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[考研] 304求調劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 304求調劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 070303一志愿西北大學學碩310找調劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 材料工程327求調劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
|
[考研] 297一志愿上交085600求調劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 328求調劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 26調劑/材料/英一數二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] 四川大學085601材料工程專碩 初試294求調劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
|
[考研] 求調劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[論文投稿]
投稿問題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
|
[考研] 0856化學工程280分求調劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|