| 查看: 2177 | 回复: 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 |
回复此楼» 猜你喜欢
我们学院常年位居 各学院 倒数第二。专业撤销的话，在编者有什么补偿？
已经有8人回复
-
公办双非，学生论文升学就业都不行。一本大学的论文升学就业没问题吧？
已经有4人回复
-
今年香江初审出结果了吗
已经有12人回复
-
迷茫毕业季
已经有7人回复
-
副教授跳槽
已经有20人回复
-
郑雅婷这种将来估计也会举报自己工作单位领导
已经有20人回复
-
还能读上博士吗
已经有13人回复
-
博士找工作真难
已经有20人回复
-
希望能25博士入学，可提前一年做科研助理
已经有7人回复
-
蹲一个6月7号江苏省青基公示
已经有30人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 精彩回帖：径向分布函数解析
已经有22人回复
- MS计算径向分布函数出现问题
已经有11人回复
- 怎样用forcite画径向分布函数
已经有14人回复
- 关于径向分布函数和配位数公式中的密度
已经有3人回复
- 对径向分布函数和配位数的不理解
已经有5人回复
- 径向分布函数与团簇结构关系
已经有5人回复
- 如何解释径向分布函数与r之间的关系？
已经有8人回复
- 径向分布函数的问题
已经有13人回复
- vmd-做径向分布函数
已经有6人回复
- forcite径向分布函数g(r)物理意义
已经有8人回复
- 请教 径向分布函数（ＲＤＦ） 问题
已经有12人回复
- 如何做径向分布函数和长程序LRO？
已经有4人回复
- 请问大家做径向分布函数的时候结果有没有趋近于1啊？
已经有7人回复
- 【求助】怎么从结果中提取径向分布函数、键长、键角
已经有17人回复
- 【讨论】陈正隆书上径向分布函数的公式是否有误？
已经有10人回复
- 【求助】应用MS计算径向分布函数
已经有15人回复
- 【求助成功】径向分布函数如何分析？
已经有14人回复
秋之愚
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 8139.8
- 红花: 1
- 帖子: 717
- 在线: 72.3小时
- 虫号: 2538185
- 注册: 2013-07-09
- 性别: GG
- 专业: 工程地质学
2楼2014-08-13 22:42:30
dummy1
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5919.2
- 散金: 32
- 红花: 1
- 帖子: 312
- 在线: 294.9小时
- 虫号: 270440
- 注册: 2006-08-06
- 专业: 理论和计算化学
3楼2014-08-14 02:45:40
秋之愚
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 8139.8
- 红花: 1
- 帖子: 717
- 在线: 72.3小时
- 虫号: 2538185
- 注册: 2013-07-09
- 性别: GG
- 专业: 工程地质学
4楼2014-08-14 19:01:25
dummy1
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5919.2
- 散金: 32
- 红花: 1
- 帖子: 312
- 在线: 294.9小时
- 虫号: 270440
- 注册: 2006-08-06
- 专业: 理论和计算化学
5楼2014-08-16 02:33:13
秋之愚
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 8139.8
- 红花: 1
- 帖子: 717
- 在线: 72.3小时
- 虫号: 2538185
- 注册: 2013-07-09
- 性别: GG
- 专业: 工程地质学
6楼2014-08-16 16:09:07
30