| 6 | 1/1 | 返回列表 |
| 查看: 1111 | 回復: 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/195
+1/46
+2/44
+1/43
+1/40
+1/39
+1/38
+1/37
+1/36
+1/34
+1/23
+1/16
+1/9
+1/5
+1/5
+1/5
+1/4
+1/3
+1/1
+1/1
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程求調劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 083000學碩274求調劑 +5 | Li李魚 2026-03-26 | 5/250 |
|
|
[考研] 0703化學求調劑,各位老師看看我!! +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 286求調劑 +8 | PolarBear11 2026-03-26 | 8/400 |
|
|
[考研] 317求調劑 +5 | 十閑wx 2026-03-24 | 5/250 |
|
|
[考研] 22408 359分調劑 +3 | Qshers 2026-03-27 | 3/150 |
|
|
[考研] 336材料求調劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 333求調劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研] 321求調劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 081200-11408-276學碩求調劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 求調劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
[考研] 化學工程085602 305分求調劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
[考研] 總分322求生物學/生化與分子/生物信息學相關調劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研] 生物學學碩,一志愿湖南大學,初試成績338 +4 | YYYYYNNNNN 2026-03-26 | 4/200 |
|
|
[考研] 打過很多競賽,085406控制工程300分,求調劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 309求調劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 化學調劑 +6 | yzysaa 2026-03-21 | 6/300 |
|
|
[考研] 調劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 335求調劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
[考研] 336化工調劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|