Discussion :

[yann_m]

Bonjour,

la question est sans doute trivial (quoique) mais je suis bloqué depuis un petit moment sans trouver de solution à mon problème...

J'ai un code qui génère un fichier de sortie a.dat, ce fichier s'organise sous forme de trois colonnes de réels x y &.
Je souhaite organiser ces données sous une autre forme dans un second fichier. Pour ne pas relancer le programme (convergence de l'ordre de la semaine), j'aimerai lire uniquement la troisième colonne de mon premier fichier à partir d'un second programme.

J'ai essayé avec le programme suivant mais par défaut je récupère la première colonne de mon fichier a.out.

Quelqu'un aurait-il une idée de l'option à ajouter ?
Merci par avance,
Yann.

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
21
22
23
24
25
26
27
28
29
        PROGRAM main
 
        IMPLICIT NONE 
 
        INTEGER   (KIND=4) :: i,j 
        REAL      (KIND=8),DIMENSION(100,100) :: mat
 
!Ouverture du fichier de entre & sortie
        OPEN(unit=1,file='a.dat')
        OPEN(unit=2,file='abis.dat')
 
!Lecture fichier entre
        DO i=1,100
                DO j=1,100
                        READ(1,*) mat(i,j)
                ENDDO
        ENDDO
 
!Reecriture sous forme d une matrice dans fichier de sortie
        DO i=1,100
                DO j=1,100
                        WRITE(2,*) mat(i,j)
                ENDDO
        ENDDO
 
        CLOSE(1)
        CLOSE(2)
 
        END PROGRAM main

[Sylvain Bergeron]

Le problème est simplement qu'en faisant :

• READ(1,*) mat(i,j)

tu demande à ton programme de lire un réel à chaque ligne. Ton programme lit donc le premier nombre qu'il trouve. Tu devras lire 3 variables à chaque ligne et n'utiliser que la 3ième valeur pour obtenir le résultat désiré.

Autres commentaires et suggestions :

• Il est préférable d'utiliser des units > 10, et non 1 et 2. C'est une bonne habitude et développer; certains compilateurs utilisent ces units par défaut pour les stdin, stdout et stderr.
• Si tu écrivais dans le fichier de sortie dès la lecture de chaque ligne, tu pourrais éviter le besoin d'un tableau et avoir un programme beaucoup plus simple.

Le programme suivant fait ce que tu veux je crois, mais en bien plus simple :

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
 
program main
 
   implicit none 
 
   integer :: i,j 
   real (kind=8) :: x,y,z
 
!Ouverture du fichier de entre & sortie
   open(unit=11,file='a.dat')
   open(unit=12,file='abis.dat')
 
!Lecture / écriture
   do i=1,100
      do j=1,100
         read(11,*) x,y,z
         write(12,*) z
      enddo
   enddo
   end program main

[yann_m]

Salut,

un double remerciement !
Le premier pour la méthode, le second pour les conseils.
J'ai en effet déjà rencontré des problèmes avec des unit > 5 sans bien savoir d'o`u cela provenait sans toutefois creuser ayant rarement plus de 5 fichiers de sortie !

Merci un grand merci à toi.

Yann.