| 查看: 2718 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
秋之愚新虫 (正式写手)
|
[求助]
径向分布函数怎么画 已有1人参与
|
|
|
求助各位大神，我在做Lennard-Jones势流体的分子动力学模拟时，能量图基本上已经画出来了，可径向分布函数的计算怎么弄都弄不出来，能有人帮帮我吗？在给出来的这个程序中switch在主程序中该怎么调用啊？？？ 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 |
回复此楼» 猜你喜欢
某些机构，以效率低为荣，以效率低作为存在感
已经有4人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有3人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
-
人气不行了
已经有11人回复
-
2026国自然函评费到账
已经有22人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有3人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
-
今天基金会出结果吗？20260819
已经有18人回复
-
我面上完蛋了
已经有14人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 怎样用forcite画径向分布函数
已经有14人回复
- 关于径向分布函数和配位数公式中的密度
已经有3人回复
- 径向分布函数与团簇结构关系
已经有5人回复
- 如何解释径向分布函数与r之间的关系？
已经有8人回复
- 径向分布函数的问题
已经有13人回复
- vmd-做径向分布函数
已经有6人回复
- forcite径向分布函数g(r)物理意义
已经有8人回复
- 请教 径向分布函数（ＲＤＦ） 问题
已经有12人回复
- 如何做径向分布函数和长程序LRO？
已经有4人回复
- 请问大家做径向分布函数的时候结果有没有趋近于1啊？
已经有7人回复
- 【求助】怎么从结果中提取径向分布函数、键长、键角
已经有17人回复
- 【讨论】陈正隆书上径向分布函数的公式是否有误？
已经有10人回复
- 【求助】应用MS计算径向分布函数
已经有15人回复
- 【求助成功】径向分布函数如何分析？
已经有14人回复
dummy1
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5841.7
- 散金: 32
- 红花: 1
- 帖子: 312
- 在线: 295.8小时
- 虫号: 270440
- 注册: 2006-08-06
- 专业: 理论和计算化学
3楼2014-08-14 02:45:40
秋之愚
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 8145.6
- 红花: 1
- 帖子: 719
- 在线: 72.3小时
- 虫号: 2538185
- 注册: 2013-07-09
- 性别: GG
- 专业: 工程地质学
4楼2014-08-14 19:01:25
dummy1
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5841.7
- 散金: 32
- 红花: 1
- 帖子: 312
- 在线: 295.8小时
- 虫号: 270440
- 注册: 2006-08-06
- 专业: 理论和计算化学
5楼2014-08-16 02:33:13
秋之愚
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 8145.6
- 红花: 1
- 帖子: 719
- 在线: 72.3小时
- 虫号: 2538185
- 注册: 2013-07-09
- 性别: GG
- 专业: 工程地质学
6楼2014-08-16 16:09:07
20