 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1397 | 回復(fù): 6 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
求助！急需macroscopic average electrostatic potential的小程序 已有1人參與
|
||||
|
想要看VASP輸出文件LOCPOT包含的體系靜電勢信息，想要可處理LOCPOT文件得到macroscopic average electrostatic potential的小程序，希望得到大家的幫助！ 或者可以實(shí)現(xiàn)高斯濾波的小程序也行啊，真心求助�。。�！ |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
dongbaojuan» 猜你喜歡
317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑
已經(jīng)有11人回復(fù)
-
332求調(diào)劑
已經(jīng)有12人回復(fù)
-
材料專碩 085600求調(diào)劑
已經(jīng)有5人回復(fù)
-
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可
已經(jīng)有14人回復(fù)
-
283求調(diào)劑（080500）
已經(jīng)有14人回復(fù)
-
材料工程專碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑 有機(jī)化學(xué)考研356分
已經(jīng)有11人回復(fù)
-
環(huán)境科學(xué)與工程334分求調(diào)劑
已經(jīng)有8人回復(fù)
-
考研調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 【求助】macroscopic average electrostatic potential
已經(jīng)有15人回復(fù)
xiaoboy185
金蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 788
- 散金: 248
- 紅花: 3
- 帖子: 702
- 在線: 472.3小時
- 蟲號: 481831
- 注冊: 2007-12-19
- 專業(yè): 函數(shù)論
4樓2014-04-25 20:04:33
vinbert
木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 6513.6
- 散金: 161
- 紅花: 9
- 帖子: 438
- 在線: 239.6小時
- 蟲號: 1031862
- 注冊: 2010-05-30
- 性別: GG
- 專業(yè): 電化學(xué)
2樓2014-04-25 09:57:58
mywai520
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 339 (大學(xué)生)
- 金幣: 10425.6
- 紅花: 45
- 帖子: 1664
- 在線: 296.9小時
- 蟲號: 265142
- 注冊: 2006-07-09
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2014-04-25 19:46:48
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時
- 蟲號: 585915
- 注冊: 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fzx2008: 金幣+5, 謝謝指導(dǎo)，謝謝分享！ 2014-04-28 19:23:37
wl7171759: 金幣+10, ★★★很有幫助, 受益了，確實(shí)自己得學(xué)著謝謝程序 2014-05-05 12:37:22
fzx2008: 金幣+5, 謝謝指導(dǎo)，謝謝分享！ 2014-04-28 19:23:37
wl7171759: 金幣+10, ★★★很有幫助, 受益了，確實(shí)自己得學(xué)著謝謝程序 2014-05-05 12:37:22
|
三樓說的很好 程式還是要學(xué)著寫，否則未?淼穆泛茈y走 implicit real*8(a-h,o-z) character*32 tmp,t dimension nt(10),r(99) logical alive allocatable ep(  ,grid(:,:,,na(!======================================================================= inquire(file='LOCPOT',exist=alive) if(alive) then write(6,*) ; write(6,*) 'LOCPOT file is found ~' open( 1,file='LOCPOT') open(11,file='LOCPOT-x.dat') open(12,file='LOCPOT-y.dat') open(13,file='LOCPOT-z.dat') do i=1,5 ; read(1,*) tmp ; end do 98 read(1,*,iostat=istat) nt ; ntype=count(nt>0) if(ntype .eq. 0) go to 98 allocate(na(ntype)) do i=1,ntype ; na(i)=nt(i) ; end do ntot=sum(na) ; allocate(ep(ntot)) do i=1,ntot ; read(1,*) tmp ; end do read(1,*) nx,ny,nz ; allocate(grid(nx,ny,nz)) read(1,*) (((grid(i,j,k),i=1,nx),j=1,ny),k=1,nz) sumgrid=sum(grid) write(6,*) " --------------- results from LOCPOT -------------------- " write(6,10) sumgrid write(6,11) nx,ny,nz write(11,'(f15.5)') (sum(grid(i,:, )/real(ny*nz),i=1,nx)write(12,'(f15.5)') (sum(grid(:,i, )/real(nx*nz),i=1,ny)write(13,'(f15.5)') (sum(grid(:,:,i))/real(nx*ny),i=1,nz) write(6,12) sumgrid/real(nx*ny*nz) write(6,*) " -------------------------------------------------------- " 10 format(" total electrostatic potential = ",f20.5) 11 format(" number of grids = ",5x,3i5) 12 format("average electrostatic potential = ",f20.5) close(1) ; deallocate(ep) end if 這是小弟寫的 fortran 處理 LOCPOT，很短 輸入檔: LOCPOT 輸出檔: LOCPOT-x.dat LOCPOT-y.dat LOCPOT-z.dat 簡單說就是沿三個方向做平面積分與平均 這只是我為了某個目的去寫 輸出結(jié)果不見得能符合你的需求 所以 還是自己寫會比問別人還快 |
5樓2014-04-25 20:12:07
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +11 | 辰沐5211314 2026-03-26 | 11/550 |
|
|---|---|---|---|---|
|
|
[考研] 277跪求調(diào)劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|
[考研] 求收留 +5 | 1943443204 2026-03-28 | 5/250 |
|
|
|
[考研] 311（085601）求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研]
 一志愿 西北大學(xué) 總分282 英語一62 求調(diào)劑
+7
|
18419759900 2026-03-25 | 8/400 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
|
[考研] 085600，材料與化工321分，求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | 研研，接電話 2026-03-24 | 7/350 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|