 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1394 | 回復(fù): 6 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
求助！急需macroscopic average electrostatic potential的小程序 已有1人參與
|
||||
|
想要看VASP輸出文件LOCPOT包含的體系靜電勢(shì)信息，想要可處理LOCPOT文件得到macroscopic average electrostatic potential的小程序，希望得到大家的幫助！ 或者可以實(shí)現(xiàn)高斯濾波的小程序也行啊，真心求助�。。�！ |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
dongbaojuan» 猜你喜歡
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑
已經(jīng)有8人回復(fù)
-
349求調(diào)劑
已經(jīng)有6人回復(fù)
-
329求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有9人回復(fù)
-
327求調(diào)劑
已經(jīng)有5人回復(fù)
-
求收留
已經(jīng)有5人回復(fù)
-
375求調(diào)劑
已經(jīng)有4人回復(fù)
-
340求調(diào)劑
已經(jīng)有6人回復(fù)
-
298求調(diào)劑
已經(jīng)有3人回復(fù)
-
085601材料工程找調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】macroscopic average electrostatic potential
已經(jīng)有15人回復(fù)
mywai520
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 339 (大學(xué)生)
- 金幣: 10425.6
- 紅花: 45
- 帖子: 1664
- 在線: 296.9小時(shí)
- 蟲號(hào): 265142
- 注冊(cè): 2006-07-09
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2014-04-25 19:46:48
vinbert
木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 6513.6
- 散金: 161
- 紅花: 9
- 帖子: 438
- 在線: 239.6小時(shí)
- 蟲號(hào): 1031862
- 注冊(cè): 2010-05-30
- 性別: GG
- 專業(yè): 電化學(xué)
2樓2014-04-25 09:57:58
xiaoboy185
金蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 788
- 散金: 248
- 紅花: 3
- 帖子: 702
- 在線: 472.3小時(shí)
- 蟲號(hào): 481831
- 注冊(cè): 2007-12-19
- 專業(yè): 函數(shù)論
4樓2014-04-25 20:04:33
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時(shí)
- 蟲號(hào): 585915
- 注冊(cè): 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ǎn)單說就是沿三個(gè)方向做平面積分與平均 這只是我為了某個(gè)目的去寫 輸出結(jié)果不見得能符合你的需求 所以 還是自己寫會(huì)比問別人還快 |
5樓2014-04-25 20:12:07
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +11 | Mr. Z 2026-03-25 | 11/550 |
|
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5 | SYA！ 2026-03-23 | 5/250 |
|
|
|
[考研]
 304材料求調(diào)劑
+4
|
鐘llll 2026-03-26 | 4/200 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研]
 085602化學(xué)工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|