| 5 | 1/1 | 返回列表 |
| 查看: 2537 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
秋之愚新蟲 (正式寫手)
|
[求助]
徑向分布函數(shù)怎么畫 已有1人參與
|
|
|
求助各位大神,我在做Lennard-Jones勢(shì)流體的分子動(dòng)力學(xué)模擬時(shí),能量圖基本上已經(jīng)畫出來了,可徑向分布函數(shù)的計(jì)算怎么弄都弄不出來,能有人幫幫我嗎?在給出來的這個(gè)程序中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 |

新蟲 (正式寫手)

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 07化學(xué)280分求調(diào)劑 +4 | 722865 2026-03-23 | 4/200 |
|
|
[考研] 081700 調(diào)劑 267分 +7 | 迷人的哈哈 2026-03-23 | 7/350 |
|
|
[考研] 材料專碩找調(diào)劑 +3 | 哈哈哈吼吼吼哈 2026-03-23 | 3/150 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項(xiàng)目,四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|