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

astringent

銅蟲 (著名寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 2210.7
帖子: 1180
在線: 354.8小時(shí)
蟲號: 900839

[交流] 【求助】高手請進(jìn)

調(diào)試正常，為什么運(yùn)行的時(shí)候出現(xiàn) segmentation fault，請大家?guī)兔纯矗@個(gè)程序哪里有問題，我調(diào)試了很長時(shí)間，總是這個(gè)錯(cuò)誤，請大家?guī)蛶兔�，不甚感激�?br /> sum = 0

   open(nin,file='TRAJ1',status='old',form &
            ='unformatted')
   open(nout,file='cswb.out',form='formatted')
   open(21,file='check')
      call readmdinit (mdtype,nin)
   loop1:  do istep=1, nstep
            call readmd  &
               (mdtype,nattot,nin,atomname,xxx,yyy,zzz,&
               cell,resname,resindex)
            open(npdb,file='newsystem.pdb',status='old')
            call readpdb (atomname,npdb,nattot,resname,resindex)
            close(npdb)
! loop over all waters
            nhb=0

         loop2: do i=1, nattot
               if ((atomname(i).eq.' O  ').and.(resname(i).eq.'WAT')) then
                  io = i
               loop3: do j=1,2
                              ih=io+j

loop4:  do iat=1,nattot
                        if((resname(iat).eq.'DC5'.or.resname(iat).eq.'DC ').and.&
                           (atomname(iat).eq.' N1 '.or.atomname(iat).eq.' N4 '.or.&
                           atomname(iat).eq.' N3 '.or.atomname(iat).eq.' O2 ')) then
                           call findhbond &
                                 (io,ih,iat,xxx,yyy,zzz,nattot,cell,ihbnum)
                           if (ihbnum==1) then
                              acc(j)=iat
                           else
                              exit loop3
                           endif
                        endif
                     enddo loop4
                  enddo loop3
                  if (acc(1)*acc(2)/=0) then
                     if (resindex(acc(1)).ne.resindex(acc(2))+1 &
                     .or. resindex(acc(1)).ne. resindex(acc(2))-1) then
                        nhb=nhb+1
                     endif
                  endif
            endif
      enddo loop2
      sum=sum+nhb
      enddo loop1
      avgnhb=sum/nstep
      write (nout,*) ihbnum

回復(fù)此樓

» 猜你喜歡

085801 總分275 本科新能源求調(diào)劑已經(jīng)有22人回復(fù)
298分 070300求調(diào)劑已經(jīng)有3人回復(fù)
材料與化工306分找調(diào)劑已經(jīng)有13人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有15人回復(fù)
260求調(diào)劑已經(jīng)有4人回復(fù)
299求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑求調(diào)劑已經(jīng)有7人回復(fù)
265求調(diào)劑已經(jīng)有15人回復(fù)
070300化學(xué)求調(diào)劑已經(jīng)有5人回復(fù)
296材料專碩求調(diào)劑已經(jīng)有13人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

重氮化高手請進(jìn) 已經(jīng)有12人回復(fù)
解譜高手請進(jìn)�。�！已經(jīng)有5人回復(fù)
氣質(zhì)高手請進(jìn)!!! 已經(jīng)有7人回復(fù)
PCR高手請進(jìn)，小女子有事相求，我剛來金幣不要嫌少啊，我只有這么多啊，全給你了已經(jīng)有20人回復(fù)
c語言高手請進(jìn)，請勿灌水，謝謝�。P(guān)于讀寫文件的問題）已經(jīng)有6人回復(fù)
水熱反應(yīng)高手請進(jìn) 已經(jīng)有17人回復(fù)
史上最難，高手請進(jìn)！已經(jīng)有9人回復(fù)
【求助】有機(jī)高手請進(jìn)-合成反應(yīng)解釋已經(jīng)有19人回復(fù)
【求助】合成高手請進(jìn) 已經(jīng)有9人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2010-11-23 17:28:39

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

astringent

銅蟲 (著名寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 2210.7
帖子: 1180
在線: 354.8小時(shí)
蟲號: 900839

問題已解決，不過感謝各位參與，沒人贈送10個(gè)金幣，呵呵。

贊一下

回復(fù)此樓

5樓2010-12-15 02:36:51

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

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

heyansheng12345

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 634.1
帖子: 71
在線: 4.9小時(shí)
蟲號: 543364

astringent(金幣+10): 2010-12-15 02:35:40

主程序和子程序之間的連接可能有問題，不很清楚，幫頂個(gè)

贊一下

回復(fù)此樓

2樓2010-12-14 13:44:36

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

zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲號: 664177

astringent(金幣+10): 2010-12-15 02:35:31

引用回帖:

Originally posted by astringent at 2010-11-23 17:28:39:
調(diào)試正常，為什么運(yùn)行的時(shí)候出現(xiàn) segmentation fault，請大家?guī)兔纯�，這個(gè)程序哪里有問題，我調(diào)試了很長時(shí)間，總是這個(gè)錯(cuò)誤，請大家?guī)蛶兔�，不甚感激�?br /> sum = 0

open(nin,file='TRAJ1',status='ol ...

這個(gè)翻譯過來，貌似叫做段錯(cuò)誤。

贊一下

回復(fù)此樓

3樓2010-12-14 13:47:15

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

cean

金蟲 (職業(yè)作家)

應(yīng)助: 23 (小學(xué)生)
金幣: 53.1
帖子: 3442
在線: 578.8小時(shí)
蟲號: 1147499

astringent(金幣+10): 2010-12-15 02:35:13

我前些天也遇到segmentation fault的問題,知道是肯定連接的不對.多變變連接的設(shè)置.

你的程序不全,沒法試.光讀不好調(diào).

贊一下

回復(fù)此樓

4樓2010-12-14 19:46:39

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

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

[考研] 296材料專碩求調(diào)劑 +12	202451007219 2026-04-02	13/650	2026-04-02 19:11 by dick_runner
[考研] 301求調(diào)劑 +4	細(xì)胞相關(guān)蛋白 2026-04-02	8/400	2026-04-02 16:30 by 271179835
[考研] 總分328生物與醫(yī)藥考數(shù)學(xué)求調(diào)劑 +3	aaadim 2026-04-02	3/150	2026-04-02 14:04 by 喬噠噠噠
[考研] 一志愿廈門大學(xué)材料工程專碩354找調(diào)劑！��！ +7	貝唄鋇鋇 2026-03-30	7/350	2026-04-02 12:19 by 1753564080
[考研] 0856初試324分求調(diào)劑 +6	想上學(xué)求調(diào) 2026-04-01	6/300	2026-04-02 11:42 by 星空星月
[考研] 289求調(diào)劑 +23	新時(shí)代材料 2026-03-27	26/1300	2026-04-02 10:29 by oooqiao
[考研] 310分求調(diào)劑 +4	成功上岸wang 2026-04-01	4/200	2026-04-01 20:35 by liu823948201
[考研] 求調(diào)劑0703 +5	周嘉堯 2026-03-31	8/400	2026-04-01 20:32 by ltltkkk
[考研] 0817化工學(xué)碩調(diào)劑 +11	努力上岸中！ 2026-03-31	11/550	2026-04-01 20:30 by 賴春艷
[考研] 0703一志愿南師大334求調(diào)劑 +4	seven7yu 2026-03-30	4/200	2026-04-01 16:10 by oooqiao
[考研] 材料0856 英一數(shù)二 323 求調(diào)劑 +9	袁sy 2026-04-01	9/450	2026-04-01 14:30 by wangjy2002
[考研] 環(huán)境工程 085701，267求調(diào)劑 +15	minht 2026-03-29	16/800	2026-04-01 10:13 by li_sujuan99
[考研] 339求調(diào)劑 +5	zjjkt 2026-03-31	5/250	2026-04-01 09:18 by JourneyLucky
[考研] 375求調(diào)劑 +7	雨夏整夜 2026-03-29	7/350	2026-03-31 18:52 by xhai2011
[考研] 生物考研337分求調(diào)劑 +4	cgxin 2026-03-30	6/300	2026-03-31 14:18 by 記事本2026
[考研] 085701求調(diào)劑初試286分 +5	secret0328 2026-03-28	5/250	2026-03-30 12:54 by fangnagu
[考研] 348求調(diào)劑 +6	小懶蟲不懶了 2026-03-28	6/300	2026-03-30 10:29 by Evan_Liu
[考研] 086000生物與醫(yī)藥調(diào)劑 +5	Feisty。 2026-03-28	9/450	2026-03-29 12:02 by longlotian
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3	llss0711 2026-03-28	6/300	2026-03-29 10:26 by llss0711
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒