Introducción a RClimDex Sarah Perkins Centre of Excellence for Climate System Science UNSW, Australia
Por qué usar R? •R es un software estadístico http://www.r-project.org/ que se puede descargar e instalar gratuitamente! •RClimDex es una interfaz gráfica (Graphical User Interface, GUI) escrita en R, que realiza un control de calidad básico a datos climáticos diarios. •Otra GUI adicional, RHtest, identifica posibles inhomogeneidades. •También hay una versión disponible en Fortran (FClimDex) •Todo el software está disponible en el sitio web de ETCCDI http://etccdi.pacificclimate.org/software.shtml
Formato de los datos para RClimDex •Archivo de texto ASCII •Columnas deben seguir la sequencia: Año, Mes, Día, PRCP, TMAX, TMIN. (NOTA: unidades de PRCP = milímetros, unidades de Temperatura = grados Celsius) •El formato descrito arriba debe ser delimitado por espacios (e.g. cada elemento separado por uno o más espacios). •Para registros de datos, datos faltantes deben estar codificados como -99.9; registros deben estar en orden de fecha calendario. Puede haber fechas faltantes.
Formato de datos para RClimDex Year month day
Precip
Tmax
Tmin
Para comenzar •Asegúrese de tener R y Rclimdex instalados ☺ •Abra R •Escriba: source(“rclimdex.r”) , O Seleccione del menu
Para comenzar •Aparecerá la ventana principal •Tanto esta ventana como R DEBEN estar abiertas mientras se calculan los índices
Para comenzar En el directorio donde guardó RClimDex, el programa generará: •Directorio Log •Guarda gráficas (plots) y archivos .csv que contienen el tipo y la localización de datos erróneos.
•Directorio de gráficas (plots) •Guarda todas las gráficas para cada índice calculado
•Directorio Indices •Guarda .csv para cada índice calculado.
•Archivo Stationameindcal.csv •Los datos luego del control de calidad que serán usados por RClimDex para calcular los índices
Usando RClimDex •RClimDex NO calculará ningún índice antes de realizar control de calidad (QC). •Algunos datos “sospechosos” son reemplazados autoáticamente como faltantes.
•Directorio “log” contiene la localización de datos fuera de rango (outliers) que QC requerirá que verifique
RClimDex y RHtest • Desarrollado y mantenido por el Climate Research Branch del Servicio Meteorológico de Canadá • Pruebas de QC de RClimDex incluyen Tmax < Tmin, -ve precipitación, pruebas de outliers •RHtest es un modelo de regresión de doble fase con una tendencia lineal para la serie completa •También se espera que los participantes hagan juicios de valor basados en sus conocimientos del tiempo y clima de la región, y en su acceso a los metadatos
Metadata confirms station relocation
Cómo funciona el control de calidad (QC) en RClimDex?
•Verifica que precipitación no sea negativa •Verifica que Tmax no sea menor que Tmin •Identifica outliers en Tmax, Tmin y DTR, definidos como valores por fuera de un valor dado por el usuario i.e. • [mean – n*std, mean+n*std] • “std” & “mean” son la desviación estándar y el promedio para el día calendario • “n” es un input del usuario
Gráficas también pueden identificar outliers Sample output from QC log file
Período frío real, o un error de digitación? Conocimiento del clima local y metadatos podrían ayudar
El papel SUYO en control de calidad (QC) •Juzgar si valores identificados como sospechosos son reales o no •Remover datos sospechosos o corregir según el caso •Es buena idea copiar el archivo original ANTES de hacer correcciones! •Una vez terminado el control de calidad vuelva a correr RClimDex para verificar si hay más datos sospechosos
Calculando indices en RClimDex •Una vez satisfecho de que todos los datos sospechosos han sido eliminados, podemos calcular algunos índices! •Calculados en datos controlados (QC) (así que asegúrese de haber guardado los datos originales antes de usar RclimDex)
Parámetros definidos por el usuario
•Le serán pedidos algunos datos importantes para calcular los índices. •Primer y último año del período base •Período de referencia para muchos índices (usualmente 30 años) •Debe estar totalmente cubierto por sus datos •Latitud y longitud de la estación •Importante para calcular las estaciones •Umbrales superiores/inferiores de Tmax/Tmin (ºC) y de precip (mm/día) •Cierta flexibilidad en los índices •Umbrales apropiados para sus datos •Utilizados para contar número de días por encima o debajo de los umbrales dados.
Seleccionando índices Escoja cuáles índices quiere calcular!
Indices entregados
•Archivos .csv (en el directorio indices) • gráficas (el el directorio plots) •Repetir para más estaciones! •Es necesario realizar un nuevo control de calidad (QC) cada vez que una nueva estación es analizada
Indices entregados
Muchas gracias por su atención!
Preguntas?
[email protected]