Objetivo n1: Gráfica de isolínea con datos repartidos en una malla irregular:

1.

Ejecute el programa: IRR_GRID.PRO

2.

El programa pregunte si quiere seguir o acabar, conteste con '' a'' para acabar la ejecución.
Hay 2 ventanas de gráfica, n 0 y n 1. Deberíamos ver la ventana n1, si no es el caso, ejecute el comando:
wshow,1

3.

Escribe en la pantalla los valores de cada punto con su posición. Para eso, ejecute el comando:
for i=0,24 do print,i,x(i),y(i),z(i)

4.

Ahora, mire la posición del cursor utilizando el programa: CURS_POS.PRO

5.

Intente ver en que posición corresponde al punto elegido, y compare el valor de la gráfica con los datos.

6.

Estudio del programa:

(a)

La copia del programa tiene unos comentarios que nos ayudaran a comprender cada paso.

(b)

Además, hay espacio para poner más comentario si se quiere.

(c)

Mire el tipo de variables utilizadas. Para eso, ejecute el comando:
help,x,y,z,xx,yy,zz

(d)

Se ve que zz es una matriz con 2 dimensiones. Esto es necesario porque los comandos contour, surface, shade_surf trabajan con matrices de 2 dimensiones.

(e)

En el programa, se especifica que los 11 niveles de isolíneas son repartidos entre:
0 y max(z) para la gráfica en la ventana n0
min(z) y max(z) para la 2nda gráfica.
La diferencia viene de los comandos
zz(where(zz eq 0))=1000
.....
zmx=max(z,min=zmn)
contour, zz,...,max_value=zmx+0.1
Esto, porque IDL asigna un valor 0 para los puntos donde no hay dato. Para asegurarse del valor de punto donde no hay datos, puede ejecutar el programa
IRR_SURF.PRO

(f)

También, observe el cambio de los valores de los ejes. En la segunda gráfica, los ejes son definidos ejecutando las líneas de comando:
xx=findgen(...
yy=findgen(...
que crean 2 vectores con lo mismo número de elemento que las dimensiones de la matriz zz.