将数组分配逻辑移至子例程的Fortran程序改造及编译报错问题

首先，咱们先解决你遇到的编译报错问题：gfortran fail with a Type mismatch in argument ‘cc’: REAL(4) to ...，这个错误非常直接——你的主程序里传递的rr是real类型的可分配数组，但子例程sub1里声明的cc是complex类型的数组，两者类型完全不匹配，编译器自然会报错。这是第一个要修正的点。

接下来咱们说核心需求：把数组的分配逻辑移到sub1里，让子例程根据实际需要的大小（比如你代码里的min_cc_sz=6）来动态分配数组，而不是主程序里硬编码rr_sz=10。要实现这个功能，有两个关键要点必须处理：

1. 必须提供子例程的显式接口 #

在Fortran中，如果子例程要操作可分配数组参数（比如在子例程里分配/释放数组），调用它的主程序必须拥有该子例程的显式接口。最规范也最省心的实现方式是用模块来封装子例程，这样编译器会自动生成显式接口，不需要手动编写接口块。

2. 统一数组类型并调整参数声明 #

既然子例程sub1里是用complex(i,i)给数组赋值，那主程序里的数组也应该改成complex类型（如果确实需要real类型，可以在赋值时做类型转换，但原逻辑里直接赋值会有问题）。同时，子例程的数组参数要声明为可分配，并指定intent属性（比如intent(out)表示子例程会分配并填充这个数组）。

改造后的完整代码示例 #

! 定义模块，封装子例程，自动提供显式接口
module subroutines_mod
    implicit none

contains

    subroutine sub1(cc)
        implicit none
        complex, dimension(:), allocatable, intent(out) :: cc
        integer, parameter :: min_cc_sz = 6
        integer :: i

        ! 先检查数组是否已分配，若已分配则释放
        if (allocated(cc)) deallocate(cc)
        ! 根据实际需要的大小分配数组
        allocate(cc(min_cc_sz))

        ! 给数组赋值
        do i = 1, min_cc_sz
            cc(i) = cmplx(i, i)  ! 用cmplx函数更规范，替代complex(i,i)
        end do
    end subroutine sub1

end module subroutines_mod

program main
    use subroutines_mod  ! 引用模块，获取sub1的显式接口
    implicit none
    complex, dimension(:), allocatable :: rr  ! 统一为complex类型

    ! 调用子例程，子例程会负责分配rr数组
    call sub1(rr)

    ! 打印数组信息和内容
    print *, "Array size: ", size(rr)
    print *, "Array content: ", rr

    ! 释放数组（也可以不写，程序结束时会自动释放，但显式释放更规范）
    if (allocated(rr)) deallocate(rr)
end program main

关键改动说明 #

• 模块封装：把sub1放到模块里，主程序通过use语句引用模块，自动获得显式接口，解决了可分配数组参数的接口问题。
• 参数调整：sub1的参数cc声明为complex, dimension(:), allocatable, intent(out)，这样子例程可以直接分配这个数组。
• 类型统一：主程序的rr改成complex类型，和子例程里的赋值逻辑匹配，消除了类型不匹配的报错。
• 动态分配：子例程里根据min_cc_sz分配数组，不再依赖主程序里的硬编码大小，实现了你想要的动态调整功能。

替代方案：不用模块的情况 #

如果你不想用模块，也可以把sub1作为内部子例程放在主程序的contains块里，这样同样会自动生成显式接口：

program main
    implicit none
    complex, dimension(:), allocatable :: rr

    call sub1(rr)

    print *, "Array size: ", size(rr)
    print *, "Array content: ", rr

    if (allocated(rr)) deallocate(rr)

contains  ! 内部子例程，自动拥有显式接口
    subroutine sub1(cc)
        implicit none
        complex, dimension(:), allocatable, intent(out) :: cc
        integer, parameter :: min_cc_sz = 6
        integer :: i

        if (allocated(cc)) deallocate(cc)
        allocate(cc(min_cc_sz))

        do i = 1, min_cc_sz
            cc(i) = cmplx(i, i)
        end do
    end subroutine sub1
end program main

为什么必须用显式接口？ #

当子例程的参数涉及可分配数组、指针、可选参数、带有intent属性的数组等情况时，编译器需要知道参数的完整类型和属性信息才能正确生成调用代码。如果没有显式接口，编译器无法验证参数是否匹配，会导致报错或者运行时的未定义行为。

备注：内容来源于stack exchange，提问作者Fravadona