应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 热点前沿问题讨论 » 求一个能砍PROCAR的程序，或者编一个程序砍PROCAR
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1187  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

甲斐之虎

银虫 (正式写手)

[求助] 求一个能砍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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
别人笑我太疯颠，我笑他人看不穿。不见五陵豪杰墓，无花无酒锄作田。
1楼 2015-04-08 21:36:01
已阅   回复此楼   关注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的回帖
查看全部 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

金幣就免了~ 大放送
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新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
赞一下(1人) 回复此楼
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的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学学硕，一志愿湖南大学，初试成绩338 +5 YYYYYNNNNN 2026-03-26 6/300 2026-03-28 17:50 by 狂炫麦当当
[考研] 305求调剂 +3 RuiFairyrui 2026-03-28 3/150 2026-03-28 17:41 by 王小欠i
[考研] 一志愿华北电力大学能动专硕，293，求调剂 +3 15537177284 2026-03-23 5/250 2026-03-28 16:11 by xxxsssccc
[考研] 283求调剂 +3 A child 2026-03-28 3/150 2026-03-28 15:41 by ms629
[考研] 347求调剂 +3 山顶见α 2026-03-25 3/150 2026-03-28 14:13 by 唐沐儿
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
[考研] 286求调剂 +12 PolarBear11 2026-03-26 12/600 2026-03-28 12:14 by zllcz
[考研] 292求调剂 +14 鹅鹅鹅额额额额� 2026-03-25 15/750 2026-03-28 08:45 by WYUMater
[考研] 266求调剂 +11 阳阳哇塞 2026-03-27 12/600 2026-03-27 17:56 by yu221
[考研] 085600材料与化工调剂 +10 A-哆啦Z梦 2026-03-23 16/800 2026-03-27 15:13 by caszguilin
[考研] 308求调剂 +7 墨墨漠 2026-03-25 7/350 2026-03-27 14:47 by 狂炫麦当当
[考研] 一志愿211，335分，0856，求调剂院校和导师 +4 倾____萧 2026-03-27 5/250 2026-03-27 11:52 by zhshch
[考研] 286求调剂 +4 lim0922 2026-03-26 4/200 2026-03-27 10:28 by guoweigw
[考研] 求调剂，一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +4 @taotao 2026-03-26 5/250 2026-03-27 08:10 by hypershenger
[考研] 271求调剂 +6 生如夏花… 2026-03-22 6/300 2026-03-26 16:48 by 张凯十八号
[考研] 一志愿武理085500机械专业总分300求调剂 +3 an10101 2026-03-24 7/350 2026-03-25 00:00 by 山鬼0-
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗？ +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 一志愿国科过程所081700，274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 材料/农业专业，07/08开头均可，过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪！？！
信息提示
关闭
请填处理意见
关闭