Interpolando dados do CMIP5 com CDO

Os arquivos de saída dos modelos globais pertencentes ao CMIP5 são extremamente diferentes entre si, além disso possui algumas peculiaridades que necessitam alguns detalhes na hora de interpolar estes dados.

O exemplo a seguir é baseado no arquivo de pressão atmosférica do Modelo HadGEM2.

A primeira coisa a fazer é descobrir quantas e quais variáveis existem em cada arquivo. Para isto utilize o seguinte comando:

> cdo sinfon arquivo.nc

Teremos como resposta algo como:

   File format: netCDF4 ZIP

    -1 : Institut     Source  Ttype                Levels  Num  Gridsize Num Dtype : Parameter name

     1 : unknown  HadGEM2-ES constant       1      1       27840      1     F32z  :   orog          

     2 : unknown  HadGEM2-ES constant     38      2              2       2     F64z  :   b_bnds        

     3 : unknown  HadGEM2-ES instant        38      2      27840       1     F32z  :   hus           

   Grid coordinates :

     1 : lonlat       > size      : dim = 27840  nx = 192  ny = 145

                        lon       : first = 0  last = 358.125  inc = 1.875  degrees_east  circular

                        lat       : first = -90  last = 90  inc = 1.25  degrees_north

                        available : xbounds ybounds

     2 : generic      > size      : dim = 2  nx = 2  ny = 0

Observe que temos três variáveis (orog, b_bnds e hus) com tamanhos diferentes, sendo que a variável b_bnds é uma variável de formato genérico (generic grid), ou seja, o cdo não consegue identificar as coordenadas geográfica da grade desta variável. Escolha a variável que você deseja interpolar, no nosso exemplo hus, e proceda da seguinte forma:

> cdo remapbil,r480x241 -selname,hus arquivo.nc output.nc

PS. Esta dica serve também para arquivos de reanálise que possuem mais de uma variável, por exemplo.

Boa Sorte e Bom Trabalho a Todos!!