24小时热门版块排行榜    

查看: 1365  |  回复: 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

专家顾问 (著名写手)

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

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的回帖
查看全部 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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 17:49 by eObJPa4gWmp6
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 9/450 2026-09-06 17:48 by myrtle
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见