print_determ_vec Subroutine

    subroutine print_determ_vec(run)
        integer, intent(in), optional :: run

        integer :: run_, iunit, i
        def_default(run_, run, 1)

        iunit = get_free_unit()
        open(iunit, file = 'determ_vecs', status = 'replace')
        associate ( rep => cs_replicas(run_))
            do i = 1, rep%determ_space_size
                write(iunit, *) rep%full_determ_vecs(1, i)
            end do
        end associate
        close(iunit)
    end subroutine print_determ_vec