±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 2330  |  »Ø¸´: 7

zyj8119

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

[½»Á÷] ¡¾×ªÌû¡¿¾¶Ïò·Ö²¼º¯Êý³ÌÐòÓë¼òµ¥ËµÃ÷ ÒÑÓÐ4È˲ÎÓë

¾¶Ïò·Ö²¼º¯Êýg(r)´ú±íÁËÇò¿ÇÄ򵀮½¾ùÊýÃܶÈ

ΪÀëÖÐÐÄ·Ö×Ó¾àÀëΪr£¬Ìå»ýΪ µÄÇò¿ÇÄÚµÄ˲ʱ·Ö×ÓÊý¡£
¾ßÌå²Î¼ûÀîÈçÉú£¬¡¶Æ½ºâºÍ·Çƽºâͳ¼ÆÁ¦Ñ§¡·¿ÆÑ§³ö°æÉ磺1995
CODE:
SUBROUTINE GR(NSWITCH)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      PARAMETER(NM=40000,PI=3.141592653589793D0,NHIS=100)
      COMMON/LCS/X0(3,-2:2*NM),X(3,-2:2*NM,5),XIN(3,-2:2*NM),
     $XX0(3,-2:2*NM),XX(3,-2:2*NM,5),XXIN(3,-2:2*NM)
      COMMON/MOLEC/LPBC(3),MOLSP,MOLSA,NBX,NBY,NBZ,NPLA,LPBCSM,NC,NN,MC
      COMMON/WALLS/HI(3,3),G(3,3),DH,AREA,VOLUME,SCM(3)
      COMMON/PBCS/HALF,PBCX,PBCY,PBCZ
        COMMON/GR_VAR/ NGR
        DIMENSION H(3,3),GG(0:NHIS),R(0:NHIS)
      EQUIVALENCE(X0(1,-2),H(1,1))
C   *****************************************************************
C      ÈçºÎÈ·¶¨·Ö×ÓÊýÃܶȣºDEN_IDEAL
C      È¡·Ö×Ó×ÜÊý×÷ΪģÄâºÐÖеÄÊýÃܶȣ¬¿É±£Ö¤²ÉÑù·Ö×Ó×ÜÊý=×Ü·Ö×ÓÊý?
C====================================================================
C        ¡¡N1=MOLSP+1
C      N2=MOLSP+NC

      DEN_IDEAL=MOLSP  

        G11=G(1,1)
      G22=G(2,2)
      G33=G(3,3)
      G12D=G(1,2)+G(2,1)
      G13D=G(1,3)+G(3,1)
      G23D=G(2,3)+G(3,2)


      IF(NSWITCH.EQ.0)THEN
          NGR=0
          DELR=HALF/NHIS
          DO I=1,NHIS
           GG(I)=0.D0
           R(I)=0.D0
          ENDDO

      ELSE IF(NSWITCH.EQ.1)THEN
         NGR=NGR+1
       DO I=1,MOLSP-1
         DO J=I+1,MOLSP
C====================================================================
C     USE PBC IN X DIRECTION:  SUITABLE FOR PBCX=1
C                              NOT GREAT PROBLEM FOR PBCX=0
C                              (THIS TIME USUALLY |DELTA X| < HALF)
C====================================================================
          XIJ=X0(1,I)-X0(1,J)
        IF(XIJ.GT.+HALF)XIJ=XIJ-PBCX
        IF(XIJ.LT.-HALF)XIJ=XIJ+PBCX
        YIJ=X0(2,I)-X0(2,J)
        IF(YIJ.GT.+HALF)YIJ=YIJ-PBCY
        IF(YIJ.LT.-HALF)YIJ=YIJ+PBCY
        ZIJ=X0(3,I)-X0(3,J)
        IF(ZIJ.GT.+HALF)ZIJ=ZIJ-PBCZ
        IF(ZIJ.LT.-HALF)ZIJ=ZIJ+PBCZ
        RSQ=XIJ*(G11*XIJ+G12D*YIJ+G13D*ZIJ)+
     $      YIJ*(G22*YIJ+G23D*ZIJ)+G33*ZIJ*ZIJ
          RRR=SQRT(RSQ)
          RRR=RRR/H(1,1)
C====================================================================
C      ÒÔÉÏÓÃÊý×éGºÍHµÄ½á¹ûÓëÏÂͬ
C      RRR=SQRT(XIJ**2+YIJ**2+ZIJ**2)
C      G11=H(1,1)**2
C====================================================================
          IF(RRR.LT.HALF)THEN
           IG=INT(RRR/DELR)
           GG(IG)=GG(IG)+2
          ENDIF
       ENDDO
         ENDDO

      ELSE IF(NSWITCH.EQ.2)THEN
        DO I=1,NHIS
           R(I)=DELR*(I+0.5D0)
        ENDDO
        DO I=1,NHIS
           VB=(4.D0/3.D0)*PI*(((I+1)**3-I**3)*(DELR**3))
           GNID=VB*DEN_IDEAL
           GG(I)=GG(I)/(NGR*MOLSP*GNID)
        ENDDO
        OPEN(UNIT=31,FILE="GR.DAT")
        DO I=1,NHIS
           WRITE(31,*)R(I),GG(I)
        ENDDO
        CLOSE(31)

        ENDIF

        RETURN
        END

»Ø¸´´ËÂ¥
ºÃºÃѧϰ£¬ÌìÌìÏòÉÏ¡£
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zyj8119

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

ÒýÓûØÌû:
Originally posted by zyj8119 at 2010-09-13 19:42:08:
¾¶Ïò·Ö²¼º¯Êýg(r)´ú±íÁËÇò¿ÇÄ򵀮½¾ùÊýÃܶÈ

ΪÀëÖÐÐÄ·Ö×Ó¾àÀëΪr£¬Ìå»ýΪ µÄÇò¿ÇÄÚµÄ˲ʱ·Ö×ÓÊý¡£
¾ßÌå²Î¼ûÀîÈçÉú£¬¡¶Æ½ºâºÍ·Çƽºâͳ¼ÆÁ¦Ñ§¡·¿Æ ...

MCÖÐÐèÒªÇóRDFÂð£¿
ºÃºÃѧϰ£¬ÌìÌìÏòÉÏ¡£
2Â¥2010-09-13 19:51:49
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

yahoohoo

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

¡ï ¡ï ¡ï ¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
lei0736(½ð±Ò+3):лл 2010-09-13 20:50:00
ÕâÑùµÄ´úÂë¿´×Ų»¹»Ã÷ÁË¡£¡£¡£¡£¡£¡£

α´úÂ룺

for (int i = 0; i < TOTN - 1; ++i)
  for (int j = i + 1; j < TOTN; ++j) {
    double dij = sqrt( pow(Pos[0]-Pos[j][0], 2) + pow(Pos[1]-Pos[j][1], 2) + pow(Pos[2]-Pos[j][2], 2));
    int kbin = func(dij); // dijËù¶ÔÓ¦µÄbinµÄÐòºÅ
    g(kbin) += 2;
  }
  // normalize
  for (int k = 0; k < NBIN; ++k)
    g(k) /= 4.0 * PI * r(k) * r(k) * dr * RHO; // r ΪµÚk¸öbinËù¶ÔÓ¦µÄ¾àÀëÖµ
3Â¥2010-09-13 20:39:36
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zyj8119

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

ÒýÓûØÌû:
Originally posted by yahoohoo at 2010-09-13 20:39:36:
ÕâÑùµÄ´úÂë¿´×Ų»¹»Ã÷ÁË¡£¡£¡£¡£¡£¡£

α´úÂ룺

for (int i = 0; i < TOTN - 1; ++i)
  for (int j = i + 1; j < TOTN; ++j) {
    double dij = sqrt( pow(Pos[0]-Pos[j][0], 2) + pow(Pos[1 ...

Õâ¸öÊÇCÓïÑԵİɣ¿
ºÃºÃѧϰ£¬ÌìÌìÏòÉÏ¡£
4Â¥2010-09-13 21:21:47
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ghcacj

ÈÙÓþ°æÖ÷ (ÖøÃûдÊÖ)

°¢³¬

ÓÅÐã°æÖ÷

¡ï ¡ï ¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
zh1987hs(½ð±Ò+2):лл 2010-09-14 08:40:12
ÒýÓûØÌû:
Originally posted by zyj8119 at 2010-09-13 21:21:47:

Õâ¸öÊÇCÓïÑԵİɣ¿

ÊÇαÂë°É
5Â¥2010-09-13 21:38:22
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zyj8119

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

ÒýÓûØÌû:
Originally posted by ghcacj at 2010-09-13 21:38:22:

ÊÇαÂë°É

µÄÈ·ÊÇαÂë¡£
ºÃºÃѧϰ£¬ÌìÌìÏòÉÏ¡£
6Â¥2010-09-14 07:56:13
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wangqingbo

Í­³æ (СÓÐÃûÆø)

¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
¿ÉÒÔÓÃÀ´Ö±½Ó¼ÆËãÂð ¶ÔÓÚÊäÈëµÄ²ÎÊýÊÇ Ê²Ã´ÒªÇó
7Â¥2010-09-16 18:41:17
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

yahoohoo

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

¡ï
ghcacj(½ð±Ò+1):лл 2010-09-17 10:01:27
ÒýÓûØÌû:
Originally posted by wangqingbo at 2010-09-16 18:41:17:
¿ÉÒÔÓÃÀ´Ö±½Ó¼ÆËãÂð ¶ÔÓÚÊäÈëµÄ²ÎÊýÊÇ Ê²Ã´ÒªÇó

±ðÈ˵ĴúÂëÄùýÀ´¾ÍÓã¬ÄãÔõôÅжϽá¹û¶Ô»¹ÊDz»¶ÔÄØ£¿
8Â¥2010-09-16 21:53:30
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ zyj8119 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 0855Çóµ÷¼Á²ÄÁÏ +11 ºìÌÒׯׯ 2026-04-04 12/600 2026-04-06 10:26 by À¶ÔÆË¼Óê
[¿¼ÑÐ] ¼ÆËã»ú11408£¬286·ÖÇóµ÷¼Á +7 ľ×ÓÄî•„ 2026-04-05 7/350 2026-04-05 19:02 by chy09050039
[¿¼ÑÐ] 282Çóµ÷¼Á +7 aaa³µÁ¾ 2026-04-02 11/550 2026-04-05 17:24 by yulian1987
[¿¼ÑÐ] 085600µ÷¼Á +9 ¶«ÕÕÕÕÕÕ 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[¿¼ÑÐ] 298·Ö 070300Çóµ÷¼Á +15 zwen03 2026-04-02 15/750 2026-04-05 12:52 by Hdyxbekcb
[¿¼ÑÐ] Ò»Ö¾Ô¸Î÷±±Å©ÁÖÐóÄÁר˶336·ÖÇóµ÷¼Á +3 5ourr 2026-04-03 3/150 2026-04-05 10:40 by JOKER0401
[¿¼ÑÐ] 11408,335·Ö£¬±¾¿Æ211£¬Çóµ÷¼Á£¬¿Éתרҵ +5 öùÀæ´óöùÓã 2026-04-03 5/250 2026-04-04 22:49 by chongya
[¿¼ÑÐ] ¿¼ÑÐÇóµ÷¼Á +3 ľÐÄÏë¼ÌÐøÉîÔì 2026-04-03 3/150 2026-04-03 21:56 by à£à£à£0119
[¿¼ÑÐ] 281Çóµ÷¼Á +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[¿¼ÑÐ] 282Çóµ÷¼Á +5 ºôÎü¶¼ÊǼõ·Ê 2026-03-31 5/250 2026-04-03 12:03 by 1753564080
[¿¼ÑÐ] ½¨»·£¬ÄÜÔ´£¬ÍÁľÀÏʦ·¹ý¿´Ò»¿´£¡£¡£¡ +5 ºÙºÙuu 2026-04-01 5/250 2026-04-03 11:47 by znian
[¿¼ÑÐ] 320Çóµ÷¼Á +3 ũҵ¹¤³ÌÓëÐÅÏ¢¼ 2026-04-03 3/150 2026-04-03 11:40 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] Ò»Ö¾Ô¸´ó¹¤Ñ§Ë¶£¬Çóµ÷¼Á +4 yub0811 2026-04-02 4/200 2026-04-02 21:36 by °ÙÁéͯ888
[¿¼ÑÐ] 318Çóµ÷¼Á +3 óÆÐÐÖÂÔ¶. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[¿¼ÑÐ] 282Çóµ÷¼Á +13 ºôÎü¶¼ÊǼõ·Ê 2026-04-01 13/650 2026-04-02 14:10 by baoball
[¿¼ÑÐ] ÍÁľ304Çóµ÷¼Á +6 ÍÃͻͻͻ£¬ 2026-03-31 7/350 2026-04-02 09:06 by coolminer
[¿¼ÑÐ] 085600 Ò»Ö¾Ô¸9 ×Ü·Ö351 Çóµ÷¼ÁѧУ +7 czhcz 2026-03-31 9/450 2026-04-01 19:24 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] Ò»Ö¾Ô¸ ÄϾ©º½¿Õº½Ìì´óѧ £¬080500²ÄÁÏ¿ÆÑ§Ó빤³Ìѧ˶ +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
[¿¼ÑÐ] ±¾2Ò»Ö¾Ô¸C9-333·Ö£¬²ÄÁÏ¿ÆÑ§Ó빤³Ì£¬Çóµ÷¼Á +9 ÉýÉý²»½µ 2026-03-31 9/450 2026-03-31 18:01 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] ¼ª´óÉúÎïѧ326·ÖÇóµ÷¼Á +3 sunnyupup 2026-03-31 3/150 2026-03-31 09:28 by longlotian
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û