GMT EN GEOFISICA.

TRABAJO FINAL 2012. VELOCIDAD DEL VIENTO - PATAGONIA

Certamen para entregar jueves 5 de julio (o antes). No hay flexibilidad en la fecha de entrega.

Se puede trabajar en este trabajo en las clases y podemos responder a preguntas. Gente que no trabajan en esta tema en las clases tienen que explicar oralmente como funciona su trabajo, y responder a preguntas, cuando entregan el trabajo y su nota dependera en eso.

En la pagina web siguiente - http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.derived.surface.html - se puede encontrar  el archivo wspd.sig995.mon.ltm.nc que contiene informacion sobre la velocidad del viento global. El archivo es en el formato NetCDF.

Para extraer datos de un archivo .nc, lo más facil es instalar NCO - http://nco.sourceforge.net/ - en ubuntu se puede hacer eso con el comando:

sudo apt-get install nco

Ahora, para crear un archivo .txt de los datos, use el comando:

ncks wspd.sig995.mon.ltm.nc > wspd.sig995.mon.ltm.txt

Para ellos que estan trabajando en los computadores de la facultad, aca esta el archivo de salida: wspd.sig995.mon.ltm.txt

En cada punto (longitud, latitud) el valor de la velocidad del viento es el promedio mensual, para las doce meses del agno (datos tomados entre 1981 y 2010). Existe un punto de datos cada 2.5 grados. El archivo generado contiene lineas de la forma:

time[1]=31 lat[52]=-40 lon[23]=57.5 wspd[18023]=7.35 m/s
 
Que tienen la longitud y latitud de la grilla, y el promedio de la velocidad del viento (m/s) para un cierto mes del agno:

time[0]=0 (promedio para enero)
time[1]=31 (promedio para febrero)
time[2]=59 (etc ...)
...
time[11]=334 (promedio para diciembre)

1) Escribe un script para graficar estos datos, para cada mes, para la region de Patagonia-Sud Pacifico en lo mejor forma posible. Recuerden comentar sus scripts. Mandeme su script, la(s) paleta(s), un .pdf de un grafico para uno de los meses, antes del limite de la fecha de entrega.
En su mapa, piensa en:
La region, ejes/bordes bien hecho, la costa Chilena, la frontera Chile-Argentina, la ciudad de Punta Arenas marcado en el mapa, representacion de topografia, una escala para la paleta usada para la velocidad, un titulo para el grafico que menciona el mes. Mas credito sera dado para cosas adicionales que pongan en sus mapas - se puede ver ejemplos cientificos ya hechos en linea por ejemplo.
En su script, piensa en:
Una opcion para elegir el mes que uno se quiere graficar, comentarios, la definicion de variables como lo de la paleta para poder modificar las opciones del script rapidamente si quieres.
Por favor, pongan sus nombres en el nombre de sus archivos.
Prefiero recibir archivos como miller_trabajo_final.sh y no 25 archivos que se llaman el mismo.

(opcional, si tiene su 12 imagenes para cada mes y quiere juntarles en un .gif, se puede usar el comando "convert" de imagemagick (ya debe ser instalado) por ejemplo)

convert -alpha off -delay 100 -loop 0 -verbose -trim +repage wspd_1.ps wspd_2.ps wspd_3.ps wspd_4.ps wspd_5.ps wspd_6.ps wspd_7.ps wspd_8.ps wspd_9.ps wspd_10.ps wspd_11.ps wspd_12.ps wspd.sig995.mon.ltm.gif

2) En la pagina web, tambien existe los valores para el promedio del viento horizontal (uwnd.sig995.mon.ltm.nc) y vertical (vwnd.sig995.mon.ltm.nc). Escriben unas oraciones sobre como se puede agregar a la mapa simbolos que representan la direccion del viento.