版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

gyldmy

铜虫 (小有名气)

1ST强帖: 1
应助: 12 (小学生)
金币: 5
红花: 1
帖子: 180
在线: 259.3小时
虫号: 1393844
注册: 2011-09-07
专业: 有机合成

[交流] 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

» 收录本帖的淘帖专辑推荐

第一性原理

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2012-07-19 13:44:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

highpeter

thx

2楼2013-01-29 15:53:33