Boucle sur write (Creation de plusieurs fichiers)

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 !

Citation:

---

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...

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:

---

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.

Citation:

---

Envoyé par 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.

---

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:

---

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

---

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

Bonjour,

Citation:

---

Envoyé par bricksters

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.

Citation:

---

Envoyé par bricksters

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:

---

1 2 3 4

plot "fichier_1.txt" replot "fichier_2.txt" ...

---

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

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.

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:

---

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

---

que

Code:

---

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.

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:

---

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 ...

---

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

Je vais essayer avec ta méthode, merci ;)

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"'