 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1874 | 回復(fù): 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文件中數(shù)據(jù)輸出的格式，即第（10+natoms）行的三個數(shù)，natoms則是整個體系中的原子的個數(shù)，得到的outdata中的數(shù)據(jù)則是要求的z方向上的nz個數(shù)據(jù)，然后你自己作圖，outdata中的數(shù)據(jù)是你的縱坐標，橫坐標是你建立模型的z方向上的長度。[ Last edited by gyldmy on 2012-8-22 at 16:29 ] |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : LOCP.f90
2012-07-19 13:34:36, 616 bytes
» 收錄本帖的淘帖專輯推薦
 第一性原理 |
» 猜你喜歡
285求調(diào)劑
已經(jīng)有6人回復(fù)
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有6人回復(fù)
-
291求調(diào)劑
已經(jīng)有6人回復(fù)
-
315分求調(diào)劑
已經(jīng)有6人回復(fù)
-
277跪求調(diào)劑
已經(jīng)有8人回復(fù)
-
295求調(diào)劑
已經(jīng)有5人回復(fù)
-
295材料工程專碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
083000學碩274求調(diào)劑
已經(jīng)有4人回復(fù)
-
275求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿211院校 344分 東北農(nóng)業(yè)大學生物學學碩，求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 《英語流行口語脫口秀 時尚會話》掃描版[PDF]
已經(jīng)有373人回復(fù)
- 人•面
已經(jīng)有8人回復(fù)
- 靜電勢為什么有正負值？
已經(jīng)有102人回復(fù)
- CTD里的峰面積單位有要求嗎？？
已經(jīng)有17人回復(fù)
- 非常異常的靜電勢圖
已經(jīng)有14人回復(fù)
- 談?wù)凪olekel做靜電勢填色等值面圖的方法及誤區(qū)
已經(jīng)有50人回復(fù)
- 今天聽以科研領(lǐng)導(dǎo)說，今年的面上項目競爭非常慘烈，我校情況不樂觀
已經(jīng)有51人回復(fù)
- 有關(guān)電勢面掃描的問題，急！�。。。。。。�！
已經(jīng)有4人回復(fù)
- 掃面電鏡對冷凍切片樣品的要求。
已經(jīng)有7人回復(fù)
- 研究生面對畢業(yè)與戀愛？你選擇哪一方？
已經(jīng)有63人回復(fù)
- 中國科大信息二處面上項目
已經(jīng)有62人回復(fù)
- 【討論】關(guān)于靜電勢圖的描述
已經(jīng)有11人回復(fù)
- 【原創(chuàng)】金屬所在ECAP方面的研究工作進展
已經(jīng)有39人回復(fù)
| thx |
2樓2013-01-29 15:53:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 gyldmy 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程求調(diào)劑 +6 | 多久上課 2026-03-27 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 265求調(diào)劑 +5 | 小木蟲085600 2026-03-27 | 5/250 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 07化學280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +9 | minht 2026-03-26 | 9/450 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
|
[考研]
 085602化學工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 22 350 本科985求調(diào)劑，求老登收留 +4 | 李軼男003 2026-03-20 | 4/200 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
|
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
|
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|