Figuras e animações

De Física Computacional
Revisão de 21h51min de 18 de abril de 2023 por Gustavobopsin (discussão | contribs)
(dif) ← Edição anterior | Revisão atual (dif) | Versão posterior → (dif)
Ir para navegação Ir para pesquisar
reset
set grid
set xlabel 'x'
set ylabel 'y'
set xrange[0:25]
set yrange[0:25]
set size ratio 1

#Inicio das posições
set term png
set output 'inicializacao.png'
plot 'posicoes.dat' w circles fill solid 0.6 border 1 lc 1 t ''

#Animações
set key out
set term gif animate delay 10
set output 'animacao1.gif'
do for [i=9900:10000] {plot 'dados_particulas.dat' index i u 3:4:(1.0) w circles fill solid 0.6 border 1 t 'tempo='.i}

set term gif animate delay 50
set output 'animacao2.gif'
do for [i=9900:10000:5] {plot 'dados_particulas.dat' index i u 3:4:(1.0) w circles fill solid 0.6 border 1 t 'tempo='.i}
reset

set grid
set xlabel 't'
set ylabel 'E'
set key out
set term png size 800,500
set title 'Energia mecânica: 124.138'
set output 'energia_tempo.png'
plot 'dados_termo.dat' u 1:2 w l lc 6 t 'Energia cinética' , \
     'dados_termo.dat' u 1:3 w l lc 7 t 'Energia potencial' , \
     'dados_termo.dat' u 1:($2+$3) w l lc 1 t 'Energia mecânica' 

unset title     
set xrange [0:300]     
set output 'energia_tempo_zoom.png'
plot 'dados_termo.dat' u 1:2 w l lc 6 t 'Energia cinética' , \
     'dados_termo.dat' u 1:3 w l lc 7 t 'Energia potencial' , \
     'dados_termo.dat' u 1:($2+$3) w l lc 1 t 'Energia mecânica' 

unset xrange
set title 'Temperatura média: 0.7237'
set ylabel 'T'
set output 'temperatura_tempo.png'
plot 'dados_termo.dat' u 1:4 w l lc 7 t ''
unset title
set xrange [0:300]
set yrange [0.1:0.8]
set output 'temperatura_tempo_zoom.png'
plot 'dados_termo.dat' u 1:4 w l lc 7 t ''

unset yrange
unset xrange
set title 'Pressão média: 0.1877'
set ylabel 'P'
set output 'pressao_tempo.png'
plot 'dados_termo.dat' u 1:5 w l lc 6 t ''
unset title
set xrange [0:300]
set output 'pressao_tempo_zoom.png'
plot 'dados_termo.dat' u 1:5 w l lc 6 t ''

unset xrange
unset yrange
set title 'Número de bins = 200'
set ylabel 'g(r)'
set xlabel 'r'
set output 'fdr.png'
plot 'dados_fdr.dat' u 1:2 w l lw 2 lc 8 t ''
reset