site
subroutine set_r_vec(this, r_vec) class(site) :: this integer, intent(in) :: r_vec(3) this%r_vec = r_vec end subroutine set_r_vec