| 查看: 1185 | 回复: 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
» 猜你喜欢
生物学学硕，一志愿湖南大学，初试成绩338
已经有6人回复
-
321求调剂
已经有7人回复
-
343求调剂
已经有4人回复
-
305求调剂
已经有3人回复
-
本科新能源科学与工程，一志愿华理能动285求调剂
已经有6人回复
-
332求调剂
已经有3人回复
-
学硕274求调剂
已经有8人回复
-
压国家一区线，求导师收留，有恩必谢！
已经有7人回复
-
081200-11408-276学硕求调剂
已经有4人回复
-
22408 359分调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1513.7
- 散金: 134
- 红花: 24
- 帖子: 702
- 在线: 326小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2015-04-10 23:01:10
★ ★ ★ ★ ★
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
5楼2015-04-11 19:40:10
