24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1200  |  回复: 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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

别人笑我太疯颠,我笑他人看不穿。不见五陵豪杰墓,无花无酒锄作田。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

★ ★ ★ ★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

討人厭的害羞圖案又出現~~ 再貼一次

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)

引用回帖:
3楼: Originally posted by jpchou at 2015-04-09 23:20:51
討人厭的害羞圖案又出現~~ 再貼一次

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 ;  ...

看似很NB的样子。
一定要精通MaterialsStudios,VASP.
4楼2015-04-10 23:01:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

甲斐之虎

银虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by jpchou at 2015-04-09 23:19:33
隨手寫了個非常陽春的 fortran
測過~ 應該是沒有問題的
如下

implicit real*8(a-h,o-z)
character*128 t,tt

open(1,file="Abnormal-PROCAR"
open(2,file="Abnormal-PROCAR.out"
...

周前輩幫忙解決大問題了!
晚生感謝周前輩一百遍,一百遍!
别人笑我太疯颠,我笑他人看不穿。不见五陵豪杰墓,无花无酒锄作田。
5楼2015-04-11 19:40:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

对马拉二胡

至尊木虫 (职业作家)

引用回帖:
3楼: Originally posted by jpchou at 2015-04-09 23:20:51
討人厭的害羞圖案又出現~~ 再貼一次

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 ;  ...

大神,如果PROCAR中只需要留下 s     py     pz     px 的数据,程序该怎么呢?万分感谢
相互学习,共同进步。
6楼2015-04-12 10:54:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

???????:
6?: Originally posted by ?????????? at 2015-04-12 10:54:33
???????PROCAR?????????? s     py     pz     px ??????????????????????л...

?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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

对马拉二胡

至尊木虫 (职业作家)

???????:
7?: Originally posted by jpchou at 2015-04-12 18:31:48
?m??????ゃ???@?N???????????N
???^~ ?@???????????????

?@????]???κ??A?????????? do loop, read, and write
?????????????????????
read(1,'(a28)') t ; write(2,'(a28)') t
????
read ...

??л?????????
相互学习,共同进步。
8楼2015-04-12 21:24:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 甲斐之虎 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 求调剂 有机化学考研356分 +3 Nadiums 2026-03-25 4/200 2026-03-29 00:50 by 544594351
[考研] 求调剂 +7 争取九点睡 2026-03-28 8/400 2026-03-28 21:07 by 争取九点睡
[考研] 本科新能源科学与工程,一志愿华理能动285求调剂 +7 AZMK 2026-03-28 11/550 2026-03-28 21:01 by xxxsssccc
[考研] 化学0703 调剂 306分 一志愿211 +4 26要上岸 2026-03-28 4/200 2026-03-28 15:30 by 1018329917
[考研] 346求调剂 一志愿070303有机化学 +3 萝卜炖青菜 2026-03-28 3/150 2026-03-28 14:11 by 唐沐儿
[考研] 0856求调剂 +11 zhn03 2026-03-25 12/600 2026-03-28 13:32 by 唐沐儿
[考研] 调剂 +3 好好读书。 2026-03-28 3/150 2026-03-28 12:04 by 王保杰33
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +9 丹青奶盖 2026-03-26 10/500 2026-03-28 07:45 by barnett0632
[考研] 328求调剂 +7 嗯滴的基本都 2026-03-27 7/350 2026-03-28 04:19 by fmesaito
[考研] 285求调剂 +4 AZMK 2026-03-27 7/350 2026-03-27 20:59 by AZMK
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 12/600 2026-03-27 14:43 by shangxh
[考研] 一志愿华东理工大学081700,初试分数271 +6 kotoko_ik 2026-03-23 7/350 2026-03-27 12:29 by 惠州彭于晏
[考研] 一志愿吉大071010,316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
[考研] 359求调剂 +4 王了个楠 2026-03-25 4/200 2026-03-27 08:43 by 不吃魚的貓
[考研] 085602 289分求调剂 +8 WWW西西弗斯 2026-03-24 8/400 2026-03-26 16:33 by 不吃魚的貓
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 10/500 2026-03-26 13:11 by 公瑾逍遥
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 各位老师您好:本人初试372分 +5 jj涌77 2026-03-25 6/300 2026-03-25 14:15 by mapenggao
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
信息提示
请填处理意见