Alterando os níveis verticais com CDO

Sugestão do Guilherme Martins!

O CDO sempre me surpreendendo!!!!

Alterando os níveis verticais com CDO.

Os níveis verticais dos modelos do CMIP5 estão em "Pa", e normalmente utilizamos em "hPa". Vamos usar o operador "setzaxis" para realizar essa alteração.

Carreguei neste tópico o arquivo "exemplo.nc" para vocês fazerem o teste.

Informações vista com "ncdump -h exemplo.nc"

double plev(plev) ;
plev:bounds = "plev_bnds" ;
plev:units = "Pa" ; <===== unidade do nível vertical
plev:axis = "Z" ;
plev:positive = "down" ;
plev:long_name = "pressure" ;
plev:standard_name = "air_pressure" ;

Informações vista com "ncdump arquivo.nc | more"

plev = 100000, 85000, 70000, 50000, 25000, 10000, 5000, 1000 ; => "Pa"

Segue a dica:

Crie um arquivo texto, e adicione as 6 linhas abaixo:

zaxistype = pressure
size = 8
name = lev
longname = pressure
units = hPa
levels = 1000 850 700 500 250 100 50 10

Salve com um nome qualquer, por exemplo, "nivel". Esse arquivo será lido pelo operador "setzaxis":

Utilize o comando abaixo:

cdo setzaxis,nivel exemplo.nc saida.nc

Esse comando altera a estrutura do seu arquivo netCDF com as novas informações em "hPa".

Para visualizar as novas alterações:

ncdump -h saida.nc

Arquivo: exemplo.nc

PS: caso o arquivo tenha dados UNDEF (que é caso), utilize o operador fillmiss ==> cdo -fillmiss -setzaxis,nivel exemplo.nc saida.nc

Boa Sorte e Bom Trabalho a Todos!