| 查看: 1288 | 回复: 5 | |||
[交流]
【求助】高手请进
|
|
调试正常，为什么运行的时候出现 segmentation fault，请大家帮忙看看，这个程序哪里有问题，我调试了很长时间，总是这个错误，请大家帮帮忙，不甚感激。 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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
国自然评审意见
已经有17人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 重氮化高手请进
已经有12人回复
- 解谱高手请进！！！
已经有5人回复
- 气质高手请进!!!
已经有7人回复
- PCR高手请进，小女子有事相求，我刚来金币不要嫌少啊，我只有这么多啊，全给你了
已经有20人回复
- c语言高手请进，请勿灌水，谢谢！（关于读写文件的问题）
已经有6人回复
- 水热反应高手请进
已经有17人回复
- 史上最难，高手请进！
已经有9人回复
- 【求助】有机高手请进-合成反应解释
已经有19人回复
- 【求助】合成高手请进
已经有9人回复
2楼2010-12-14 13:44:36
3楼2010-12-14 13:47:15
4楼2010-12-14 19:46:39
5楼2010-12-15 02:36:51
6楼2015-04-27 15:22:45
投票:
10