sitesym_replace_d_matrix_band Subroutine

  subroutine sitesym_replace_d_matrix_band()
    !==================================================================!
    use w90_parameters, only: num_wann

    implicit none

    !write(stdout,"(a)") '-- sitesym_replace_d_matrix_band --'
    !write(stdout,"(a)") 'd_matrix_band is replaced by d_matrix_wann'
    deallocate (d_matrix_band)
    allocate (d_matrix_band(num_wann, num_wann, nsymmetry, nkptirr))
    d_matrix_band = d_matrix_wann

    return
  end subroutine sitesym_replace_d_matrix_band