释放模块中声明的变量时出错

Question

我的 Fortran 90 代码出现问题，涉及解除分配在模块中声明的数组，然后在子例程中分配和初始化。在我的程序中，我在 real*8, dimension(:), allocatable :: test 之类的模块中声明了一堆数组。然后，在一个初始化子程序中，我 use 模块，分配变量 allocate(test(8))，并用test = 0.d0初始化。

在此之后，我可以print*, test 并获得适当的输出：0.E+0 0.E+0 0.E+0 0.E+0 0.E+0 0.E+0 0.E+0 0.E+0。此外，调用allocated(test) 会返回.TRUE.。尽管如此，还是出了点问题。直接在调用allocate() 之后，调用sizeof(test) 返回0，调用deallocate(test) 会引发以下错误：

lib-4422 : UNRECOVERABLE library error 
  A DEALLOCATE statement argument points to a portion of the
 original allocation.
  Original size in bytes for the argument is 512
  Current size in bytes for the argument is 0
Segmentation fault

这一切都发生在一个更大的代码中，我在整个代码中都使用了这些数组，没有任何错误。我只是在寻找内存泄漏时试图在程序结束时释放内存时才注意到这个问题。我试图制作一个仅执行上述操作的简单程序（即，在模块中声明，在子程序中分配和初始化，然后打印数组并在同一个子程序中释放它）。然而，这个简单的代码可以正常工作并且没有错误地运行。因此，我对在较大代码的上下文中可能导致其行为不端的原因感到非常困惑。此外，如果在我较大的代码中，我将声明行从模块移动到子例程，一切都会正常运行。

任何建议将不胜感激！提前致谢，

~BCL

Answer 1

在处理不是单个子例程内部的数组时，我总是在声明它们时使用 pointer 属性而不是属性。在实践中，两者的工作方式相同。例如：

module Mod1
contains
  subroutine sub1(Array)
    !*** Declaration of Variables
    implicit none
    !*** External Variables - Arguments
    real(kind=8), pointer, intent(inout) :: Array(:)
    !*** Internal Variables
    real(kind=8), allocatable            :: InternalArr(:)
    !*** Memory allocation
    allocate(Array(1:8))
    allocate(InternalArr(9:10))
    !*** End of Subroutine
  end subroutine sub1
end module Mod1

调用sub1 后，您的数组将在代码的任何位置保持其维度。您可以根据需要deallocate 和allocate。

在使用pointer 之前，我遇到了一些类似您描述的问题，有人建议我改用pointer。从那以后，我已经非常习惯以这种方式声明数组（如果它们打算保留在我一直使用allocatable 的同一个子例程或模块上），并且它一直被证明是控制数组的好方法。

我真的希望这会有所帮助...这是我的第一个答案 :) 编辑：我更正了示例代码中的一个小错误