Discussion :

[djocin]

Hello

J'ai de nouveau un petit soucis que je n'arrive pas à résoudre, peut-être sauriez-vous m'aider.

Dans une boucle, je souhaite écrire dans un fichier un tableau de dimension 1 sur une même ligne, puis à l'itération suivante je voudrais passer à la ligne en dessous pour écrire les nouvelles valeurs sur une même ligne, ainsi de suite.

Voilà ce que j'ai écrit mais qui ne fonctionne pas :

Code :

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

1
2
3
4
5
6
7
 
do i=1,p
open(1,nomfich,formatted,unknown)
call bidule(RB)
write(1,'(f31.30)')RB
close(1)
enddo

où RB est un vecteur, que je veux donc écrire sur une ligne différente à chaque itération.

Si vous avez la solution ça me soulagerait !

[kango]

bonjour,

c'est probablement la question qui est la plus fréquemment posée mais:

ça veut dire quoi "ça ne fonctionne pas" ?

là par exemple dans ton cas le spécificateur de format ne devrait afficher qu'un réel de la forme:

Code :

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

.123456789012345678901234567890

donc pas de nombres négatifs, ni de nombres supérieurs à 1.

si ton vecteurs contient 3 nombres:

Code :

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

fmt='(3(f9.6,x)'

par exemple

s'il contient un nombre inconnu de nombre, mais que tu sais borner (par exemple pas plus de 100):

Code :

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

fmt='(100(f9.6,x)'

si en plus tu ne sais pas le borner:

Code :

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

1
2
3
4
do i=1,n
    write(unit,'(f9.6,x)',advance='no')RB(i)
end do
write(unit,*)

[djocin]

Merci de la réponse.

J'ai vu que la question était récurrente, mais aucune solution proposée n'a fonctionné (cad que ça ne m'affiche pas la valeurs comme je veux).

Par contre petite précision sur ma boucle : elle est faite sur j, mais les coeff de RB ne sont pas les RB(j). Le call bidule me stocke la valeur de RB.

J'ai déja essayé avec advance='no', ou avec un $ à la fin du format, mais sans résultat.

En fait ça fonctionne lorsque j'écris mon format sous la forme 8(f32.31), cad en mettant explicitement que j'ai 8 termes dans mon vecteurs. Mais je voudrais pouvoir changer ce nombre de termes, l'idéal serait de pouvoir écrire n(f32.31), le nombre n étant obtenu par un read*,n.

Une petite idée?

[kango]

ok, dis nous comment tu déclares RB:

par exemple:

Code :

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

real,dimension(1:100):: RB

dis nous aussi avec un exemple comment tu veux que ton fichier texte apparaisse.

ensuite, montre nous comment est ta boucle

on peut pas deviner ce que tu ne nous dis pas, même si c'est évident pour toi.

bref, pose une question PRECISE et tu auras une réponse PRECISE.

merci

[Sylvain Bergeron]

Plusieurs choses à dire :

• Si tu open le fichier dans la boucle, tu dois utiliser la close "append" pour déplacer le pointeur à la fin et ajouter une ligne. La solution classique est plutôt d'ouvrir le fichier avant la boucle et de le fermer après la boucle.
• Il est fortement conseillé d'utiliser un unit > 10, et non 1. Les petits units sont souvent prédéfinis par le compilateur.
• Pour écrire le vecteur sur une ligne, le plus simple est d'avoir un seul write écrivant tout le vecteur d'un coup. Le format peut prévoir un nombre d'élément supérieur à la taille du vecteur. Ex "(99999f13.6)"

Si la solution 99999 ne te convient pas, tu peut utiliser la technique suivante. Le format peut être une variable caractère. Tu peux donc construire le format à l'exécution. Il est possible "d'écrire" dans une variable caractère plutôt que dans un fichier avec un write. Par exemple :

Code :

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

1
2
3
4
5
6
7
8
9
 
character (40) :: MonFormat
...
<< saisie de n >>
...
write(MonFormat,'(a,i0,a)') '(',n,'f13.6)'
...
write(11,MonFormat) RB
...

Note: Il est déconseillé d'utiliser la technique write(11,'(<n>f13.6)') n,RB car elle est propriétaire (DEC --> Compaq --> Intel)

[djocin]

je crois que je vais utilisé la technique du 99999 devant le format, c'est la plus simple

merci