Los datos de reanálisis son datos meteorológicos que describen el estado de la atmósfera en una malla regular de puntos a diferentes alturas que cubren todo el planeta. Estos datos se obtienen a partir del análisis de los datos meteorológicos históricos disponibles y la simulación con modelos meteorológicos de manera que se genera una malla regular de datos. Se puede obtener más información de este tipo de datos desde los enlaces disponibles en su entrada en Wikipedia.
Para representar los mapas utilizo el programa Grid Analysis and Display System, GRADS.
El primer paso es descargar los datos de reanálisis NCEP/NCAR. Para ello debemos registrarnos en el servidor de datos de UCAR y navegar hasta la página de los datos de reanálisis. En este ejemplo he descargado los datos de septiembre de 1996; de entre las posibilidades que ofrecen estos reanálisis he utilizado los datos “pgb.f00 – 2.5 degree analyses on 17 pressure levels, plus other surfaces and layers”. Estos datos tienen cobertura global y una resolución espacial de 2.5 grados (demasiado gruesa para estudiar fenómenos locales) en 17 niveles de presión en la vertical.
El procedimiento que viene a continuación se ha realizado sobre Linux (Ubuntu 10.04) pero debería poder hacerse también en Windows aunque nunca he probado GRADS sobre Windows.
El primer paso es descomprimir el archivo ya que el fichero de reanálisis mensual es un agrupamiento de ficheros diarios.
tar -xvf A09542-199609.pgb.f00
de manera que se producen una serie de ficheros con los datos cada 6 horas
pgb.f0096090100 pgb.f0096090106 pgb.f0096090112 pgb.f0096090118 ...............
Una vez descomprimidos los datos deben crearse los ficheros .ctl e .idx necesarios para que GRADS pueda representar las distintas variables contenidas en los archivos. Para ello es necesario el script grib2ctl.pl de Wesley Ebisuzaki.
grib2ctl.pl pgb.f0096090100 > pgb.f0096090100.ctl gribmap -i pgb.f0096090100.ctl -0 <-- Esto es un cero
Ahora ya se puede abrir GRADS y comenzar con los mapas:
grads (para arrancar el programa)
Dentro de GRADS veamos un ejemplo sencillo)
open pgb.f0096090100.ctl (abre el fichero ctl) q file (muestra información acerca del fichero y sus variables) set z 1 (selecciona el primer nivel de presión, 1000 hPa) set gxout shaded (para pintar zonas de color en lugar de isolíneas) d tmpprs (temperatura al nivel de presión seleccionado) set gxout contour (ponemos los contornos de los niveles) d tmpprs
Mapa de temperatura al nivel de 1000 hPa, datos de reanálisis del 31 de septiembre de 1996
A partir de aquí se puede ir experimentando con las distintas opciones de GRADS.
