| 5 | 1/1 | 返回列表 |
| 查看: 2536 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
秋之愚新蟲 (正式寫手)
|
[求助]
徑向分布函數(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 |

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

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

![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|
百度網(wǎng)盤 |
360云盤 |
千易網(wǎng)盤 |
華為網(wǎng)盤
在新窗口頁面中打開自己喜歡的網(wǎng)盤網(wǎng)站,將文件上傳后,然后將下載鏈接復(fù)制到帖子內(nèi)容中就可以了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +7 | 紙魚ly 2026-03-21 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5 | 葵梓衛(wèi)隊 2026-03-18 | 7/350 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|