版塊導航: 正在加載中...

調(diào)劑小程序

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (6119)
>考研 (1487)
>導師招生 (1255)
>文獻求助 (464)
>蟲友互識 (274)
>考博 (244)
>碩博家園 (133)
>博后之家 (111)
>基金申請 (108)
>招聘信息布告欄 (89)
>論文投稿 (85)
>休閑灌水 (68)
>第一性原理 (65)
>教師之家 (64)
>公派出國 (41)
>外文書籍求助 (32)

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

chk0521

金蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 863.1
散金: 34
紅花: 1
帖子: 45
在線: 84.3小時
蟲號: 494798
注冊: 2008-01-14
性別: GG
專業(yè): 半導體光電子器件

[交流] 【求助成功】Band－procar 和Dos－procar 程序執(zhí)行出錯已有8人參與

我從版內(nèi)下載的分析PROCAR的程序Band－procar.f 和Dos－procar.f 通過g77 編譯成可執(zhí)行文件。但是執(zhí)行后出現(xiàn)segmentation fault (core dumped)錯誤提示。這是為什么�。课疑暇W(wǎng)查了一下這表示指針使用出錯。我的編程水平太差實在查不出問題在哪。麻煩哪位高人指點一下，十分感謝！以下是源程序
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 ]

回復此樓

» 猜你喜歡

0854電子信息求調(diào)劑已經(jīng)有7人回復
0854人工智能方向招收調(diào)劑已經(jīng)有3人回復
調(diào)劑已經(jīng)有4人回復
一志愿北化315 求調(diào)劑已經(jīng)有3人回復
有機合成求助已經(jīng)有5人回復
材料學碩333求調(diào)劑已經(jīng)有3人回復
289求調(diào)劑已經(jīng)有7人回復
收08調(diào)劑生已經(jīng)有6人回復
材料專碩找調(diào)劑已經(jīng)有5人回復
085602 289分求調(diào)劑已經(jīng)有5人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

二硫蘇糖醇加入還原性膠束中的具體操作方法和加入量。已經(jīng)有10人回復
請問高分子和有機藥物合成對身體影響大嗎已經(jīng)有8人回復
求助：巖石薄片巖礦分析和XRD得到的各礦物成分及含量有較大差異，怎么解釋？已經(jīng)有6人回復
循環(huán)伏安循環(huán)極化和動電位再活化曲線有什么分別？已經(jīng)有9人回復
水文局和海事局，知道的來說一下已經(jīng)有13人回復
得不到純產(chǎn)物，哪出錯了，親們幫忙看看啊已經(jīng)有3人回復
一個二倍體植物的全基因組測序程序及花費？已經(jīng)有18人回復
用c語言編程怎來表達時滯！謝謝！已經(jīng)有6人回復
請問學位論文查重會不會和自己的文章重復呀？已經(jīng)有11人回復
Fortran的subroutine和function有什么區(qū)別，是在于輸入輸出的特點嗎？已經(jīng)有13人回復
關于CASTEP 中計算時出錯已經(jīng)有6人回復
出錯 warning#2 #180 已經(jīng)有7人回復
MS一打開就出錯，任何任務都不能運行已經(jīng)有12人回復
DOS圖中怎么判斷非金屬性，半金屬性，金屬性已經(jīng)有7人回復
請教：orca 中用CCSD(T)計算單點能出錯--提示內(nèi)存不足該怎么辦？已經(jīng)有10人回復
MS軟件老出錯，求指導已經(jīng)有3人回復
制備丙烯酸乳液，十二烷基磺酸鈉和十二烷基苯磺酸鈉能否通用，大俠們幫幫忙啊~~~ 已經(jīng)有8人回復

1樓 2010-06-21 18:59:33

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

repent03

金蟲 (小有名氣)

1ST強帖: 1
應助: 3 (幼兒園)
金幣: 1339.3
散金: 6
帖子: 87
在線: 40小時
蟲號: 1068065
注冊: 2010-08-02
專業(yè): 半導體材料

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

據(jù)本身前經(jīng)驗只要修改parameter里面的數(shù)字增大即可

贊一下

回復此樓

9樓2012-03-22 00:45:22

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

查看全部 9 個回答

lvjian8596

木蟲 (小有名氣)

1ST強帖: 1
應助: 0 (幼兒園)
貴賓: 0.014
金幣: 3531.7
紅花: 2
帖子: 296
在線: 209.1小時
蟲號: 545009
注冊: 2008-04-13
性別: GG
專業(yè): 凝聚態(tài)物性I:結構、力學和

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+1):感謝交流 2010-06-24 08:50:43

segmentation fault (core dumped)出現(xiàn)這個問題很有可能是數(shù)組越界！
程序本身沒有問題，但是當你的數(shù)據(jù)過大時可能出現(xiàn)問題！
比如一個數(shù)組a（100），但是你有200個數(shù)據(jù)讀進來，這時可能出現(xiàn)這個問題！

贊一下(2人)

回復此樓

2樓2010-06-24 08:44:53

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

chk0521

金蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 863.1
散金: 34
紅花: 1
帖子: 45
在線: 84.3小時
蟲號: 494798
注冊: 2008-01-14
性別: GG
專業(yè): 半導體光電子器件

zzy870720z:如果問題解決，幫忙把【求助】改為【求助成功】，謝謝合作，如果在整理出解決辦法，給予獎勵，呵呵 2010-06-24 17:32:28

謝謝斑竹！確實是數(shù)組越界。已經(jīng)解決了:）

[ Last edited by chk0521 on 2010-6-25 at 16:37 ]

贊一下(1人)

回復此樓

3樓2010-06-24 15:51:51

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

calos818

木蟲 (著名寫手)

新人菜鳥

應助: 11 (小學生)
金幣: 3125.7
散金: 77
紅花: 3
帖子: 1201
在線: 730.8小時
蟲號: 698507
注冊: 2009-02-08
專業(yè): 凝聚態(tài)物性 II ：電子結構

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流

引用回帖:

Originally posted by chk0521 at 2010-06-24 15:51:51:
謝謝斑竹！確實是數(shù)組越界。已經(jīng)解決了:）

[ Last edited by chk0521 on 2010-6-25 at 16:37 ]

請問你做了哪些修改的，麻煩你提示下，謝謝

贊一下(1人)

回復此樓

4樓2010-10-24 18:23:05

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

查看全部 9 個回答

普通表情龍兔虎貓高級回復 (可上傳附件)

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

[考研] 0703化學調(diào)劑，求導師收 +6	天天好運來上岸?/a> 2026-03-24	6/300	2026-03-24 17:19 by dick_runner
[考研] 291求調(diào)劑 +3	HanBeiNingZC 2026-03-24	3/150	2026-03-24 16:34 by barlinike
[考研] 086003食品工程求調(diào)劑 +4	淼淼111 2026-03-24	4/200	2026-03-24 11:53 by 544594351
[考研] 279分求調(diào)劑一志愿211 +18	chaojifeixia 2026-03-19	20/1000	2026-03-24 10:34 by dolphin_ycj
[考研] 一志愿河北工業(yè)大學0817化工278分求調(diào)劑 +7	jhybd 2026-03-23	12/600	2026-03-24 09:03 by jhybd
[考研] 284求調(diào)劑 +10	Zhao anqi 2026-03-22	10/500	2026-03-24 00:08 by Equinoxhua
[考研] 335分 \| 材料與化工專碩 \| GPA 4.07 \| 有科研經(jīng)歷 +4	cccchenso 2026-03-23	4/200	2026-03-23 23:00 by 徐ckkk
[考研] 材料專碩英一數(shù)二306 +8	z1z2z3879 2026-03-18	8/400	2026-03-23 20:49 by baobaoye
[考研] 生物學一志愿985，分數(shù)349求調(diào)劑 +6	zxts12 2026-03-21	9/450	2026-03-23 18:37 by macy2011
[考研] 一志愿東華大學化學070300，求調(diào)劑 +7	2117205181 2026-03-21	8/400	2026-03-22 22:55 by chixmc
[考研] 308求調(diào)劑 +3	墨墨漠 2026-03-21	3/150	2026-03-22 16:54 by i_cooler
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學員8dgXkO
[考研] 302求調(diào)劑 +12	呼呼呼。。。。 2026-03-17	12/600	2026-03-21 17:29 by ColorlessPI
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 0817 化學工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

chk0521

[交流] 【求助成功】Band－procar 和Dos－procar 程序執(zhí)行出錯 已有8人參與

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

repent03

lvjian8596

chk0521

calos818

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[交流] 【求助成功】Band－procar 和Dos－procar 程序執(zhí)行出錯已有8人參與

» 本主題相關價值貼推薦，對您同樣有幫助: