LINUX GMT - EJEMPLO 1
De
esta página web - http://www.esrl.noaa.gov/psd/data/gridded/data.cmap.html
se puede obtener el archivo precip.mon.ltm.nc
que contiene información sobre la precipitación global. A cada punto
(longitud, latitud) los valores de precipitación son el promedio diario
para los doce meses del año (datos tomados entre los años 1979 - 2000).
Existe un punto de datos cada 2.5 grados latitud y longitud.
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 aptitude install nco
Ahora, para crear un archivo .txt de los datos, use el comando:
ncks precip.mon.ltm.nc
> precip.mon.ltm.txt
Acá esta el archivo generado - precip.mon.ltm.txt
- contiene líneas de la forma
time[11]=8016 lat[71]=-88.75
lon[139]=348.75 precip[124411]=0.473636 mm/day
Que tienen la longitud y latitud de la grilla, y el promedio de la
precipitación (mm por día) para cada mes del año:
time[0]=0 hours (promedio para enero)
time[1]=744 hours (promedio para febrero)
time[2]=1416 hours (etc ...)
time[3]=2160 hours
time[4]=2880 hours
time[5]=3624 hours
time[6]=4344 hours
time[7]=5088 hours
time[8]=5832 hours
time[9]=6552 hours
time[10]=7296 hours
time[11]=8016 hours (promedio para diciembre)
SCRIPT 1 PARA GRAFICAR LOS DATOS:
Script - ejemploA1.sh
Paleta de colores - GMT_meanprecip.cpt
Instrucciones: Bajar el script y la paleta de colores a un directorio que contiene los datos. Cambiar el variable ${topografia} en el script para dar el camino a tu topografia topo_6.2.img. Correr el script.
Se puede obtener la topografía aquí: http://mttmllr.com/topo.htm
Ejemplo imagen, por la región de sudamérica:
SCRIPT 2 PARA GRAFICAR LOS DATOS:
Script - ejemploA2.sh
Instrucciones:
Bajar el script a un directorio que contiene los datos. Cambiar los
variables ${topografia}, ${paleta_topo} y ${paleta} en el script para dar los caminos a tus archivos. Correr el script.
Se puede obtener la topografía aquí: http://mttmllr.com/topo.htm
Ejemplo imagen por la región de sudamérica (mes 7):
