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
Partager