Introducción al GrADS
Contents
1. correcto Sin embargo para otros fines calculos podemos dejar variando tantas dimensiones como queramos Ejemplos Un corte vertical dejando la latitud constante gt set lat 40 gt set lev 1000 100 gt set lon 90 50 gt sett1 Obtener una marcha temporal en un punto y nivel fijos gt set lat 34 gt set lon 58 gt set lev 1000 gt sett1 10 Obtener el valor de una variable en un punto y tiempo determinados Agregar al caso anterior gt sett2 En este caso al usar display no generamos un gr fico sino que el valor aparece en la pantalla donde ingresamos los comandos La excepci n a la regla de las dimensiones es la posibilidad de mostrar animaciones Si dejamos 3 dimensiones variando y una de ellas es el tiempo obtendremos una animaci n de gr ficos Si queremos volver a como estaban seteadas las dimensiones en el momento que el grads abre el ctl utilizamos el comando reset Otro comando que realiza esta operaci n es el comando reinit pero a diferencia de reset este comando cierra todos los archivos de datos y tenemos que abrirlos nuevamente Superposici n de gr ficos En grads hasta que no pongamos el comando c clear la pantalla no se borra y todos los gr ficos que realicemos se ir n superponiendo en el orden mismo orden que los realicemos Superponer en 500 hPa el gr fico de altura geopotencial y el de temperatura Expresiones El GrADS permite realizar operaciones matem ticas con las2. Cuentas y m s cuentas Algunos ejemplos de cuentas que se pueden hacer en GrADS no son todos gt d mag u v Calcula la magnitud del vector u v gt d hcurl u v Calcula la vorticidad relativa gt d hdivg u v Calcula la divergencia de u v gt d pow var exp Calcula varexP gt d var lev 500 lev var 1000 Hace la resta del valor de var en 500 el valor de var en 1000 gt d var t 2 var t 1 Idem anterior pero en el tiempo gt d ave var t 1 t 10 Grafica el promedio de var entre los tiempos 1 y 10 gt d sqrt u u v v Magnitud del viento usando la raiz cuadrada Funciones de GrADS sin cos tan y sus inversas min max minloc maxloc aave abs asum cdiff const exp fndlvl hcurl hdivg log log10 mag maskout scorr skip smth9 sart vint Hay m s Maskout La funci n maskout reemplaza los valores de una variable por el valor de dato faltante cuando se cumple una condici n que nosotros impongamos gt a maskout variable condicion Condicion es una variable o una expresion cuando es lt 0 los valores de variable se reemplazan por el valor de dato faltante y por ende no son graficados si hacemos un display En este caso el resultado del maskout lo guardamos en una variable que se llama a a maskout t t 273 a tiene los valores de t cuando t es mayor que 273 para el resto de los puntos tiene el valor de dato faltante Podemos usar esto para no mostrar los valores de las
3. DS que el GrADS conozca ese valor para que no lo interprete dset gdas pibal pgban y4 m2 d2 h2 ae al de ja senal a Nana aX e que queremos graficar pan 2 999E 20 produced by grib2ctl v0 9 12 5p32I dtype grib 3 options yrev ydef 181 linear 90 000000 1 xdef 360 linear 0 000000 1 00000 tdef 10 linear 00Z15Dec02 O6hr zdef 26 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10 vars 84 no4LFTXsfc 0 132 1 0 surface Best 4 layer lifted index K Qu informaci n necesita el GrADS T tulo Simplemente por prolijidad dset gdas pibal pgban y4 m2 d2 h2 index ana idx Do a 0 9 12 5p32l dde pen 3 options yrev ydef 181 linear 90 000000 1 xdef 360 linear 0 000000 1 00000 tdef 10 linear 00Z15Dec02 O6hr zdef 26 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 7050 30 20 10 vars 84 no4LFTXsftc 0 132 1 0 surface Best 4 layer lifted index K Qu informaci n necesita el GrADS dla del tipo de atos dset gdas pibal pgban y4 m2 d2 h2 index ana idx options template undef 9 999E 20 na ra 100z pgrbof00 y grib2ctl vO 9 12 5p321 ydef 181 lea 90 000000 1 xdef 360 linear 0 000000 1 00000 tdef 10 linear 00Z15Dec02 O6hr zdef 26 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 7050 30 20 10 vars 84 no4LFTXsftc 0 132 1 0 surface B
4. Introducci n al GrADS Curso Intensivo en Asimilaci n de datos Buenos Aires 2008 GrADS Grid Analysis and Display System es un software gratuito que sirve para facilitar el acceso manipulaci n y visualizaci n de datos vinculados con las ciencias de la tierra Esta presentaci n constituye una introducci n b sica al manejo del programa as como una discusi n de sus principales caracter sticas Se discutir n los siguientes aspectos Inicio del programa Tipos de datos que podemos visualizar en GrADS El archivo de descripci n de datos CTL Comandos b sicos Introducci n a los scripts Comentarios finales Inicio del programa En la versi n Linux iniciamos el programa desde una terminal tipeando el comando grads en algunas m quinas gradsnc gt grads Esto despliega una ventana adicional donde se muestran los gr ficos Display los comandos ser n ingresados por la terminal e En la ventana de comandos terminal GrADS pregunta si la ventana gr fica abrir en modo apaisado landscape o vertical portrait El modo apaisado tiene una dimensi n de 11x8 5 pulgadas y el modo vertical una de 8 5x11 El tama o de la ventana gr fica puede ser modificado a posteriori solo en Linux Modo landscape Modo portrait El modo a seleccionar depende de lo que se quiera graficar Tipos de archivos que se pueden graficar en GrADS El GrADS puede leer los siguientes formatos Binario gen
5. ariable que antes pero en una regi n cercana a Buenos Aires gt set lat 37 32 la menor siempre primero gt Set lon 60 55 En este caso le decimos al grads considere solo las latitudes entre 37 y 32 y las longitudes entre 60 y 55 Si graficamos nuevamente la variable el gr fico solo mostrar la regi n seleccionada Hasta ahora solo graficamos el primer tiempo del archivo porque el grads cuando abre un archivo autom ticamente se situa en el primer tiempo y hasta que no le indiquemos otra cosa grafica las variables en el primer tiempo Para modificar el tempo podemos hacer lo siguiente gt sett2 Le indica al grads que se pare en el segundo tiempo del archivo de datos Si queremos ir a una fecha en particular podemos utilizar la opci n time gt set time 12719dec2002 la fecha va en el formato HHzDDMMMYYYY MMM son las 3 primeras letras del mes en ingl s Para obtener informaci n sobre en que tiempo estamos actualmente podemos usar el comando gt q time Este comando nos devolver la fecha o rango de fechas como veremos m s adelante que corresponden a las dimensiones actuales Graficar la altura geopotencial para el tiempo 3 del archivo De manera an loga podemos cambiar el nivel en la vertical gt set lev 500 Con este comando le decimos al GrADS que lo que queremos graficar es el valor de las variables en el nivel de 500 hPa El GrADS por defecto se posiciona en primer nivel del archivo de datos en nu
6. ario que se genere un ctl Comandos b sicos e C mo abrir un archivo Para abrir un archivo utilizando un ctl usamos el comando open gt open ruta_a_mi_archivo archivo ctl Si existe alg n problema el GrADS mostrar un mensaje de error Para obtener informaci n de lo que contiene el archivo podemos utilizar el comando q file o query file Este comando nos devolver informaci n sobre el tama o de la ret cula en la que est n disponibles los datos y nos dar una lista de las variables que se encuentra en el archivo Esta descripci n es generada a partir del ctl es decir que toda la informaci n all contenida la podemos obtener directamente mirando el ctl A continuaci n vamos a graficar una variable Para graficar utilizamos el comando display d Por ejemplo du Va a graficar la componente zonal del viento en contornos Sobre todo el dominio disponible En nuestro caso el dominio disponible es la regi n adyacente a Sudam rica El comando c clear permite borrar la pantalla totalmente Una de las principales ventajas del GrADS es que trata los datos como variables de 4 dimensiones x y z t Eso nos permite graficar series temporales cortes verticales cartas en planta etc Para poder hacer todas estas cosas debemos aprender a controlar las dimensiones dimension enviroment El comando set es el que nos permite hacer esto Por ejemplo supongamos que queremos dibujar la misma v
7. delante Por qu utilizar scripts Para ahorrar tiempo GrADS no tiene una interface gr fica y no permite deshacer los cambios que realizamos sobre un gr fico Un gr fico particular puede llevar un n mero grande de comandos Si los ingresamos a mano podemos equivocarnos o podemos querer m s tarde repetir el gr fico con un ligero cambio con lo cual deber amos ingresar todos los comandos de nuevo Haciendo un script reducimos la posibilidad de equivocarnos y podemos repetir un mismo gr fico realizando peque as modificaciones a nuestro script original Por eso es mejor desde el momento que uno decide que gr fico hacer comenzar a armarlo desde un script y no ingresando comandos por la pantalla Un script es un archivo ascii se puede editar con el vi gedit nedit o similares En cada rengl n se coloca un comando de GrADS en el mismo orden que lo har amos si estuvi ramos ingresando los comandos por pantalla En general a los scripts de GrADS por convenci n se les asigna la extensi n gs pero cualquier extensi n puede ser utilizada Para ejecutar un script de este tipo debemos usar el comando exec gt exec mi_script gs Ejecuta mi_script gs es importante que le aclaremos al GrADS la ruta completa al lugar donde est el script Ejemplo de un script reinit open pepe ctl sett10 set lev 500 d skip u 3 v printim viento500hPa png white Existe otro tipo de scripts que son los que permiten realizar calculos m
8. direcci n x o y en algunas versiones se puede tambi n calcular las diferencias centradas en z y t Ejemplo Como calcular la derivada del geopotencial respecto de x dx cdiff lon x 111000 cos lat 3 14159 180 Esto calcula el dx en metros dz cdiff Z x Esto calcula el dz dzdx dz dx Finalmente este es el valor de la derivada de z con respecto a x Notar que lat y lon son variables que pueden o no estar en el ctl pero que independientemente de que est n o no el GrADS las tiene definidas las calcula a partir de lo que dice el ctl y podemos utilizarlas en los c lculos Ejemplo C mo calcular amos el par metro de Coriolis Tipos de gr ficos en GrADS Adem s de poder operar con las variables y cambiar las dimensiones que queremos graficar GrADS nos permite realizar diferentes tipos de gr ficos El tipo de gr fico se selecciona con el comando set gxout tipo_de grafico donde tipo_de grafico puede ser alguno de los siguientes Para campos 2D Shaded sombreado diferentes colores para los diferentes rangos de la variable Contour contornos grafico de contornos isolineas de la variable seleccionada Grfill shaded grid boxes Fgrid shaded grid boxes y valores Stream lineas de corriente Vector vectores Barb barbas Para gr ficos 1D bar grafico de barras errbar grafico con barras de error line grafico de lineas linefill grafico de lineas pero rellenando con colores en
9. emplo el primer Qu informaci n necesita el GrADS dset gdas pibal pgban y4 m2 d2 h2 index ana idx options template undef 9 999E 20 title gdas1 t00z pgrbof00 produced by grib2ctl v0 9 12 5p321 dtype grib 3 options yrev ydef 181 linear 90 000000 1 xdef 360 linear 0 000000 1 tdef 10 linear 0OZ15Dec n mero es el n mero de niveles verticales disponibles para esa variable en este caso O indica un solo nivel Al final de cada variable se encuentra una descripci n de la misma La lista 000 culmina con la sentencia 06hr ENDVARS zdef 26 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 noz Xstc 0 132 1 O surface Best 4 layer lifted index K Los archivos binarios generados con fortran requieren de un CTL para saber cual es el nombre de las variables y la ubicaci n espacial a la cual corresponden los datos almacenados en el archivo Los archivos GRIB poseen dentro del archivo toda esta informaci n pero para poder visualizarlos correctamente en el GrADS es necesario construir un archivo CTL y generar un archivo indice que permitan que GrADS interprete este formato Algunos archivos NetCDF pueden ser abiertos directamente con el GrADS sin necesidad de generar un archivo ctl ya que tambi n contienen toda la informaci n necesaria sobre los datos guardados en el archivo Sin embargo en algunos casos particulares es igualmente neces
10. erado con fortran por ejemplo GRIB GRldded Binary es un estandar de la OMM para la transmici n de datos en puntos de ret cula NetCDF SDS Otros formato muy utilizado para guardar las salidas de los modelos num ricos y los an lisis NetCDF y GRIB son archivos que contienen Metadata o informaci n sobre lo que contienen Los archivos binarios por lo general no y tienen que venir acompa ados de un archivo que contenga informaci n sobre lo que contienen Para saber que informaci n est contenida en un archivo el GrADS utiliza un archivo de descripci n de datos Data descriptor file que es un archivo ASCII que se puede editar f cilmente y que utiliza por convenci n la extensi n ctl El archivo de descripci n de datos El nombre del archivo de Qu informaci n necesita el GrADS datos con la ruta incluida o dset gdas pibal pgban y4 m2 d2 h2 el simbolo para indicar que index anatdz se encuentra en la misma options template carpeta que el ctl undef 9 999E 20 title gdas1 t00z pgrbof00 produced by grib2ctl vO0 9 12 5p321 dtype grib 3 options yrev ydef 181 linear 90 000000 1 xdef 360 linear 0 000000 1 00000 tdef 10 linear 00Z15Dec02 O6hr zdef 26 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 7050 30 20 10 vars 84 no4LFTXsftc 0 132 1 0 surface Best 4 layer lifted index K El archivo indice necesario unicamente cuando
11. est 4 layer lifted index K Esta opci n se puede poner junto con la primera e indica que el orden lo los puntos Qu informaci n necesita el GrADS dset gdas pibal pgban y4 m2 d2 h2 index ana idx options template undef 9 999E 20 title gdas1 t00z pgrot00 Maca by grib2ctl vO en latitud es al rev s que como el GrADS los lee habitualmente del 360 linear o 000000 1 00000 tdef 10 linear 00Z15Dec02 06hr zdef 26 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 7050 30 20 10 vars 84 no4LFTXsftc 0 132 1 0 surface Best 4 layer lifted index K Dimensi n espacial de las matrices Cuantos puntos en x y en y Qu informaci n necesita el GrADS tiene Como es el dset gdas pibal pgbanl y4 m2 d2 h2 espaciamiento entre puntos index ana idx y en el caso de ser lineal options template como en el ejemplo cual es undef 9 999E 20 el valor inicial y el title gdas1 t00z pgrbof00 incremento o produced by grib2ctl vO 9 12 5p32 espaciamiento dtype sl 3 tdef 10 linear 00Z15Dec02 O6hr zdef 26 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 705030 20 10 vars 84 no4LFTXsftc 0 132 1 0 surface Best 4 layer lifted index K Dimensi n temporal Cuantos tiempos contiene el archivo cual es la fecha Qu informaci n necesita el GrADS HHzDDMMMYYYY del primer tiempo y con que frecuenc
12. estro ejemplo 1000 hPa Si el nivel elegido no est disponible GrADS buscar el nivel m s cercano Otra forma de cambiar el nivel es usar el comando gt setz 2 Este comando cambia el valor del nivel al segundo nivel para saber a que nivel corresponde debemos ver el ctl y ver el segundo valor de la lista de niveles Esto pude ser til para hacer c lculos Graficar una variable en el nivel de 500 hPa El comando gt q dims Nos da informaci n sobre como est n configuradas las dimensiones en un momento dado es decir el nivel niveles tiempo tiempos latitudes longitudes etc que estamos considerando en un momento dado En grads todas las dimensiones pueden dejarse fijas o variables Una dimensi n es fija cuando le asignamos un nico valor gt set lev 1000 fija el valor del nivel en 1000 Y una dimensi n es variable cuando le asignamos un rango de valores gt set lev 1000 500 le dice al grads que hagamos lo que hagamos considere todos los niveles entre 1000 y 500 hPa El comando q dims tambi n nos indica que dimensiones est n fijas y cuales est n como variables Importante GrADS solo hace gr ficos en 2D y 1D es decir que al momento de hacer un gr fico solo pueden estar variando 2 dimensiones puede ser cualquier combinaci n pero no pueden ser m s de 2 Si tenemos m s de 2 dimensiones seteadas como variables y pedimos que el GrADS grafique obtendremos una animaci n siempre y cuando set loopdim est
13. ia temporal est n disponibles los datos dset gdas pibal pgban y4 m2 d2 h2 index ana idx options template undef 9 999E 20 title gdas1 t00z pgrbof00 produced by grib2ctl v0 9 12 5p321 dtype grib 3 options yrev ydef La linear 90 000000 1 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 7050 30 20 10 vars 84 no4LFTXsftc 0 132 1 0 surface Best 4 layer lifted index K Qu informaci n necesita el GrADS dset gdas pibal pgban y4 m2 d2 h2 index ana idx options template undef 9 999E 20 title gdas1 t00z pgrbof00 produced by grib2ctl v0 9 12 5p32 dtype grib 3 options yrev ydef 181 linear 90 000000 1 xdef 360 linear 0 000000 1 00000 tdef 10 linear 0OZ15Dec02 O6hr Niveles verticales Cuantos niveles verticales son la opcion levels en este caso se utiliza cuando los niveles verticales no est n equiespaciados a continuaci n debemos especificar los niveles de presi n o altura disponibles en el archivo de datos La opci n levels tambi n se puede usar para las dimensiones x e y 1000 975 950 o 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 vars 84 no4LFTXsftc 0 132 1 0 surface Best 4 layer lifted index K Vars indica el numero de variables presentes en el archivo de datos A continuaci n viene una lista con el nombre de cada variable y algunos n meros que indican propiedades de las mismas Por ej
14. n rengl n abajo Agregar barra de colores Si hicimos un gr fico en shaded podemos agregar una barra de colores en donde se indiquen los intervalos de la variable y que color corresponde a cada intervalo Para hacer esto tipeamos gt run colorbar una vez que tengamos realizado el gr fico Cambiar el tipo de mapa GrADS tiene 3 mapas que vienen con la distribuci n los mapas se llaman hires medres y lowres gt set mpdset mapa Cambiar el tipo de proyecci n GrADS trabaja con m ltiples proyecciones por defecto utiliza la mercator Para cambiar por ejemplo a proyecci n polar se puede usar el comando mproj gt set mprojsps Estereografica polar para el HS Guardar el grafico en un archivo Hay al menos 2 maneras Utilizando el comando enable print gt print gt disable print gt enable print mi_figura gmf gt grafico lo que quiero gt print gt disable print El formato gmf me permite guardar muchos gr ficos en un mismo archivo esto quiere decir que puedo hacer varios prints de varias figuras distintas antes de hacer un disable print y cerrar definitivamente el archivo El formato gmf puede ser visualizado con un programa gratuito que viene para Windows o bien convertido a otros formatos con las herramientas gxgif gxeps etc que vienen con la distribuci n del GrADS El comando printim es la segunda opci n y la m s sencilla para generar figuras que no requieran mucha calidad gt Genero una grafic
15. o gt printim mi_grafico gif gif x800 y600 white Este comando genera un archivo llamado mi_grafico gif en formato gif con 600 puntos en y 800 en x y en fondo blanco El formato de salida por defecto del printim es png gt printim mi_grafico png Guarda el gr fico que estamos viendo en un archivo png Como generar varios gr ficos en la misma Pagina Supongamos que queremos hacer una figura con 2 paneles en la misma p gina entonces podemos usar el comando vpage para graficar en 2 sectores de la p gina separados Vpage crea una p gina virtual MODO LANDSCAPE dentro de la p gina total gt set vpage 0 5 5 0 8 5 Le dice que la p gina virtual se extiende entre O y 5 5 en X y entre O y 8 5 en Y 8 5 pulgadas ai Y Cuando realicemos un gr fico Y este aparecer en la regi n roja del display el resto quedar libre para que realicemos otro gr fico 11 pulgadas X Cuando lo guardemos a una Origen A imagen ambos ir n al mismo archivo 5 5 pulgadas Antes de definir una nueva p gina virtual debemos decirle que ignore la anterior utilizando el comando set vpage off Introducci n a los scripts As como ingresamos las instrucciones a mano por la pantalla GrADS tiene la posibilidad de leer una lista de instrucciones de un archivo de texto y ejecutarlas en el mismo orden Al mismo tiempo GrADS posee su propio lenguaje que permite realizar c lculos m s complejos como veremos m s a
16. s complejos La descripci n de estos se puede encontrar en la secci n GrADS scripting language de la p gina de GrADS Las herramientas que provee este lenguaje se asemejan a las de un lenguaje de programaci n tipo fortran Estos scripts se ejecutan con el comando run Para escribir estos scripts se debe seguir una serie de normas diferentes a las necesarias para generar los scripts que se ejecutan con exec Comentarios finales Recursos de ayuda en el uso de GrADS Indice de funciones http grads iges org grads gadoc gadocindex html Guia del usuario hitp grads iges org grads gadoc users html Manual en PDF ftp grads iges org grads beta doc gadoc151 paf Guia r pida imprescindible http grads iges org grads gadoc reference_card paf Guia para scripts http grads iges org grads gadoc reference card scl paf Donde bajar el software http grads iges org grads downloads html
17. se ngbanl Yy4Y m2 d29 trabaja con los datos GRIB Qu informaci n necesita el GrADS contiene informaci n de pta donde est cada variable en undef 9 999E 20 el archivo GRIB title gdas1 t00z pgrof00 produced by grib2ctl v0 9 12 5p32I dtype grib 3 options yrev ydef 181 linear 90 000000 1 xdef 360 linear 0 000000 1 00000 tdef 10 linear 0O0Z15Dec02 06hr zdef 26 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10 vars 84 no4LFTXsfc 0 132 1 0 surface Best 4 layer lifted index K Opciones adicionales sobre el formato de los datos o A Template indica que los datos Qu informaci n necesita el GrADS est n distribuidos en varios archivos generalmente dset A y4 m2 d2 h2 s ci identificados con el nombre de la fecha Sequential big endian little title dEl 100Z ao produced by grib2ctl v0 9 12 5p321 endian son otras opciones dtype grib 3 posibles vinculadas con el options yrev formato de los archivos ydef 181 linear 90 000000 1 binarios xdef 360 linear 0 000000 1 00000 tdef 10 linear 00Z15Dec02 06hr zdef 26 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 7050 30 20 10 vars 84 no4LFTXsftc 0 132 1 0 surface Best 4 layer lifted index K En los archivos de datos por lo general se asigna un valor a los datos faltantes Es necesario Qu informaci n necesita el GrA
18. tornos se dibujan cada 2 Por ejemplo cada 2 K en el caso de la temperatura o cada 2 hPa en el caso de la presi n gt set ccolor 2 Fijamos el color de las isolineas en el valor 2 gt set ccolor rainbow En este caso cada contorno se dibuja con un color distinto y el orden es el que se muestra en la siguiente diapositiva Los colores que vienen por defecto en el GrADS y el n mero asociado a cada uno The GrADS Default colors Col Description Sample R G B O background 0 00 black by default 1 foreground 255 255 255 white by default 2 red 250 60 60 3 green 02200 4 dark blue 30 60 255 5 light blue 0 200 200 6 magenta 240 0 130 7 yellow 230 220 50 8 orange 240 130 40 9 purple 160 0 200 10 yellow green 160 230 50 11 medium blue O 160 255 12 dark yellow 230 175 45 13 aqua 0210140 14 dark purple 130 0 220 15 gray 170170170 The GrADS Default Rainbow Sequence Col Description Sample R G B 9 purple 160 0 200 14 dark purple 1100220 4 dark blue 30 60 255 11 medium blue O 160 255 5 light blue 0 200 200 13 aqua 0210140 3 green 0220 O 10 yellow green 160 230 50 7 yellow 230 220 50 12 dark yellow 230 175 45 8 orange 240 130 40 2 red 250 60 60 6 magenta 240 O 130 Agregar un t tulo Una vez que hicimos un gr fico podemos agregar un t tulo utilizando el comando gt draw title Grafico de geopotencial en 500 hPa Si el t tulo es muy largo podemos usar 1 para indicarle al GrADS que tiene que seguir en u
19. tre 2 l neas scatter scatter plot una variable en funci n de la otra stat muestra estad sticas de la variable en cuesti n Ejemplo Grafico 2D de viento gt resel gt set lev 500 gt set gxout vector gt d uv Grafica los vientos en 500 hPa Un vector por cada punto de ret cula Cuando los puntos de ret cula son muchos este gr fico se puede ver muy mal Para eso tenemos la funci n skip La funci n skip se queda con un dato cada cierto n mero de puntos de ret cula y pone dato faltante en los dem s u2 skip u 3 u2 est tiene un dato en uno de cada 3 puntos en x e y de los que tenia u gt set gxout vector gt d skip u 3 v Nos muestra un campo de vectores con menos flechas No es necesario aplicar el skip a u y a v ya que si u es dato faltante el grads ya no puede reconstruir el vector en ese punto y por ende no mostrar la flecha aunque el valor de v est disponible Skip permite saltear diferente cantidad de puntos en x y en y u2 skip u 3 2 Se queda con 1 de cada 3 puntos en x y 1 de cada 2 en y Esto es til cuando la resoluci n en x y en y son muy diferentes o bien cuando el eje y en realidad representa la vertical M s comandos tiles para graficar Elegir los contornos sombreados en un gr fico gt set clevs valor1 valor2 valor3 Nos permite elegir que contornos queremos dibujar antes de hacer un grafico de contour o shaded gt set cint 2 Fijamos un intervalo los con
20. variables Por ejemplo si queremos dibujar la temperatura en C pero la tenemos en K entonces debemos restarle 273 16 gt di 273 16 realiza la operaci n deseada y la grafica Tambi n podemos definir nuevas variables que sean el resultado de alg n c lculo Para eso contamos con la funci n define gt define tc t 273 16 Esto crea una nueva variable tc que no estaba en el archivo original y le asigna el valor de t 273 16 tc tiene las dimensiones que ten a el grads en el momento de hacer la cuenta gt set lat 90 0 gt set lon 100 50 gt setit1 gt set lev 1000 gt define tc t 273 16 Define una variable tc que solo esta definida en el nivel de 1000 en el tiempo 1 y en la regi n seleccionada Si en lugar de eso hacemos gt set lat 90 O gt set lon 100 50 gt sett1 10 gt set lev 1000 100 gt define tc t 273 16 tc estar definida para los primeros 10 tiempos del archivo y para los niveles entre 1000 y 100 En las versiones m s nuevas no es necesario poner define y la variable tc se puede crear simplemente usando el siguiente comando gt tc t 273 16 Para plotear la temperatura en F gt d t 273 16 9 5 32 Tambi n podemos realizar operaciones que involucren 2 o m s variables gt densidad lev 287 T Calcula la densidad del aire seco Lev es una variable que tiene la informaci n del nivel que en el caso de los datos nuestros es el valor de la presi n a la cual nos encontramos
21. variables que estan por debajo de la topografia a maskout variable z topo Z es el valor de la altura geopotencial y topo es el valor de la topografia esto puede variar de un set de datos a otro Cuando la altura geopotencial es menor que la altura de la topografia entonces reemplazamos por el valor de dato faltante Para hacer un ejemplo de esto en nuestro caso necesitamos ver los datos de topografia que est n en un archivo de datos diferente Como combinar la informaci n de dos archivos distintos Primero abrimos los 2 archivos ejecutamos el comando open 2 veces 1 para cada archivo de datos Cada archivo se identifica con un n mero de acuerdo con el orden en que los fuimos abriendo Por ejemplo para graficar la variable u del archivo 1 el primero que abrimos gt d u Para graficar la variable topo en el segundo archivo que abrimos gt d topo 2 Importante GrADS solo nos permite hacer cuentas que involucren variables del archivo 1 y 2 juntas si las ret culas de datos son iguales podr a ser que una fuera m s chica que la otra pero los puntos de la m s chica deber an coincidir con los puntos de la mayor gt set lev 850 gt d maskout t z topo 2 Combino variables del archivo 1 y 2 en una misma cuenta usando el n para identificar de que archivo proviene la informaci n La funci n cdiff y el c lculo de derivadas La funci n cdiff calcula las diferencias centradas de una variable con respecto a la
Download Pdf Manuals
Related Search