Discussion :

[jal13]

Bonjour à tous

J'aurais une question par rapport à la manipulation création de fichier.
En fait j'utilise un logiciel qui simule la saturation dans le temps dans un sol. Le fichier de sortie (plot.out) est monté comme suit

VARIABLES = "X", "Y", "NmoleF", "So", "Sw"
ZONE T = "3.000 ", I = 41 , J = 41, F=POINT
0.0000 20.00 -10.00 0.0000 0.2000
0.5000 20.00 -10.00 0.0000 0.2000
1.000 20.00 -10.00 0.0000 0.2000
1.500 20.00 -10.00 0.0000 0.2000
..
ZONE T = "3.001 ", I = 41 , J = 41, F=POINT
0.0000 20.00 -10.00 0.0000 0.2000
0.5000 20.00 -10.00 0.0000 0.2000
1.000 20.00 -10.00 0.0000 0.2000
1.500 20.00 -10.00 0.0000 0.2000

où x et y sont les coordonnées du maillage.

La 5ème colonne (saturation) doit être extraite pour faire des manipulations.

J'utilise donc le code fortran comp.f pour faire ceci.
Cependant dans mon code je n'arrive pas à conserver les titres Zone T=...
De plus, selon les simulations le nombre de zone varie, donc a moins d'augmenter manuellement les commandes suivantes dans le fichier comp.f, je ne vois pas comment faire

do k=1,nz
do i=1,nx
nd = (k-1)*nx + i
read(90,*) x(nd), z(nd), nmoleF(nd), So(nd), Sw(nd)
enddo
enddo

write(91,4)(x(i),z(i),(xmoln(i,k), k=1,3), i=1,nn)
4 format(1X,5e15.5)


Est ce que vous auriez une solution svp

Merci de votre aide

NB les fichiers sont en pieces jointes

Bonne journée

[genteur slayer]

l'idée c'est de parcourir plusieur fois ton fichier texte:

Code :

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

1
2
3
4
5
6
7
8
9
10
 
character(len=4) ::toto
...
open(unit=40,file="plot.out",status="old",position="rewind",iostat=ios)
NZONE=0
do while(ios==0)
  read(unit=40,iostat=ios,fmt=*)toto
  if (toto=="ZONE") NZONE=NZONE+1
end do
...

à partir de là, tu connais le nombre de zone!
tu peux donc réserver les tableaux avec la bonne taille, rembobiner le fichier et le lkire à nouveau avec toutes les informations.

[Sylvain Bergeron]

Autre piste : lire une ligne à la fois sous forme de chaîne de caractères, examiner le contenu, puis la traiter.

Ça donne quelque chose comme :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
character (100) Ligne
...
nd = 1
do
   read(90,'(a)',iostat = k) Ligne
   if (k /= 0) exit
   if (Ligne(1:4) /= 'ZONE') then
      read(Ligne,*) x(nd), z(nd), nmoleF(nd), So(nd), Sw(nd)
      nd = nd + 1
   endif
enddo

Entre la méthode proposée par genteur slayer et celle-ci, il n'y en a pas une de meilleur a priori. C'est vraiment une question de contexte.

[jal13]

Merci beaucoup je vais essayer ca de suite

bonne journée à vous