在两级子例程中传递假定形状的数组 (Fortran 90)

Question

我在 Fortran 90 中调用具有假定形状数组的连续子例程时遇到了问题。更具体地说，我调用了两个级别的子例程，将假定形状数组作为参数传递，但最终该数组丢失了。为了演示它，可以按照下面的代码。

  program main

  INTERFACE
     subroutine sub1(x)
       real, dimension(:):: x
       real C
     end subroutine sub1

     subroutine sub2(x)
       real, dimension(:):: x
       real C
     end subroutine sub2
  END INTERFACE

  real, dimension(:), allocatable:: x

  allocate(x(1:10)) ! First executable command in main
  x(1) = 5.
  call sub1(x)
  write(*,*) 'result = ',x(1)
  deallocate(x)
  end program main

  subroutine sub1(x) ! The first subroutine
  real, dimension(:):: x
  real C
  call sub2(x)
  end subroutine sub1

  subroutine sub2(x) ! The second subroutine
  real, dimension(:):: x
  real C
  C2=x(1)
  end subroutine sub2

很快，main 分配 x 然后调用 sub1(x)。然后 sub1 调用 sub2(x)。这意味着分配的数组被传递给一个子程序，该子程序将它传递给另一个子程序。我希望 sub2 中的数组与我在 main 中创建的数组相同，但没有。使用 gdb 作为工具来探索它，我得到了这个：

1) 在 main 中，就在调用 sub1 之前，数组 x 被完美定义：

(gdb) p x
$1 = (5, 0, 0, 0, 0, 0, 0, 0, 0, 0)

2) 在 sub1 中，就在调用 sub2 之前，x 也被很好地定义：

(gdb) p x
$2 = (5, 0, 0, 0, 0, 0, 0, 0, 0, 0)

3) 然而，在 sub2 内部，x 有一个意想不到的值，甚至它的维度也是绝对错误的：

(gdb) p x
$3 = ()
(gdb) 什么是 x
type = REAL(4) (0:-1)

因此，x 已成功从 main 传递到 sub1，但没有从 sub1 传递到 sub2。我一直在使用 Intel Fortran 和 gfortran，结果相同。

我已经为此苦苦挣扎了很长时间。任何帮助将不胜感激。
G.奥利维拉。

Answer 1

使用假定形状的虚拟参数需要显式接口。 在您的主程序中，您为两个子例程提供了显式接口，但这些接口不会传播到子例程本身。即使您已将所有代码放入一个源文件中，子例程也会被编译为单独的单元。
这意味着 sub1 没有可用于 sub2 的显式接口，因此使用隐式接口，其中参数 x 被假定为实数标量。

只需将两个子例程放在一个模块中并将use 该模块放在主程序中，就可以避免所有这些，自动使显式接口可用。这样就不用自己提供接口了，容易出错。

作为旁注，我建议在所有代码中使用隐式 none。