site
subroutine set_k_vec(this, k_vec) class(site) :: this integer, intent(in) :: k_vec(3) this%k_vec = k_vec end subroutine set_k_vec