应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助！急需macroscopic average electrostatic potential的小程序
71/1返回列表
查看: 1382  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wl7171759

新虫 (小有名气)

[求助] 求助！急需macroscopic average electrostatic potential的小程序 已有1人参与

想要看VASP输出文件LOCPOT包含的体系静电势信息，想要可处理LOCPOT文件得到macroscopic average electrostatic potential的小程序，希望得到大家的帮助！
或者可以实现高斯滤波的小程序也行啊，真心求助啊！！！
回复此楼

» 本帖已获得的红花（最新10朵）

dongbaojuan

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-04-01 12:05:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vinbert

木虫 (正式写手)
LZ，你找到macrscopic average的程序了么？能不能给我发一份，谢谢先。我邮箱是vinberg.lee@mail.dlut.edu.cn
赞一下 回复此楼
欲望是一切的源动力
2楼2014-04-25 09:57:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)
自己写吧，这个现在是个流行的话题。有程序的人也不会给出的。相信自己！！！！！
赞一下(1人) 回复此楼
3楼2014-04-25 19:46:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoboy185

金虫 (正式写手)
好多人都问这个问题i
记得有个小程序可以通过LOCPOT求出xy面的平均

得到这个xy的平均之后  可以直接用一些数据处理软件得到一个线性的平均，也就是那个macrscopic average
比如说origin
楼主可以试试看
赞一下 回复此楼
4楼2014-04-25 20:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
簡單說就是沿三個方向做平面積分與平均

這只是我為了某個目的去寫
輸出結果不見得能符合你的需求
所以
還是自己寫會比問別人還快
赞一下(2人) 回复此楼
5楼2014-04-25 20:12:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dongbaojuan

木虫 (小有名气)
送红花一朵
楼主得到这个处理程序了吗？得到的话，能不能发我一份啊？不胜感激~
邮箱：bjdong12s@imr.ac.cn
赞一下 回复此楼
6楼2015-01-13 18:48:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tsglss

铁杆木虫 (著名写手)
学习

发自小木虫Android客户端
回复此楼
做好自己
7楼2018-02-07 16:35:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wl7171759 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +3 西南交通专材355 2026-03-19 3/150 2026-03-20 09:34 by 每天只摆一小会
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 本人考085602 化学工程 专硕 +17 不知道叫什么！ 2026-03-15 19/950 2026-03-19 15:06 by 尽舜尧1
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
信息提示
关闭
请填处理意见
关闭