 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 26 | 1/1 | 返回列表 |
| 查看: 4271 | 回復: 25 | |||||
| 【獎勵】 本帖被評價17次，作者qli0218增加金幣 13.25 個 | |||||
[資源]
【分享】DOSCAR數(shù)據(jù)整理小程序
|
|||||
|
我在作DOS圖的時候一點點導數(shù)據(jù)確實很麻煩，于是自己編了個小程序，在我筆記本里的F90編的，和大家分享一下，要是有什么建議盡管提，我在Polish一下。 我算的是一個原子的DOS是分1000行的，取自己感興趣的一組或多組軌道分析，自動將0點調(diào)節(jié)為費米能級。輸出的文件就可以直接導入到Origion里作圖了，我暫時沒用到自旋，所以就沒考慮，也沒計算，兩列里直接取一列了。 希望大家指點 程序： PROGRAM DOSCAR IMPLICIT NONE REAL::R,E,E1 INTEGER::J,I,N,N1,NUM REAL,DIMENSION(0:18)::X OPEN(1,FILE='DOSCAR',STATUS='OLD',ACCESS='SEQUENTIAL',FORM='FORMATTED',ACTION='READ') OPEN(10,FILE='DOSE.DAT',STATUS='REPLACE',ACCESS='SEQUENTIAL',FORM='FORMATTED') READ(1,'(I4)'),J DO I=1,4 READ(1,*) END DO READ(1,'(42X,F11.8)'),E WRITE(*,'("系統(tǒng)中有元素",I2,"個，費米能級為",F9.5)')J,E !WRITE(10,*),J,E DO I=1,1001 READ(1,*) END DO PRINT*,'請輸入要分析的原子序數(shù)：' READ(*,*),E1 PRINT*,'請輸入要分析的軌道：（s軌道（1）；p軌道（2）；d軌道（4）及合軌道（數(shù)字加和））' READ*,N1 N=1 DO NUM=0 DO J=1,1000 READ(1,'(4X,F7.3,18(2X,E10.4E2))'),(X(I),I=0,18) X(0)=X(0)-E IF(N==E1) THEN SELECT CASE(N1) CASE (1) !S WRITE(10,'(1X,F7.3,1(2X,F7.4))')X(0),X(1) CASE (2) !P WRITE(10,'(1X,F7.3,3(2X,F7.4))')X(0),(X(I),I=3,8,2) CASE (3) !S+P WRITE(10,'(1X,F7.3,4(2X,F7.4))')X(0),(X(I),I=1,8,2) CASE (4) !D WRITE(10,'(1X,F7.3,5(2X,F7.4))')X(0),(X(I),I=9,18,2) CASE (5) !S+D WRITE(10,'(1X,F7.3,6(2X,F7.4))')X(0),X(1),(X(I),I=9,18,2) CASE (6) !P+D WRITE(10,'(1X,F7.3,8(2X,F7.4))')X(0),(X(I),I=3,18,2) CASE (7) !S+P+D WRITE(10,'(1X,F7.3,9(2X,F7.4))')X(0),(X(I),I=1,18,2) END SELECT NUM=1 END IF END DO IF(NUM==0) READ(1,*) IF(NUM==1) EXIT N=N+1 END DO CLOSE(1) CLOSE(10) END [ Last edited by mazuju028 on 2010-11-28 at 22:52 ] |
回復此樓» 收錄本帖的淘帖專輯推薦
 vasp | 程序與腳本 | 計算-vasp |
» 猜你喜歡
求助
已經(jīng)有4人回復
-
0805材料320求調(diào)劑
已經(jīng)有3人回復
-
297求調(diào)劑
已經(jīng)有11人回復
-
268求調(diào)劑
已經(jīng)有9人回復
-
一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有10人回復
-
278求調(diào)劑
已經(jīng)有7人回復
-
26考研一志愿中國石油大學(華東)305分求調(diào)劑
已經(jīng)有5人回復
-
280求調(diào)劑
已經(jīng)有9人回復
-
279求調(diào)劑
已經(jīng)有5人回復
-
286求調(diào)劑
已經(jīng)有5人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 博士論文緒論整理處理發(fā)一個綜述不涉及學術(shù)道德吧？
已經(jīng)有4人回復
- band dos 小程序的一些疑問？？？
已經(jīng)有19人回復
- 畢業(yè)論文整理成一篇論文，與畢業(yè)論文重復率較高？會算抄襲嗎？
已經(jīng)有10人回復
- 【活動】Show 腳本，贏獎勵
已經(jīng)有23人回復
- [01-19]精挑細選綠色軟件推薦下載[清清心整理]【已搜索無重復】
已經(jīng)有8人回復
- 【分享】處理DOSCAR文件的腳本和程序
已經(jīng)有19人回復
自己坐沙發(fā)了，歡迎大家提意見，Copy完了別忘評分，嘿嘿 
2樓2009-03-24 07:00:18
3樓2009-03-24 15:07:01
4樓2009-03-24 22:32:52
5樓2009-03-25 08:01:24
6樓2009-03-27 21:22:37
7樓2009-03-29 11:11:31
9樓2009-03-30 16:03:20
10樓2009-03-30 20:08:59
11樓2009-03-31 23:01:32
12樓2009-05-19 20:45:35
guohuazhong
至尊木蟲 (職業(yè)作家)
- 1ST強帖: 1
- 應助: 52 (初中生)
- 貴賓: 0.032
- 金幣: 35131.5
- 帖子: 3924
- 在線: 885.6小時
- 蟲號: 550098
16樓2009-08-13 16:04:39
17樓2010-02-06 06:06:13
18樓2010-02-10 16:27:59
強烈支持~ 19樓2010-11-22 10:43:29
20樓2010-11-28 20:30:10
22樓2010-11-29 10:29:02
23樓2011-05-19 08:40:46
24樓2011-09-21 17:29:41
25樓2012-02-29 15:10:44
簡單回復
wuchenwf8樓
2009-03-29 19:39
回復
dalongmao13樓
2009-05-19 20:48
回復
好
wuli814樓
2009-05-19 21:37
回復
hczheng89915樓
2009-08-13 05:12
回復
leo肖21樓
2010-11-29 10:10
回復
zplou26樓
2015-01-06 13:10
回復
五星好評
相關版塊跳轉(zhuǎn)
我要訂閱樓主 qli0218 的主題更新
| 26 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +5 | 嘉年新程 2026-03-15 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 【考研調(diào)劑】化學專業(yè) 281分，一志愿四川大學，誠心求調(diào)劑 +9 | 吃吃吃才有意義 2026-03-19 | 9/450 |
|
|
|
[考研] 初試 317 +3 | 半拉月丙 2026-03-20 | 3/150 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 材料學學碩080502 337求調(diào)劑-一志愿華中科技大學 +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
|
[考研] 一志愿中國石油大學（華東） 本科齊魯工業(yè)大學 +3 | 石能偉 2026-03-17 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
|
[考研] 一志愿華中科技大學，080502，354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 081700化工學碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 收復試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
|
[考研] 東南大學364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|