Discussion :

[malgachie]

Bonjour,

Je cherche à écrire un tableau dans un fichier avec description du format. Cela fonctionne très bien, mais j'ai observé qu'au-delà d'une certaine taille par ligne, mon output se divise en deux lignes (comme si un nombre de caractères maximum était atteint).

Ex. de programme (structure similaire au programme initial) :

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
 
	PROGRAM output_probl
	IMPLICIT NONE
	INTEGER, PARAMETER :: alen=110
	DOUBLE PRECISION, DIMENSION(alen)::a
	INTEGER :: i,j
 
	OPEN(1,file='output.dat')
2	format(I5,1x,110(e14.6e3,1x))
 
	DO i=1,15
	 DO j=1,alen
	  a(j)=(i*j**2)*0.0123456789
	 ENDDO
	 write(1,2)i,(a(j),j=1,alen)
	ENDDO
 
	END PROGRAM output_probl

A noter que l'output se passe très bien si je change le format pour l'écriture, par ex en passant de 14 à 8 caractères - ce qui n'est évidemment pas une solution satisfaisante. Tout se passe donc comme si une taille maximale de caractères par ligne était atteinte.
Un détail pas trop compliqué doit m'échapper, mais je cherche sans succès depuis quelque temps.

Quelqu'un voit-il une solution à mon problème (sans changer le format d'output) ?

Merci beaucoup pour votre aide.

m.

[FR119492]

Salut!
Je ne suis pas certain que ma réponse résoudra ton problème, mais il y a une chose que beaucoup de gens ignorent, à savoir que, pour des raisons historiques, l'utilisation de numéros d'unités inférieurs à 10 peut être "dangereuse". Alors, essaie de modifier ton programme avec Open(11,... etc.
Jean-Marc Blanc

[malgachie]

Bonjour Jean-Marc,

Merci pour ta réponse.
En fait j'ai utilisé ces références par simplicité pour l'exemple (et effectivement j'aurais plutôt dû utiliser 20 et 30 par exemple), mais ce n'est pas la source du problème (qui persiste si l'on change cela, j'ai vérifié).
Ce n'est donc malheureusement pas cela le problème...

m.

[Sylvain Bergeron]

J'ai légèrement modifié ton programme comme suit :

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
	PROGRAM output_probl
	IMPLICIT NONE
	INTEGER, PARAMETER :: alen=11000
	DOUBLE PRECISION, DIMENSION(alen)::a
	INTEGER :: i,j
	character (40) ::MonFmt
 
	OPEN(11,file='output.dat')
	write(MonFmt,'(a,i0,a)')'(I5,1x,',alen,'(e14.6e3,1x))'
 
	DO i=1,15
	 DO j=1,alen
	  a(j)=(i*j**2)*0.0123456789
	 ENDDO
	 write(11,MonFmt)i,(a(j),j=1,alen)
	ENDDO
 
	END PROGRAM output_probl

Et je ne suis pas capable de reproduire ton problème. En fait, mon éditeur fait des siennes si alen est grand, mais c'est uniquement un problème d'éditeur incapable d'afficher un fichier très large. Il affiche en insérant des « word-wraps » à partir de la colonne 4096.

Peut-être as-tu aussi un simple problème d'affichage.

[malgachie]

Bonjour Sylvain,

J'ai trouvé où était le problème, et tu tenais en effet le bon bout.
Il n'y avait aucun souci avec l'écriture de l'array, et comme tu le suggérais c'était un problème d'affichage. Je pensais avoir un souci à ce niveau car mon post-processing plantait, mais c'était lié à une autre erreur, très discrète, et je cherchais erronément du côté d'un problème d'écriture.

Merci pour l'aide !

m.

[Sylvain Bergeron]

Alors, tu me dois une bière !