| 6 | 1/1 | 返回列表 |
| 查看: 2533 | 回復(fù): 5 | ||
秋之愚新蟲 (正式寫手)
|
[求助]
徑向分布函數(shù)怎么畫 已有1人參與
|
|
求助各位大神,我在做Lennard-Jones勢流體的分子動力學(xué)模擬時,能量圖基本上已經(jīng)畫出來了,可徑向分布函數(shù)的計算怎么弄都弄不出來,能有人幫幫我嗎?在給出來的這個程序中switch在主程序中該怎么調(diào)用??? subroutine gr(switch,box,density,rx,ry,rz,g) dimension rx(108),ry(108),rz(108),g(125) double precision rx,ry,rz,g double precision delg,pi,box,r22,rr,density,nid integer ngr,ig,nhis,switch nhis=125 delg = box/(2*nhis) pi=3.1415926 ! if (switch.eq.0)then ngr=0 !delg=box/(2*nhis) do i=1,nhis g(i)=0 enddo ! elseif (switch.eq.1)then ngr=ngr+1 do i=1,107 do j=i+1,108 xr=rx(i)-rx(j) xr=xr-box*nint(xr/box) yr=ry(i)-ry(j) yr=yr-box*nint(yr/box) zr=rz(i)-rz(j) zr=zr-box*nint(zr/box) r22=xr*xr+yr*yr+zr*zr rr=sqrt(r22) if(rr.lt.box/2) then ig=int(rr/delg) g(ig)=g(ig)+2 endif enddo enddo ! elseif(switch.eq.2) then do i=1,nhis r=delg*(i+0.5) vb=((i+1)**3-i**3)*delg**3 nid=(4/3)*pi*vb*density g(i)=g(i)/(ngr*108.0*nid) write(*,*) nid,r,g(i) enddo ! endif return end |

新蟲 (正式寫手)

木蟲 (正式寫手)
新蟲 (正式寫手)

木蟲 (正式寫手)
新蟲 (正式寫手)

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 299求調(diào)劑 +4 | shxchem 2026-03-20 | 5/250 |
|
|
[考研] 317求調(diào)劑 +8 | 申子申申 2026-03-19 | 13/650 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|