| 5 | 1/1 | 返回列表 |
| 查看: 1114 | 回復: 5 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[交流]
【求助】高手請進
|
|||
|
調試正常,為什么運行的時候出現(xiàn) segmentation fault,請大家?guī)兔纯,這個程序哪里有問題,我調試了很長時間,總是這個錯誤,請大家?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 |
» 搶金幣啦!回帖就可以得到:
+1/286
+1/82
+1/78
+1/50
+1/38
+2/36
+1/35
+1/33
+1/14
+1/12
+1/9
+1/9
+2/8
+1/6
+1/4
+1/4
+1/2
+1/2
+1/1
+1/1
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[碩博家園] 招收生物學/細胞生物學調劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|---|---|---|---|---|
|
[碩博家園] 求調劑 有機化學考研356分 +3 | Nadiums 2026-03-25 | 4/200 |
|
|
[考研] 305求調劑 +7 | RuiFairyrui 2026-03-28 | 7/350 |
|
|
[考研] 085701求調劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
[考研] 085600 286分 材料求調劑 +7 | 麻辣魷魚 2026-03-27 | 8/400 |
|
|
[考研] 339求調劑,想調回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 291求調劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 考研調劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 07化學280分求調劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 085600,材料與化工321分,求調劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 求調劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 289求調劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 086000生物與醫(yī)藥292求調劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
[考研] 環(huán)境專碩324分求調劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 一志愿天津大學339材料與化工求調劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
[考研] 調劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 070300,一志愿北航320求調劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|