24СʱÈÈÃŰæ¿éÅÅÐаñ    

²é¿´: 1286  |  »Ø¸´: 5
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

astringent

Í­³æ (ÖøÃûдÊÖ)


[½»Á÷] ¡¾ÇóÖú¡¿¸ßÊÖÇë½ø

µ÷ÊÔÕý³££¬ÎªÊ²Ã´ÔËÐеÄʱºò³öÏÖ 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
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

astringent

Í­³æ (ÖøÃûдÊÖ)


ÎÊÌâÒѽâ¾ö£¬²»¹ý¸Ðл¸÷λ²ÎÓ룬ûÈËÔùËÍ10¸ö½ð±Ò£¬ºÇºÇ¡£
5Â¥2010-12-15 02:36:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 6 ¸ö»Ø´ð

heyansheng12345

Òø³æ (СÓÐÃûÆø)


astringent(½ð±Ò+10): 2010-12-15 02:35:40
Ö÷³ÌÐòºÍ×Ó³ÌÐòÖ®¼äµÄÁ¬½Ó¿ÉÄÜÓÐÎÊÌ⣬²»ºÜÇå³þ£¬°ï¶¥¸ö
2Â¥2010-12-14 13:44:36
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zyj8119

ľ³æ (ÖøÃûдÊÖ)


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µÄ»ØÌû

cean

½ð³æ (Ö°Òµ×÷¼Ò)


astringent(½ð±Ò+10): 2010-12-15 02:35:13
ÎÒǰЩÌìÒ²Óöµ½segmentation faultµÄÎÊÌâ,ÖªµÀÊǿ϶¨Á¬½ÓµÄ²»¶Ô.¶à±ä±äÁ¬½ÓµÄÉèÖÃ.

ÄãµÄ³ÌÐò²»È«,û·¨ÊÔ.¹â¶Á²»ºÃµ÷.
4Â¥2010-12-14 19:46:39
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[»ù½ðÉêÇë] ¹ú×ÔÈ»ÃæÉϸ´ÅÌ~»¶Ó­ÌÖÂÛ (½ð±Ò+15) +13 ÇçÌì¼ÓÓÍ 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[»ù½ðÉêÇë] ÓÐûÓÐÈÔûÊÕµ½ÐÅÏ¢µÄ +3 µÂÉÐÖÐÐÐ 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[»ù½ðÉêÇë] »ù½ð²»ÖУ¬¹²Ãã +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by »ðÐdz¬ÈËxi
[»ù½ðÉêÇë] 2026ÄêÒ¶ÆóËï»ù½ð +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[»ù½ðÉêÇë] »ù½ðϵͳʲôÄÚÈÝҲûÓÐ 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[»ù½ðÉêÇë] Ôõô²é°¡ +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[»ù½ðÉêÇë] ÄÄλ¸ßÈËÖÐÁË£¬°Ñ²éѯµ½µÄ½ØÍ¼Ìù³öÀ´ÈÃÎÒ¿´¿´£¬ÈÃÎÒ³¤³¤¼ûʶ +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[»ù½ðÉêÇë] ÃæÉϺÏ×÷µ¥Î»¸ÇÕ +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[»ù½ðÉêÇë] ÉêÇëɾ³ý±¾Ìû +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by Äþ¾²ÖÂÔ¶sy
[»ù½ðÉêÇë] »ù½ðδÖУ¬ÕâÖÖ´ð¸´ÊÇÄ£°åÂ𣿠+5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[»ù½ðÉêÇë] Ôõô¿´Çà»ùÖÐÁËûÓа¡ +5 Ò¶¾Å΢ 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[»ù½ðÉêÇë] Ϊʲô ¹ú¼Ê(µØÇø)ºÏ×÷Óë½»Á÷ÏîÄ¿ ûÓзŰñ£¿ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by ±±¾©À³Òð±à¼­
[»ù½ðÉêÇë] ÎÒ²»Àí½â£¡ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[»ù½ðÉêÇë] ³öÀ´ÁË +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by ±¦±´³æ×Ó
[»ù½ðÉêÇë] Ϊʲô¹ú×ÔÈ»²»ÄÜÖ±½Ó¹«²¼ +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[»ù½ðÉêÇë] ¹ú¼ÊºÏ×÷¿É²éÁË£¬ÖÐÁËÃæÉÏ (EPI+1)(½ð±Ò+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[»ù½ðÉêÇë] ¹úºÏ¿É²éÁË +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[Ö°³¡ÈËÉú] ѧÉúײ¼û¸¨µ¼Ô±ËÍÍâÂô£¬µÚ¶þÌìÈ«°à¶¼³ÁĬÁË +3 ¾¨ÓãÈÚ½ð_Õã½­_É 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[»ù½ðÉêÇë] Å£À´£¡Ã×À´£¡ÃæÀ´£¡ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[½Ìʦ֮¼Ò] Ìø²ÛºóÔÚÑÐÏîÄ¿Ôõô°ì£¿ +5 ¼òµ¥»¯xn 2026-08-22 10/500 2026-08-23 12:38 by ¼òµ¥»¯xn
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û