Discussion :

[Le_Boucher]

Bonjour à tous,

J'ai fais une simulation de transport d'humidité dans un materiau poreux par fortran, et j'ai imprimé mes résultats.
Je cherche à réaliser une animation, j'entends par la une succession d'image, qui montrerai l'évolution de mon humidité en fonction du temps, dans l'espace.

J'ai plusieurs blocs de données rangée de la manière suivante :

Code :

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

1
2
3
4
5
6
7
#bloc     i
zk  xk  yk
...
 
#bloc     i+1
zk  xk  yk
...

Une fois sur GNUPLOT, j'execute tout d'abord ce script :

Code :

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

1
2
3
4
set term gif animate
set output 'animation.gif'
imax=23
i=0

Puis ensuite, j'execute celui ci : (Vous l'aurez compris, mon fichier avec les données s'appelle "DATA.dat" )

Code :

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

1
2
3
4
5
6
7
set pm3d
set xrange [0:0.02]
set yrange [0:0.01]
splot "DATA.dat" with pm3d at s ; index i
pause 0.2
i=i+1
if (i<=imax) reread

il m'affiche :

Code :

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

1
2
3
gnuplot> splot "DATA.xls" with pm3d at s ; index i
                                         ^
"datagnu.gnu", line 4: All points y value undefined

(le chapeau sous le point virgule)


J'ai littéralement débuté aujourd'hui en GNUPLOT, et je suis assez fier de moi d'avoir réussi à imprimer certains résultats de la manière que je voulais, mais je n'arrive pas à réaliser d'animation!
Merci d'avance de vous pencher sur la question, et veuillez pardonner mes fautes d'orthographes je ne suis pas très bon !

[Flodelarab]

Bonjour

Une ligne blanche : rupture dans la courbe.
Deux lignes blanches : autre bloc.

Ici, ton fichier d'entrée ne doit pas être le premier mais le second.

Code :

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

1
2
3
4
5
6
7
#bloc     i
zk  xk  yk
...
 
#bloc     i+1
zk  xk  yk
...

Code :

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

1
2
3
4
5
6
7
8
#bloc     i
zk  xk  yk
...
 
 
#bloc     i+1
zk  xk  yk
...

De plus, tu n'as pas respecté la syntaxe du mot "index". Après le point virgule, point de salut.

[Le_Boucher]

Merci d'avoir répondu!

Pour le saut de ligne, justement je cherchais à faire une rupture, ce que je veux c'est qu'il me défile les différentes surfaces une par une, parce que c'est une surface qui évolue au cours du temps.

pour "index" du coup je sépare avec une virgule?

[Flodelarab]

Ce n'est pas la question.
À ta place, je ferais:

Code :

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

splot "DATA.dat" index i

[Le_Boucher]

Finalement j'ai réussi ! voici les scripts utilisés !

dans mon code fortran je fais :

Code :

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

    CALL EXECUTE_COMMAND_LINE( "gnuplot 'format.gnu'")

dans format.gnu je fais :

Code :

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

1
2
3
4
5
set term gif animate
set output 'animation.gif'
imax=22
i=0
load 'datagnu.gnu'

dans datagnu.gnu je fais :

Code :

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

1
2
3
4
5
6
7
8
9
set xrange [0:0.02]
set yrange [0:0.01]
set zrange [0:0.606]
set dgrid3d 30,30
set hidden3d
splot "DATA.dat" index i
pause 1
i=i+1
if (i<=imax) reread

Et j'obtiens bien un gif animé cohérent avec ma simulation !! Le voici ci joint !



Sachant que mes données sont bien écrit bloc par bloc et séparé par deux saut de ligne comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
#BLOC         0
0.200000000E-01	0.800000000E-02	0.600000000
0.200000000E-01	0.900000000E-02	0.600000000
0.200000000E-01	0.100000000E-01	0.600000000
...
 
 #BLOC           1
 0.00000000    	 0.00000000    	0.443667971
 0.00000000    	0.100000000E-02	0.446617914
...

En tout cas merci pour les réponses!!! À bientôt !