Format écriture et "advance='no'"

Bonjour

Voila, je suis en train d'écrire une subroutine me permettant d'exporter des résultats au format VTk.

Le format est un format ASCII dont la structure est la suivante:

DATASET RECTILINEAR_GRID
DIMENSIONS nx ny nz
X_COORDINATES nx [dataType]
x0 x1...x(n-1)
Y_COORDINATES ny [dataType]
y0 y1...y(n-1)
Z_COORDINATES nz [dataType]
z0 z1...z(n-1)



Voici un extrait de mon code:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

write(file,'(a)')"# vtk DataFile Version 4.2" write(file,'(a)')"Fichier vtk" write(file,'(a)')"ASCII" write(file,'(a)')"DATASET RECTILINEAR_GRID" write(file,'(a,3i4)')"DIMENSIONS",nx,ny,nz write(file,'(a,i4)')"X_COORDINATES",nx write(file,'(f5.3)')(R(I),I=1,nx) write(file,'(a,i4)')"Y_COORDINATES",ny write(file,"(f5.3)")(Z(I),I=1,ny) write(file,'(a,i4)')"Z_COORDINATES",nz write(file,"(f5.3)")(1._8,I=1,nz) write(file,'(3a)')"SCALARS","TEST","type float" do I=1,nx do J=1,ny do K=1,nz write(file,'(f10.2)')Conc_sol(I,J,K) enddo enddo enddo

---

Mon problème vient des lignes

Code:

---

write(file,"(f5.3)")(Z(I),I=1,ny)

---

Je voudrais éviter le retour charriot mais conserver le controle sur le format d'écriture. Or, si j'utilise "advance='no'", j'ai un retour à la ligne.

Si j'abandonne l'idée du format et utilise une étoile *, j'ai un retour à la ligne même avec le advance no...

Une idée

Je ne sais pas pour advance=no. Je ne l'utilise presque jamais.

Il y a quelques techniques antérieures à l'introduction du advance=no. Ce sont les suivantes :

Format plus grand

Rien n'oblige le format à prévoir le nombre exacte de variables de la liste. Si le format est trop court, il y a boucle et changement de ligne (c'est ce que tu observes). Si le format est trop long, l'écriture arrête quand un masque du format n'a pas de variable correspondante (par exemple, write(,'(i2,"=",i0)')12 imprimera 12=

Donc, tu peux faire write(file,"(9999f5.3)")(Z(I),I=1,ny)

Cette technique requiert évidemment que tu aies une idée du nombre maximale d'élément à imprimer.

Pour les puristes

character (20) MonFmt
...
write(MonFmt,'("(",i0,"f5.3)")') ny
write(file,MonFmt)(Z(I),I=1,ny)

À éviter

Les formats variables de type <n> qui sont des extensions.

Bonjour,

Pour complèter la réponse de Sylvain (la solution "pour puristes" est parfaite; j'y adhère à 100%):
- Le 'advance=no' fonctionne sans problème mais il ne peut pas être utilisé conjointement à un format libre '*'
- Une autre technique pour supprimer le retour à la ligne est de mettre un '$' en fin de format:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

program ecriture   implicit none   character(len=10) :: phrase1="premier" character(len=10) :: phrase2="deuxieme" integer :: nombre=7   write(*,*) "Utilisation du $ pour supprimer le retour a la ligne:" write(*,'(a,i2,$)') trim(phrase1),nombre write(*,*) trim(phrase2)   write(*,*) 'Utilisation de advance="no" pour supprimer le retour a la ligne:' write(*,'(a,i2)',advance="no") trim(phrase1),nombre write(*,*) trim(phrase2)   end program ecriture

---

Bonne continuation.

Merci beaucoup.

Bonsoir

Citation:

---

- Une autre technique pour supprimer le retour à la ligne est de mettre un '$' en fin de format:

---

Le problème est que l'utilisation du $ ne fait pas partie de la norme Fortran 90 et n'est donc pas portable.
Ainsi, sous linux avec gfortran, la compilation de ton code me renvoie :

Code:

---

1 2 3 4

write(*,'(a,i2,$)') trim(phrase1),nombre 1 Error: Extension: $ descriptor at (1)

---

Je vais peut-être passer pour un intégriste (j'ai l'habitude ;) ) mais je ne pense pas qu'il faille s'affranchir de la norme, car des années plus tard, on peut se retrouver obligé de reprendre tout un code ou poser de gros problèmes à quelqu'un qui doit le reprendre.

Citation:

---

Envoyé par Grame

...Je vais peut-être passer pour un intégriste (j'ai l'habitude ;) ) mais je ne pense pas qu'il faille s'affranchir de la norme, car des années plus tard, on peut se retrouver obligé de reprendre tout un code ou poser de gros problèmes à quelqu'un qui doit le reprendre.

---

Je seconde !

Est-ce qu'il y a plus frustrant que de devoir comprendre le programme d'un autre et le corriger uniquement parce qu'il contient des extensions non portable ? Si on est obligé parce qu'il n'y a pas de solution dans la norme, ça peut être acceptable (la situation était fréquente avec les normes F66 et F77), mais quand une solution existe dans la norme, c'est inacceptable.

Bonjour,

Entièrement d'accord avec le fait qu'il faut éviter de sortir de la norme (en ce qui concerne la programmation en tout cas...).
Je ne savais pas que le '$' n'était pas dans la norme Fortran 90; j'emploie habituellement advance="no" .

Citation:

---

Ainsi, sous linux avec gfortran, la compilation de ton code me renvoie :

Code:

---

1 2 3 4

write(*,'(a,i2,$)') trim(phrase1),nombre 1 Error: Extension: $ descriptor at (1)

---

---

Ca c'est très surprenant... chez moi, sous linux avec gfortran (version 4.1.2), ça compile sans problème...

Bonjour

Citation:

---

Je ne savais pas que le '$' n'était pas dans la norme Fortran 90; j'emploie habituellement advance="no" .

Ca c'est très surprenant... chez moi, sous linux avec gfortran (version 4.1.2), ça compile sans problème...

---

Non, c'est parce que j'utilise les options:

Code:

---

1 2	-pedantic -std=f95

---

qui permettent justement de détecter ce qui n'est pas standard.

Code:

---

1 2	man gfortran (ou mieux : info gfortran)

---

est ton ami.

N : je suis en version 4.5.1 également sous linux