| 查看: 343 | 回复: 2 | ||
[求助]
算一种材料的动力学模拟，老师给我发的一个程序一直不能用，忘高人能够指点一二 已有1人参与
|
|
错误如这图片，请高人指点发生什么事情了 程序如下 program main implicit none character*15 :: charac integer :: i,j,k,status,temp,step,nband_shift,flagd integer :: read_bands_total,read_bands_h,read_bands_e real(8) :: Efermi integer, allocatable :: itera_steps(  real(8), allocatable :: energy(:, open(unit=99,file='ReadOUTCAR.inp') read(99,*) step ! Number of MD steps read(99,*) nband_shift ! nband_shift+1 is the first band which is output read(99,*) read_bands_h ! Number of hole energy bands to output read(99,*) read_bands_e ! Number of electron energy bands to output read(99,*) Efermi ! Fermi energy close(99) read_bands_total=read_bands_h+read_bands_e allocate(itera_steps(step),energy(read_bands_total,step)) open(18,file='iteration-steps.dat') do i=1,step read(18,*) flagd,itera_steps(i) end do close(18) open(32,file='OUTCAR') read(32,*,IOSTAT=status) charac j=1 k=0 do while(status==0) charac=adjustl(charac) if(index(charac,'E-fermi')==1) then k=k+1 if(k>itera_steps(j)) then write(*,*) 'error!!' stop else if(k==itera_steps(j)) then do i=1,4+nband_shift read(32,*) end do do k=1,read_bands_total read(32,*) temp, energy(k,j) end do j=j+1 k=0 end if end if if (j>step) exit read(32,*,IOSTAT=status) charac end do close(32) open(32,FILE='hole-result.dat') do i=1,step do j=1,read_bands_h write(32,'(F12.6)',advance='no') energy(j,i)+Efermi end do write(32,*) end do close(32) open(32,FILE='elec-result.dat') do i=1,step do j=read_bands_h+1,read_bands_total write(32,'(F12.6)',advance='no') energy(j,i)+Efermi end do write(32,*) end do close(32) deallocate(itera_steps,energy) end program  4Q{N`3WHZS9B}EU(~1ZXH4T.png |
回复此楼» 猜你喜欢
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
如何用Image J 软件计算发泡材料的泡孔尺寸和泡孔密度？跪求指导！！！
已经有16人回复
- 材料计算经典：从基本原理到材料性质的《计算材料学》（英文）
已经有503人回复
- 从传统到低维——分享关于Bi2Te3体系热电性能模拟计算的的文章
已经有69人回复
- 在FLUENT的材料设置中，空气设置如果设为理想气体，那么大家觉得粘性设置还有用么？
已经有21人回复
- 材料模拟与计算 VASP软件培训班，好资源跟大家分享！
已经有14人回复
- 如何模拟材料加工过程(例如，铸造，锻造和焊接)中金属材料内部组织结构的变化？
已经有9人回复
- 随着温度变化的材料属性值怎么算
已经有8人回复
- Abaqus软件可以算孔隙介质材料的受力问题吗？如果可以，孔隙有液体作用时还可以算吗？
已经有10人回复
- 请问大面积复制自己上一个发表的文章的前言、材料方法的文字算是抄袭吗
已经有17人回复
- 测量CO2在材料中吸附量有一个单位是cm3(STP)/g。这个单位的计算应该采用那一种算法?
已经有7人回复
- 用VASP做出了一种材料的声子谱，想计算出它的热导率这个要如何从中提取信息呢？
已经有12人回复
- [原创]一种基于能带计算结果的载流子有效质量计算方法
已经有174人回复
- 我们用VASP计算材料的力学性能，是绝对零度的。有意义吗？正解者有红花相送哦
已经有8人回复
- 当体系中包含两种材料时，如何用VASP计算其中一种材料的功函？
已经有11人回复
weiwei43
至尊木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 18316.1
- 红花: 6
- 帖子: 737
- 在线: 94.7小时
- 虫号: 250282
- 注册: 2006-05-13
- 专业: 机械动力学
2楼2015-06-11 09:14:38
weiwei43
至尊木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 18316.1
- 红花: 6
- 帖子: 737
- 在线: 94.7小时
- 虫号: 250282
- 注册: 2006-05-13
- 专业: 机械动力学
3楼2015-06-11 09:17:22
