FORTRAN: deklaration in sub-sub ?

[ogx]

Hey, ich hab eine Frage zu einem Fortran Code.

Ich programmiere gerade eine Subroutine als Teil eines Umfangreichen Programms. Jetzt habe ich folgende Anordnung welche mir Probleme bereitet:

Hauptgrogramm(geht mich nix an)

-----------

subroutine MeineSub(..)

!hier wird eine Variable x festgelegt welche in einer so zu

!sagen sub-sub die grösse einer Matrix festlegt.

implicit none

integer x

parameter (x=5)

...

call subroutine Subsub(x)

...

end

------------

subroutine Subsub(x)

implicit none

integer x

real*8 aa(x,x)

...

end

---------------------

Jetzt habe ich das problem, dass in der subsub zwar die grösse der matrix richtig dimensioniert wird, aber ich mit ihr nichts berechnen kann, nicht mal null setzten ist möglich. im debugger erhalte ich nur "floating point" einträge.

Ist diese anordnung überhaupt möglich, also eine deklaration ein einer untersub ohne dass die matrix in der übergeordneten subrotine deklariert wird?? Wäre blöd weil ich sie nur in der SubSub brauche und nicht reinübergeben wollte.

wenn ich ganz normal durch eine zahl dimensioniere funktioniert alles. und der debugger sagt mir auch dass x in der subsub die richtige zahl enthällt.

Aber ich muss das programm leider variablen-abhängig halten.

Ich hoffe jemand kann mir hefen!

ogx

[casio]

Dynamische Arrays zu generieren, die nicht übergeben werden, dass ist nicht möglich, da eine Konstante in diesem Fall erwartet wird!

Ich hätte folgendes anzubieten, vielleicht hilft es dir weiter:

Ich würde die Passage mit (integer x, parameter (x=5)) in ein Include file auslagern, dann könntest du in abhängigkeit von x folgendermaßen arbeiten:

Include file inc.for

C

      INTEGER*4

     * X

C

      PARAMETER(X=5)


Dein Programm:

      program test

C

      IMPLICIT NONE

C

      INCLUDE 'INC.FOR'      

C      

      CALL SUB1()

C      

      end



      subroutine sub1()

C

      IMPLICIT NONE 

C

      INCLUDE 'INC.FOR' 

C

      INTEGER*4

     * I, J

C

      REAL*8

     * AA(X,X)

C

      DO I=1, X 

       DO J=1, X

        AA(I,J)=0.0d0

       ENDDO

      ENDDO

C     

      end

Bearbeitet 17. März 2010 von casio