program test0202
      use convertharmonics
      implicit none

      integer(i16b) :: indet(2), detplus, detminus
      integer(i16b), allocatable :: zbasislist(:, :)
      real(rk), allocatable :: zcoefs(:, :)
      real(rk) :: zoefs
      integer :: pos
      !allocate(zbasislist(10, 2))
      !allocate(zcoefs(10, 2))

      !pos =13
      !write(*, *) sign_m(pos, detplus, detminus)
      !sign_m tested all correct
     
      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      !Need to change neact in convertharmonics
      !1e tests
      !indet(1:2) = (/21_i16b, 0_i16b/)!Correct
      !indet(1:2) = (/0_i16b, 21_i16b/)!Correct
      !indet(1:2) = (/0_i16b, 512_i16b/)!Correct for 1 e from 1 - 512
      !indet(1:2) = (/512_i16b, 0_i16b/)!Correct

      !2e tests
      !indet(1:2) = (/3_i16b, 0_i16b/)!Correct all m=0  
      !indet(1:2) = (/0_i16b, 3_i16b/)! Correct
      indet(1:2) = (/0_i16b, 24_i16b/)
      !indet(1:2) = (/0_i16b, 20_i16b/)
      call Y2Z_det(indet, 1._rk,  zbasislist, zcoefs, 0)

      deallocate(zbasislist)
      deallocate(zcoefs)

end program test0202
