subroutine init_sites_sujun(this)
! order of the lattice sites
! 4 7 10
! 3 6 9
! 1 2 5 8
!
class(sujun) :: this
this%sites(1) = site(1, 4, [2,4,8,10])
this%sites(2) = site(2, 4, [1,3,5,7])
this%sites(3) = site(3, 4, [2,4,6,8])
this%sites(4) = site(4, 4, [1,3,7,9])
this%sites(5) = site(5, 4, [2,6,8,10])
this%sites(6) = site(6, 4, [3,5,7,9])
this%sites(7) = site(7, 4, [2,4,6,10])
this%sites(8) = site(8, 4, [1,3,5,9])
this%sites(9) = site(9, 4, [4,6,8,10])
this%sites(10) = site(10,4, [1,5,7,9])
end subroutine init_sites_sujun