动态数组排名

Question

我有一个模块中的几个数组变量是动态的，稍后分配给模块外的两个子例程之一。但是，在一个子程序中，我希望数组是一维的，而在另一个子程序中，我希望它是二维的。

原则上我希望在模块中这样的事情，但我不相信这是可能的声明区域？：

if (option1) then 
    real (kind=8), allocatable :: arr1(:) 
else 
    real (kind=8), allocatable :: arr1(:,:) 
endif 

是否有可能分配申报办法有维是动态的？

编辑1：我这样做的原因是我添加一个新的子程序到现有的代码库，但我想向后兼容。 arr1仅由两个独立的子程序使用，主程序完全不使用它。下面是该理念的一些更完整的代码：

program myprog 
use inputs 

call read_inputs 

if (option1) then 
    call do1 
else 
    call do2 
endif 

contains 
    subroutine read_inputs 
     use inputs 
     use mymod 
     !!!read from file .logical. option1, integers N1, N2 

     !allocate arrays 
     if (option1) then 

     else 

     endif 
    end subroutine read_inputs 

    subroutine do1 
     use inputs 
     use mymod 

     allocate(arr1(N1)) 

     !do stuff with arr1 
    end subroutine do1 

    subroutine do2 
     use inputs 
     use mymod 

     allocate(arr1(N1,N2)) 

     !do stuff with arr1 
    end subroutine do2 

end program 

module inputs 
    logical :: option1 
    integer :: N1, N2 

end module inputs 

module mymod 
    use inputs 
    !!!!can I do something here to make the rank of arr1 dynamic? I don't think the following will work 
    if (option1) 
     real (kind=8), allocatable :: arr1(:) 
    else 
     real (kind=8), allocatable :: arr1(:,:) 
    endif 
end module mymod 

我可能只是在mymod，ARR1两个独立的变量，arr1_new。我只是希望避免这种情况。

Answer 1

我认为“古代”的方式做这样的事情是通过第一个元素，而不是整个阵列和阵列的单独大小：

program dyn_array 
    implicit none 
    integer :: a(2, 3) 
    integer :: i 
    call set1d(a(1,1), size(a)) 
    do i = 1, 3 
     write(*, '(2I4)') a(:,i) 
    end do 
    contains 
     subroutine set1d(array, s) 
      implicit none 
      integer, intent(in) :: s 
      integer, intent(out) :: array(s) 
      integer :: i 
      do i = 1, s 
       array(i) = 3 * i 
      end do 
     end subroutine set1d 
end program dyn_array 

Answer 2

“你能传递一个二维数组需要一维数组并获得正确大小的子程序？“

您可以使用重塑，但如果您的代码是依靠编译器来帮助，然后2D到1D是dicey。你可以在之前和之后使用RESHAPE ...或者你可以有2个例程，我们可以调用set1d和set2d。 然后在模块中，您可以选择您想要使用的模型。你可以有整数，浮点数，复数，字节。

您将调用数组$集（数组，S）

MODULE Stuff 
    PUBLIC :: Array$Set 

    PRIVATE 
    INTERFACE Array$Set 
    MODULE PROCEDURE Set1d_Float, Set1D_Double, set2D_Float, Set2D_Double 
END INTERFACE Array$Set 

CONTAINS 

SUBROUTINE Set1D_Float(Array,S)... 
!$DIR ATRIBUTES ASUUME_ALIGND:64    :: Array 
REAL, DIMENSION(:,:), CONTIGUOUS, INTENT(INOUT) :: Array 
REAL, DIMENSION(:,:),    INTENT(IN ) :: S 
REAL, DIMENSION(2)        :: Shapez 
... 
Shapez = Shape(Array) 
DO I = 1, Shapez(1) 
    DO J = 1, Shapez(2) 
... 
END SUBROUTINE Set1D_Float 

END MODULE Stuff 

---

对于示例： 如果（选项1），然后 实（KIND = 8），可分配:: ARR1（:) 其他 实（KIND = 8），可分配:: ARR1（:, :) ENDIF

我建议这样的：

!DIR ATTRIBUTES ALIGN:64  :: Arr1 
REAL, DIMENSION(:), ALLOCATABLE :: Arr1 
... 


if (option1) then 
    ALLOCATE(Arr1(<#>)) 
else 
    ALLOCATE(Arr1(<#>*<#2>)) 
    RESHAPE(Arr1, SHAPE=/(#1,#2)) !Check the syntax 
endif 

CALL Array$Set(Arr1,s) !It'll find the right one... 

!... at the bottom ... 
IF(ALLOCATED(Arr1)) DEALLOCATE(Arr1) 
END PROGRAM