| 查看: 1864 | 回复: 1 | ||||
[交流]
vasp 扫面静电势（electrostatic potential）线 已有1人参与
|
|
本来在也在论坛里看到了一个处理LOCPOT文件的fortran语言，但是编辑不通过，所以自己写了一个！在自己的机子上是可以运行的，拿来分享一下！下面贴出程序语言： parameter(nx=60,ny=60,nz=60, natoms=3) real pot(nx*ny*nz),pot2(nx,ny,nz),pot3(nz) open（1,file="LOCPOT"） open（2,file="outdata"） do i=1,7 read(1,*) enddo do i=1,natoms+3 read(1,*) enddo do i=1,nx*ny*nz/5 read(1,*)pot((i-1)*5+1),pot((i-1)*5+2),pot((i-1)*5+3),pot((i-1)*5+4),pot((i-1)*5+5) enddo a=0 do h=1,nx do k=1,ny do l=1,nz a=a+1 pot2(h,k,l)=pot(a) enddo enddo enddo pot3(nz)=0 do i=1,nz do j=1,nx do k=1,ny pot3(i)=pot3(i)+pot2(k,j,i)/(nx*ny/1.0d0) enddo enddo write(2,*)pot3(i) enddo close(1) close(2) end 修改的话，只需要改红字的部分，其中nx,ny,nz的值分别是你的LOCPOT文件中数据输出的格式，即第（10+natoms）行的三个数，natoms则是整个体系中的原子的个数，得到的outdata中的数据则是要求的z方向上的nz个数据，然后你自己作图，outdata中的数据是你的纵坐标，横坐标是你建立模型的z方向上的长度。[ Last edited by gyldmy on 2012-8-22 at 16:29 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : LOCP.f90
2012-07-19 13:34:36, 616 bytes
» 收录本帖的淘帖专辑推荐
 第一性原理 |
» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
《英语流行口语脱口秀 时尚会话》扫描版[PDF]
已经有373人回复
- 人•面
已经有8人回复
- 静电势为什么有正负值？
已经有102人回复
- CTD里的峰面积单位有要求吗？？
已经有17人回复
- 非常异常的静电势图
已经有14人回复
- 谈谈Molekel做静电势填色等值面图的方法及误区
已经有50人回复
- 今天听以科研领导说，今年的面上项目竞争非常惨烈，我校情况不乐观
已经有51人回复
- 有关电势面扫描的问题，急！！！！！！！！！
已经有4人回复
- 扫面电镜对冷冻切片样品的要求。
已经有7人回复
- 研究生面对毕业与恋爱？你选择哪一方？
已经有63人回复
- 中国科大信息二处面上项目
已经有62人回复
- 【讨论】关于静电势图的描述
已经有11人回复
- 【原创】金属所在ECAP方面的研究工作进展
已经有39人回复
| thx |
2楼2013-01-29 15:53:33
