| 查看: 1349 | 回复: 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
» 猜你喜欢
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有19人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
★ ★ ★ ★ ★
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.5小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2015-04-10 23:01:10
5楼2015-04-11 19:40:10
对马拉二胡
至尊木虫 (职业作家)
- 应助: 27 (小学生)
- 金币: 27355.6
- 散金: 244
- 红花: 13
- 帖子: 4251
- 在线: 943.2小时
- 虫号: 2471734
- 注册: 2013-05-19
- 性别: GG
- 专业: 无机非金属材料
6楼2015-04-12 10:54:33
|
?m??????ゃ???@?N???????????N ???^~ ?@??????????????? ?@????]???κ??A?????????? do loop, read, and write ????????????????????? read(1,'(a28)') t ; write(2,'(a28)') t ???? read(1,'(a66)') t ; write(2,'(a66)') t ???????????? a28 ( a66 ) ?@????????????? ??? input file (PROCAR) ?????? ion s py pz px dxy dyz dz2 dxz dx2 f-3 f-2 f-1 f0 f1 f2 f3 tot ^ ?@??? 66 ?????? format ?? a66 ??? ion s py pz px dxy dyz dz2 dxz dx2 f-3 f-2 f-1 f0 f1 f2 f3 tot ^ ??????????@???????????????????????λ??????? |
7楼2015-04-12 18:31:48
对马拉二胡
至尊木虫 (职业作家)
- 应助: 27 (小学生)
- 金币: 27355.6
- 散金: 244
- 红花: 13
- 帖子: 4251
- 在线: 943.2小时
- 虫号: 2471734
- 注册: 2013-05-19
- 性别: GG
- 专业: 无机非金属材料
8楼2015-04-12 21:24:14
