| 查看: 4297 | 回复: 25 | |||||
| 【奖励】 本帖被评价17次，作者qli0218增加金币 13.25 个 | |||||
[资源]
【分享】DOSCAR数据整理小程序
|
|||||
|
我在作DOS图的时候一点点导数据确实很麻烦，于是自己编了个小程序，在我笔记本里的F90编的，和大家分享一下，要是有什么建议尽管提，我在Polish一下。 我算的是一个原子的DOS是分1000行的，取自己感兴趣的一组或多组轨道分析，自动将0点调节为费米能级。输出的文件就可以直接导入到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(*,'("系统中有元素",I2,"个，费米能级为",F9.5)')J,E !WRITE(10,*),J,E DO I=1,1001 READ(1,*) END DO PRINT*,'请输入要分析的原子序数：' READ(*,*),E1 PRINT*,'请输入要分析的轨道：（s轨道（1）；p轨道（2）；d轨道（4）及合轨道（数字加和））' 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 |
» 猜你喜欢
291分调剂
已经有9人回复
-
调剂求收留
已经有34人回复
-
291 求调剂
已经有38人回复
-
22408 312求调剂
已经有17人回复
-
一志愿华中农业071010，320求调剂
已经有6人回复
-
290调剂生物0860
已经有41人回复
-
291求调剂
已经有3人回复
-
211本科材料化工求调剂
已经有23人回复
-
山东省基金2026
已经有9人回复
-
药学求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 博士论文绪论整理处理发一个综述不涉及学术道德吧？
已经有4人回复
- band dos 小程序的一些疑问？？？
已经有19人回复
- 毕业论文整理成一篇论文，与毕业论文重复率较高？会算抄袭吗？
已经有10人回复
- 【活动】Show 脚本，赢奖励
已经有23人回复
- [01-19]精挑细选绿色软件推荐下载[清清心整理]【已搜索无重复】
已经有8人回复
- 【分享】处理DOSCAR文件的脚本和程序
已经有19人回复
自己坐沙发了，欢迎大家提意见，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
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
回复
五星好评
20