[交流] MS能带处理方法程序

        program msbs
        implicit none
        integer::i,j,k,n,l
        real*8,dimension(:,:),allocatable::d
        character(len=1)::a
        character(len=30)::a1
        write(*,'(15x,"该软件为计算处理ms band structure小程序")')
        write(*,'("使用方法：")')
        write(*,'("1、在ms中画出band structure；")')
        write(*,'("2、将band structure图export为.csv格式文件；")')
        write(*,'("3、将得到的文件用exel打开，设置数据格式为15位小数；")')
        write(*,'("4、再将此文件另存为.txt文件。")')
        write(*,'("于是即可利用此程序处理band structure数据")')
        write(*,'("得到dat文件可直接导入到origin中作图")')
        write(*,'("注意：文件名由数字和字母组成，不区分大小写，最长30
        1字节")')
        write(*,*)
        write(*,'("--------------------")')
        write(*,*)
        write(*,'("请输入分析的文件名：（不用写后缀）")')
        read*,a1
        open(1,file=trim(a1)//'.txt',status='old',access='sequential',
        1form='formatted',action='read')
        open(2,file=trim(a1)//'.dat',status='replace',access='sequential',
        1form='formatted')
        i=0
        n=0
        do
        read(1,'(a1)',iostat=l)a
        if (l/=0) exit
        if (a=="1") n=n+1
        i=i+1
        end do
        allocate(d(i/n-1,n+1))
        rewind(1)
        do j=1,i/n-1
        read(1,'(f17.15,f21.15)')d(j,1),d(j,2)
        end do
        k=3
        read(1,*)       
        do j=i/n+1,i
        if (mod(j,i/n)==0) then
        k=k+1
        read(1,*)       
        else
        read(1,'(a30)')a1
        backspace(1)
        if (a1(20:20)=="-") then
        read(1,'(17x,f21.15)')d(mod(j,i/n),k)
        else
        read(1,'(17x,f20.15)')d(mod(j,i/n),k)
        end if
        end if
        end do
        do j=2,n+1
        if (maxval(d(:,j))*minval(d(:,j+1))<=0) then
        print*,abs(minval(d(:,j+1))-maxval(d(:,j)))
        exit
        end if
        end do
        k=j
        do j=1,i/n-1
        write(2,'(2x,11f9.4)')d(j,1),d(j,k-4:k+5)
        end do
        deallocate(d)
        close(1)
        close(2)
        end program msbs

2楼: Originally posted by 小红豆 at 2018-12-18 22:17:27
费米能级只有一个办

3楼: originally posted by qli0218 at 2018-12-20 11:22:19
只有一个“办”是什么意思？我就想知道ms里面画出来的能带图是不是已经把费米能级移动到了0 ev。谢谢！...