从大小数组中分配数组

Question

假设我有一个这样声明的维度数组：

integer, dimension(5) :: dims
dims = (/ 5, 6, 7, 8, 9 /)

我怎样才能最简单地使用这个数组来分配另一个数组，使用dims 的元素来指定各个维度的大小？假设第二个数组是这样声明的：

real, dimension(:,:,:,:,:), allocatable :: dset

为什么我不能执行以下操作？

allocate (dset(dims(:)))

当我尝试编译时得到Error: Rank mismatch in array reference。

我意识到我可以做到：

allocate (dset(dims(1), dims(2), dims(3), dims(4), dims(5)))

但我正在寻找更容易扩展的东西。

Answer 1

由于编译器给出的确切原因，您不能写allocate(dset(dims(:)))：它们有不同的等级。如果您要打印shape(dims)，您会在屏幕上看到5；如果您要打印shape(dset)，您会看到0 0 0 0 0 或5 6 7 8 9（取决于您是否分配了它）； dims 是 rank-1 数组，dset 是 rank-5 数组。

因此，分配可分配对象的唯一方法几乎是通过您认为不可扩展的显式方法。