Вы хотите использовать размещаемые массивы:
ли объявлять размеры массива перед любым другим кодом? Например, я написал следующий упрощенный пример кода:
PROGRAM mytest
IMPLICIT NONE
INTEGER :: i, j, k, mysum
! Let array c be a k-by-k**2 array
! Determine k within the program by some means...for example,
mysum=0
DO i=1, 3
mysum=mysum+1
END DO
k=mysum
REAL, DIMENSION(k, k**2) :: c
WRITE(*,*) "k=", k
WRITE(*,*) "k**2=", k**2
WRITE(*,*)
DO i=1,size(c,1)
WRITE(*,"(100(3X,F3.1))") (c(i,j), j=1,size(c,2))
END DO
END PROGRAM mytest
Дело в том, что я хочу создать массив.c
этоk
-по-k**2
по размеру иk
определяется только другими вычислениями в коде;k
неизвестно в самом начале.
Но приведенный выше код дает мне следующее сообщение об ошибке во время компиляции:
mytest.f90:13.31:
REAL, DIMENSION(k, k**2) :: c
1
Error: Unexpected data declaration statement at (1)
где строка 13 в моем коде это строка, где я, наконец, объявляюc
: REAL, DIMENSION(k, k**2) :: c
.
С другой стороны, если я вместо этого объявитьk
и указать его размеры заранее,
PROGRAM mytest
IMPLICIT NONE
INTEGER :: i, j, k, mysum
REAL, DIMENSION(3,9) :: c
! Let array c be a k-by-k**2 array
! Determine k within the program by some means...for example,
mysum=0
DO i=1, 3
mysum=mysum+1
END DO
k=mysum
WRITE(*,*) "k=", k
WRITE(*,*) "k**2=", k**2
WRITE(*,*)
DO i=1,size(c,1)
WRITE(*,"(100(3X,F3.1))") (c(i,j), j=1,size(c,2))
END DO
END PROGRAM mytest
Я получаю правильный вывод:
k= 3
k**2= 9
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
Но так как я не знаюk
заранее я не могу сделать именно это в моем реальном коде. Есть ли способ «объявить» массивc
первоначально, без указания его размеров, а затем позже укажите размеры, как только значениеk
известен?