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 ]

回复此楼