| 查看: 4263 | 回复: 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 |
» 猜你喜欢
0856材料求调剂
已经有11人回复
-
298求调剂
已经有7人回复
-
311求调剂
已经有7人回复
-
材料学调剂
已经有4人回复
-
面上模板改不了页边距吧？
已经有3人回复
-
高分子化学与物理调剂
已经有5人回复
-
版面费该交吗
已经有18人回复
-
285求调剂
已经有5人回复
-
材料类求调剂
已经有4人回复
-
材料调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 博士论文绪论整理处理发一个综述不涉及学术道德吧？
已经有4人回复
- band dos 小程序的一些疑问？？？
已经有19人回复
- 毕业论文整理成一篇论文，与毕业论文重复率较高？会算抄袭吗？
已经有10人回复
- 【活动】Show 脚本，赢奖励
已经有23人回复
- [01-19]精挑细选绿色软件推荐下载[清清心整理]【已搜索无重复】
已经有8人回复
- 【分享】处理DOSCAR文件的脚本和程序
已经有19人回复
3楼2009-03-24 15:07:01
自己坐沙发了，欢迎大家提意见，Copy完了别忘评分，嘿嘿 2楼2009-03-24 07:00:18
6楼2009-03-27 21:22:37
12楼2009-05-19 20:45:35
