Discussion :

[Dreyliciouss]

Bonjour,

Je m'adresse à vous tous sur ce forum car je n'ai pas trouvé la réponse à ma question.
En effet, je voudrais savoir s'il est possible d'allouer un tableau dont les dimensions sont lues dans un fichier binaire et de retourner ce tableau en sortie.

Je sais qu'allouer un tableau et lui attribuer ses dimensions dynamiquement en les récupérant dans un fichier sont possible mais je sais que l'instruction ALLOCATABLE n'est pas compatible avec l'instruction INTENT.

Est ce que quelqu'un a déjà réalisé quelque chose de similaire?

Merci d'avance pour vos réponses.

[genteur slayer]

tu peux passer par un module:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
module toto
  real(kind=8),dimension(:),allocatable :: titi
end module toto

puis dans une autre routine:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
subroutine litfichier
use toto
implicit none
 
...
open(unit=15,.....)
  read(15)taille
...
  if (allocated(titi)) deallocate(titi)
  allocate(titi(taille))
...
  read(15)titi
...
  return
end subroutine litfichier

ensuite tu n'a qu'à charger le module toto dans la/les routine qui utilise le tableau titi !!!

[Dreyliciouss]

Merci pour cette réponse
J'ai donc suivi votre méthode qui m'a l'air vraiment bien. Mais j'avoue que j'ai du mal à voir comment déclarer le tableau et dimensions dans chaque fonction.

[Dreyliciouss]

J'ai essayé d'inclure tes suggestions dans mon code mais j'obtiens les erreurs suivante :

1) Number of subscripts declared and used do not match (PROJ)
2) Illegal structure component reference.

PROJ correspond à titi dans ton exemple.
La deuxième erreur intervient lorsque j'essai de lire dans le fichier (je récupère ce que je lis dans le tableau PROJ avec l'instruction READ).

Je ne comprend pas du tout la première erreur. Auriez-vous une idée de ce à quoi elle peut correspondre?

[Sylvain Bergeron]

1) Number of subscripts declared and used do not match (PROJ)

C'est un problème de dimension. Tu fais référence à un tableau avec un nombre d'indice différent du nombre de dimensions déclaré. Par exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
real, dimension(:,:) :: Tableau ! 2 dim
...
Tableau(3) = ... ! 1 dim

[Dreyliciouss]

Nikel! L'erreur a disparue par contre, à la dernière ligne de code lorsque j'essai d'allouer les dimensions a et b lues précédemment dans le fichier, il me dit "ALLOCATE argument must be ALLOCATABLE or POINTER. (PROJ)". Cela vient sans doute de la déclaration de a et b.

Mais quand je déclare a et b avec la ligne suivante :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

INTEGER, ALLOCATABLE :: a,b

Il me dit :

ALLOCATABLE or POINTER array-spec must be deferred shape. (A)
ALLOCATABLE or POINTER array-spec must be deferred shape. (A)

et l'erreur précédente est toujours présente. Je ne sais pas comment les déclarer correctement.