±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 597  |  »Ø¸´: 3
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

astringent

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

[½»Á÷] ¡¾ÇóÖú¡¿ÎªÊ²Ã´°ÑÑ­»·±äÁ¿µÄÖÕÖµ¸Ä±ä£¬±àÒë¾Í³ö´íÁË£¿ ÒÑÓÐ3È˲ÎÓë

Ô­À´µÄÑ­»·±äÁ¿µÄÖÕÖµÊÇ2000£¬¿ÉÊÇÐèÒª¶à×öÒ»¶Î£¬ÐèÒªÔö¼Óµ½5000£¬¿ÉÊDZàÒëʱ¾Í³öÏÖÏÂÃæµÄ´íÎóÁË£º
ld: rel32 out of range in start in start from /usr/lib/crt1.10.5.o
collect2: ld returned 1 exit status
ÊÇʲôԭÒòÄØ£¿
ÎҵijÌÐòÖÐÓëÑ­»·±äÁ¿ÓйصIJ¿·ÖÊÇ£º
integer,parapemer::nstep=5000  (nstepΪѭ»·±äÁ¿µÄÖÕÖµ)
real*8 xcm(nattot,nstep),ycm(nattot,nstep),zcm(nattot,nstep)
real*8 msd(nstep)
ntime=0
        msd=0.0d0
        do i=1,nstep-1,10
           do j=i+1,nstep
              m=j-i
              ntime(m)=ntime(m)+1
              do k=1,nwat
                 xdiff = xcm(p(k),j) - xcm(p(k),i)
                 ydiff = ycm(p(k),j) - ycm(p(k),i)
                 zdiff = zcm(p(k),j) - zcm(p(k),i)
                 msd(m) =msd(m)+xdiff*xdiff+ydiff*ydiff+zdiff*zdiff
               enddo
            enddo
          enddo
          do i=1,nstep-1
             msd(i)=msd(i)/(dble(nwat)*dble(ntime(i)))
             write(nout,*) dble(i)*dt,msd(i)
          enddo
»Ø¸´´ËÂ¥
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

astringent

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

ÒýÓûØÌû:
Originally posted by snoopyzhao at 2011-01-29 10:14:20:
ÄãÕâÀdble, ntime µÄ dimension ÊǶàÉÙ£¿£¿

¿ÉÄÜÊÇÕ⼸¸öxcm(nattot,nstep),ycm(nattot,nstep),zcm(nattot,nstep)Ì«´ó£¬ÆäÖÐnattot=28000£¬nstep=5000
ÇëÎʸÃÈçºÎ½â¾ö°¡£¿
3Â¥2011-01-29 16:07:04
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 4 ¸ö»Ø´ð

snoopyzhao

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

¡ï ¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
ÓàÔó³É(½ð±Ò+1): лл²ÎÓëÓ¦Öú£¡ 2011-02-20 15:50:59
ÄãÕâÀdble, ntime µÄ dimension ÊǶàÉÙ£¿£¿
2Â¥2011-01-29 10:14:20
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

scatcj

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

¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
¸ù¾ÝÊý×é´óС£¬¶¯Ì¬·ÖÅä
Troubleisafriend!
4Â¥2011-02-22 23:50:47
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤£¨0856£©304ÇóBÇøµ÷¼Á +3 Çñgl 2026-03-30 7/350 2026-03-30 20:25 by 12022246319
[¿¼ÑÐ] 285Çóµ÷¼Á +4 AZMK 2026-03-30 7/350 2026-03-30 20:24 by laoshidan
[¿¼ÑÐ] ¸÷λÀÏʦºÃ£¬ÎÒµÄһ־ԸΪ±±¾©¿Æ¼¼´óѧ085601²ÄÁÏר˶ +10 Koxui 2026-03-28 10/500 2026-03-30 19:33 by Ô´_2020
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄÏ¿ª´óѧ0710ÉúÎïѧ359Çóµ÷¼Á +5 ÍÃÍÃÍÃ111223314 2026-03-29 7/350 2026-03-30 18:29 by ÍÃÍÃÍÃ111223314
[¿¼ÑÐ] 328Çóµ÷¼Á +8 àŵεĻù±¾¶¼ 2026-03-27 8/400 2026-03-30 17:20 by Wang200018
[¿¼ÑÐ] ¿¼Ñе÷¼Á +5 Sanmu-124 2026-03-26 5/250 2026-03-30 13:36 by chemdavid
[»ù½ðÉêÇë] ÃæÉÏ5BÄÜÉÏ»áÂ𣿠+7 redcom 2026-03-29 7/350 2026-03-30 12:44 by °®¸ãÑо¿µÄСº¢
[¿¼ÑÐ] »·¾³¿ÆÑ§Ó빤³Ì334·ÖÇóµ÷¼Á +6 ÍõÒ»Ò»ÒÀÒÀ 2026-03-30 8/400 2026-03-30 11:52 by yjolah
[¿¼ÑÐ] 296Çóµ÷¼Á +10 ±Ë°¶t 2026-03-29 10/500 2026-03-30 10:50 by ̽123
[¿¼ÑÐ] 085600£¬²ÄÁÏÓ뻯¹¤321·ÖÇóµ÷¼Á +10 ´ó²öС×Ó 2026-03-28 10/500 2026-03-29 23:35 by ·ÉÐÐÈÕ¼ÇÎ÷
[¿¼ÑÐ] 290Çóµ÷¼Á +3 dfffsar 2026-03-29 3/150 2026-03-29 22:38 by ëë밢Ī2
[¿¼²©] 26É격×Ô¼ö +6 whh869393 2026-03-24 6/300 2026-03-29 21:05 by nxgogo
[¿¼ÑÐ] 311Çóµ÷¼Á +6 ¶¬Ê®Èý 2026-03-24 6/300 2026-03-29 20:45 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] Ò»Ö¾Ô¸Î人Àí¹¤£¬×Ü·Ö321£¬Ó¢Ò»Êý¶þ£¬ÇóÀÏʦÊÕÁô¡£ +11 nnnnnnn5 2026-03-25 11/550 2026-03-29 20:42 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] ¿¼Ñе÷¼Á +7 СÀ¯ÐÂ±Ê 2026-03-29 7/350 2026-03-29 19:00 by ѧԱ8dgXkO
[¿¼ÑÐ] 343Çóµ÷¼Á +6 °®î¿°í 2026-03-29 6/300 2026-03-29 12:00 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] Çóµ÷¼Á +7 ÕùÈ¡¾Åµã˯ 2026-03-28 8/400 2026-03-28 21:07 by ÕùÈ¡¾Åµã˯
[¿¼ÑÐ] 347Çóµ÷¼Á +3 ɽ¶¥¼û¦Á 2026-03-25 3/150 2026-03-28 14:13 by ÌÆãå¶ù
[¿¼ÑÐ] 292Çóµ÷¼Á +4 ÇóÇóÁËÊÕÏÂÎÒ°É£ 2026-03-26 4/200 2026-03-27 10:37 by zhshch
[¿¼ÑÐ] 292Çóµ÷¼Á +4 ¶ì¶ì¶ì¶î¶î¶î¶î¶ 2026-03-24 4/200 2026-03-24 16:41 by peike
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û