print_hamiltonian Subroutine

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

        HElement_t(dp), allocatable :: full_H(:, :)
        integer :: iunit

        call calc_determin_hamil_full(full_H, rep)
        iunit = get_free_unit()
        open(iunit, file = 'semistoch-hamil', status = 'replace')
        call print_matrix(full_H, iunit)
        close(iunit)

    end subroutine print_hamiltonian