24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1193  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 甲斐之虎 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 300求调剂,材料科学英一数二 +8 leaflight 2026-03-24 8/400 2026-03-29 01:31 by fmesaito
[考研] 348求调剂 +3 小懒虫不懒了 2026-03-28 3/150 2026-03-29 00:39 by 544594351
[考研] 085600,专业课化工原理,321分求调剂 +4 大馋小子 2026-03-28 4/200 2026-03-29 00:32 by 544594351
[考研] 0856求调剂 +6 楒桉 2026-03-28 6/300 2026-03-29 00:31 by 544594351
[考研] 315求调剂 +4 akie... 2026-03-28 5/250 2026-03-28 21:05 by zhq0425
[考研] 083000学硕274求调剂 +8 Li李鱼 2026-03-26 8/400 2026-03-28 20:33 by 加油向未来啊
[考研] 312,生物学求调剂 +3 小译同学abc 2026-03-28 3/150 2026-03-28 15:32 by 落睿可思
[考研] 化学0703 调剂 306分 一志愿211 +4 26要上岸 2026-03-28 4/200 2026-03-28 15:30 by 1018329917
[考研] 求化学调剂 +4 wulanna 2026-03-28 4/200 2026-03-28 13:37 by 唐沐儿
[考研] 085602 307分 求调剂 +7 不知道叫什么! 2026-03-26 7/350 2026-03-28 09:57 by 神马都不懂
[考研] 材料与化工考研调剂 +17 孅華 2026-03-22 17/850 2026-03-28 08:35 by WYUMater
[考研] 291求调剂 +7 孅華 2026-03-22 7/350 2026-03-28 04:02 by fmesaito
[考研] 化学308分求调剂 +8 你好明天你好 2026-03-23 9/450 2026-03-27 14:01 by 杨光于青云
[考研] 求调剂 +6 林之夕 2026-03-24 6/300 2026-03-27 08:38 by hypershenger
[考研] 求调剂 +8 Auroracx 2026-03-22 8/400 2026-03-26 19:55 by 不吃魚的貓
[考研] 085602化学工程求调剂。 +4 平乐乐乐 2026-03-26 4/200 2026-03-26 17:57 by fmesaito
[考研] 总分293求调剂 +6 加一一九 2026-03-25 8/400 2026-03-26 13:30 by yujianx
[考研] 机械学硕总分317求调剂!!!! +4 Acaciad 2026-03-25 4/200 2026-03-25 19:59 by hanserlol
[考研] 282求调剂 +3 wcq131415 2026-03-24 3/150 2026-03-25 12:16 by userper
[考研] 求调剂一志愿武汉理工大学材料工程(085601) +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
信息提示
请填处理意见