Linux GMT. Trabajo Final, 2011.

Graficando datos geofísicos

Fecha de recepción: al final del semestre (antes del 6 Oct)

En este trabajo, se puede mandar respuestas en la forma de texto (.txt, .cpt, .sh etc.) o en la forma de .pdf
Con el comando ps2pdf  se puede facilmente convertir un archivo .ps a un archivo .pdf


El siguiente archivo - velocidad_intersismica.txt - contiene para 53 estaciones de remedición de la red SAGA (GFZ-Potsdam) los vectores de velocidad intersísmica (datos de GPS), en milímetros al año con respecto al ITRF2005, en las columnas 4 (este-Ve05), 5 (norte-Vn05) y 6 (vertical-Vu05) con sus errores respectivos (columnas 7-9). La velocidad con respecto a un Sudamérica fijo esta en las columnas 12 (este-Ve-sam) y 13 (norte-Vn-sam).

Su trabajo: Grafique esta información (en 1 gráfico) para la región de Chile entre Concepción y Puerto Montt.

Cosas que considerar son:
(i) Una buena basemap,
(ii) Cuidades de Concepción, Valdivia y Puerto Montt marcadas,
(iii) Costa de Chile, y borde Chile-Argentina marcados,
(iv) Topografía, con una buena paleta,
(v) Vectores para representar la velocidad horizontal (norte, este) de los sitios de GPS,
(vi) Una escala para los vectores, con unidades.

Cosas adicionales son:

(i) Una manera de representar la velocidad vertical.
(ii) Encontrar una manera de representar los errores en las velocidades horizontales,
(iii) Dar una opción en el script, para poder elegir si quieres graficar columnas 4 y 5, o columnas 12 y 13, para representar el campo de velocidades,
(iv) Representar el error en la velocidad vertical (creo que eso es casi imposible).

Mandenme todos los scripts, paletas, comentarios que tienen (en formato de texto) el gráfico(s) (en formato de pdf). Es importante que ustedes pongan sus nombres en el nombre de sus archivos para que sabemos de quien son. También es necesario hacer una "entrevista" cuando entregas su trabajo para sabemos que entienden el trabajo que hicieron. Se puede hacer esto por teléfono.