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

568286346

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

[ÇóÖú] Êý×鳬³öÎÊÌâ

¸÷λ³æÓÒ£¬±¾ÈËÔÚÓÃfortran±à³ÌÐò´¦ÀíÎÊÌâʱÓöµ½ÁËÊý×鳬³öµÄÎÊÌ⣨³¬³öÁËfortranĬÈÏÔÊÐíµÄ×î´óÖµ£©£¬ÇëÎÊÓÐʲô·½·¨¿ÉÒÔ½â¾öÕâ¸öÎÊÌâÂ𣿠ÒòΪÎÒÔÚ´¦ÀíÊý¾ÝʱÐèÒª¿ªºÜ´óµÄÊý×飬ËùÒÔ»¹Íû´ó¼ÒÌṩµãÒâ¼û°¡¡£
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

¼á³Ö²ÅÄÜʵÏÖ×Ô¼ºµÄÀíÏ룡
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
»ØÌûÖ§³Ö ( ÏÔʾ֧³Ö¶È×î¸ßµÄǰ 50 Ãû )

pippi6

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

¹¤³ÌºÍ¿ÆÑ§ÊýÖµ¼ÆËã×Éѯ

²éÁËÒ»ÏÂlzµÄÇóÖú¼Ç¼£¬http://muchong.com/bbs/space.php?uid=1785423&view=thread£¬·¢ÁË9¸öÇóÖúÌû£¬Ö»ÓÐÒ»¸ö¶ÒÏÖÁË¡£ÎÒ¿´ÁË¿´£¬´ó²¿·Ö½â´ð¶¼Ó¦¸ÃÊÇÓаïÖúµÄ¡£lzƽ³£×÷Ϊ¿É¼ûÒ»°ß¡£

» ±¾Ìû¸½¼þ×ÊÔ´Áбí

  • »¶Ó­¼à¶½ºÍ·´À¡£ºÐ¡Ä¾³æ½öÌṩ½»Á÷ƽ̨£¬²»¶Ô¸ÃÄÚÈݸºÔð¡£
    ±¾ÄÚÈÝÓÉÓû§×ÔÖ÷·¢²¼£¬Èç¹ûÆäÄÚÈÝÉæ¼°µ½ÖªÊ¶²úȨÎÊÌ⣬ÆäÔðÈÎÔÚÓÚÓû§±¾ÈË£¬Èç¶Ô°æÈ¨ÓÐÒìÒ飬ÇëÁªÏµÓÊÏ䣺xiaomuchong@tal.com
  • ¸½¼þ 1 : ¿ìÕÕ1.jpg
  • 2013-07-01 16:31:25, 110.01 K
12Â¥2013-07-01 16:36:08
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
ÆÕͨ»ØÌû

pippi6

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

¹¤³ÌºÍ¿ÆÑ§ÊýÖµ¼ÆËã×Éѯ

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
ÄãÒª¶à´óµÄÊý×é°¡£¿Ëµ¸öÊýÁ¿¼¶¡£
2Â¥2013-06-30 18:06:31
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

568286346

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

ÒýÓûØÌû:
2Â¥: Originally posted by pippi6 at 2013-06-30 18:06:31
ÄãÒª¶à´óµÄÊý×é°¡£¿Ëµ¸öÊýÁ¿¼¶¡£

±ÈÈç˵n_conf=200000,n_mol=196,Êý×éÉèÊÇx£¨n_conf£¬n_mol£©
¼á³Ö²ÅÄÜʵÏÖ×Ô¼ºµÄÀíÏ룡
3Â¥2013-06-30 18:25:11
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

pippi6

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

¹¤³ÌºÍ¿ÆÑ§ÊýÖµ¼ÆËã×Éѯ

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

ÒýÓûØÌû:
3Â¥: Originally posted by 568286346 at 2013-06-30 18:25:11
±ÈÈç˵n_conf=200000,n_mol=196,Êý×éÉèÊÇx£¨n_conf£¬n_mol£©...

ÎÒÊÔÁË£¬Ã»ÎÊÌâ°¡¡£
  integer, parameter :: n_conf=200000,n_mol=196
  real*4 :: ddx(n_conf,n_mol)
4Â¥2013-06-30 18:39:21
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

568286346

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

ÒýÓûØÌû:
4Â¥: Originally posted by pippi6 at 2013-06-30 18:39:21
ÎÒÊÔÁË£¬Ã»ÎÊÌâ°¡¡£
  integer, parameter :: n_conf=200000,n_mol=196
  real*4 :: ddx(n_conf,n_mol)...

program main
        implicit none
        integer   nn_conf, nn_mol, n_conf, n_mol1,n_mol2
        parameter (nn_conf=200000, nn_mol=196)
        real*8  len,b
        integer a, i, j, m, n, t,num2(nn_conf,nn_mol)
        parameter(m=4,n=2)
        character*4 name1
        character*20 filename
        real*8 xc_1(nn_conf,nn_mol), yc_1(nn_conf,nn_mol), zc_1(nn_conf,nn_mol)
        real*8 xcc(nn_mol),ycc(nn_mol),zcc(nn_mol)
        real*8 xcc_1, ycc_1, zcc_1, xcc_2, ycc_2, zcc_2
        real*8 xc_c, yc_c, zc_c
        real*8 dx_c, dy_c, dz_c
        real*8 num_ch4(nn_conf),num_co2(nn_conf), d_ch4(nn_conf),d_co2(nn_conf), r_c

    write(*,18)
18        format(/'     ###################################################################')
    write(*,19)
19        format('     #####################################################################')
        write(*,20)
20  format('    ###                                                                 ###')
        write(*,21)
21  format('   ###    MEAN SQUARED DISPLACEMENT --- Software Tools for Deffusion     ###')
    write(*,22)
22  format('   ##                                                                     ##')
        write(*,30)
30  format('   ##                       Version 1.1 December 2012                     ##')
    write(*,23)
23  format('   ##                                                                     ##')
        write(*,24)
24  format('   ##             Copyright <c> Zhen Yang & Guobing 2012-2022             ##')
    write(*,25)
25  format('   ##                                                                     ##')
        write(*,26)
26  format('   ###                         All Rights Reserved                       ###')
        write(*,27)
27  format('    ###                                                                 ###')
    write(*,28)
28        format('     #####################################################################')
    write(*,29)
29        format('      ###################################################################')





        write(*,11)
11        format(/'Enter the filename of archive (.arc):',$)
    read(*,*) filename
    open(1, file=filename)

        write(*,12)
12        format(/'Enter the filename of output (.txt):'$)
    read(*,*) filename
    open(2, file=filename)

        write(*,13)
13        format(/,'Enter the number of configurations:', $)
    read(*,*) n_conf
     
        write(*,14)
14        format(/,'Enter the number of methane molecules:', $)
    read(*,*) n_mol1

    write(*,15)
15        format(/,'Enter the number of carbon dioxide molecules:', $)
    read(*,*) n_mol2

        write(*,16)
16        format(/,'Enter the length of simulation box:', $)
    read(*,*) len

        write(*,38)
38        format(/,'Enter the time of dump(ps):', $)
    read(*,*) b

        write(*,*) ' Calculating ...... '

  do i =1 , n_conf
   
        do j = 1, n_mol1 + n_mol2

             read(1, *) a, name1, xc_1(i,j), yc_1(i,j), zc_1(i,j),num2(i,j)

    end do

  end do

  do j = 1, n_mol1 + n_mol2

    xcc(j) = xc_1(1,j)
   
    ycc(j) = yc_1(1,j)

        zcc(j) = zc_1(1,j)
  
  end do

  do i = 1 , n_conf-1
   
        do j = 1 , n_mol1 + n_mol2

      dx_c = xc_1(i+1,j) - xcc(j)
       
          dy_c = yc_1(i+1,j) - ycc(j)
       
          dz_c = zc_1(i+1,j) - zcc(j)

          dx_c = dx_c - anint(dx_c/len)*len

          dy_c = dy_c - anint(dy_c/len)*len

          dz_c = dz_c - anint(dz_c/len)*len

          xcc(j) = xc_1(i+1,j)
          
          ycc(j) = yc_1(i+1,j)

      zcc(j) = zc_1(i+1,j)
          
          xc_1(i+1,j) = xc_1(i,j) + dx_c

          yc_1(i+1,j) = yc_1(i,j) + dy_c

          zc_1(i+1,j) = zc_1(i,j) + dz_c
   
        end do
  
  end do
         
  do t = 1, n_conf/2

      num_ch4(t) = 0.0

      num_co2(t) = 0.0

  end do
  
  do i=1, n_conf

    do j=1, n_mol1 + n_mol2

      xcc_1 = xc_1(i,j)

          ycc_1 = yc_1(i,j)

          zcc_1 = zc_1(i,j)

          do t=1,n_conf/2

            if ( i + t <= n_conf ) then
             
                  xcc_2 = xc_1(i+t,j)

                  ycc_2 = yc_1(i+t,j)

                  zcc_2 = zc_1(i+t,j)

                  xc_c = xcc_1 - xcc_2
                  
                  yc_c = ycc_1 - ycc_2
                  
                  zc_c = zcc_1 - zcc_2
                  
                  r_c = xc_c**2 + yc_c**2 + zc_c**2
       
          if(num2(i,j)==1)then

                    num_ch4(t) = num_ch4(t) + r_c
         
                  elseif(num2(i,j)==3)then

                    num_co2(t) = num_co2(t) + r_c
        
                  end if

                else
                 
                  goto 10
            
                end if
       
      end do

10   continue
  
    end do
  
  end do

  do t=1, n_conf/2
   
        d_ch4(t) = real(num_ch4(t))/(real(n_conf-t)*real(n_mol1))

    d_co2(t) = real(num_co2(t))/(real(n_conf-t)*real(n_mol2))

        write(2,"(f12.6,2f14.6)" b*t, d_ch4(t)*0.01,d_co2(t)*0.01
  
  end do

end

¶ÔÓÚÉÏÃæµÄ³ÌÐò£¬Èç¹û¾Í°´ÕÕÉèÖõÄ×î´óÖµÊäÈëµÄ»°£¬±àÒëµÄʱºò»á³öÏÖwarning£¬Õâ¸öʱºò¸Ã²Éȡʲô°ì·¨ÄØ£¿
¼á³Ö²ÅÄÜʵÏÖ×Ô¼ºµÄÀíÏ룡
5Â¥2013-06-30 19:00:56
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

pippi6

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

¹¤³ÌºÍ¿ÆÑ§ÊýÖµ¼ÆËã×Éѯ

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

ÒýÓûØÌû:
5Â¥: Originally posted by 568286346 at 2013-06-30 19:00:56
program main
        implicit none
        integer   nn_conf, nn_mol, n_conf, n_mol1,n_mol2
        parameter (nn_conf=200000, nn_mol=196)
        real*8  len,b
        integer a, i, j, m, n, t,num2(nn_conf,nn_mol)
        parameter(m ...

ÄãÊÇ˵ÀàËÆÏÂÃæµÄwarning£¿
..\files\code_release\slab.exe : warning LNK4084: total image size 1106726912 exceeds max (268435456); image may not run

ºöÂÔ¾ÍÊÇ¡£Äã²»ÐèÒªËùÈκÎÊÂÇé¡£Æäʵ£¬Ç°Á½Ìì¸ÕÌÖÂÛ¹ýµÄ¡£
http://muchong.com/bbs/viewthread.php?tid=6059476&fpage=1
6Â¥2013-06-30 19:16:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

568286346

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

ÒýÓûØÌû:
6Â¥: Originally posted by pippi6 at 2013-06-30 19:16:51
ÄãÊÇ˵ÀàËÆÏÂÃæµÄwarning£¿
..\files\code_release\slab.exe : warning LNK4084: total image size 1106726912 exceeds max (268435456); image may not run

ºöÂÔ¾ÍÊÇ¡£Äã²»ÐèÒªËùÈκÎÊÂÇé¡£Æäʵ£¬Ç°Á½Ìì¸ÕÌÖÂÛ ...

ÎÒÉÏ´ÎÔËÐйý£¬Ìáʾarrays exceeded¡£
¼á³Ö²ÅÄÜʵÏÖ×Ô¼ºµÄÀíÏ룡
7Â¥2013-06-30 19:33:22
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

pippi6

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

¹¤³ÌºÍ¿ÆÑ§ÊýÖµ¼ÆËã×Éѯ

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï
568286346: ½ð±Ò+3, ¡ïÓаïÖú 2013-07-01 17:33:47
ÒýÓûØÌû:
7Â¥: Originally posted by 568286346 at 2013-06-30 19:33:22
ÎÒÉÏ´ÎÔËÐйý£¬Ìáʾarrays exceeded¡£...

ÎÒÔÚCVF6.6 Xp ÏÂÔËÐйý£¬Ã»ÎÊÌâ¡£ ÎÒÖ»ÓÐ3G¡£ ÄãÊÇʲôϵͳ°¡£¿ ¶øÇÒ£¬×¼È·µÄ´íÎóÐÅÏ¢ÊÇʲô°¡
8Â¥2013-06-30 19:51:37
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

568286346

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

ÄÚÈÝÒÑɾ³ý
¼á³Ö²ÅÄÜʵÏÖ×Ô¼ºµÄÀíÏ룡
9Â¥2013-06-30 19:57:32
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ѧ³ædj

½ð³æ (ÕýʽдÊÖ)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
Èç¹ûÐèÒªºÜ´óµÄÊý×飬ÄãÊÔÒ»ÊÔ¶¨ÒåÒ»¸ö¿É±ä´óСµÄÊý×é
ÇîÔò¶ÀÉÆÆäÊ£¬´ïÔò¼æ¼ÃÌìÏÂ
10Â¥2013-07-01 08:43:25
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ 568286346 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +5 Ò»ÑùYWY 2026-04-06 5/250 2026-04-06 13:08 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 0854Çóµ÷¼Á +4 assdll 2026-04-05 4/200 2026-04-06 12:29 by ÖзÉÔº¿Õ¹ÜѧԺÑ
[¿¼ÑÐ] µ÷¼Á +9 ²»·ê´º 2026-04-05 10/500 2026-04-06 11:58 by ¿¼ÑÐѧУÕеãÈË
[¿¼ÑÐ] 301Çóµ÷¼Á +7 ϸ°ûÏà¹Øµ°°× 2026-04-03 7/350 2026-04-06 11:47 by lijunpoly
[¿¼ÑÐ] 086000ÉúÎïÓëÒ½Ò©298µ÷¼ÁÇóÖú +9 ÔªÔªÇàÇà 2026-03-31 12/600 2026-04-05 21:03 by ѧԱ8dgXkO
[¿¼ÑÐ] 11408,335·Ö£¬±¾¿Æ211£¬Çóµ÷¼Á£¬¿Éתרҵ +5 öùÀæ´óöùÓã 2026-04-03 5/250 2026-04-04 22:49 by chongya
[¿¼ÑÐ] ¿¼Ñе÷¼Á +6 15615482637 2026-04-04 6/300 2026-04-04 22:43 by yu221
[¿¼ÑÐ] Çóµ÷¼Á +4 êɹ¦? 2026-04-03 4/200 2026-04-04 21:58 by hemengdong
[¿¼ÑÐ] Ò»Ö¾Ô¸ÉϺ£´óѧÉúÎïѧ346 +3 ÉϺ£´óѧ346µ÷¼Á 2026-04-03 3/150 2026-04-04 20:20 by dongzh2009
[¿¼ÑÐ] 331Çóµ÷¼Á +3 niby 2026-04-02 3/150 2026-04-04 19:56 by À¶ÔÆË¼Óê
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©¿Æ¼¼´óѧ²ÄÁϹ¤³Ì085601£¬Çóµ÷¼Á +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[¿¼ÑÐ] Çóµ÷¼Á +3 wos666 2026-04-03 3/150 2026-04-04 05:16 by gswylq
[¿¼ÑÐ] 282Çóµ÷¼Á +20 ycy1201 2026-04-01 22/1100 2026-04-04 00:42 by userper
[¿¼ÑÐ] 085601Ò»Ö¾Ô¸±±Àí325·ÖÇóµ÷¼Á +6 ÕÒµ÷¼Á£¬£¬ 2026-04-02 6/300 2026-04-03 22:20 by –¹Æ?
[¿¼ÑÐ] Çóµ÷¼Á +8 akdhjs 2026-04-03 8/400 2026-04-03 18:17 by ´÷άING
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÖÐũҵ071010£¬×Ü·Ö320Çóµ÷¼Á +7 À§À§À§À§À¤À¤ 2026-04-02 7/350 2026-04-03 17:26 by Yuena_Wang
[¿¼ÑÐ] ÄÜÔ´¶¯Á¦ µ÷¼Á +3 ²»ÆÆ²»Á¢0 2026-04-02 3/150 2026-04-02 12:46 by ffffjjjj
[¿¼ÑÐ] 262Çóµ÷¼Á +9 ÀøÖ¾Ò»¶¨·¢ÎÄÕ 2026-03-31 10/500 2026-04-01 12:22 by sunshine0013
[¿¼ÑÐ] 326Çóµ÷¼Á +4 áÌáÌ×Ð 2026-03-31 4/200 2026-04-01 09:58 by ÎҵĴ¬Îҵĺ£
[¿¼ÑÐ] 370Çóµ÷¼Á +3 080700µ÷¼Á 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û