24小时热门版块排行榜    

查看: 1117  |  回复: 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的回帖

对马拉二胡

至尊木虫 (职业作家)

引用回帖:
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的回帖
查看全部 8 个回答

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的回帖
信息提示
请填处理意见