plot Module Reference

Functions/Subroutines
subroutine	band_plot ()
	Output band dispersion. More...
subroutine	fermi_plot ()
	Output FermiSurfer file. More...

Function/Subroutine Documentation

band_plot()

subroutine plot::band_plot

Output band dispersion.

Definition at line 32 of file plot.F90.

    !
    use constant, only : htr2ev
    use k_point, only : kvec, nk
    use kohn_sham, only : nbnd, eval, ef
    !
    integer :: fo = 20, ibnd, ik
    real(8) :: xk(nk), dkvec(3), dxk
    !
    write(*,*) "  Output band.dat"
    open(fo, file="band.dat")
    !
    xk(1) = 0.0d0
    do ik = 2, nk
       dkvec(1:3) = kvec(1:3,ik) - kvec(1:3,ik-1)
       dxk = sqrt(dot_product(dkvec(1:3), dkvec(1:3)))
       xk(ik) = xk(ik-1) + dxk
    end do
    !
    write(fo,*) "# k energy[eV]" 
    do ibnd = 1, nbnd
       do ik = 1, nk
          write(fo,*) xk(ik), (eval(ibnd, ik)-ef)*htr2ev
       end do
       write(fo,*) ""
    end do
    !
    close(fo)
    !

References kohn_sham::ef, kohn_sham::eval, constant::htr2ev, k_point::kvec, kohn_sham::nbnd, and k_point::nk.

Referenced by pwdft().

Here is the caller graph for this function:

fermi_plot()

subroutine plot::fermi_plot

Output FermiSurfer file.

Definition at line 65 of file plot.F90.

    !
    use atm_spec, only : bvec
    use k_point, only : kgrd, nk, kvec
    use kohn_sham, only : eval, nbnd, evec, ef
    use gvec, only : g_wf
    use constant, only : htr2ev
    !
    integer :: ik, ibnd, fo = 20, ii, ipw
    real(8) :: kgv(g_wf%npw,3), avf(nbnd,nk), vf(3)
    !
    do ik = 1, nk
       do ipw = 1, g_wf%npw
          kgv(ipw,1:3) = kvec(1:3,ik) + matmul(bvec(1:3,1:3), dble(g_wf%mill(1:3,g_wf%map(ipw))))
       end do
       do ibnd = 1, nbnd
          do ii = 1, 3
             vf(ii) = dble(sum(conjg(evec(1:g_wf%npw,ibnd,ik)) &
             &                     * evec(1:g_wf%npw,ibnd,ik)  &
             &                      * kgv(1:g_wf%npw,ii     )  ))
          end do
          avf(ibnd,ik) = sqrt(dot_product(vf(1:3),vf(1:3)))
       end do
    end do
    !
    write(*,*) "  Output vf.frmsf"
    open(fo, file = "vf.frmsf")
    write(fo,*) kgrd(3), kgrd(2), kgrd(1)
    write(fo,*) 1
    write(fo,*) nbnd
    write(fo,*) real(bvec(1:3,3))
    write(fo,*) real(bvec(1:3,2))
    write(fo,*) real(bvec(1:3,1))
    do ibnd = 1, nbnd
       do ik = 1, nk
          write(fo,*) real((eval(ibnd,ik) - ef)*htr2ev)
       end do
    end do
    do ibnd = 1, nbnd
       do ik = 1, nk
          write(fo,*) real(avf(ibnd,ik))
       end do
    end do
    close(fo)    
    !

References atm_spec::bvec, kohn_sham::ef, kohn_sham::eval, kohn_sham::evec, gvec::g_wf, constant::htr2ev, k_point::kgrd, k_point::kvec, kohn_sham::nbnd, and k_point::nk.

Referenced by pwdft().

Here is the caller graph for this function: