| 查看: 1513 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求助！急需macroscopic average electrostatic potential的小程序 已有1人参与
|
|||
|
想要看VASP输出文件LOCPOT包含的体系静电势信息，想要可处理LOCPOT文件得到macroscopic average electrostatic potential的小程序，希望得到大家的帮助！ 或者可以实现高斯滤波的小程序也行啊，真心求助啊！！！ |
回复此楼» 本帖已获得的红花（最新10朵）
dongbaojuan» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
-
是这周出结果还是下周出结果？
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有49人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有5人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有7人回复
-
filecode
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
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
