24小时热门版块排行榜

返回列表

astringent

铜虫 (著名写手)

应助: 16 (小学生)
金币: 2210.7
帖子: 1180
在线: 354.8小时
虫号: 900839

[交流] 【求助】高手请进

调试正常，为什么运行的时候出现 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

回复此楼

» 猜你喜欢

不自信的我已经有5人回复
磺酰氟产物，毕不了业了！已经有4人回复
论文终于录用啦！满足毕业条件了已经有16人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

重氮化高手请进已经有12人回复
解谱高手请进！！！已经有5人回复
气质高手请进!!! 已经有7人回复
PCR高手请进，小女子有事相求，我刚来金币不要嫌少啊，我只有这么多啊，全给你了已经有20人回复
c语言高手请进，请勿灌水，谢谢！（关于读写文件的问题）已经有6人回复
水热反应高手请进已经有17人回复
史上最难，高手请进！已经有9人回复
【求助】有机高手请进-合成反应解释已经有19人回复
【求助】合成高手请进已经有9人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

heyansheng12345

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 634.1
帖子: 71
在线: 4.9小时
虫号: 543364

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

主程序和子程序之间的连接可能有问题，不很清楚，帮顶个

赞一下

回复此楼

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
帖子: 2936
在线: 1329.4小时
虫号: 664177

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

引用回帖:

Originally posted by astringent at 2010-11-23 17:28:39:
调试正常，为什么运行的时候出现 segmentation fault，请大家帮忙看看，这个程序哪里有问题，我调试了很长时间，总是这个错误，请大家帮帮忙，不甚感激。
sum = 0

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

这个翻译过来，貌似叫做段错误。

赞一下

回复此楼

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖