RClimDex (1
Contents
1. Lista de ndices Clim ticos Formato de los Datos de entrada Definici n de los ndices C lculo de Los l mites y periodo base de los ndices de temperatura Preguntas Frecuentes sobre R para Windows 1 Introducci n ClimDex es un programa basado en Microsoft Excel que proporciona un paquete computacional f cil de usar para el c lculo de ndices de extremos clim ticos para monitorear y detectar cambio clim tico Fue desarrollado por Byron Gleason del National Climate Data Centre NCDC de NOAA y ha sido usado en talleres CCI CLIV AR sobre ndices clim ticos desde el 2001 El objetivo original fue el de colocar ClimDex en un ambiente que no dependa de un sistema operativo determinado Fue muy natural usar R como nuestra plataforma dado que R es gratis pero al mismo tiempo es un software muy robusto y poderoso para gr ficos y an lisis estad stico Puede ser ejecutado tanto en ambientes Windows como Unix En el 2003 se descubri que el m todo utilizado por ClimDex y otros programas para calcular ndices de temperatura basados en percentiles resultaba en inhomogeneidad en las series de ndices Una soluci n a este problema requiere un procedimiento bootstrap que se hace imposible de implementar en un ambiente Excel Esto ha hecho a n m s urgente la necesidad de desarrollar este paquete basado en R RClimDex 1 0 est dise ado para proporcionar una interfase amigable para calcular ndices de extremos clim ticos2. 2004 Avoiding inhomogeneity in percentile based indices of temperature extremes J Climate submitted 21 Anexo E Preguntas y Respuesta Frecuentes sobre R para Windows A Adobe Acr obat Docun 22
3. La m xima temperatura m xima diaria cada mes es entonces TX max 7 x4 7 TNx Sea Tn la temperatura m nima diaria en el mes k periodo j La m xima temperatura m nima diaria cada mes es entonces TN xs max T ny 8 TXn Sea Txwy la temperatura m xima diaria en el mes k periodo j La m nima temperatura m xima diaria cada mes es entonces TXny min Tx4 9 TNn Sea Tn la temperatura m nima diaria en el mes k periodo j La m nima temperatura m nima diaria en cada mes es entonces TNnyg min Tny 10 Tn10p Sea Tni la temperatura m nima diaria en el d a ien el periodo j y sea TninlO el d a calendario del percentil 10 centrado en una ventana de 5 d as calculado usando el m todo del Anexo D El porcentaje del tiempo es determinado donde Tni lt TninlO0 11 Tx10p Sea Tx la temperatura m xima diaria en el d a ien el periodo j y sea Txinl0 el d a calendario del percentil 10 centrado en una ventana de 5 d as calculado usando el m todo del Anexo D El porcentaje del tiempo es determinado donde Txi lt Txinl0 12 Tn90p Sea Tn la temperatura m nima diaria en el d a ien el periodo j y sea Tnin90 el d a calendario del percentil 90 centrado en una ventana de 5 d as calculado usando el m todo del Anexo D El porcentaje del tiempo es determinado donde Tni gt Tnin90 17 13 Tx90p Sea Tx la temperatura m xima diaria en el d a en el periodo j y sea T7xin90
4. el d a calendario del percentil 90 centrado en una ventana de 5 d as calculado usando el m todo del Anexo D El porcentaje del tiempo es determinado donde Txi gt Txin90 14 WSDI Sea Tx la temperatura m xima diaria en el d a en el periodo j y sea Txin90 el d a calendario del percentil 90 centrado en una ventana de 5 d as calculado usando el m todo del Anexo D Entonces el n mero de d as por periodo es sumado donde en intervalos de por lo menos 6 d as consecutivos Txi gt Txin90 15 CSDI Sea Tn la temperatura m nima diaria en el d a en el periodo j y sea TxinlO el d a calendario del percentil 10 centrado en una ventana de 5 d as calculado usando el m todo del Anexo D Entonces el n mero de d as por periodo es sumado donde en intervalos de por lo menos 6 d as consecutivos Tni lt Tninl0 16 DTR Sean Tx y Tny las temperaturas diarias m ximas y m nimas respectivamente en el d a en el periodo j Si I representa el n mero de d as en j entonces I AR Txi Tni DTR I 17 RX1day Sea RR el total diario de precipitaci n en el d a en el periodo j Entonces los valores m ximos de 1 d a para el periodo j son Rxlday max RR 18 Rx5day Sea RRy la cantidad de precipitaci n para el intervalo de cinco d as terminando en k periodo j Entonces los valores m ximos de 5 d as para el periodo j son Rx3day max RR4 18 19 SDII Sea RR la cantidad di
5. en un a o no hay m s de 15 d as faltantes No se calcular un valor anual si cualquier mes de datos esta faltando Para ndices con umbral el umbral se calcula si por lo menos el 70 de los datos est n presentes Para indicadores de duraci n de temporada marcados con un una temporada puede continuar en el a o siguiente y se la cuenta en el a o en que la temporada termina e g una temporada fr a CSDD en el Hemisferio Norte comenzando el 31 Diciembre 2000 y terminando el 6 Enero 2001 se cuenta dentro del n mero total de temporadas fr as en el 2001 1 FDO Sea Tn la temperatura m nima diaria en el d a en el periodo j Cuente el n mero de d as cuando Tn lt 0 C 2 SU25 Sea Tx la temperatura m xima diaria en el d a i periodo j Cuente el n mero de d as cuando Txi gt 25 C 3 IDO Sea Tx la temperatura m xima diaria en el d a en el periodo j Cuente el n mero de d as cuando Txi lt 0 C 4 TR20 Sea Tn la temperatura m nima diaria en el d a ien el periodo j Cuente el n mero de d as cuando Tni gt 20 C 5 GSL Sea Ty la temperatura media en el d a ien el periodo j Cuente el n mero de d as entre la primera ocurrencia de por lo menos 6 d as consecutivos con Ti gt 5 C Y la primera ocurrencia despu s de 1 Julio 1 Enero en HS de por lo menos 6 d as consecutivos con Ti lt 5 C 16 6 TXx Sea Txw la temperatura m xima diaria en el mes k periodo j
6. hay necesidad de llenar el par metro Station name or code pues este es usado luego Despu s de establecer el par metro n haga clic en OK para continuar 1 Data QC 21946 may take 5 minutes click OK to continue En PC s m s lentas este proceso puede tomar varios minutos Ventanas aparecer n si valores no razonables son encontrados Por ejemplo cuando la temperatura m nima diaria es mayor que la temperatura m xima diaria el siguiente mensaje aparece rms RClimDex NS a G Errors in temperature please check the log file Data QC Indices Calculation Exit Si hay valores negativos en los valores de precipitaci n diaria aparte de los valores faltantes que son codificados como 99 9 aparece el siguiente mensaje RClimDex ES a G Errors in prep please check the log file Data QC Indices Calculation Exit Si hay valores extremos aparece la siguiente ventana y Outliers found please check the log file X myXPfetccd RClimDexo0g 21946_tepstdQC csw Una ventana aparece una vez que el CC es completado Al mismo tiempo cuatro archivos Excel 21946tempQC csv 21946prcpQC csv 21946tepstdQC csv y 21946indcal csv son creados en un subdirectorio llamado log Los primeros dos archivos contienen informaci n sobre valores no razonables para temperatura y precipitaci n El tercer archivo marca todos los valores extremos posibles en temperatura d
7. menos 6 d as consecutivos en que TN lt 10th percentil Diferencia media mensual entre TX y TN M ximo mensual de precipitaci n en 1 d a M ximo mensual consecutivos de precipitaci n en 5 d as Precipitaci n anual total dividida para el n mero de d as h medos definidos por PRCP gt 1 0mm en un a o 13 UNIDAD D as D as D as D as E C C C D as D as D as D as C Mm d a R10 R20 Rnn CDD CWD R95p R99p PRCPTOT Number of heavy precipitation days N mero de d as con precipitaci n intensa Number of very heavy precipitation days N mero de d as con precipitaci n muy intensa Number of days above nnmm N mero de d as sobre nn mm Consecutive dry days D as secos consecutivos Consecutive wet days D as h medos consecutivos Very wet days D as muy h medos Extremely wet days D as extremadamente secos Annual total wet day precipitation Precipitaci n total anual en los d as h medos N mero de d as en un a o en que PRCP gt 10mm N mero de d as en un a o en que PRCP gt 20mm N mero de d as en un a o en que PRCP gt nn mm nn es un par metro definido por el usuario N mero m ximo de d as consecutivos con RR lt 1mm N mero m ximo de d as consecutivos con RR gt 1mm Precipitaci n anual total en que RR gt 95 percentil Precipitaci n anual total en que RR gt 99 percentil Precipitaci n anual total en los
8. que ser realizado de nuevo en el nuevo set de datos antes que estos cambios se reflejen en el c lculo de los ndices 3 3 C lculo de ndices RClimDex es capaz de calcular todos los 27 ndices b sicos enumerados en el Ap ndice A Sin embargo los usuarios pueden calcular solamente aquellos ndices que quieran Despu s de seleccionar ndices Calculation del men principal se pide al usuario que establezca algunos par metros para el c lculo de los ndices La ventana Set Parameter Values permite al usuario introducir el primer y ultimo a o del periodo base para el c lculo de los umbrales la latitud de la estaci n Hemisferio Sur es negativa para determinar en que hemisferio est localizada la estaci n un umbral de precipitaci n diaria definido por el usuario P en mm para calcular el n mero de d as en la precipitaci n diaria excede este umbral el indicador Rnn y 4 umbrales definidos por el usuario El User defined Upper Limit of Day High permite el c lculo del n mero de d as cuando la temperatura m xima diaria ha excedido este umbral El User defined Lower Limit of Day High permite el c lculo del n mero de d as cuando la temperatura m xima diaria est debajo de este valor El User defined Upper Limit of Day Low permite el c lculo del n mero de d as cuando la temperatura m nima diaria ha excedido este umbral El User defined Lower Limit of Day Low permite el c lculo del n
9. some demos help for on line help or help start for a HTML browser interface to help Type q to quit R gt 1 R 1 8 1 A Language and Environment Una vez que el c digo fuente es cargado exitosamente aparecer el men principal de RClimDex RClimDex 1 0 Load Data and Run QC Indices Calculation Exit 3 2 Cargar Datos y ejecutar Control de Calidad El control de calidad CC de los datos es un prerrequisito para el c lculo de los ndices El CC de RClimDex desarrolla el siguiente procedimiento 1 Reemplaza todos los datos faltantes actualmente codificados como 99 9 en un formato interno que reconoce R i e NA no disponible y 2 Remplaza todos los valores no razonables por NA Estos valores incluyen a cantidades de precipitaci n diarias menores que cero y b Temperatura m xima diaria menor que temperatura m nima diaria Adicionalmente CC tambi n identifica valores extremos outliers en temperaturas diarias m ximas y m nimas Los valores extremos son valores diarios que se encuentran fuera de una regi n definida por el usuario Actualmente esta regi n se define como n veces la desviaci n est ndar del valor del d a esto es media n std media n std Donde std representa la desviaci n est ndar para el d a y n es una entrada del usuario Seleccione Load Data and Run QC del Men RClimDex para abrir una ventana como se muestra abajo Esto permite al usuario selecciona
10. Calcula todos los 27 ndices b sicos recomendados por el Equipo de Expertos de CCI CLIVAR para Climate Change Detection Monitoring and ndices ETCCDMI as como tambi n otros ndices de temperatura y precipitaci n con l mites definidos por el usuario Los 27 ndices b sicos incluyen casi todos los ndices calculados por ClimDex Versi n 1 3 Esta versi n de RClimDex ha sido desarrollada bajo R 1 84 Debe poder ser ejecutada con R 1 84 o una versi n posterior Uno de los principales objetivos de construir ndices de extremos clim ticos es para usarlos para estudios de monitoreo y detecci n de cambios clim ticos Esto requiere que los ndices sean homogeneizados La homogenizaci n de datos ha sido planificada pero no implementada en esta versi n Actualmente RClimDex incluye solamente un procedimiento simple de control de calidad de datos que fue proporcionado por ClimDex Como en ClimDex se requiere que los datos tengan un control de calidad antes de que los ndices sean calculados Este manual de usuario proporciona instrucciones paso a paso sobre 1 La instalaci n de R y el establecimiento del ambiente de usuario 2 Control de Calidad de datos clim ticos diarios 3 C lculo de los 27 ndices b sicos 2 Instalaci n y ejecuci n de R R es un lenguaje y ambiente para c lculos y gr ficos estad sticos Es una implementaci n GNU del lenguaje S desarrollado por John Chambers y colegas de los Laboratorios Bell ant
11. RClimDex 1 0 Manual del Usuario Por Xuebin Zhang y Feng Yang Climate Research Branch Environment Canada Downsview Ontario Canada Versi n en Espa ol por Jos Luis Santos Centro Internacional Para la Investigaci n del Fen meno de El Ni o CUFEN Septiembre 20 2004 Reconocimientos RClimDex es desarrollado y mantenido por Xuebin Zhang y Feng Yang del Departamento de Investigaci n Clim tica del Servicio Meteorol gico de Canad Su desarrollo inicial fue financiado por la Agencia Canadiense de Desarrollo Internacional a trav s del Proyecto de Cambio Clim tico de la Cooperaci n Canad China C5 Lisa Alexander Francis Zwiers Byron Gleason David Stephenson Albert Klan Tank Mark New Lucie Vincent y Tom Peterson hicieron importantes contribuciones al desarrollo y prueba del paquete Versiones anteriores de RClimDex han sido usadas en Talleres de CCI CLIVAR ETCCDMI en Cape Town Africa del Sur Mayo 31 Junio 4 2004 y en Maceio Brasil Agosto 9 14 2004 Los Instructores y Participantes de estos talleres proporcionaron sugerencias muy valiosas para el mejoramiento de RClimDex TABLA DE CONTENIDO 1 Introducci n 2 Instalaci n y ejecuci n de R 2 1 C mo instalar R 2 2 C mo ejecutar R 3 C mo usar RClimDex 3 1 Cargando RClimDex 3 2 Control de Calidad de Datos 3 3 C lculo de Indices 4 Problemas conocidos 5 Reporte de problemas Anexo A Anexo B Anexo C Anexo D Anexo E
12. ake more than 5 minutes to compute all the indices Please be patient you will be informed once computations are done L ml Una ventana aparecer una vez que los ndices seleccionados han sido calculados Calculation Done Indices calculation completed Plots are in X myXP etced RClimDex plots La serie de ndices resultantes son almacenados en un subdirectorio titulado ndices en formato Excel Los archivos de ndices tienen nombres 21946_XXX cvs donde XXX representa el nombre del ndice Columnas de datos se separan por coma Para prop sitos de visualizaci n graficamos series anuales junto con tendencias calculadas por regresi n lineal de m nimos cuadrados l nea s lida y regresi n lineal con ponderamientos locales l nea punteada Estad sticas del ajuste de curvas lineal se muestran en los gr ficos Estos gr ficos se guardan en un subdirectorio llamado plots en formato JPEG Los nombres de archivos de los gr ficos siguen la misma regla excepto que cvs se cambia a jpg Seleccione Indices Calculation del men principal para calcular ndices adicionales seleccione Data QC y repita el proceso mencionado arriba Seleccione Exit si todos los c lculos requeridos han sido completados 11 4 Problemas conocidos Hay un problema conocido en esta versi n y versiones anteriores de RClimDex El programa se detendr si el primer a o de datos disponibles es el mismo que
13. aria de precipitaci n en d as h medos w RR gt Imm en el periodo j Si W representa el n mero de d as h medos en j entonces 20 R10 Sea RR la cantidad diaria de precipitaci n en el d a en el periodo j Cuente el n mero de d as donde RRi gt 10mm 21 R20 Sea RR la cantidad diaria de precipitaci n en el d a ien el periodo j Cuente el n mero de d as donde RRi gt 20mm 22 Rnn Sea RR la cantidad diaria de precipitaci n en el d a ien el periodo j Si nnrepresenta cualquier valor razonable de precipitaci n diaria entonces cuente el n mero de d as donde RRi gt nnmm 23 CDD Sea RR la cantidad diaria de precipitaci n en el d a en el periodo j Cuente el m s grande n mero de d as consecutivos donde RRi lt Imm 24 CWD Sea RR la cantidad diaria de precipitaci n en el d a en el periodo j Cuente el m s grande n mero de d as consecutivos donde RRi gt Imm 25 R95pTOT Sea RRw la cantidad diaria de precipitaci n en un d a h medo w RR 21 0mm en el periodo j y sea RRw95 el percentil 95 de precipitaci n en los d as h medos en el periodo 1961 1990 Si W representa el n mero de d as h medos en el periodo entonces 19 W ROS pj X RR where RRwj gt RRwn95 w 1 26 R99p Sea RR la cantidad diaria de precipitaci n en un d a h medo w RR 21 0mm en el periodo j y sea RRwn99 el percentil 99 de precipitaci n en los d as h medos en el per
14. d as h medos RR gt 1mm 14 D as D as D as D as mm ANEXO B Formato de entrada de datos Todos los archivos de datos que son le dos o escritos est n en formato de list formatted La nica excepci n es el primer archivo de datos que es procesado en el paso de Quality Control Este archivo de entrada de datos tiene varios requerimientos 1 Archivo de texto ASCII Columnas en las siguientes secuencias A o Mes D a PRCP TMAX TMIN NOTA unidades de PRCP mil metros y unidades de Temperatura grados Celsius 3 El formato descrito arriba debe ser delimitado por espacios e g cada elemento separado por uno o m s espacios 4 Para los registros de datos los datos faltantes deben se codificados como 99 9 los registros de datos deben estar en orden cronol gico Se permite fechas faltantes Ejemplo de formato de datos para el archivo de datos inicial e g usados en el paso de Quality Control 1901 1 1 99 9 30 1 6 8 1901 2 799 E E NS 1901 1 3 III AMO 1901 1 4 99 9 1 Rel 1901 1 7 9991 80 804 15 ANEXO C Definici n de los ndices A continuaci n se incluyen definiciones para los indicadores enumerados en el Ap ndice A Por razones pr cticas en esta versi n del software no todos los ndices son calculados en una base mensual ndices mensuales son calculados si en un mes hay no m s de 3 d as faltantes mientras que los valores anuales son calculados si
15. el primer a o del periodo base Esto es causado por algunos c lculos que requieren datos m s all del l mite del periodo base El c lculo de los ndices de temperatura basados en percentiles es un ejemplo Una manera de evitar este problema es a adir un registro extra para el d a con valores marcados como perdidos justo antes del inicio del periodo base Por ejemplo si el periodo base es 1961 1990 y los datos tambi n empiezan en 1961 uno puede a adir 1960 12 31 99 9 99 9 99 9 como la primera l nea del archivo de entrada 5 Reporte de Problemas Por favor reporte cualquier problema error a yang fengOec gc ca junto con los mensajes de error y los datos que se usaron para el c lculo de los ndices Esto ser de mucha ayuda para producir una mejor versi n en un futuro cercano Tambi n apreciar amos sus sugerencias para una mayor mejor a 12 FDO SU25 IDO TR20 GSL TXx TNx TXn TNn TN10p TX10p TN90p TX90p WSDI CSDI DTR RX1day Rx5day SDII Nombre del Indicador Frost days D as de heladas Summer days D as de verano Ice days D as de hielo Tropical nights Noches tropicales Growing season Length Duraci n de la estaci n de cultivo Max Tmax Max Tmin Min Tmax Min Tmin Cool nights Noches fr as Cool days D as fr os Warm nights Noches calientes Warm days D as calientes Warm spell duration indicador Indicador de la duraci n de per
16. es AT amp T ahora Tecnolog as Lucent S plus proporciona una implementaci n comercial del lenguaje S 2 1 Como instalar R RClimDex requiere el paquete base de R y la interfase gr fica de usuario TclITk La instalaci n de R requiere un procedimiento muy simple 1 Conectarse al sitio Web del proyecto R que es http www r project org y 2 Seguir los enlaces para descargar la versi n m s reciente de R para el sistema operativo de su computador desde cualquier sitio espejo de CRAN Para Microsoft Windows 95 98 2000 y XP descargar el programa de Windows setup Ejecute ese programa y R se instalar autom ticamente en su computador con un atajo a R en su desktop El TclTk se incluye en la instalaci n default de R 1 9 0 o versiones posteriores Podr a ser necesario que lo instale en forma separada si usted est ejecutando una versi n anterior de R Para Linux descargue los binarios precompilados apropiados y siga las instrucciones para instalar R Para otros sistemas Unix puede que tenga que descargar el c digo fuente y compilarlo usted mismo 2 2 Como ejecutar R Bajo el ambiente Windows haga un doble clic en el icono de R en su desktop o ejec telo a trav s del men Start de Windows Esto usualmente lo lleva a la interfase de usuario de R Para algunos computadores puede que usted necesite primero establecer una variable ambiente llamada HOME Para m s detalles vea las Preguntas frecuentes FAQ de W
17. es posible hacer una estimaci n exacta de los umbrales debido a incertidumbres en el muestreo Para proporcionar estimados temporalmente consistentes de la tasa de excedencia a lo largo del periodo base y fuera del periodo base adaptamos el siguiente procedimiento Zhang et al 2004 para estimar la tasa de excedencia para el periodo base a El periodo base de 30 a os se divide en un a o fuera de base el a o para el cual la excedencia se va a estimar y un periodo base que consiste en los restantes 29 a os a partir de los cuales se estimar n los umbrales b Se construye un bloque de datos de 30 a os usando el set de datos del periodo base de 29 a os y a adiendo un a o adicional de datos a partir del periodo base 1 e uno de los a os en el periodo base es repetido Este bloque de 30 a os construido es usado para estimar los umbrales c El a o fuera de base es entonces comparado con los umbrales y se obtiene la tasa de excedencia para el a o fuera de base d Pasos b y c se repiten 28 veces repitiendo cada uno de los 28 a os en base para as construir el bloque de 30 a os e El ndice final para el a o fuera de base es obtenido promediando los 29 estimados obtenidos en los pasos b c y d Referencia Hyndman R J y Y Fan 1996 Sample quantiles in statistical packages The American Statistician 50 361 367 Zhang X G Hegerl F W Zwiers y J Kenyon
18. iaria junto con las fechas en que aquellos valores ocurrieron El ltimo archivo contiene los datos una vez hecho el control de calidad y que ser n usados para el c lculo de los ndices Note que en este archivo solamente los valores faltantes y no razonables son reemplazados por NA posibles valores extremos que sean marcados NO son cambiados Para una f cil visualizaci n 4 archivos PDF que contienen gr ficos de las series de tiempo valores faltantes graficados como puntos rojos de precipitaci n diaria y temperaturas diarias m ximas y m nimas tambi n son guardados en log En este punto el usuario puede revisar los datos en el archivo 21946tepstdQC csv para determinar si cualquier valor marcado como un valor extremo es en realidad un valor extremo El archivo 21946indcal csv puede ser modificado bajo Windows usando Excel y cualquier editor bajo Unix si es necesario tomar cualquier acci n Despu s de completar este paso el usuario puede hacer clic en OK en la siguiente ventana para proceder con el c lculo de los ndices m i If you have checked data Xx myXP etccd RClimDex 21946indcal csw click OK to continue Note que los ndices son calculados a partir de los datos que se han hecho el CC El archivo de datos de entrada original no es alterado De manera que si un usuario decide modificar el archivo de datos original para corregir algunos de los datos con problemas el procedimiento Load Data y Run QC tiene
19. indows Ap ndice E si tiene usted problemas Bajo un ambiente unix simplemente ejecute R para que entre a la consola R Salga de R tipeando q en la consola R tanto bajo Windows como Unix Bajo Windows tambi n puede hacer clic en el men File y luego Exit 3 Como usar RClimDex 3 1 Cargando RClimDex En el prompt de la consola de R gt digite source rclimdex r Esto cargar RClimDex en el ambiente R Puede que tenga que incluir el path completo antes del nombre de archivo rclimdex r O en el caso que su computador este conectado a Internet puede descargar la versi n m s reciente del sitio Web ETCCDMI digitando source http cccma seos uvic ca ETCCDMV RClimDex rclimdex r Bajo Windows RClimDex puede tambi n ser cargado del men Escoja File del men RGui y luego seleccione Source R code Esto traer una nueva ventana dentro de la cual usted puede seleccionar el c digo fuente rclimdex r desde el directorio donde el programa fue guardado RGui File Edit Misc Packages Windows Help E de 03 The R Foundation for Statistical Computing 2003 11 21 ISBN 3 900051 00 3 re and comes with ABSOLUTELY NO WARRANTY o redistribute it under certain conditions or l cence for distribution details ive project with many contributors Sdrs for more information and on how to cite R in publications citation Type demo for
20. iodo 1961 1990 Si W representa el n mero de d as h medos en el periodo entonces W R99 p Y RR where RRw gt RRun99 w 1 27 PRCPTOT Sea RR la cantidad diaria de precipitaci n en el d a en el periodo j Si I representa el n mero de d as en j entonces I PRCPTOT 9_ RR i 1 20 ANEXO D Estimaci n de umbrales y c lculo del periodo base para los ndices de temperatura Estimaci n emp rica de cuantil El cuantil de una distribuci n se define como Q p F p inf x FG gt p 1 lt p lt l Donde F x es la funci n de distribuci n Sea X X m la estad stica de orden de X gt X i e valores ordenados de X y sea p la definici n del cuantil ith Los cuantiles muestrales pueden ser generalmente escritos como 0 p 1 7 X X jaa gt Hyndman y Fan 1996 sugieren una f rmula para obtener estimados medios un biased del cuantil haciendo int p n gt 1 p 3 y haciendo y p n 1 p 3 j donde int u es el entero m s grande que no sea mayor que u El cuantil emp rico es igual al valor en la muestra m s grande o al m s peque o cuando lt o j gt n respectivamente Esto es estimados de cuantiles correspondientes a p lt n 1 son fijados al valor m s peque o de la muestra y aquellos correspondientes a p gt n n 1 son igualados al valor m s grande de la muestra Procedimiento Bootstrap para la estimaci n de la tasa de excedencia para el periodo base No
21. iodos calientes Cold spell duration indicator indicador de la duraci n de periodos fr os Diurnal temperature range rango diurno de temperatura Max 1 day precipitation amount Cantidad M xima de precipitaci n en un d a Max 5 day precipitation amount Cantidad M xima de precipitaci n en 5 d as Simple daily intensity index ndice simple de intensidad diaria ANEXO A Lista de Los ndices Clim ticos b sicos de ETCCDMI Definici n N mero de d as en un a o cuando TN m nimo diario lt 0 C N mero de d as en un a o cuando TX m ximo diario gt 25 C N mero de d as en un a o cuando TX m ximo diario lt 0 C N mero de d as en un a o cuando TN m nimo diario gt 20 C Anual 1st Ene a 31 Dic en HN 1 Julio a 30 Junio en HS cuenta entre el primer periodo de por lo menos 6 d as con TG gt 5 C y primer periodo despu s de Julio 1 Enero 1 en HS de 6 d as con TG lt 5 C Valor mensual m ximo de temperatura m xima diaria Valor mensual m ximo de temperatura m nima diaria Valor mensual m nimo de temperatura m xima diaria Valor mensual m nimo de temperatura m nima diaria Porcentaje de d as cuando TN lt 10th percentil Porcentaje de d as cuando TX lt 10th percentil Porcentaje de d as cuando TN gt 90th percentil Porcentaje de d as cuando TX gt 90th percentil Contaje anual de d as con por lo menos 6 d as consecutivos en que TX gt 90th percentil Contaje anual de d as con por lo
22. mero de d as cuando la temperatura m nima diaria esta debajo de este l mite Estos ndices son llamados SUmm FDmm TRmm IDmm donde mm corresponde a valores definidos por el usuario Este paso incluye cierto procesamiento de datos por lo que tomar algunos segundos para terminar iD gt Ju Q i User defined parameters for Indices Calculation First year of base period ETS Last year of base period ECT Latitude of this station location A Longitude of this station location AS User defined upper threshold of daily maximum temperature 2 User defined lower threshold of daily maximum temperature A User defined upper threshold of daily minimum temperature EN User defined lower threshold of daily minimum temperature A User defined daily precipitation threshold A oK CANCEL y Una vez que este paso ha sido completado una ventana aparecer para permitir al usuario seleccionar los ndices que desea calcular Todos los ndices son seleccionados por default Quite la marca de los ndices que no se necesiten luego haga clic en OK pare realizar el c lculo Dependiendo de los ndices seleccionados este procedimiento puede tomar un momento 10 Check desired indices SU25 FDO TR20 IDO User Defined SU FD TR ID GSL growing season length TXx TXn TNx TNn TX10p TX90p TN10p TN90p WSDI CSDI DTR Rx1day Rx5day SDII R10mm R20mm Ennmm CDD CWD R95p R99p PRCPTOT OK CANCEL It may t
23. r cargar el archivo de datos a partir del cual los ndices ser n calculados History G te ESA E My Documents My Computer File name Files of type TEXT Files txt y My Network R El nombre de archivo debe se de la forma stationname txt Los valores en el archivo deber an ser en el formato descrito en el Ap ndice B En este men usamos datos de una estaci n cuyos datos fueron almacenados en el archivo ASCH 21946 txt para prop sitos de demostraci n Una ventana como se muestra abajo aparecer una vez que los datos para la estaci n 21946 sean cargados exitosamente i Data 21946 loaded click OK to continue Mensajes de error aparecer n en la consola de R si este paso no ha sido completado exitosamente Esto usualmente es causado por datos de entrada en formato equivocado Por favor compare su formato con nuestros datos de ejemplo si usted ve tales mensajes Valores no razonables son identificados autom ticamente pero la identificaci n de valores extremos en datos de temperatura necesita input del usuario Set Parameters for Data OC Station name or code Criteria number of Standard Deviation la El valor default para n es 3 Criterio en la ventana Set Parameters for Data QC pero este n mero puede ser cambiado por el usuario Como un valor de 3 puede marcar un n mero de valores muy grande los usuarios pueden desear comenzar estableciendo este valor a 4 No
Download Pdf Manuals
Related Search