| 6 | 1/1 | 返回列表 |
| 查看: 1110 | 回復(fù): 5 | |||
astringent銅蟲(chóng) (著名寫(xiě)手)
|
[交流]
【求助】高手請(qǐng)進(jìn)
|
|
調(diào)試正常,為什么運(yùn)行的時(shí)候出現(xiàn) segmentation fault,請(qǐng)大家?guī)兔纯,這個(gè)程序哪里有問(wèn)題,我調(diào)試了很長(zhǎng)時(shí)間,總是這個(gè)錯(cuò)誤,請(qǐng)大家?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/85
+1/75
+1/74
+2/44
+2/40
+2/38
+1/38
+1/33
+1/31
+1/19
+1/14
+1/10
+1/7
+1/5
+1/5
+1/4
+1/3
+1/1
+1/1
+1/1
銀蟲(chóng) (小有名氣)
金蟲(chóng) (職業(yè)作家)
銅蟲(chóng) (著名寫(xiě)手)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +4 | 阿姨復(fù)古不過(guò) 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +5 | 哇呼哼呼哼 2026-03-27 | 7/350 |
|
|
[考研] 070300求調(diào)劑306分 +3 | 26要上岸 2026-03-27 | 3/150 |
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
[考研] 303求調(diào)劑 +7 | 安憶靈 2026-03-22 | 8/400 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +3 | IceGuo 2026-03-26 | 4/200 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 5/250 |
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書(shū) 2026-03-23 | 4/200 |
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣(mài)魚(yú) 2026-03-26 | 3/150 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|