Discussion :

[manaliac]

Bonjour,

je suis entrain de coder un programme qui lit un fichier texte ou sont stockées des valeurs par colonne de la manière suivante :

x y z p1 p2
-8.509950340E-002 8.100000024E-002 0.00E+000 0.00E+000 0.00E+000
-8.420702815E-002 8.145157248E-002 0.00E+000 0.00E+000 0.00E+000
-8.331453800E-002 8.190311491E-002 0.00E+000 0.00E+000 0.00E+000
-8.242201060E-002 8.235459030E-002 0.00E+000 0.00E+000 0.00E+000
-8.152943850E-002 8.280596882E-002 0.00E+000 0.00E+000 0.00E+000
-8.063680679E-002 8.325722814E-002 0.00E+000 0.00E+000 0.00E+000
-7.974408567E-002 8.370833099E-002 0.00E+000 0.00E+000 0.00E+000

Mon nombre de lignes vient d'une information lue plus haut dans le fichier et qui est de : Nbi*Nbj*Nbk.
Je veux qu'il assimile chaque colonne lue à une variable tridimentionnelle de taille (Nbi,Nbj,Nbk), mais je n'y arrive pas, il ne lit pas le fichier!

Est ce que quelqu'un aurait une idée??

[genteur slayer]

et bien tu ouvre ton fichier, c'est un fichier texte donc cette commande devrai suffire:

Code :

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

1
2
 
open(unit=60,file="fichier",status="old",position="rewind",ierr=ios)

60est le numéro d'unité, si tu l'utilise déjà, prends-en un autre, fichier est à remplacer par le nom de ton fichier (éventuellement avec les répertoire) old signifie que le fichier existe déjà et oblige à renvoyé une erreur si il ne le trouve pas rewind est optionnel mais force à commencer au début du fichier
ios revoie le code d'érreur (0 si touot s'est bien passé)

du coup, ensuite:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
if (ios/=0) then 
  print*,"erreur à l'ouverture du fichier"
else !si il n'y a pas d'erreur
  read(60,*)NLigne  !on lit le nombre de ligne
  read(60,*)           !on saute une ligne qui contient les commentaires
  allocate (x(NLigne),y(NLigne),z(NLigne),p1(NLigne),p2(NLigne))
  do i=1,NLigne
    read(60,*) x(i),y(i),z(i),p1(i),p2(i)
  end do
  close(60)
end if

pour que cela marche bien, il faut déclarer tes variables:

Code :

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

1
2
 
real(kind=8),dimension(:),allocatable :: x,y,z,p1,p2

il faut selement pensé à déallouer tes tableau quand tu n'en a plus besoin

Code :

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

1
2
 
deallocate(x,y,z,p1,p2)

une fois que tu as lu et stocké x,y,z,p1,p2, il devient facile de faire des changement d'indices pour obtenir des tableau à plusieur dimensions... (note: en mémoire quelque soit le nombre d'indice, tout est stocké de manière linéaire!!)

[manaliac]

Bonjour, en fait, j'ai déjà réussi à faire cela, et ça marche très bien ! Merci, ça confirme la justesse de mon programme !

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
19
20
	READ (7,'(a3,i3,a4,i1,a4,i3)') char1, Nbi, char2, Nbj, char3, Nbk
			PRINT *, ' --> Lecture de la sourface d integration:',Nbi,Nbj,Nbk
     !! Allocation en memoire pour les variables.....avant leur lecture!!!!!
	 !!*********************************************************************** 
	       d_end = Nbi*Nbj*Nbk
	        allocate(  X_Centers(d_end),   stat=stats)
			allocate(  Y_Centers(d_end),   stat=stats)
			allocate(  Z_Centers(d_end),   stat=stats)
			allocate(  ReTfP(d_end),   stat=stats)
			allocate(  ImTfP(d_end),   stat=stats)
 
 
     !!****************************************************
	 !!*******************
			READ (7,*)
			READ (7,*)
			do il=1, d_end 
	 		READ (7,*) X_Centers(il),Y_Centers(il ),Z_Centers(il ),ReTfP(il ),ImTfP(il )
 
			enddo

Est ce que tu aurais une idée par contre, de comment transformer un tableau unidimensionnel en tridimensionnel??
par exemple je stocke tout dans x(Nbi*Nbj*Nbk) avec Nbi, Nbj,Nbk, les trois dimensions recherchée, et je veux passer à un X_nouveau(Nbi, Nbj, Nbk).

Je n'arrive pas trop à trouver la correspondance d'indices !

Merci bien

Manal

[genteur slayer]

alors cela dépendra de l'ordre dans lequel tes données sont rangé, mais, mettons un rangement par le dernier indice puis le second et enfin le premier
tab_tri(1,1,1) = tab_uni(1)
tab_tri(1,1,2) = tab_uni(2)
... etc...
tab_tri(i,j,k) = tab_uni(n)
...

il ne rest qu'à trouver la fonction n=f(i,j,k)

alors tu peux faire une relation du style:

Code :

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

1
2
3
4
5
 
forall(i=1:Nbi,j=1:Nbj,k=1:Nbk)
  n=k+(i-1)*Nbi+(j-1)*Nbj
  tab_tri(i,j,k)=tab_uni(n)
end forall

en fonction de l'ordre dans lequel sont rangé tes donné dans le tableau tab_uni, il te faudra changer la ligne 'n=k+(i-1)*Nbi+(j-1)*Nbj' pour avoir une bonne corespondance