 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 9 | 1/1 | 返回列表 |
| 查看: 2952 | 回復: 8 | |||
[交流]
【求助成功】Band－procar 和Dos－procar 程序執(zhí)行出錯 已有8人參與
|
|
我從版內下載的分析PROCAR的程序Band－procar.f 和Dos－procar.f 通過g77 編譯成可執(zhí)行文件。但是執(zhí)行后出現(xiàn)segmentation fault (core dumped)錯誤提示。這是為什么��？我上網查了一下這表示指針使用出錯。我的編程水平太差實在查不出問題在哪。麻煩哪位高人指點一下，十分感謝！以下是源程序 implicit real*8(a-h,o-z) parameter (nbd = 200) parameter (nkd = 500) parameter (nxd = 300) parameter (natmd = 20) dimension a(3,3),b(3,3),c(3),e(nkd,nbd),sk(nkd,3) dimension xx(nxd) ,wei(nkd) dimension dump(20),oc(nkd,nbd,natmd,4) open(7,file='PROCAR',form='FORMATTED',status='OLD') pi = 3.141592654 read(7,103) dump write(*,*) 'Spin polarized calculation? (no=1,yes=2):' read (*,*) ispin if ((ispin.ne.1).and.(ispin.ne.2)) then write(*,*) ' INPUT ERROR, ispin must equal to 1 or 2 ' stop endif c write(*,*) 'Enter # of interval (npoints) and division (ndiv):' c read (*,*) npoints,ndiv open(9,file='KPOINTS',form='FORMATTED',status='OLD') read(9,100) temp read(9,*) ndiv write(*,*) 'Enter the range of energy to plot:' read (*,*) er1,er2 emin=min(er1,er2) emax=max(er1,er2) write(*,*) 'Enter the value of fermi energy:' read(*,*) ef if (ispin.eq.1) then open(11,file='band.dat') elseif (ispin.eq.2) then open(11,file='band-up.dat') open(12,file='band-dn.dat') endif open(8,file='POSCAR',form='FORMATTED',status='OLD') read(8,100) temp c write(6,100) temp 100 format(20a4) read (8,*) aa c WRITE(6,*) aa c c *** read lattice constant from POSCAR** c do i=1,3 read (8,*) (a(i,j),j=1,3) c WRITE(6,500) (a(i,j),j=1,3) 500 format (3f12.8) enddo do i=1,3 do j=1,3 a(i,j)=aa*a(i,j) enddo c WRITE(6,500) (a(i,j),j=1,3) enddo c c *** read lattice vector from POSCAR*** c volume=a(1,1)*a(2,2)*a(3,3)+a(1,2)*a(2,3)*a(3,1) &+a(1,3)*a(2,1)*a(3,2)-a(1,1)*a(2,3)*a(3,2) &-a(1,2)*a(2,1)*a(3,3)-a(1,3)*a(2,2)*a(3,1) do i=1,3 if (i .eq. 1) then j=2 k=3 else if (i .eq. 2) then j=3 k=1 else j=1 k=2 endif c(1)=a(j,2)*a(k,3)-a(j,3)*a(k,2) c(2)=a(j,3)*a(k,1)-a(j,1)*a(k,3) c(3)=a(j,1)*a(k,2)-a(j,2)*a(k,1) do j=1,3 b(i,j)=2*pi*c(j)/volume c WRITE (6,*) b(i,j) enddo enddo do 9000 isp=1,ispin read(7,104) nk,nband,nion do 1000 k = 1,nk read(7,103) dump read(7,105) kp,(sk(k,j),j=1,3),wei(k) c write(6,105) kp,(sk(k,j),j=1,3),wei(k) read(7,103) dump do nb = 1,nband read(7,106) nb1,e(k,nb),occ c write(6,106) nb1,e(k,nb),occ read(7,103) dump read(7,103) dump c write(6,*) 'nion=',nion niont = nion +1 if (nion .eq. 1) niont = 1 do ion = 1,niont read(7,107) (oc(k,nb,ion,j),j=1,4) c write(6,107) (oc(k,nb,ion,j),j=1,4) enddo read(7,103) dump c write(6,103) dump enddo 1000 continue weight = 0.0 do k = 1, nk weight = weight + wei(k) enddo do k = 1,nk wei(k) = wei(k) / weight enddo 101 format(10x,f9.5) 102 format(f10.5) 103 format(20a4) 104 format(16x,i3,20x,i5,19x,i4) 105 format(10x,i3,5x,3f11.8,13x,f11.8) 106 format(4x,i4,9x,f14.8,7x,f12.8) 107 format(3x,4f7.3) c c *** find reciprocal lattice vector *** xx(1) = 0.0 nn = 1 do k = 1,nk-1 dkx=(sk(k+1,1)-sk(k,1))*b(1,1) + (sk(k+1,2)-sk(k,2))*b(2,1) & + (sk(k+1,3)- sk(k,3))*b(3,1) dky=(sk(k+1,1)-sk(k,1))*b(1,2) + (sk(k+1,2)-sk(k,2))*b(2,2) & + (sk(k+1,3)- sk(k,3))*b(3,2) dkz=(sk(k+1,1)-sk(k,1))*b(1,3) + (sk(k+1,2)-sk(k,2))*b(2,3) & + (sk(k+1,3)- sk(k,3))*b(3,3) del = sqrt ( dkx**2 + dky**2 + dkz**2 ) nn = nn +1 xx(nn) = xx(nn-1) + del enddo do n=1,nband if (mod(n,2).ne.0) then do k=1,nk ee = e(k,n) - ef if ( ee .gt. emax ) ee = emax if ( ee .lt. emin ) ee = emin write (10+isp,300) xx(k),ee enddo elseif (mod(n,2).eq.0) then do i=nk,1,-1 ee = e(i,n) - ef if ( ee .gt. emax ) ee = emax if ( ee .lt. emin ) ee = emin write (10+isp,300) xx(i),ee enddo endif enddo 300 format (f12.8,2x,f12.8) if (mod(nband,2) .ne. 0) then write (10+isp,300) xx(nk),emin write (10+isp,300) xx(1),emin else write (10+isp,300) xx(1),emin endif c c *** write xx-ee *** c npoints=nk/ndiv do n=2,npoints kk=(n-1)*ndiv write (10+isp,300) xx(kk),emin write (10+isp,300) xx(kk),emax write (10+isp,300) xx(kk),emin enddo write (10+isp,300) xx(nk),emin write (10+isp,300) xx(nk),emax write (10+isp,300) xx(1),emax write (10+isp,300) xx(1),emin zero=0.0 write (10+isp,300) xx(1),zero write (10+isp,300) xx(nk),zero 9000 continue stop end [ Last edited by chk0521 on 2010-6-25 at 16:36 ] |
回復此樓» 猜你喜歡
0817 化學工程 299分求調劑 有科研經歷 有二區(qū)文章
已經有4人回復
-
0854可跨調劑，一作一項核心論文五項專利，省、國級證書40+數一英一287
已經有5人回復
-
070300化學319求調劑
已經有6人回復
-
一志愿985，本科211，0817化學工程與技術319求調劑
已經有6人回復
-
一志愿西南交大，求調劑
已經有3人回復
-
286求調劑
已經有6人回復
-
焦慮
已經有11人回復
-
331求調劑（0703有機化學
已經有7人回復
-
298-一志愿中國農業(yè)大學-求調劑
已經有5人回復
-
297求調劑
已經有5人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 二硫蘇糖醇加入還原性膠束中的具體操作方法和加入量。
已經有10人回復
- 請問高分子和有機藥物合成對身體影響大嗎
已經有8人回復
- 求助：巖石薄片巖礦分析和XRD得到的各礦物成分及含量有較大差異，怎么解釋？
已經有6人回復
- 循環(huán)伏安 循環(huán)極化和動電位再活化曲線有什么分別？
已經有9人回復
- 水文局和海事局，知道的來說一下
已經有13人回復
- 得不到純產物，哪出錯了，親們幫忙看看啊
已經有3人回復
- 一個二倍體植物的全基因組測序程序及花費？
已經有18人回復
- 用c語言編程怎來表達時滯！謝謝！
已經有6人回復
- 請問學位論文查重會不會和自己的文章重復呀？
已經有11人回復
- Fortran的subroutine和function有什么區(qū)別，是在于輸入輸出的特點嗎？
已經有13人回復
- 關于CASTEP 中計算時出錯
已經有6人回復
- 出錯 warning#2 #180
已經有7人回復
- MS一打開就出錯，任何任務都不能運行
已經有12人回復
- DOS圖中怎么判斷非金屬性，半金屬性，金屬性
已經有7人回復
- 請教：orca 中用CCSD(T)計算單點能出錯--提示內存不足 該怎么辦？
已經有10人回復
- MS軟件老出錯，求指導
已經有3人回復
- 制備丙烯酸乳液，十二烷基磺酸鈉和十二烷基苯磺酸鈉能否通用，大俠們幫幫忙啊~~~
已經有8人回復
lvjian8596
木蟲 (小有名氣)
- 1ST強帖: 1
- 應助: 0 (幼兒園)
- 貴賓: 0.014
- 金幣: 3531.7
- 紅花: 2
- 帖子: 296
- 在線: 209.1小時
- 蟲號: 545009
- 注冊: 2008-04-13
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
2樓2010-06-24 08:44:53
3樓2010-06-24 15:51:51
calos818
木蟲 (著名寫手)
新人菜鳥
- 應助: 11 (小學生)
- 金幣: 3125.7
- 散金: 77
- 紅花: 3
- 帖子: 1201
- 在線: 730.8小時
- 蟲號: 698507
- 注冊: 2009-02-08
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
4樓2010-10-24 18:23:05
huangyc
木蟲 (正式寫手)
- 應助: 2 (幼兒園)
- 金幣: 2242.1
- 散金: 12
- 紅花: 2
- 帖子: 411
- 在線: 312.5小時
- 蟲號: 624241
- 注冊: 2008-10-12
- 專業(yè): 理論和計算化學
5樓2010-10-25 22:09:24
6樓2011-09-21 00:19:44
luding861011
木蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 1923.2
- 紅花: 2
- 帖子: 22
- 在線: 1207.1小時
- 蟲號: 1193967
- 注冊: 2011-01-21
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
7樓2011-09-21 11:21:52
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
你好，關于Band－procar.f 和Dos－procar.f ,請問你是怎么解決“segmentation fault (core dumped)這個問題的，即數組越界： http://www.gaoyang168.com/bbs/viewthread.php?tid=2158806 非常感謝而你！ lyyqdu |
8樓2012-03-21 11:07:28
repent03
金蟲 (小有名氣)
- 1ST強帖: 1
- 應助: 3 (幼兒園)
- 金幣: 1339.3
- 散金: 6
- 帖子: 87
- 在線: 40小時
- 蟲號: 1068065
- 注冊: 2010-08-02
- 專業(yè): 半導體材料
9樓2012-03-22 00:45:22
相關版塊跳轉
我要訂閱樓主 chk0521 的主題更新
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[教師之家] 焦慮 +8 | 水冰月月野兔 2026-03-13 | 11/550 |
|
|---|---|---|---|---|
|
|
[考研] 0703化學求調劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 268求調劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
|
[考研] 277調劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 機械專碩325，尋找調劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 一志愿211 0703方向310分求調劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 070300化學學碩求調劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 中科院材料273求調劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研]
 機械專碩調劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[考研] 22408總分284求調劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 材料工程327求調劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
|
[考研] 學碩285求調劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 266求調劑 +4 | 學員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 304求調劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 求材料調劑 085600英一數二總分302 前三科235 精通機器學習 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] 310求調劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 求調劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 290求調劑 +7 | ADT 2026-03-12 | 7/350 |
|