| ²é¿´: 1706 | »Ø¸´: 12 | ||
568286346ľ³æ (СÓÐÃûÆø)
|
[ÇóÖú]
Êý×鳬³öÎÊÌâ
|
|
| ¸÷λ³æÓÒ£¬±¾ÈËÔÚÓÃfortran±à³ÌÐò´¦ÀíÎÊÌâʱÓöµ½ÁËÊý×鳬³öµÄÎÊÌ⣨³¬³öÁËfortranĬÈÏÔÊÐíµÄ×î´óÖµ£©£¬ÇëÎÊÓÐʲô·½·¨¿ÉÒÔ½â¾öÕâ¸öÎÊÌâÂ𣿠ÒòΪÎÒÔÚ´¦ÀíÊý¾ÝʱÐèÒª¿ªºÜ´óµÄÊý×飬ËùÒÔ»¹Íû´ó¼ÒÌṩµãÒâ¼û°¡¡£ |
» ²ÂÄãϲ»¶
»¯¹¤Çóµ÷¼Á£¡
ÒѾÓÐ9È˻ظ´
301Çóµ÷¼Á
ÒѾÓÐ13È˻ظ´
308Çóµ÷¼Á
ÒѾÓÐ13È˻ظ´
²ÄÁϹ¤³Ì310ר˶µ÷¼Á
ÒѾÓÐ15È˻ظ´
304Çóµ÷¼Á
ÒѾÓÐ9È˻ظ´
081700£¬311£¬Çóµ÷¼Á
ÒѾÓÐ16È˻ظ´
080500Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
312Çóµ÷¼Á
ÒѾÓÐ10È˻ظ´
070300»¯Ñ§279Çóµ÷¼Á
ÒѾÓÐ19È˻ظ´
Çóµ÷¼Á Ò»Ö¾Ô¸Î÷ÄϽ»Í¨´óѧ085701»·¾³¹¤³Ì 282·Ö
ÒѾÓÐ9È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
ÒѾÅźÃÐòµÄÊý×éµÄÇóºÍÎÊÌâ
ÒѾÓÐ13È˻ظ´
ÇóÖúÊý×éÔ½½çÎÊÌ⣡
ÒѾÓÐ8È˻ظ´
sshÉÏÔËÐгÌÐò³öÏÖsegmentation fault
ÒѾÓÐ6È˻ظ´
ÇóÖú£º¹ØÓÚfortran½¨Êý×éÎÊÌ⣬лл£¡
ÒѾÓÐ4È˻ظ´
VC++Öк¯Êý·µ»ØÊý×éÖ¸Õë»òÕß´øÖ¸ÕëµÄ½á¹¹ÌåµÄ±àÒ뷽ʽÊÇ·ñ¿ÉÈ¡£¿
ÒѾÓÐ6È˻ظ´
Êý×鳬³öÏÞÖÆÔõô°ì£¿
ÒѾÓÐ3È˻ظ´
¹ØÓÚcÓïÑÔµÄÒ»¸öСÎÊÌâ ×øµÈ»Ø¸´°¡ ¼±°¡ ÊǹØÓÚ½¨Á¢Êý×éµÄ
ÒѾÓÐ9È˻ظ´
¡¾ÇóÖú¡¿C# ÖÐÈçºÎÉùÃ÷½á¹¹ÌåÊý×é²¢¸³Öµ
ÒѾÓÐ13È˻ظ´
¡¾ÇóÖú¡¿DPD³ÌÐò±àд±ß½çÈçºÎ´¦Àí£¿
ÒѾÓÐ12È˻ظ´
¡¾ÇóÖú¡¿Ò»¸öÇóÊý×é¾ø¶ÔÖµµÄСÎÊÌâ
ÒѾÓÐ7È˻ظ´
¡¾ÇóÖú¡¿C++³ÌÐòËùÐèÄÚ´æÔ¶Ô¶´óÓÚÎïÀíÄÚ´æÔõô°ì£¿
ÒѾÓÐ8È˻ظ´

pippi6
Ìú¸Ëľ³æ (ÖøÃûдÊÖ)
¹¤³ÌºÍ¿ÆÑ§ÊýÖµ¼ÆËã×Éѯ
- Ó¦Öú: 413 (˶ʿ)
- ¹ó±ö: 0.002
- ½ð±Ò: 7116.5
- É¢½ð: 15
- ºì»¨: 63
- Ìû×Ó: 1639
- ÔÚÏß: 798.9Сʱ
- ³æºÅ: 2469437
- ×¢²á: 2013-05-14
- רҵ: ¼ÆËãÊýѧÓë¿ÆÑ§¹¤³Ì¼ÆËã
| ²éÁËÒ»ÏÂ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
pippi6
Ìú¸Ëľ³æ (ÖøÃûдÊÖ)
¹¤³ÌºÍ¿ÆÑ§ÊýÖµ¼ÆËã×Éѯ
- Ó¦Öú: 413 (˶ʿ)
- ¹ó±ö: 0.002
- ½ð±Ò: 7116.5
- É¢½ð: 15
- ºì»¨: 63
- Ìû×Ó: 1639
- ÔÚÏß: 798.9Сʱ
- ³æºÅ: 2469437
- ×¢²á: 2013-05-14
- רҵ: ¼ÆËãÊýѧÓë¿ÆÑ§¹¤³Ì¼ÆËã
2Â¥2013-06-30 18:06:31
568286346
ľ³æ (СÓÐÃûÆø)
- Ó¦Öú: 31 (СѧÉú)
- ½ð±Ò: 789
- Ìû×Ó: 288
- ÔÚÏß: 153.3Сʱ
- ³æºÅ: 1785423
- ×¢²á: 2012-04-29
- רҵ: ÀíÂۺͼÆË㻯ѧ

3Â¥2013-06-30 18:25:11
pippi6
Ìú¸Ëľ³æ (ÖøÃûдÊÖ)
¹¤³ÌºÍ¿ÆÑ§ÊýÖµ¼ÆËã×Éѯ
- Ó¦Öú: 413 (˶ʿ)
- ¹ó±ö: 0.002
- ½ð±Ò: 7116.5
- É¢½ð: 15
- ºì»¨: 63
- Ìû×Ó: 1639
- ÔÚÏß: 798.9Сʱ
- ³æºÅ: 2469437
- ×¢²á: 2013-05-14
- רҵ: ¼ÆËãÊýѧÓë¿ÆÑ§¹¤³Ì¼ÆËã
4Â¥2013-06-30 18:39:21
568286346
ľ³æ (СÓÐÃûÆø)
- Ó¦Öú: 31 (СѧÉú)
- ½ð±Ò: 789
- Ìû×Ó: 288
- ÔÚÏß: 153.3Сʱ
- ³æºÅ: 1785423
- ×¢²á: 2012-04-29
- רҵ: ÀíÂۺͼÆË㻯ѧ
|
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.01end do end ¶ÔÓÚÉÏÃæµÄ³ÌÐò£¬Èç¹û¾Í°´ÕÕÉèÖõÄ×î´óÖµÊäÈëµÄ»°£¬±àÒëµÄʱºò»á³öÏÖwarning£¬Õâ¸öʱºò¸Ã²Éȡʲô°ì·¨ÄØ£¿ |

5Â¥2013-06-30 19:00:56
pippi6
Ìú¸Ëľ³æ (ÖøÃûдÊÖ)
¹¤³ÌºÍ¿ÆÑ§ÊýÖµ¼ÆËã×Éѯ
- Ó¦Öú: 413 (˶ʿ)
- ¹ó±ö: 0.002
- ½ð±Ò: 7116.5
- É¢½ð: 15
- ºì»¨: 63
- Ìû×Ó: 1639
- ÔÚÏß: 798.9Сʱ
- ³æºÅ: 2469437
- ×¢²á: 2013-05-14
- רҵ: ¼ÆËãÊýѧÓë¿ÆÑ§¹¤³Ì¼ÆËã
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
|
ÄãÊÇ˵ÀàËÆÏÂÃæµÄ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
568286346
ľ³æ (СÓÐÃûÆø)
- Ó¦Öú: 31 (СѧÉú)
- ½ð±Ò: 789
- Ìû×Ó: 288
- ÔÚÏß: 153.3Сʱ
- ³æºÅ: 1785423
- ×¢²á: 2012-04-29
- רҵ: ÀíÂۺͼÆË㻯ѧ

7Â¥2013-06-30 19:33:22
pippi6
Ìú¸Ëľ³æ (ÖøÃûдÊÖ)
¹¤³ÌºÍ¿ÆÑ§ÊýÖµ¼ÆËã×Éѯ
- Ó¦Öú: 413 (˶ʿ)
- ¹ó±ö: 0.002
- ½ð±Ò: 7116.5
- É¢½ð: 15
- ºì»¨: 63
- Ìû×Ó: 1639
- ÔÚÏß: 798.9Сʱ
- ³æºÅ: 2469437
- ×¢²á: 2013-05-14
- רҵ: ¼ÆËãÊýѧÓë¿ÆÑ§¹¤³Ì¼ÆËã
8Â¥2013-06-30 19:51:37
568286346
ľ³æ (СÓÐÃûÆø)
- Ó¦Öú: 31 (СѧÉú)
- ½ð±Ò: 789
- Ìû×Ó: 288
- ÔÚÏß: 153.3Сʱ
- ³æºÅ: 1785423
- ×¢²á: 2012-04-29
- רҵ: ÀíÂۺͼÆË㻯ѧ

9Â¥2013-06-30 19:57:32
ѧ³ædj
½ð³æ (ÕýʽдÊÖ)
- Ó¦Öú: 18 (СѧÉú)
- ½ð±Ò: 892.4
- É¢½ð: 63
- Ìû×Ó: 688
- ÔÚÏß: 89.9Сʱ
- ³æºÅ: 1829966
- ×¢²á: 2012-05-22
- ÐÔ±ð: GG
- רҵ: Ó¦ÓõØÇòÎïÀíѧ

10Â¥2013-07-01 08:43:25














»Ø¸´´ËÂ¥
b*t, d_ch4(t)*0.01,d_co2(t)*0.01