Preguntas frecuentes sobre ZIVIS (Blog de Zivis)

plot de una trayectoria

20070430140200-gnuplot.jpg

El primer paso para ver una trayectoria es encontrar un fichero lgv3t_conTr_C0_xxxxx_5 que este terminado pero no lo hayais mandado aun. Un truco es detener comunicaciones mientras se enreda con esto para asegurarse que el fichero no lo acaba enviando boinc y lo perdemos a mitad del juego.

El segundo paso es conseguir que el formato pase a ser ascii. El fichero es una secuencia de vectores (X,Y,Z,v) en floats (de 4 bytes, por tanto 16 bytes por dato) y ademas tiene un vector cero cada 1000 pasos. Un truco en unix, aprovechando que la linea es de 16 bytes, es usar directamente od -An -tf4 lgv3t_conTr_C0_30-06-27-15-84_0_5 | grep -v 0.000 Ya lo siento Razorblade, pero desde windows tendreis que buscaros algun otro comando o escribir un programita de pasar binario a texto.

El fichero texto producido se puede poner directamente en gnuplot. De hecho en unix no se necesita ni usar un fichero de texto intermedio. Podemos llamar a gnuplot y dar por ejemplo la orden

gnuplot> splot ’< od -An -tf4 resultados/lgv3t_conTr_C0_30-06-27-15-84_0_5 | grep -v 0.000’ using 1:2:3

y esa linea hace todo el trabajo! Ademas, es posible rotar la imagen con el raton.

30/04/2007 14:02 #.

Comentarios » Ir a formulario

gravatar.com
Autor: Razorblade

"desde windows tendreis que buscaros algun otro comando"

Bueno, para cosillas así se puede tirar de Cygwin, ¿no?

Fecha: 30/04/2007 16:46.


gravatar.com
Autor: Razorblade

Para el que lo quiera probar bajo Windows, he subido los archivos necesarios a: http://rapidshare.com/files/28754334/Cygwin.rar

El procedimiento es similar, y tendréis que hacerlo desde el Símbolo del Sistema (lo ideal es que copiéis el archivo lgv3t_conTr_C0_xxxxx_5 a la carpeta Cygwin para no liarnos con las rutas):

od -An -tf4 lgv3t_conTr_C0_xxxxx_5

Fecha: 30/04/2007 19:46.


gravatar.com
Autor: Razorblade

grep -v 0.000 >> salida.dat

El archivo salida.dat es el que se puede abrir con GNUPLOT:
gnuplot> splot salida.dat using 1:2:3


PD: la primera línea de este mensaje es continuación de la última línea del anterior, y entre medio hay que poner la barra vertical (parece que el blog la toma como fin del mensaje).

Fecha: 30/04/2007 20:26.


gravatar.com
Autor: Razorblade

Disculpad, lo he vuelto a subir porque se me olvidó meter un archivo: http://rapidshare.com/files/28810702/Cygwin.rar

Ah, y un par de correcciones:
-No hace falta poner dos >, con uno basta: grep -v 0.000 > salida.dat
-Se me olvidaron las comillas en el nombre del fichero para gnuplot: splot 'salida.dat' using 1:2:3

Fecha: 01/05/2007 02:16.


gravatar.com
Autor: José Luis

Habrá casos en que os salgan trayectorias con una pinta "extraña", que no se parezcan a la que ha subido Alejandro.

Se trata de trayectorias que:

a) Han escapado en seguida, y no ha dado tiempo a que den ni una vuelta a TJ-II.

b) Han rebotado varias veces en espejos magnéticos (allí donde el campo magnético es más intenso, cerca de una bobina) y aunque hayan durado más tiempo, tampoco completan una vuelta entera.

En esos casos, para ver bien la imagen y hacerse una idea de lo que ha pasado, conviene hacer, antes del splot,

set xrange [-2:2]
set yrange [-2:2]
set zrange [-0.5:0.5]

Dibujar con líneas en lugar de con puntos

splot 'salida.dat' using 1:2:3 w l

también ayuda.

Fecha: 01/05/2007 12:30.


gravatar.com
Autor: Alejandro Rivero

Gracias Razorblade y Jose Luis!

Lo del "espejo magnetico" esta en una de las trasparencias de la conferencia, ¿no? A ver si la subimos a la web mañana.

Fecha: 01/05/2007 19:33.


Añadir un comentario




No será mostrado.




BIFI - CIEMAT - Ay. de Zaragoza / © Ayuntamiento de Zaragoza / webmunicipal@zaragoza.es / Mapa Web / Miembro W3c / XHTML 1.0 / CSS 2.0 / Accesibilidad