No Fortran 90, as dimensões do array precisam ser declaradas previament

Question

Jun 28, 2011, 06:00 PM

No Fortran 90, as dimensões do array precisam ser declaradas previament

É necessário declarar as dimensões da matriz antes de qualquer outro código? Por exemplo, escrevi o seguinte código de exemplo simplificado:

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

O ponto que estou tentando destacar é que gostaria de criar uma matrizc isso ék-pork**2 em tamanho ek é determinado apenas por outros cálculos dentro do código;k não é conhecido no início.

Mas, o código acima me fornece a seguinte mensagem de erro em tempo de compilação:

mytest.f90:13.31:

  REAL, DIMENSION(k, k**2) :: c
                               1
Error: Unexpected data declaration statement at (1)

onde a linha 13 do meu código é a linha em que finalmente declaroc: REAL, DIMENSION(k, k**2) :: c.

Por outro lado, se eu declarark e especifique suas dimensões antecipadamente,

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

Recebo a saída correta:

 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

Mas, desde que eu não seik de antemão, não posso fazer exatamente isso no meu código real. Existe alguma maneira de "declarar" a matrizc inicialmente, sem especificar suas dimensões, e depois especifique as dimensões assim que o valor dek é conhecido

leaveComments