将分配的数组从子例程传递到fortran中的主程序；使用还是模块？接口？

问题描述

我想在子程序中分配一个数组，然后在主程序中使用该数组并将其传递给其他子程序。过去（F77？）传递可以在一个通用块中完成，但如今，最受欢迎的过程似乎是使用模块。当我尝试执行此操作时，如代码示例中所示，编译器会告诉我

Rank mismatch in argument ‘f’ at (1) (scalar and rank-1)

显然，主程序认为'f'是一个标量：但是，我读这段代码意味着我已经在子程序内部和主程序中将其声明为一维数组。我想念什么？

我已经尝试了各种变体，例如将变量声明为模块的一部分，但是我想不出什么使编译过程无错误（并且有些产生了更多的错误；-(）。任何见解都值得赞赏。

          module subs
        contains
        subroutine makef(f)
        end subroutine makef
      end module subs
c-----------------------------------------------------------------------
      program work

      use subs
      implicit none
        real,allocatable :: f(:)

      call makef(f)

      write (*,*) f
      stop
      end
c---------------------------------------------------------------------
      subroutine makef(f)
      implicit none

      real,allocatable,intent(out) :: f(:)
      integer :: i
      integer :: is

      is=10
      allocate(f(-is:is))

      do i=-is,is
        f(i)=i
      end do
      return
      end subroutine makef

解决方法

Fortran中的模块与其他语言中的头文件不同，后者仅提供有关在其他地方定义的内容的信息。存在“延迟定义”（子模块）的概念，但是在这种情况下，模块应该说出关于子例程的所有内容，而不仅仅是尝试指出其存在。

在问题的示例中，我们有：主程序；具有模块过程subs的模块makef；外部子例程makef。

主程序使用模块subs及其过程makef，因此在主程序中对makef的引用是该模块过程，而不是外部子例程makef

模块子例程makef的参数f不包含声明语句，从而使其成为隐式声明的标量/外部函数。这是编译器的消息。在模块中使用implicit none，就像在这里的主程序和外部子例程中一样。

子例程的整个定义应放在模块中：

module subs
  implicit none
contains
  subroutine makef(f)
    real,allocatable,intent(out) :: f(:)
    integer :: i
    integer :: is

    is=10
    allocate(f(-is:is))

    do i=-is,is
      f(i)=i
    end do
  end subroutine makef
end module subs

或者，如果确实要引用外部过程的后续实现，则接口模块可以在模块中使用，而无需声明子例程本身。在这种情况下，仍然有必要指定完整的接口：

module subs
  implicit none

! An interface block to give an explicit interface to the external subroutine makef
  interface
     subroutine makef(f)
       implicit none
       real,intent(out) :: f(:)
     end subroutine makef
  end interface
end module subs

在这种情况下，不要使用接口块。

,

您仅将第一行和最后一行的副本放入模块中。那行不通。您必须将整个子例程移入模块。