| 5 | 1/1 | 返回列表 |
| 查看: 1180 | 回復(fù): 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金幣,求大俠出手相助 |

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗: +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 ^ ??????????@???????????????????????λ??????? |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗: +20 |
|
隨手寫了個非常陽春的 fortran 測過~ 應(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)驗: +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 |
金蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +3 | 丟掉懶惰 2026-03-27 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 13/650 |
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +5 | Li李魚 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +14 | 荷包蛋hyj 2026-03-26 | 14/700 |
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
[考研]
|
孅華 2026-03-22 | 16/800 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +8 | 北道巷 2026-03-24 | 8/400 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
[考研] 機械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當牛 2026-03-24 | 4/200 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運來來來 2026-03-22 | 5/250 |
|