| 查看: 1116 | 回复: 7 | ||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | ||
[求助]
求一个能砍PROCAR的程序,或者编一个程序砍PROCAR
|
||
|
如题,我计算出的PROCAR是带有f轨道的,但是我不希望f轨道存在,想砍掉所有的f轨道,附件中是希望被砍掉的轨道,格式为: ion s py pz px dxy dyz dz2 dxz dx2 f-3 f-2 f-1 f0 f1 f2 f3 tot 希望砍为: ion s py pz px dxy dyz dz2 dxz dx2 tot 在此特奉上100金币,求大侠出手相助 |
» 本帖附件资源列表
-
欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com - 附件 1 : Abnormal-PROCAR.7z
2015-04-08 21:35:48, 3.91 M
» 猜你喜欢
职称评审没过,求安慰
已经有41人回复
回收溶剂求助
已经有7人回复
硝基苯如何除去
已经有3人回复
A期刊撤稿
已经有4人回复
垃圾破二本职称评审标准
已经有17人回复
投稿Elsevier的Neoplasia杂志,到最后选publishing options时页面空白,不能完成投稿
已经有22人回复
EST投稿状态问题
已经有7人回复
毕业后当辅导员了,天天各种学生超烦
已经有4人回复
求助文献
已经有3人回复
三无产品还有机会吗
已经有6人回复
» 本主题相关价值贴推荐,对您同样有帮助:
PROCAR里面的相因子究竟是如何对应到波函数系数的
已经有3人回复
用dos_procar处理dos数据不准
已经有3人回复
PROCAR里数据是什么意思
已经有8人回复
dos-procar.f90
已经有94人回复
如何把16位dos程序变成64位程序
已经有8人回复
vasp.5.3.3 的PROCAR 和PARCHG结果相差很大
已经有8人回复
从PROCAR中提取态密度的shell代码
已经有42人回复
【求助成功】Band-procar 和Dos-procar 程序执行出错
已经有8人回复
【求助】procar文件如何提取dxy,dz2,dy等的LDOS
已经有6人回复
【求助】请教一下PROCAR的格式
已经有11人回复


8楼2015-04-12 21:24:14
★ ★ ★ ★ ★
fzx2008: 金币+5, 多谢出手 2015-04-10 08:27:51
fzx2008: 金币+5, 多谢出手 2015-04-10 08:27:51
|
隨手寫了個非常陽春的 fortran 測過~ 應該是沒有問題的 如下 implicit real*8(a-h,o-z) character*128 t,tt open(1,file="Abnormal-PROCAR" ![]() open(2,file="Abnormal-PROCAR.out" ![]() read(1,'(a28)') t ; write(2,'(a28)') t do i=1,2 read(1,'(a66)') t ; write(2,'(a66)') t do nk=1,30 read(1,'(a )') t ; write(2,'(a )') t read(1,'(a75)') t ; write(2,'(a75)') t do nb=1,100 read(1,'(a )') t ; write(2,'(a )') t read(1,'(a50)') t ; write(2,'(a50)') t read(1,'(a )') t ; write(2,'(a )') t read(1,'(a66,49x,a7)') t,tt ; write(2,'(a66,a7)') t,tt do natom=1,32 read(1,'(a66,49x,a7)') t,tt ; write(2,'(a66,a7)') t,tt end do read(1,'(a66,49x,a7)') t,tt ; write(2,'(a66,a7)') t,tt read(1,'(a )') t ; write(2,'(a )') t read(1,'(a66)') t ; write(2,'(a66)') t do natom=1,32*2 read(1,'(a66)') t ; write(2,'(a66)') t end do end do read(1,'(a )') t ; write(2,'(a )') t end do end do end 金幣就免了~ 大放送 |
» 本帖已获得的红花(最新10朵)
2楼2015-04-09 23:19:33
|
討人厭的害羞圖案又出現~~ 再貼一次 implicit real*8(a-h,o-z) character*128 t,tt open(1,file="Abnormal-PROCAR" ) open(2,file="Abnormal-PROCAR.out" ) read(1,'(a28)') t ; write(2,'(a28)') t do i=1,2 read(1,'(a66)') t ; write(2,'(a66)') t do nk=1,30 read(1,'(a )') t ; write(2,'(a )') t read(1,'(a75)') t ; write(2,'(a75)') t do nb=1,100 read(1,'(a )') t ; write(2,'(a )') t read(1,'(a50)') t ; write(2,'(a50)') t read(1,'(a )') t ; write(2,'(a )') t read(1,'(a66,49x,a7)') t,tt ; write(2,'(a66,a7)') t,tt do natom=1,32 read(1,'(a66,49x,a7)') t,tt ; write(2,'(a66,a7)') t,tt end do read(1,'(a66,49x,a7)') t,tt ; write(2,'(a66,a7)') t,tt read(1,'(a )') t ; write(2,'(a )') t read(1,'(a66)') t ; write(2,'(a66)') t do natom=1,32*2 read(1,'(a66)') t ; write(2,'(a66)') t end do end do read(1,'(a )') t ; write(2,'(a )') t end do end do end |
3楼2015-04-09 23:20:51
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1513.7
- 散金: 134
- 红花: 24
- 帖子: 702
- 在线: 326小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II :电子结构

4楼2015-04-10 23:01:10













回复此楼

甲斐之虎