print_basis Subroutine

    subroutine print_basis(rep)
        type(core_space_t), intent(in) :: rep

        integer :: iunit, i, nI(nel)

        iunit = get_free_unit()
        open(iunit, file = 'semistoch-basis', status = 'replace')

        do i = 1, rep%determ_space_size
            call decode_bit_det(nI, rep%core_space(:, i))
            call write_det(iunit, nI, .true.)
        end do
        close(iunit)

    end subroutine print_basis