| 8 | 1/1 | 返回列表 |
| 查看: 1178 | 回復(fù): 7 | ||||
[求助]
求一個(gè)能砍PROCAR的程序,或者編一個(gè)程序砍PROCAR
|
|
如題,我計(jì)算出的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金幣,求大俠出手相助 |

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
隨手寫了個(gè)非常陽春的 fortran 測(cè)過~ 應(yīng)該是沒有問題的 如下 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 金幣就免了~ 大放送 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
討人厭的害羞圖案又出現(xiàn)~~ 再貼一次 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 |
金蟲 (正式寫手)


至尊木蟲 (職業(yè)作家)

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
?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 ^ ??????????@???????????????????????λ??????? |
至尊木蟲 (職業(yè)作家)

| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 322求調(diào)劑 +3 | 舊吢 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +5 | 嗯嗯嗯嗯2 2026-03-27 | 5/250 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過A區(qū)國(guó)家線 +9 | 哈哈157349 2026-03-21 | 9/450 |
|
|
[考研] 269專碩求調(diào)劑 +9 | 金恩貝 2026-03-21 | 9/450 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +3 | 辰沐5211314 2026-03-26 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +4 | XQTJZ 2026-03-26 | 4/200 |
|
|
[考研] 333求調(diào)劑 +7 | 87639 2026-03-21 | 12/600 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +4 | 李軼男003 2026-03-20 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[有機(jī)交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
[考研] 求助 +5 | 夢(mèng)里的無言 2026-03-21 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|