eliocamp / tesis

Repositorio para mi tesis de licenciatura.

Repositorio para poner mi tesis de licenciatura. Nada demasiado intersante por ahora.

Datos

Por ahora hay algunos datos (en la carpeta DATA) pero git me cortó el chorro gratuito por pasarme del giga. Así que la idea es, por ahora, subirlos a esta carpeta de Mega.

tesis's Issues

Corridas

Hacer corridas sacando los modelos de suelo y hielo.

• Sin hielo
• Sin suelo (ni hielo)

Composición de campos

Hacer composición de campos de geopotencial en 700 y en 200/300 para casos de onda 3 'activos'. Considerando 'activos' como el percentil 66 de cada mes para la amplitud media de La onda en los niveles entre 700 y 100.

Convergencia de humedad

Buscar variable de vigencia de humedad en varias bajas. Debería tener relación con las correlación negativa con la pp que encontré a nivel estación.

Problemas en stat_fill_contour()

En ncep.mean septiembre anomalía de temperatura en 500 y 300.

ERA

Usar datos de ERA. Para probar la sensibilidad de los resultados al reanalisis utilizado y porque es un mejor producto que NCEP.

Hosting para los datos

Originalmente (#4) la idea era usar git-lfs para hostear los datos en github, pero tiene un límite de 1Gb gratuito, así que tengo que conseguir otro sistema. Opciones:

• Hostearlo yo en mi "server".
  Está bueno porque tengo todo el control sobre los datos y no tengo que preocuparme porque se suba. Además, no tengo prácticamente límite de espacio. Lo malo es que dependo de tenerlo siempre corriendo y quizás no es una opción a largo plazo. Además bajar los datos seguramente es un parto con la pésima velocidad de upload que tengo.
• Torrent
  Sería otra forma de hostearlo yo pero que reduce un poco las preocupaciones técnicas y de uptime. El principal problema es que no es dinámico y no es fácil agregar nuevos archivos o modificar los exsitentes.
• Mega
  Es lo que estoy viendo ahora. Permite 50Gb gratis. Lo malo es que hay que subirlos a mano o usar un cliente de sincronización que no tengo ganas de instalar en la compu.
• Figshare
  No lo conozco.
• Otras.
  ¡Consultar!

Composiciones/regresiones speedy

Faltan hacer composiciones o regresiones para las distintas corridas de speedy.

• Control
• CLIM
• NOICE
• NOLAND

Usar Rds en vez de .Rda o .csv

Salvar los datos en .Rds (con compresión gzip, ver benchmark) ya que es como un .Rda pero que permite cargar y asignar el resultado a una variable. Evitar .csv porque, si bien son fáciles de manejar, son muy pesados para usar en el repo.

• Cambiar NCEP a Rds
• Cambiar Speedy a Rds
• Subir archivos al repo (ver git-lfs)
• Crear una función que pase de .Rds a .csv para replicabilidad.

Ver acá.

Corrida con sst zonal

Corrida de speedy usando condiciones de sst climatológicas y zonalmente simétricas

• Cambiar el archivo de las condiciones de contorno para que la SST sea sólo la componente zonal.
• Correr el modelo.
• Ver qué onda.

Correlación espacial QS3 y GH

• Usar R^2 en vez de la correlación.
• Sombrear a partir de un valor.

URLs para referencias que son página web (como metR)

Presentaciones

Explorar cómo hacer presentaciones en powerpoint directamente desde R. En 4fe2bda está hecho con ioslides, pero mejor aún sería poder mandarlo directo a Power Point, así es compatible con lo que hacen los demás.

Función para convertir longitud

Hacer una función que convierta la longiutd entre 0:360 y -180:180.

Flujos de actividad de onda

Pasar el script para calcular los flujos a R y calcularlos.

Hacer un outline de la tesis

Agregar xwrap

Todas las figuras que lo requieran deben tener xwrap = c(0, 360).

Referencia para SVD

Necesito una referencia bibliográfica para justificar que componentes principales utilizando SVD es más estable que utilizando autovalores de XX'.

¿Qué es la fase?

Aprender a interpretar la fase... es un issue un poco open ended.

Arreglar el script de interpolación

El script para interpolar NCEP a la grilla de SPEEDY no anda bien y manda cualquier fruta. El problema es el orden en el que guarda los datos, seguramente, que se hizo quilombo al usar foreach().
Hay que arreglar eso. Quizás volviendo a un for() de toda la vida que no corre en paralelo pero al menos funciona 😞

Interacción entre ondas, SAM, ENSO

Ver qué relación hay entre las ondas 1-2 que acompañan la onda 3 y con el SAM y la señal de El Niño 3.4.

Nombres de variables

Decidir cómo nombrar a las variables en el texto y en los gráficos.

Labels para los contornos

Crear un geom_contourlabel() para marcar los contours.

Obtener datos

Datos para usar

• GH
• U
• V
• SLP
• PP
• OLR
• T diaria en 850hPa
• v diaria en 850hPa
• u diaria en 300hPa
• v diaria en 300hPa

Calcular:

• $\nu f + \zeta$
• u'v' en 300hPa
• v'T' en 850hPa

Normalización

Normalizar la amplitud con respecto a la amplitud media del mes??

filled contour

Hacer contornos llenos que funcionen de verdad.

geom_contourlabel() cuando hay pocos contornos.

Mejorar ReadNC

Hacer que la función sea genérica y descubra automáticamente cómo convertir la coordenada temporal. Esto va a ahorrar muchos dolores de cabeza.

Fechas usadas para composición

Tabla (o gráfico) con la fechas usadas para cada composición.

Validación de las corridas

Hacer una pequeña validación para mostrar que las corridas están bien hechas.

• Mapa de $\sigma$ para un mes y cada corrida (ver que lo que tiene que ser 0, sea 0).
• Mostrar una serie temporal en distintos puntos (uno sobre el océano, otro sobre el suelo, otro sobre el hielo).

Regresiones

• Agregar las regresiones de gh en amplitud y fase para la corrida control.
• Usar los regresores normalizados.
• Probar cómo se ve haciendo la regresión simple en vez de múltiple.
• Proyectar también en mercator (u otra hemisférica).

Análisis climatológico

Muchas cosas y muy abierto...

Poner el mapa en el gráfico 13

El gráfico de x = fase, y = amplitud, color = selección.

Quizás hacer un mapa lat, lon (fase) y tamaño = amplitud, color = selección. Con un alpha para que se vea dónde hay overplotting.

Flechas

Usar geom_segment() para las flechas no sirve cuando se usan proyecciones ya que se deforman (por ejemplo usando mercator en altas latitudes). Hay que hacer algo mejor.

Explicar Antártida

Aclarar en las figuras con la antártida que el coso blanco es eso.

Gran figura

La gran figura:

Para cada mes, mostrar la onda teórica y el campo observado para NCEP y CONTROL. Agregar una indicación de la actividad de SAM y ENSO.

Ciclo anual speedy

Agregar el ciclo anual de la onda 3 para la corrida control y la climatológica.

Hacer corrida con sst climatológica

• Modificar los archivos de anomalías para que sean 0 Se cambia un parámetro y usa climatológica.
• Correr el modelo
• Fijarse que esté todo bien.
• Hacer análisis de onda 3

Arreglar presentación

Está medio rota la presentación.

Regresiones

Agregar regresiones de campos en cada índice

• Geopotencial (300 HPa)
• OLR
• SST

Ondas teóricas

Figuras de las ondas teóricas (o sea, lo que da Fourier) para cada mes usado en las composiciones, junto con el promedio.

Cómo mostrar la amplitud de onda 3

La media de la amplitud o la amplitud de la media.

En el primer caso, obtengo una climatología de la amplitud independientemente de la fase. Es compatible con el índice que uso pero no da una climatología correcta de la onda 3 estacionaria.

Como mínimo tiene que haber una discusión sobre este dilema y luego elegir una u otra forma y hacerlo consistente.

Corregir el código para que siga la guía de estilo

El código como está ahora es un desastre y no sigue la guía de estilo que elaboré. Tengo que ir y limpiarlo asegurándome de no romper nada.