NCEP Climate Forecast System Reanalysis (CFSR)
ds093.0: 6-hourly Products
 Reanalysis  1980-2010(*)

http://rda.ucar.edu/datasets/ds093.0/index.html
Copia de los archivos para uso interno del CIMA
Datos Disponibles via FTP y DODS(**)

FTP: ftp://usuario:password@datos.cima.fcen.uba.ar
DODS: http://gds.cima.fcen.uba.ar:9090/dods/



Resolucion horizontal Sudamerica (SA): 0.5° x 0.5°  -  15°N a 65°S  y  30°W a 85°W       ( 110 x 160 Longitud/Latitud)
Global (GL):         2.5° x 2.5°  -  90°N a 90°S  y  0°E a 357.5°E       ( 144 x   73 Longitud/Latitud)

Resolucion vertical
(vbles de Nivel)
32 niveles de presión (mbar):
1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550
  500 450 400 350 300 250 225 200 175 150 125 100 70 50 30 20 10

Período 31 Años : 1980/01/01 a 2010/12/31  

Frecuencia  4 campos sinopticos por día : 0 6 12 18 h
 Frecuencia del Pronostico:   0 6 12 18 h

 

VARIABLES DE SUPERFICIE

Nombre

codigo CFSR

 Codigo GRIB

Producto

Unidades

Albedo
3%217-0.2-1:0.19.1
ALBDOsfc
1h AV
(0-1)
Clear sky downward longwave flux
3%217-0.2-1:0.5.196

var05196sfc

1h AV
W/m^2 s
Clear sky upware longwave flux
3%217-0.2-1:0.5.195

var05195sfc

1h AV
W/m^2 s
Clear sky downward solar flux
3%217-0.2-1:0.4.196

var04196sfc

1h AV
W/m^2 s
Clear sky upware solar flux
3%217-0.2-1:0.4.198

var04198sfc

1h AV
W/m^2 s
Convective precipitation
3%217-0.2-1:0.1.10
ACPCPsfc
6h AC
kg m^2
Dewpoint temperature
3%217-0.2-1:0.0.6

DPT

AN
K
Direct evaporation fron bare soil
3%217-0.2-1:0.3.198

var23198sfc

1h AV
W/m^2
Downward longwave radiation flux
3%217-0.2-1:0.5.192
var05192sfc
1h AV
W/m^2
Downward shortwaveradiation flux
3%217-0.2-1:0.4.192
var04192sfc
1h AV
W/m^2
Large scale precipitation (non conec.)
3%217-0.2-1:0.1.9
NCPCP
6h AC
kg/m^2
Latent heat flux
3%217-0.2-1:0.0.10
LHTFLsfc
1h AV
W/m^2
Maximun temperature
3%217-0.2-1:0.0.4
TMAX2m
6h PE
K
Minimun temperature
3%217-0.2-1:0.0.5
TMIN2m
6h PE
K
Precipitable water
3%217-0.2-1:0.1.3
PWAT
AN
kg/m^2
Pressure reduced Mean Sea Level
3%217-0.2-1:0.3.1
PRMSLmsl
AN
Pa
Sensible heat flux
3%217-0.2-1:0.0.11
SHTFLsfc
1h AV
W/m^2
Soil Moisture Content
3%217-0.2-1:2.0.3
SOILM

1h FC

kg/m^2

Soil type

3%217-0.2-1:2.3.0

SOTYP

1h FC

kg/m^2

Total precipitation

3%217-0.2-1:0.1.8

APCP

6h AC

kg/m^2

Upward longwave radiation flux
3%217-0.2-1:0.5.193
var05193sfc
1h AV
W/m^2
Upward shortwaveradiation flux
3%217-0.2-1:0.4.193
var04193sfc
1h AV
W/m^2
Volumetric Soil moisture content
3%217-0.2-1:2.0.192
var20192
1h FC
 

Water Runoff

3%217-0.2-1:2.0.5

WATRsfc

6h AC

kg/m^2

En el CIMA están disponibles, dependiendo de cada variable en particular, estos Productos:
(AN) Analisis - (FC) Forecast  - (AC) Acumulated - (AV) Average -(PE) Period

 

VARIABLES A NIVELES DE PRESION (AN:analisis)

Nombre

codigo CFSR

 nombre GRIB

Unidades

Geopotential height

3%217-0.2-1:0.3.5

HGT

gmp

Temperature

3%217-0.2-1:0.0.0

TMP

K

U wind component (Eastward)

3%217-0.2-1:0.2.2

UGRD

m s**-1

V wind component (Northward)

3%217-0.2-1:0.2.3

VGRD

m s**-1

Specific humidity

3%217-0.2-1:0.1.0
SPFH

Kg Kg**-1

Vertical velocity

3%217-0.2-1:0.2.8

VVEL

Pa s**-1

 

 

(**)DODS, Catálogo de archivos remoto; es posible navegar por sus directorios, permite inspeccionar la estructura y metadatos de los archivos almacenados. También abrir archivos a distancia (servidor OpenDAP).

Ejemplo:   Abrir un archivo con Grads (albedo)

Obtener la URL del archivo a graficar  en el DODS

Buscar en http://gds.cima.fcen.uba.ar:9090/dods/         la variable a graficar

Buscar en Reanalisis  > Sudamerica  >> Superficie >>> ALBDO >>>> info

Copiar la dirección indicada en Data URL

http://gds.cima.fcen.uba.ar:9090/dods/Reanalisis/Sudamerica/Superficie/ALBDO

Para abrir este archivo con Grads: en lugar de usar   “open”  poner “sdfopen”

>grads

ga-> sdfopen http://gds.cima.fcen.uba.ar:9090/dods/Reanalisis/Sudamerica/Superficie/ALBDO

ga-> d albdosfc

 

Nota: Usar Grads version 2.0.a6  o superior,  para manejar archivos grib2 y opendap.

Ejemplo:   Abrir un archivo con Matlab v2008/2009(tmax)

Obtener la URL del archivo a graficar  en el DODS

En matlab,  al abrir el archivo, hay que agregar el nombre de la variable y las dimensiones de los arreglos que la contienen:

http://gds.cima.fcen.uba.ar:9090/dods/

Buscar en Sudamerica  > Superficie >> TMAX  >>> info

URL:      http://gds.cima.fcen.uba.ar:9090/dods/Sudamerica/Superficie/TMAX

            Buscar:   Variables....  tmax2m: Grid...

                 tmax2m:  Array of 32 bit Reals [time = 0..45290][lat = 0..160][lon = 0..110]

   Entonces, para abrir el archivo es:

x=loaddap('

http://gds.cima.fcen.uba.ar:9090/dods/Reanalisis/Sudamerica/Superficie/TMAX?tmax2m[0:0][0:160][0:110]

')

 

Hay que poner el URL, ?, el nombre de la variable, Y  los límites de las dimensiones a recuperar en el orden indicado (tpo, lat,lon).

 

NOTA, a partir de la versión 12a, Matlab trae incluido el toolbox para abrir archivos via OpenDAP,  para versiones anteriores hay que bajar  loaddap:

Matlab loaddap Tool 3.7.2, 31 May 2012 Matlab loaddap Tool 3.7.2, 31 May 2012
 http://www.opendap.org/matlab-loaddap

Seguir las instrucciones, Y  leer el readme del loaddap, (explica como configurar matlab para que funcione):
 
-------------
Be sure to add loaddap's directory to your matlabpath! For example, you could add
 
     addpath /usr/local/bin
 
to startup.m.  The order of the paths' in the matlabpath variable is important!
However you add the paths, /usr/local/bin should appear before /usr/local/share/loaddap.
If the order is reversed you will get a error in matlab about mismatched parentheses. 
Here's what matlabpath looks like when it's set correctly:
 
<code>
     >> matlabpath
 
                MATLABPATH
 
        /usr/local/bin
        /Applications/MATLAB71/toolbox/matlab/general
..
.
</code>

Ejemplo:   Abrir un archivo con Matlab v2012a  (tmax)

(gentileza Alfredo Rolla)

Buscar el URL del archivo en el GDS (Idem anterior)

Para abrir el archivo:
===========================================
%Abrir un archivo en gds y ver dimensiones , 
%recortar una variable y graficar un tiempo de tmax

f='
http://gds.cima.fcen.uba.ar:9090/dods/Reanalisis/Sudamerica/Superficie/TMAX
';

ncid=netcdf.open(f);

[name,xtype,dimids,natts]=netcdf.inqVar(ncid,0);
name
[name,xtype,dimids,natts]=netcdf.inqVar(ncid,1);
name
[name,xtype,dimids,natts]=netcdf.inqVar(ncid,2);
name
[name,xtype,dimids,natts]=netcdf.inqVar(ncid,3);
name
%format longG;
Times=netcdf.getVar(ncid,0);
Lats=netcdf.getVar(ncid,1);
Lons=netcdf.getVar(ncid,2);

%recortamos el archivo para ver un tiempo solamente
tmax=netcdf.getVar(ncid,3,[1,1,1],[110,160,1]);
tmax=squeeze(tmax);
%imagesc(flipdim(tmax,x));figure(gcf);
imagesc(flipdim(tmax',1)-273.15);figure(gcf);
colorbar();

resultado