format d'écriture dans un fichier à partir d'une variable

Salut à tous,

C'est une question très simple que je me pose mais je n'ai pas trouvé de sujet relatif sur le forum pour y répondre.

En fait, j'aimerais écrire (en fortran 90) un tableau A de ni x nj éléments dans un fichier *.dat.
Je compte obtenir dans le fichier exemple1.dat les éléments de A sous forme de ni lignes et nj colonnes.
Jusque là, ça va...

Code:

---

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

(implicit none)   integer :: i,j integer :: ni,nj real (kind=8),dimension(:,:),allocatable :: A   . . .   open(30,file="exemple1.dat",status='unknown',form='formatted')   do i=1,ni write(30,1120) (A(i,j),j=1,nj) enddo   close(30)

---

Dans un premier temps, imaginons que A comprenne 20 x 20 éléments, le tableau en sortie devra donc avoir 20 lignes et 20 colonnes. Pour ce faire, j'écris :

Code:

---

1 2	1120 format (20(e12.5,1x))

---

Maintenant, comme je traite des tableaux qui n'ont pas forcément les mêmes dimensions (toujours en utilisant le même programme), comment dois-je faire pour rendre ce bout de code plus généraliste ?

J'avais essayé :

Code:

---

1 2	1120 format (ni(e12.5,1x))

---

... mais cela renvoie une erreur :

Code:

---

1 2 3 4

exemple.f90(15): error #6186: This character is not valid in format list. [n] 1120 format (ni(e12.5,1x)) -------------^

---

Pouvez-vous m'aider svp ? :roll:

Solution #1

Le format peut prévoir plus de plage que la liste de variable n'en exige.

Code:

---

1 2 3 4 5

do i=1,ni write(30,1120) (A(i,j),j=1,nj) enddo 1120 format(9999(e12.5,1x))

---

Solution #2

Générer le format.

Code:

---

1 2 3 4 5 6 7

character (20) :: MonFmt ... write(MonFmt,'(a,i0,a)') '(',nj,'(e12.5,1x))' do i=1,ni write(30,MonFmt) (A(i,j),j=1,nj) enddo

---

Solution #3

Tu utilises la clause advance='no' du write et tu écris élément par élément.

Solution #4

NON RECOMMANDÉE. Marchera avec certains compilateurs (Intel) mais c'est non portable.

Code:

---

1 2 3 4 5

do i=1,ni write(30,1120) (A(i,j),j=1,nj) enddo 1120 format(<nj>(e12.5,1x))

---

Merci beaucoup Sylvain, effectivement, il n'y avait pas à chercher très loin :ptdr:

Je vais tout de suite tester tout ça. Merci encore!