【獎勵】本帖被評價17次，作者qli0218增加金幣 13.25 個

qli0218

木蟲 (小有名氣)

[資源] 【分享】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 ]

回復此樓

1樓 2009-03-23 14:44:00

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

qli0218

木蟲 (小有名氣)

自己坐沙發(fā)了，歡迎大家提意見，Copy完了別忘評分，嘿嘿

贊一下

回復此樓

2樓2009-03-24 07:00:18

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hdh912

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1180.9
帖子: 842
在線: 72.9小時
蟲號: 387166

★★★ 三星級,支持鼓勵

雖然不用VASP 鼓勵感謝下

贊一下

回復此樓

3樓2009-03-24 15:07:01

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

muzi1984

木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 12083.2
帖子: 226
在線: 111.6小時
蟲號: 600181

支持一下

回復此樓

4樓2009-03-24 22:32:52

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

veryman

木蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 6793.4
帖子: 479
在線: 206.6小時
蟲號: 77254

★★★★★ 五星級,優(yōu)秀推薦

強烈支持~

回復此樓

5樓2009-03-25 08:01:24

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hakuna

木蟲 (知名作家)

★★★★★ 五星級,優(yōu)秀推薦

*****鼓勵原創(chuàng)*****

贊一下

回復此樓

6樓2009-03-27 21:22:37

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cbFeng200881

金蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 262.8
帖子: 2772
在線: 699.4小時
蟲號: 524662

★★★ 三星級,支持鼓勵

split_dos...

回復此樓

7樓2009-03-29 11:11:31

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zdhlover

榮譽版主 (職業(yè)作家)

★★★★★ 五星級,優(yōu)秀推薦

支持原創(chuàng)

回復此樓

9樓2009-03-30 16:03:20

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zsf888

銅蟲 (小有名氣)

好好好

回復此樓

10樓2009-03-30 20:08:59

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hpli517

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 12312.5
帖子: 158
在線: 117.6小時
蟲號: 600347

強烈支持

回復此樓

11樓2009-03-31 23:01:32

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

405963512

銅蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1127.1
帖子: 301
在線: 9.9小時
蟲號: 726885

你好

引用回帖:

Originally posted by qli0218 at 2009-3-23 14:44:
我在作DOS圖的時候一點點導數(shù)據(jù)確實很麻煩，于是自己編了個小程序，在我筆記本里的F90編的，和大家分享一下，要是有什么建議盡管提，我在Polish一下。
我算的是一個原子的DOS是分1000行的，取自己感興趣的一組或 ...

處理dos確實很麻煩，請問如何用呢？謝謝！

贊一下

回復此樓

12樓2009-05-19 20:45:35

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖