Discussion :

[diddix]

Bonjour,
Je dois écrire un code en fortran qui traite des données récupérées dans un fichier texte.
Je suis débutant et ne sais pas vraiment me servir des formats.
Mes données sont représentées comme ceci:

Code :

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

1
2
3
4
5
6
7
 
x1            y1           z1
x2            y2           z3
x3            y3           z3
.              .            .
.              .            .
.              .            .

Cependant certains xi, yi ou zi peuvent être des blancs (absence de donnée). Je ne peux donc pas utiliser l'instruction:

read(étiquette_du_fichier,*) xi, yi, zi

De plus les espace entre les colonnes x y et z est aléatoire.
Le nombre de colonnes est fixe, la premiere numerote les lignes.

Exemple typique:

Code :

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

1
2
3
4
 
1               6  23     8
2  1               9      78
3  15           8  16

A partir de ceci, j'aimerais pouvoir récupérer 3 tableaux x, y, z et t, correspondant aux 4 colonnes, tels que:

x(1)=b; x(2)=1; x(3)=15

y(1)=6; y(2)=b; y(3)=8

z(1)=23; z(2)=9; z(3)=16

t(1)=8; t(2)=78; t(4)=b

b étant un sentinelle quelconque signifiant l'absence de valeur.

Quelqu'un connaîtrait-il une solution?

Merci d'avance.

[genteur slayer]

alors... si les valeurs sont rangée n'importe comment, comment fait-tu pour savoir qui est le banc dans l'histoire??? je veux dire:

Code :

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

1
2
3
 
 0          8
      0     8

on peux voir ces lignes comme:
x1=0 y1=b z1=8
x2=b y2=0 z2=8
mais aussi:
x1=0 y1=8 z1=b
x2=0 y2=8 z2=b
....
après, pour ces deux lignes, un 'read(X,*)x,y fonctionne parfaitement... il se fout de savoir où sont les blanc.....

[diddix]

En fait je me suis trompe le nombre de colonnes est fixe.
Aucune colonne n'est vide.
La premiere colonne numerote les lignes.

[genteur slayer]

si tu fais du fortran90:

Code :

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

1
2
3
4
5
6
 
open(unit=15,file="data.dat",status="old",iostat=ios)
  do while(ios==0)
    read(15,*)i,x(i),y(i),z(i)
  end do
close(15)

devrais faire la blague!
il te faut déclarer: deux entier (i et ios) et trois tableaux assez grand...

pour trouver dynamiquement la taille des tableaux:

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
17
18
 
integer i,ios
real(kind=8),dimension(:),allocatable x,y,z
!
open(unit=15,file="data.dat",status="old",iostat=ios)
  do while(ios==0)
    read(15,*)i
  end do
  allocate(x(i),y(i),z(i))
  rewind(15)
  do while(ios==0)
    read(15,*)i,x(i),y(i),z(i)
  end do
close(15)
!
..... traitement de tes valeurs.....
!
deallocate(x,y,z)

en fait, vu que tes collonne ne sont pas bien rangée, il ne faut pas utiliser les format qui sont fixe mais plutôt le '*' qui sera variable!!