Discussion :

[bricksters]

Bonjour,
j'ai un petit problème sur les entrées sortie (que je maitrise pas totalement).
A l'intérieur d'une boucle qui calcul des vecteurs, je voudrais enregistrer le vecteur obtenu dans un fichier, puis lors de la boucle suivante, le nouveau vecteur dans un nouveau fichier jusqu'a n !

do j=1,10
f=j
U(1)=fx
********
Calcul du vecteur en fonction de fx
x=....
********

do k=1,n
write(j, fmt=*) k, x(k)
enddo

if (x_max>100) exit
enddo

Au final, j'obtiens des fichiers fort.1 ... fort.10 avec les composantes du vecteurs comme je veux.
Néanmoins, je trouve pas ça "très propre" car je controle pas le nom du fichier et surtout je n'ai pas besoin d'utiliser open et close...

[yogitetradim]

tu peux par exemple, ouvrir toutes tes unités dans premier temps, en donnant un nouveau nom de fichier a chaque fois puis faire ces écritures...et enfin fermer.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
do k=1,n
write(ficname,'(a3,i0)')"fic",n
open(j, file=ficname)
end do
 
do k=1,n
write(j, fmt=*) k, x(k)
enddo
 
do k=1,n
close(j)
end do

je penses que ceci devrait correspondre à ta demande, si je ne me trompes pas.

[Sylvain Bergeron]

tu peux par exemple, ouvrir toutes tes unités dans premier temps, en donnant un nouveau nom de fichier a chaque fois puis faire ces écritures...et enfin fermer.

Il serait plus « joli » de déplacer l'ouverture et la fermeture dans la boucle, donc de n'avoir qu'un fichier ouvert à la fois.

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
 
do j=1,10
   f=j
   U(1)=fx
   ...
!  Calcul du vecteur en fonction de fx
   x = ....
   ...
!  Ouverture du fichier
   write(NomFichier,'(a,i0,a)')"fic",n,'.txt'
   open(11, file=NomFichier) ! il est inutile d'utiliser un unit variable
 
!  Ecriture
   do k=1,n
      write(11, fmt=*) k, x(k)
   enddo
 
   close(11)
   if (x_max>100) exit
enddo

[bricksters]

Merci Sylvain Bergeron, j'obtiens exactement ce que je voulais. Parfait.

Juste un chose pour comprendre, on utilises toujours l'unité n°11 mais pourtant plusieurs fichiers sont crées ? Je pensais que le unit=x où x représente un fichier et un seul...

Autrement, savez vous s'il serait possible d'utiliser gnuplot pour tracer les fichiers obtenu sur une même courbe ?

Car je fais ensuite call system('gnuplot graph_fx.dat')
mais dans graph_fx.dat je suis obligé d'ecrire le nom de tout mes fichiers un à un, c'est moyen !
J'ai vu cette solution
c = c + 1
#... some gnuplot commands ...
if(c<100) reread

mais ensuite impossible d'appeler mon fichier avec fic'c'.txt

[bobbyboy]

Bonjour,

Je pensais que le unit=x où x représente un fichier et un seul...

C'est vrai mais uniquement durant la durée d'ouverture du fichier (instruction OPEN). A la fin de la boucle, l'instruction CLOSE(11) libère le numéro, tu peux donc le ré-utiliser à l'itération suivante pour ouvrir un autre fichier.

Autrement, savez vous s'il serait possible d'utiliser gnuplot pour tracer les fichiers obtenu sur une même courbe ?

Car je fais ensuite call system('gnuplot graph_fx.dat')
mais dans graph_fx.dat je suis obligé d'ecrire le nom de tout mes fichiers un à un, c'est moyen !

Je n'ai pas bien compris ce que tu voulais faire. Si tu veux tracer tous tes résultats sur un même graphe à la fin de ta boucle, tu peux créer un fichier "trace.plt" avec les instruction gnuplot :

Code trace.plt :

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

1
2
3
4
 
plot "fichier_1.txt"
replot "fichier_2.txt"
...

Puis l'appeler comme tu fais avec call system('gnuplot trace.plt').

[bricksters]

Ok pour la boucle sur le write, je n'y avais pas pensé !

Justement pour le fichier plot trace.plt, si j'avais beaucoup de fichier, je serais obliger de taper replot "fichier_[1 jusqu'a n].txt", je ne pourrais pas faire une boucle sur le replot ?

Comme ceci :

do i=1,n
replot 'fichier_'i'.txt'
enddo

J'éviterai de devoir ecrire une nouvelle ligne pour chaque fichier.

[Ehouarn]

Bonjour,

Le plus simple serait d'écrire le fichier 'graph_fx.dat' dans ton programme fortran (avant le "call system(gnuplot graph_fx.dat") puisque tu connais dans ce dernier le nombre (et les noms) des fichiers de données à utiliser.

Au passage, je recommanderai plutôt d'utiliser la syntaxe:

Code :

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

1
2
 
plot "fichier_1.txt","fichier_2.txt","fichier_3.txt"

que

Code :

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

1
2
3
4
 
plot "fichier_1.txt"
replot "fichier_2.txt"
replot "fichier_3.txt"

Car dans ce dernier cas, gnuplot fait trois plot: un premier avec fichier_1.txt, qui est immédiatement écrasé par un deuxième avec fichier_1.txt et fichier_2.txt, qui est lui-même écrasé par un troisième plot utilisant les trois fichiers.

[Sylvain Bergeron]

Un (léger) avantage de Fortran sur plusieurs langage est que le délimiteur de chaîne de caractère peut être " ou ', à condition d'utiliser le même aux 2 bouts de la chaîne. Ça évite souvent le recours aux """"...

Dans ton cas, tu peux faire :

Code :

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

1
2
3
4
5
6
7
8
 
   character (1000) Ligne
   ...
   Ligne = 'plot "' // NomFichier(1) // '"'
   do k= 2, n
      Ligne = trim(Ligne) // ',"' // NomFichier(k) // '"'
   enddo
   ...

[bricksters]

Ah j'étais en train de retaper mon message !

Je vais essayer avec ta méthode, merci

[Sylvain Bergeron]

Je ne l'ai pas écrit initialement, mais tu peux combiner les writes dans la construction :

write(Ligne,'(a,a,i,a)') trim(Ligne), ', "Fic',k,'.txt"'