版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (6964)
>考研 (2002)
>導師招生 (1633)
>文獻求助 (519)
>考博 (316)
>蟲友互識 (307)
>碩博家園 (168)
>休閑灌水 (154)
>博后之家 (129)
>基金申請 (118)
>招聘信息布告欄 (114)
>論文投稿 (99)
>第一性原理 (79)
>教師之家 (74)
>公派出國 (50)
>外文書籍求助 (32)

返回列表

gyldmy

銅蟲 (小有名氣)

1ST強帖: 1
應助: 12 (小學生)
金幣: 5
紅花: 1
帖子: 180
在線: 259.3小時
蟲號: 1393844
注冊: 2011-09-07
專業(yè): 有機合成

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

回復此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內容負責。
本內容由用戶自主發(fā)布，如果其內容涉及到知識產權問題，其責任在于用戶本人，如對版權有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : LOCP.f90

2012-07-19 13:34:36, 616 bytes

» 收錄本帖的淘帖專輯推薦

第一性原理

» 猜你喜歡

考研調劑已經有4人回復
281求調劑已經有4人回復
0805 316求調劑已經有6人回復
085601求調劑總分293英一數二已經有3人回復
08工學調劑已經有17人回復
340求調劑已經有4人回復
311求調劑已經有3人回復
食品專碩一志愿雙一流 328 已經有4人回復
材料調劑已經有6人回復
300分，材料，求調劑，英一數二已經有5人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

highpeter

thx

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