| 8 | 1/1 | 返回列表 |
| 查看: 1176 | 回復(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è)非常陽(yáng)春的 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ā)表 | |
|---|---|---|---|---|
|
[考研] 324求調(diào)劑 +4 | wysyjs25 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績(jī)338 +4 | YYYYYNNNNN 2026-03-26 | 4/200 |
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
[考研] 292求調(diào)劑 +9 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 10/500 |
|
|
[考研] 考研調(diào)劑 +8 | 小蠟新筆 2026-03-26 | 8/400 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 3/150 |
|
|
[考研] 311求調(diào)劑 +4 | 勇敢的小吳 2026-03-20 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +3 | 你好明天你好 2026-03-23 | 3/150 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|