| 查看: 1515 | 回复: 6 | |||
[求助]
求助！急需macroscopic average electrostatic potential的小程序 已有1人参与
|
|
想要看VASP输出文件LOCPOT包含的体系静电势信息，想要可处理LOCPOT文件得到macroscopic average electrostatic potential的小程序，希望得到大家的帮助！ 或者可以实现高斯滤波的小程序也行啊，真心求助啊！！！ |
回复此楼» 本帖已获得的红花（最新10朵）
dongbaojuan» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
-
是这周出结果还是下周出结果？
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
vinbert
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 6513.6
- 散金: 161
- 红花: 9
- 帖子: 438
- 在线: 239.6小时
- 虫号: 1031862
- 注册: 2010-05-30
- 性别: GG
- 专业: 电化学
2楼2014-04-25 09:57:58
mywai520
铁杆木虫 (著名写手)
- 应助: 339 (大学生)
- 金币: 10425.6
- 红花: 45
- 帖子: 1664
- 在线: 296.9小时
- 虫号: 265142
- 注册: 2006-07-09
- 专业: 凝聚态物性 II ：电子结构
3楼2014-04-25 19:46:48
xiaoboy185
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 788
- 散金: 248
- 红花: 3
- 帖子: 702
- 在线: 472.3小时
- 虫号: 481831
- 注册: 2007-12-19
- 专业: 函数论
4楼2014-04-25 20:04:33
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fzx2008: 金币+5, 谢谢指导，谢谢分享！ 2014-04-28 19:23:37
wl7171759: 金币+10, ★★★很有帮助, 受益了，确实自己得学着谢谢程序 2014-05-05 12:37:22
fzx2008: 金币+5, 谢谢指导，谢谢分享！ 2014-04-28 19:23:37
wl7171759: 金币+10, ★★★很有帮助, 受益了，确实自己得学着谢谢程序 2014-05-05 12:37:22
|
三樓說的很好 程式還是要學著寫，否則未?淼穆泛茈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 簡單說就是沿三個方向做平面積分與平均 這只是我為了某個目的去寫 輸出結果不見得能符合你的需求 所以 還是自己寫會比問別人還快 |
5楼2014-04-25 20:12:07
dongbaojuan
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 2364.6
- 散金: 1100
- 红花: 6
- 帖子: 253
- 在线: 161.5小时
- 虫号: 1504764
- 注册: 2011-11-22
- 专业: 金属功能材料
6楼2015-01-13 18:48:31
tsglss
铁杆木虫 (著名写手)
- 应助: 41 (小学生)
- 金币: 6014.2
- 红花: 19
- 帖子: 1875
- 在线: 940小时
- 虫号: 321641
- 注册: 2007-03-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2018-02-07 16:35:57
