 |
INFORMACIÓN IMPORTANTE ...
La entrega de esta tarea debería ser en un documento en el formato
.txt (o .sh, (o .m)), con su
nombre
en el
nombre del archivo (por ejemplo, m_miller_gmt_tarea2_preg1.sh). Para
entregar la tarea, mandar
el
archivo por email a ambos profesor y ayudante. Fecha de entrega:
viernes 17 de mayo a las 23:59
horas.
Importante:
Deben comentar sus documentos.
Esta tarea está basada sobre los mismos
archivos de datos de la tarea 1:
Tarea 2, Pregunta 1:
MAULE-aftershocks-catalogue-lange-etal2012.dat
El cual es un catálogo de sismicidad de las réplicas del terremoto de
Maule, 2010, con 20189 réplicas
registradas entre las 2010-03-15 y las 2010-09-30.
Escribe un script para obtener el momento sísmico acumulativo para
cada día en el archivo (es decir, la suma de los momentos sísmicos por
día de todos los eventos desde el punto de inicio del catálogo,
hasta el día de interés).
Un ejemplo de momento acumulativo esta en este imagen:
http://www.iris.edu/hq/gallery/d/1526-2/CumulativeMomentMag.jpg
La relación entre el magnitud de momento, Mw, y el momento sísmico, M0, (en Nm) es:
|
y en términos de awk, para convertir la columna ($15) de Mw al M0 para un evento individual:
awk '{print 10^(1.5*$15 + 9.1)}'
Use julday para convertir la fecha al formato: año jdy.
Cuando corres el script; debes obtener algo así:
2010 074 2.579e+15
2010 075 1.161e+16
2010 076 4.709e+16
etc ...
Guarde los datos en un archivo llamado data.txt. Se puede usar el script de octave plot_data_t2.m para graficar los datos:
octave plot_data_t2.m
Tarea 2, Pregunta 2:
sst.mnmean.1996.1997.1998.txt.tgz
Es un catálogo (global)
de la temperatura de la superficie del mar (sst) - el promedio mensual,
entre los meses de enero de
1996 y de diciembre de 1998.
Escribe un script que trabaja con 5 argumentos (año-mes; lat1; lat2;
lon1; lon2) para aislar los
datos de la temperatura de la superficie del mar para un cierto mes,
dentro de un cierto rango de
latitudes y longitudes, y que escribe las columnas de latitud,
longitud, temperatura (x, y, z)
dentro de un nuevo archivo con un nombre asociado con el año-mes que
corresponde.
Use el comando case para relacionar el tiempo dado en la segunda columna con el año-mes que representa.
This document was generated using the LaTeX2HTML translator Version 2008 (1.71)
Copyright © 1993, 1994, 1995, 1996,
Nikos Drakos, Computer Based Learning Unit, University of Leeds.
Copyright © 1997, 1998, 1999,
Ross Moore,
Mathematics Department, Macquarie University, Sydney.
The command line arguments were:
latex2html -split 0 2013_gmt_tarea2.tex
The translation was initiated by matt on 2013-05-06