Cristiano Prestrelo

Criando Ensemble mean dos dados do CMIP5

2016-12-01T14:47:00.000-02:00

Trabalhar com dados do CMIP5 pode ser um pouco confuso pois nem sempre eles apresentam estruturas iguais. No meu caso eu tenho um conjunto de arquivos de três modelos (CanESM2, CSIRO e HadGEM2) do CMIP5 e estes dados estão dispostos em arquivos de 5-10 anos, ou seja, uma série de 100 anos temos 10-20 arquivos.

Para juntar todos os arquivos em um só eu utilizei o NCO (Netcdf Operator) através do comando ncrcat. O script abaixo irá juntar todos os arquivos, de cada modelo, em apenas um arquivo.

#!/bin/bash

#########################################################

# O nome dos meu arquivos tem a seguinte estrutura:

# tos_Omon_CanESM2_esmHistorical_r1i1p1_199001-200012.nc

# $var $modelos $xpt

#########################################################

dir_inp=( '/dados/data/ESM-PP/SST/Historical/' ) # Diretório

var=( 'tos_Omon' ) # Variável

modelos=( 'CanESM2' 'CSIRO' 'HadGEM2' ) # Pasta/Modelos

xpt=( 'esmHistorical_r1i1p1' ) # Ensemble

# loop para contar o número de arquivos referente a cada modelo

for mdl in ${modelos[@]}; do

# numero de arquivos neste ensemble member

num_arq=$( ls ${dir_inp}${var}_${mdl}_${xpt}*.nc | wc -w )

echo "" $mdl

# se houver apenas um arquivo continue para o próximo loop

if [ ${num_arq} -le 1 ]; then

echo "Há apenas um arquivo em" ${dir_inp}.

continue

echo "Existem" ${num_arq} "arquivos no" ${dir_inp}

# Data inicial de cada arquivo (o nome do primeiro arquivo inclui a data inicial)

yyyymm_in=$( ls ${dir_inp}${var}_${mdl}_${xpt}*.nc | sed -n '1p' | cut -d '_' -f 6 | cut -d '-' -f 1 )

# Data final de cada arquivo ( o nome do último arquivo inclui a data final)

yyyymm_fi=$( ls ${dir_inp}${var}_${mdl}_${xpt}*.nc | sed -n "${fl_nbr}p" | cut -d '_' -f 6 | cut -d '-' -f 2 )

echo $yyyymm_in, $yyyymm_fi

## Juntando todos os arquivs em apenas um só arquivo

ncrcat -O ${dir_inp}${var}_${mdl}_${xpt}*.nc ${dir_inp}${var}_${xpt}_${mdl}_${yyyymm_in}-${yyyymm_fi}

done

exit 0

Este script calculou em primeiro lugar os pesos para cada célula da grade do modelo, antes de usar os pesos para realizar a interpolação cúbica da grade nativa para a nova grade regular. Eu segui as mesmas etapas para todos os meus arquivos e reuniu os comandos em um script que irá varrer todos os arquivos no diretório definido. Os arquivos do CMIP5 estão em formato de Netcdf porém, o espaçamento de grade destes arquivos são diferentes. Portanto, para calcular uma grade retangular genérica para todos os dados para poder comparar os arquivos e executar operações com matrizes a solução foi utilizar o cdo para gerar esta grade retangular, com latitudes de -90 a 90 e longitude de 0 a 360.

Para isto execute estes comandos:

#/bin/sh
modelos=( 'CanESM2' 'CSIRO' 'HadGEM2' ) # Pastas para cada modelo
dir_inp=( '/dados/data/ESM-PP/SST/Historical/$modelos/' )
num_arq=$( ls ${dir_inp}*.nc | wc -w )
echo " Foram encontrados " num_arq " no diretório "

modelos=$( ls *.nc )

for model in ${models[@]}; do
echo "Degradando o modelo: " $model
nome=$( ls $model | cut -d '.' -f 1 )

cdo genbil,r360x180 ${dir_inp}${var}_${xpt}_${mdl}_${yyyymm_in}-${yyyymm_fi} ${dir_inp}${var}_wgts_${xpt}_${mdl}_${yyyymm_in}-${yyyymm_fi}
cdo remap,r360x180,${dir_inp}${var}_wgts_${xpt}_${mdl}_${yyyymm_in}-${yyyymm_fi} ${dir_inp}${var}_${xpt}_${mdl}_${yyyymm_in}-${yyyymm_fi} ${dir_inp}${var}_regular_${xpt}_${mdl}_${yyyymm_in}-${yyyymm_fi}
done

Em seguida, era preciso calcular as médias anuais de TSM em todos os arquivos, pois isso torna um pouco mais fácil de usar comandos do CDO depois. Por isso, foi feito um loop sobre todos os arquivos (nova grade) calculando a média anual. Uma vez que todas as grades nativas são diferentes, e alguns modelos têm valores fill_value ou mascaras isso pode gerar um problema quando você for calcular o ensemble mean. Para não ter problema o programa faz o loop apenas sobre os valores presente em cada ponto de grade e executa a média. Por isso, é necessário definir o intervalo de valores aceitáveis no arquivo que é, do ponto de congelamento da água da superfície do mar até um valor ridiculamente alto. O comando 'setvrange' do cdo cuida disso para nós antes de calcular as médias anuais.

#/bin/sh
modelNumbers=$( ls *.nc | wc -w )
echo "Encontrados " ${modelNumbers}
rm -f YRLYAVG/*.nc
modelos=$( ls *.nc )

for model in ${modelos[@]}; do
echo "Calculando media anual: " $model

name=$( ls $model | cut -d '.' -f 1 )

cdo -setvrange,271.45,320 $model $model.tmp
cdo yearmean -selyear,1950/2005 -selvar,tos $model.tmp ${name}_yrly.nc
rm $model.tmp

cdo sinfo ${name}_yrly.nc
mv ${name}_yrly.nc YRLYAVG/.

done

cd YRLYAVG

Este simples comando calcula o ensemble mean de todos os arquivos de médias anuais cada um representando um cenário do modelo.

cdo ensmean *.nc ensemble_media_historico.nc

Em seguida, repetimos todos esses comandos tanto para as projeções futuras RCP8.5 quanto para os arquivos históricos. Como preparação final desses arquivos antes da análise eu queria estimar a média de longo prazo do ensemble mean. Primeiro, é preciso calcular a climatologia e foi escolhido 1961-1990 como climatologia de referência.

cdo timselavg,30 -seldate,19610101,19900101 ensemble_media_historico.nc ensemble_media_1961_1990.nc
cdo sinfov ensemble_media_1961_1990.nc

Estes comandos calcula a climatologia das simulações históricas dos modelos. Em seguida, repete os comandos para os períodos de 2006-2025 e 2025-2050 para obter as projeções futuras.

Boa Sorte e Bom Trabalho a Todos!!

Formato de Dados Meteorológicos

2016-10-23T14:22:00.000-02:00

Existe alguns tipos de dados que são facilmente encontrados no dia a dia dos estudos Meteorológicos, entre eles destacam-se: GRIB, netCDF e HDF. Duas características importante comuns a estes dados é que eles são portáteis e auto-descritivos. Portáteis pois podem ser manipulados em qualquer máquina (AMD, Intel, ...) e em qualquer sistema (Windows, Linux, ...). Auto-descritivos pois podem ser examinados e lidos por qualquer software (apropriado) sem que o usuário saiba detalhes estruturais do arquivo. Além disso, informações adicionais sobre os dados, chamada de “metadado (dado sobre dado)” pode ser incluída no arquivo. Metadados típicos podem incluir informações textuais sobre o conteúdo e sobre as unidades de cada variável (ex., “Temperatura” e "C") ou informações numéricas que descrevem as coordenadas (p.e., tempo, nível de altura, latitude e longitude) que se aplicam as variáveis.

O formato destes dados tem evoluído ao longo do tempo com o objetivo de atender as necessidades da comunidade científica e civil. Por isso, alguns formatos possuem várias versões anteriores. Infelizmente para as pessoas que utilizam estes dados, as versões mais recentes nem sempre são compatíveis com as versões anteriores, apesar de fazer parte do mesmo formato de dado (e com mesma nomenclatura). Esta característica tem gerado confusão e críticas por parte dos usuários, pois a cada versão lançada exige uma adaptação necessária.

Os dados mais conhecidos são:

GRIB1: GRIdded Binary (Edition 1), World Meteorological Organization

GRIB2: GRIdded Binary (Edition 2), World Meteorological Organization

netCDF3: Network Common Data Form, (Version 3.x), Unidata (UCAR/NCAR)

netCDF4: Network Common Data Format, (Version 4.x), Unidata (UCAR/NCAR

HDF4: Hierarchical Data Format, (Version 4.x), NCSA/NASA

HDF4-EOS2: HDF4-Earth Obseving System, (Version 2; georeferenced data)

HDF5: Hierarchical Data Format, (Version 5.x), NCSA/NASA

HDF5-EOS5: HDF5-Earth Obseving System, (Version 5; georeferenced data)

GeoTIFF: Georeferenced raster imagery

Existe uma pequena diferença entre netCDF/HDF/HDF-EOS e GRIB, no caso do primeiro conjunto de dados estamos falando de formatos de arquivos de dados enquanto que no caso do GRIB é um formato específico de gravação de dados. Os tipos de arquivos netCDF/HDF/HDF-EOS são formatos de arquivos que possuem regras sobre o conteúdo de cada arquivo. Por exemplo, no caso dos arquivos NetCDF uma regra existente é que todos os nomes das variáveis devem ser exclusivas (isto é uma regra). No caso do HDF ele permite que multiplas variáveis possuam mesmo nome, porém as variáveis devem estar em diferentes “grupos” de variáveis. Esses “grupos” podem ser complicados, mas de forma simples pode ser comparado ao sistema de diretório do Linux (ou pastas do Windows). No caso dos dados tipo GRIB, cada registro do GRIB-1 (aka, 'message') contém informações para duas dimensões horizontais (p.e., Latitude e Longitude) para um tempo e um nível. Já os arquivos GRIB-2 permite que cada “registro” contenha várias grades e níveis para cada tempo. Portanto, um conjunto de “registros GRIB” é chamado de “arquivo GRIB”. Diferente dos arquivos netCDF/HDF/HDF-EOS os arquivos GRIB não possuem regras que determinem a ordem deste conjunto de “registros GRIB” (p.e., estes registros podem estar em ordem cronológica aleatório).

Um exemplo da estrutura de grupos que contém um arquivo HDF.

NetCDFType

Segue uma lista de exemplos de como ler e escrever estes tipos de dados utilizando o NCL.

Neste link você pode ver outro conjunto de scripts NCL para ler alguns dados específicos.

Mais detalhes:

netCDF
HDF
GRIB
Binary
Arquivo Texto (ASCII)

Fonte: NCAR.

Boa Sorte e Bom Trabalho a Todos!

Espiral da Temperatura Global

2016-10-08T22:23:00.002-03:00

Recentemente o Pesquisador Ed Hawkins do National Centre for Atmospheric Science (NCAS) publicou em sua página do twitter uma animação espiral do aumento da Temperatura Global e esta animação já compõem várias apresentações de vários pesquisadores em diversos eventos. Para quem deseja utiliza-la em suas aulas, apresentações e etc estou disponibilizando os links originais.

Espiral do aumento da Temperatura Média Global. Clique para ver em tamanho maior.

A animação em forma de espiral representa mudança da temperatura média global de uma forma visualmente atraente e simples. O ritmo da mudança é imediatamente óbvio, especialmente ao longo das últimas décadas. A relação entre as temperaturas globais atuais e os limites alvo discutidos internacionalmente também são claros sem muita interpretação complexa necessária.

Os dados são do HadCRUT4.4 em representam a Anomalia da Temperatura de Janeiro de 1850-Março de 2016, em relação a média de 1850-1900 e pode ser acessado aqui.

Características que você pode observar:

1877-1878: El Nino intenso responsável pelo aumento da temperatura;
1880-1910: pequeno resfriamento, parcialmente devido a erupções vulcânicas;
1910-1940: aquecimento, em parte devido à recuperação de erupções vulcânicas, pequeno aumento na atividade solar e variabilidade natural;
1950-1970: temperaturas praticamente inalteradas com o resfriamento através dos aerossóis de sulfato mascarando o aquecimento dos gases do efeito estufa;
1980-hoje: forte aquecimento, com aumento maior das temperaturas em 1998 e 2016 devido a eventos fortes de El Nino.

Porque inicia em 1850?

Porque é o início do conjunto de dados do HadCRUT4, e não temos dados de temperatura, mais antigos, suficientes para construir de forma confiável a temperatura média global.

O que as cores significam?

As cores representam o tempo. Roxo para os primeiros anos, passando pelo azul, do verde para o amarelo para a maioria dos últimos anos. A escala de cores utilizado é chamado de "viridis" e as ilustrações foram feitas em MATLAB.

Além dessa animação Ed também disponibilizou uma animação espiral da variação do Gelo Marinho do Ártico.

Volume diário do gelo marinho do Ártico estimado pela reconstrução PIOMAS de 1979-presente, produzindo uma espiral da variação do volume de gelo marinho.

Por fim, segue a animação da variação da concentração de CO2 na atmosfera medido na estação de Mauna Loa no Hawaii.

Boa Sorte e Bom Trabalho a Todos!!

Quanto vale a Meteorologia?

2016-07-12T21:35:00.000-03:00

Quanto vale $$ uma Previsão de tempo ou clima? Quanto vale um mapa de vento? Quanto vale uma consultoria em Meteorologia?

Recentemente fui convidado para ministrar uma palestra sobre o Mercado de Trabalho de Meteorologia, curioso que este sempre foi um assunto que me chamou muita a atenção, por motivos óbvios. Neste texto eu pretendo ter a pretensão de trazer uma reflexão para todos nós Profissionais/Estudantes de Meteorologia sobre o mercado de Meteorologia.

Fiz uma busca na internet, mas não encontrei informações sobre as oportunidades de emprego/empresas disponíveis atualmente no Brasil e fazer uma estimativa destas oportunidades poderia ficar muito fora da realidade. Portanto, este texto não trata das oportunidades/vagas/quantidades, mas do mercado em si, suas características, empregado, empregador e ensino.

Todo mercado de trabalho tem suas características peculiares influenciadas pelo "sentimento" do empregado e do empregador, mas principalmente pelo perfil do formador do profissional. Alguns setores são predominantemente privado outros apresentam uma proporção um pouco mais dividida e outros são efetivamente Públicos. Onde nós estamos?

Atualmente o mercado de trabalho em Meteorologia oferece oportunidades de emprego no setor Público e no setor Privado, estando no setor Público as maiores oportunidades. Além disso, todos nós sabemos que o ensino de Meteorologia nível médio/superior é exclusivamente Público. Por estes motivos identificamos que a primeira ou a principal ou até mesmo a escolha mais segura dos Profissionais de Meteorologia está no setor Público. Entretanto, tenho observado que a grande maioria dos funcionários públicos ou dos pretendentes ao setor nos cargos de Meteorologia demonstram uma nítida falta de conhecimento do que é ser um servidor público. Entendo e aceito que pode não ser a maioria e eu possa estar enganado, mas este não é o ponto. Desta forma, é válido falarmos um pouco sobre o que é ser ou o que deveríamos levar em consideração ao escolhermos seguir uma carreira no serviço público:

Antes de mais nada permitam-me fazer uma observação: definitivamente somos o país do contraditório. Nos últimos anos, observamos o surgimento de uma febre na busca por concursos públicos. Na esteira desta febre surgiu-se um mercado de cursinhos e publicações especializadas no assunto, criou-se até um adjetivo para qualificar aqueles que almejam uma vaga no serviço público, CONCURSEIROS. O contraditório está no fato deste povo todo escolher representantes, presidentes principalmente, que pregam o estado mínimo, ou seja, que não querem tantas vagas de concursos assim. Uma enorme contradição, mas o tema aqui não será política.

Antes de decidirmos seguir ou não uma carreira no serviço público deveríamos (Meteorologistas ou não) ser capazes de responder duas perguntas simples:

Qual a motivação de um profissional escolher o serviço público?

R: Carreira? Ambiente de Trabalho? A qualidade/excelência das Instituições? Salário?...ESTABILIDADE?

Qual a sua opinião sobre a carreira e principalmente sobre qualidade do serviço público no Brasil hoje?

R: Funciona?, O atendimento é satisfatório? Tem qualidade? É ético?

É possível com essas motivações (Carreira, Salário, ESTABILIDADE e etc) resolver os atuais problemas (serviço de baixa qualidade) do serviço público?

E se mudarmos as motivações do servidor público para: Vocação, Ideologia, Patriotismo, vontade de ajudar ao próximo e principalmente VONTADE de trabalhar para melhorar o SERVIÇO Público? Pois bem, a decisão de seguir no setor Público não pode ser somente pela carreira, salário, estabilidade e etc, mas sim principalmente pela ideologia. Garanto que em algum momento das nossas carreiras encontramos funcionários públicos infelizes com seu trabalho e vemos está infelicidade refletindo nas suas atividades diárias ou até mesmo na falta delas.

Uma pesquisa recente mostrou que o motivo de quase 15 milhões de jovens brasileiros escolherem fazer um concurso não é vocação, ideologia ou patriotismo e sim ESTABILIDADE, mas batalhar por Estabilidade não é pouco para quem é inteligente e determinado?

Será que alguém passaria em um concurso público sem dedicação e inteligência/conhecimento?

Não seria um anseio conservador e de pouca ambição batalhar por: Estabilidade?

Se todo mundo decidisse ser Servidor Público, teríamos 200 milhões de vagas? Afinal de contas alguém precisa pagar imposto.

Existe a tal ESTABILIDADE?

Em 2013 os EUA mandaram para casa sem salários mais de 800 mil funcionários públicos;

Na Grécia mais de 100 mil foram demitidos.

A Estabilidade NÃO existe! Precisamos aprender a viver com essa realidade.

Portanto, o problema identificado mostra que temos Profissionais de Meteorologia que escolhem seguir no serviço público sem ter perfil e as consequências destas escolha têm sido desastrosas para Meteorologia tanto no setor operacional quanto (principalmente) no ensino.

Raciocine um pouco. Como foi sua formação dentro dos cursos técnicos ou das Universidades de Meteorologia? Seus professores aconselharam, seu departamento deu condições para você seguir no setor Privado? Se o único caminho lhe indicado foi o serviço Público você foi orientado ou preparado para tal? Alguém lhe falou que você poderia criar uma empresa de Meteorologia ou trabalhar como Consultor? Vejo com muita preocupação essa falta de informação durante a formação dos nossos profissionais de Meteorologia.

Porque isso seria uma preocupação? Existe um concurso aberto para Professor do curso de Meteorologia da UFAL (07/2016) com duas vagas onde o requisito para inscrição é ser Graduado em Meteorologia e temos aproximadamente 30 candidatos (Doutores) inscritos. São no mínimo 30 Doutores "desempregados". Em tempos de crise econômica nós sabemos que não existe previsão do surgimento de novas vagas para novos concursos, estas vagas da UFAL são provenientes de aposentadoria.

Existem na minha visão duas consequência graves deste "encaminhamento" dos Profissionais de Meteorologia para o setor público, que são eles: a visão de que os serviços de Meteorologia são gratuitos e empecilho do desenvolvimento do setor privado de Meteorologia do Brasil. Estas consequência influenciam fortemente o consumidor dos serviços de Meteorologia.

Para mim o maior dos problemas gerado por essa cultura é a nossa Previsão de Tempo. A Previsão de Tempo no Brasil NÃO tem valor (R$) tanto o setor público (CPTEC, INMET, etc) quanto no privado (Climatempo, SOMAR, etc) a Previsão de Tempo é distribuída gratuitamente em seus portais ou produtos.

Percebam que nas ementas dos cursos de Meteorologia espalhados pelo Brasil a principal "obrigação" é formar Profissionais de Meteorologia capazes de fazer uma Previsão de Tempo, ou seja, a principal habilidade de um Meteorologista Graduado é fazer Previsão do Tempo. Habilidade essa distribuída gratuitamente. Agora como dizer para um estudante de Meteorologia que existe "vida" fora do serviço público se ao mesmo tempo eu digo para ele que a principal habilidade que ele tem é oferecida de forma gratuita pela comunidade? Imagina informar a um Engenheiro Civil que a construção de prédios é gratuita? Ou a um arquiteto que os projetos de casas, prédios e etc é distribuída de forma gratuita?

Diante deste cenário de crise e incertezas eu pergunto. Estamos preparados para ser um consultor ou empreendedor de Meteorologia? Quanto custa uma Previsão de tempo ou clima? Como eu calculo o preço de uma simulação climática ou um mapa de vento? Isto é ensinado nas salas de aula? Chutar um valor entre 10 e 100 mil é saudável para o setor privado? Os preços praticados pelas empresas privadas estão dentro da realidade?

Diante de um cenário de crise econômica e de poucas oportunidade de concursos me preocupa muito a existência destes diversos fatores que prejudicam o crescimento do setor privado em Meteorologia no Brasil. Sabemos que as instituições Públicas disponibilizam de forma gratuita para população outros produtos além da Previsão de Tempo. Além disso, há relatos de que destas instituições concorrendo com as instituições privadas por projetos e por espaço no mercado privado.

As Universidade continuam formando Meteorologistas os cursos técnicos continuam formando técnicos em Meteorologia as pós-graduação continuam formando Mestres e Doutores em Meteorologia sabemos que não teremos concursos públicos nos próximos anos, sabemos que as instituições de fomento (CNPQ, CAPES e as Fundações de Amparo a Pesquisa) não terão Bolsas suficientes para manter tantos Meteorologistas e que cada vez mais criamos mecanismos para enfraquecer o setor privado.

Não passou da hora de discutirmos o fortalecimento do setor privado de Meteorologia? Não passou da hora de utilizarmos o setor Público de Meteorologia para fortalecer o setor privado e não o contrário? Não passou da hora de prepararmos nossos Profissionais de Meteorologia para o setor privado ou para criar o setor privado de Meteorologia do Brasil? Não passou da hora dos cursos de Meteorologia do Brasil terem disciplinas voltadas a administração e empreendedorismo? Vale a ressalva o curso de Ciências Atmosféricas da UNIFEI é uma excelente exceção neste ponto.

Precisamos pensar urgentemente no mercado de trabalho brasileiro de meteorologia. As instituições Públicas precisam encontrar seu melhor papel neste mercado. As empresas de meteorologia precisam valorizar seus produtos e sua mão de obra. Os meteorologista precisam saber que existe outras opções além da pós. As instituições de ensino precisam preparar seus estudantes para o setor privado. Para aqueles que acham que criar um mercado de trabalho, criar seu próprio emprego empreender em Meteorologia é muito difícil eu lhes digo que fazer uma pós viver de Bolsa esperar um concurso público não é nada fácil.

Para quem tem dúvida se Meteorologia dar dinheiro existe uma estimativa de que o setor privado de Meteorologia rendeu no ano de 2014 algo próximo de 30 milhões de reais no ano, detalhe este setor está concentrado na região sudeste. Além disso, para quem não sabe a IBM comprou a Weather Channel no ano passado por 2 bilhões de dólares.

Boa Sorte e Bom Trabalho a Todos!

Definição de El nino e La Nina

2016-03-17T11:29:00.001-03:00

O Termo El Niño é utilizado desde os anos 1600 pelos pescadores da costa do Equador e Peru para se referir a corrente quente do Oceano Pacífico que aparecia próximo do Natal e durava alguns meses. O termo Niño em espanhol durante a época do Natal é utilizado para se referir ao menino Jesus, todavia esta corrente quente próximo da costa diminuía a quantidade de alimento dos peixes obrigando os pescadores a se afastarem cada vez mais da costa para encontrar o pescado ou muitas vezes utilizar este tempo para reparar seus equipamentos ou ficar com suas famílias. No entanto, em alguns casos este período de correntes mais quentes perduravam até Maio ou Junho o que provocava uma ruptura na temporada de pesca, por isso por muitas vezes o El Niño era chamado de "menino travesso".

A definição do fenômeno El Niño evoluiu muito ao longo do tempo e no início provocou muita confusão sobre seu real significado. No ano de 1997 o pesquisador Kevin E. Trenberth do NCAR (National Center for Atmospheric Research) escreveu um artigo a pedido do CLIVAR (Climate and Ocean Variabilit, Predictability and Change) do WCRP (Wolrd Climate Research Programme) para tentar definir tal fenômeno. Na época foi feito uma consulta com a comunidade científica e uma tentativa (bem sucedida por sinal) de resumir as principais características sabidas na época sobre o El Niño foi feita. Este ano é emblemático pois se trata do ano em que ocorreu o mais intenso El Nino já registrado (iremos falar dele em breve link).

Ao longo dos anos as pesquisas científicas mostraram que, além do El Niño existia um período em que as condições meteorológicas observada eram inversa ao padrão observado durante eventos de El Niño ao qual foi chamada de La Niña, anos depois as pesquisas mostraram que estes padrões influenciavam a atmosfera o que levou a criação do termo ENOS (El Niño Oscilação Sul) do inglês ENSO (El Nino South Oscillation). Portanto, para definirmos o fenômeno precisamos ao mesmo tempo definirmos a sua fase oposta ou La Niña e o ENOS.

Definição do El Niño, La Niña e ENOS

O El Niño é caracterizado pelo enfraquecimento de larga escala dos ventos alísios e aquecimento das camadas superficiais da região leste e central do oceano Pacífico. Além disso, durante um evento de El Niño a pressão atmosférica sobre a região do Pacífico central é mais baixa. Portanto, ao descobrirmos que durante um evento de El Niño observamos mudanças na Temperatura da Superfície do Mar (TSM), nos ventos Alísios e na Pressão ao Nível Médio do Mar (PNMM) chegamos a conclusão que se trata de fenômeno de interação entre o oceano e a atmosfera que conhecemos como ENOS.

Antes de continuarmos vamos abrir aspas ("") para falarmos um pouco sobre estas três características anteriormente citadas (Ventos Alísios, TSM e PNM).

Ventos Alísios

O termo trade winds (Ventos de Comércio, tradução livre) é utilizado para se referir aos ventos utilizados pelos navegadores e seus navios a vela na época das grandes navegações, este ventos foram muito importantes para a expansão do império europeu para as Américas possibilitando criar rotas da Europa para as Américas e regiões do oceano Pacífico. Os ventos alísios são ventos que ocorrem na região tropical durante todo o ano. Os alísios são resultado do fluxo associado ao cinturão de alta pressão subtropical (altas semi-permanentes), localizados nas latitudes de 30ºN e 30ºS (Figura 1), que apresentam orientação de sudeste no Hemisfério Sul e de Nordeste no Hemisfério Norte (representado pelas setas vermelhas na parte superior da Figura 1) a confluência destes dois fluxos na região equatorial produz um fluxo predominantemente de leste localizados próximos da superfície e a convergência de massa produz a Zona de Convergência Intertropical ZCIT (região equatorial da Figura 1).

Figura 1: Representação esquemática horizontal (parte superior) dos cinturões de alta pressão e dos ventos Alísios e vertical (parte inferior) de parte da circulação geral.

2. Temperatura da Superfície do Mar - TSM

A Temperatura da Superfície do Mar medida através de Boias, Navios e após 1979 Satélites são importantes para determinar as condições de El Niño e La Niña. Para estudarmos o fenômeno El Niño nós meteorologistas calculamos Índices a partir da anomalia TSM o mais conhecido deles é o Oceanic Niño Index (ONI). O ONI se refere uma série temporal trimestral da anomalia de TSM sobre as regiões que nós identificamos o El Niño Figura 2.

Figura 2: Regiões referente ao tipo de El Niño identificado.

3. Pressão ao Nível Médio do Mar - PNMM

A Pressão ao Nível Médio do Mar é uma variável extremamente importante na identificação do fenômeno El Niño. Durante a ocorrência do El Niño observa-se a existência do que nós Meteorologistas chamamos de "gangorra barométrica" está gangorra está relacionada a uma diferença na PNMM em duas regiões do Pacífico; uma sobre a Ilha do Tahiti (17ºS-149ºW) e outra sobre o norte da Austrália (Darwin: 12ºS-130ºE) Figura 3. A diferença entre a PNMM medida na estação meteorológica do Tahiti e a PNMM medida na estação de Darwin produz um Índice chamado de Southern Oscillation Index (SOI) esta diferença começou a ser documentada pelos pesquisadores Sir Gilber Thomas Walker e Bliss e registrada em seus artigos publicados em 1930, 1932 e 1937. O Índice SOI é uma medida das flutuações de grande escala na pressão do ar que ocorrem entre o Pacífico tropical oeste e leste durante os episódios de El Niño e La Niña (Figura 4). A fase negativa do SOI representa a pressão do ar abaixo do normal sobre o Tahiti e pressão do ar acima do normal sobre Darwin. Períodos prolongados de valores negativos (positivos) do SOI coincidem com TSM positivas (negativas) em todo o leste do Pacífico tropical características típicas observadas durante episódios de El Niño (La Niña).

Figura 3: Localização das estações Meteorológicas do Tahiti e Darwin.

Figura 4: Valores mensais do SOI. A Linha preta refere-se a uma média trimestral do SOI.

Portanto, o termo El Niño refere-se ao aquecimento anômalo da TSM observado nas regiões destacadas na Figura 2. No entanto, a fase positiva do ENOS refere-se à interação do oceano com atmosfera do clima em larga escala ligados a um aquecimento periódico das TSM nas regiões do Pacífico central e leste e mudanças na intensidade dos ventos Alísios.

Por outro lado, o termo La Niña, evento frio do ENOS ou fase fria do ENOS são usados para se referir aos períodos em que as anomalias da Temperatura da Superfície do Leste e Centro do oceano Pacífico são negativas e os padrões dos ventos Alísios estão mais fortes. O termo "El Viejo" e "anti-El Niño" também foram termos utilizados para se referir a fase fria do ENOS. No entanto, estes termos são utilizados com menos frequência e o termo La Niña é o mais comum.

Uma última característica a ser observada refere-se a Circulação de Walker. Na Figura 1 podemos observar uma representação simplificada do fluxo de ar sobre a região do Equador durante o inverno (verão) do hemisfério norte (sul). Este padrão de ventos associados a circulação de Walker é observado durante as condições médias (normais) conhecido como fase neutra do ENOS. As características importantes observadas são: o ramos ascendentes da circulação Walker intenso (mais forte) em todo o continente Marítimo (Austrália, Indonésia e etc) e mais fraco sobre o oste da África e no norte da América do Sul. As regiões onde observam-se movimento ascendentes (descendentes) estão associadas a formação (inibição) da precipitação para o trimestre Dezembro-Janeiro-Fevereiro.

Figura 5: Representação simplificada da Circulação de Walker durante a fase Neutra do ENOS durante o Verão do Hemisfério Sul (DJF). Fonte ENSO-blog.

Durante a fase quente do ENOS o padrão da circulação de Walker sofre mudanças significativas Figura 6. As cores mais quente na Figura 6 representa as anomalias positivas de TSM associadas ao fenômeno El Niño. Observe associado as anomalias de TSM está deslocam dos ramos ascendentes (descendentes) localizados sobre o continente Marítimo e sobre a América do Sul (Pacífico central e Índico oeste) associados a circulação de Walker. Além disso, observa-se o enfraquecimento dos ventos Alísios de leste próximo da superfície. Esta mudanças no padrão de circulação também provoca mudança no padrão de precipitação que se tornam mais intensas (fracas ou inexistentes) sobre o Pacífico central e norte da África (América do Sul e Continente Marítimo).

Figura 6: Representação simplificada da Circulação de Walker durante a fase Quente do ENOS durante o Verão do Hemisfério Sul (DJF). Fonte ENSO-blog.

Na fase fria do ENOS período em que as águas do Pacífico oeste (continente Marítimo) são ainda mais quente em comparação a fase neutra do ENOS (Figura 5) e as águas do Pacífico leste são ainda mais frias. O padrão da circulação de Walker observado durante a fase fria do ENOS é semelhante ao padrão observado durante a fase neutra do ENOS porém mais intensificado. Observa-se que durante esta fase o ar úmido e quente se eleva ainda mais sobre o continente marítimo e sobre a América do Sul elevando a precipitação a valores acima da média. No leste do oceano Pacífico, onde a TSM é mais fria que a média, o ramo descendente da circulação de Walker mais intenso contribui com a redução maior dos totais pluviométricos da região que normalmente já são baixos.

Figura 7: Representação simplificada da Circulação de Walker durante a fase Fria do ENOS durante o Verão do Hemisfério Sul (DJF). Fonte ENSO-blog.

Este é o primeiro de uma série de post que irei fazer sobre El Nino, La Nina e ENOS. Os próximos post irão tratar dos impactos de ENOS sobre a variabilidade climática do Globo e América do Sul, a previsão do El Nino e as condições do atual El Niño. Em breve!

Referências Bibliográficas:

http://www.cpc.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.shtml

https://www.ncdc.noaa.gov/teleconnections/enso/indicators/soi/

https://en.wikipedia.org/wiki/Trade_winds

https://en.wikipedia.org/wiki/Walker_circulation

http://www.rmets.org/publications/classic-papers

http://www.nature.com/news/2009/090923/full/461455a.html

http://biophysics.sbg.ac.at/atmo/elnino.htm

https://www.climate.gov/news-features/department/enso-blog

http://www.pmel.noaa.gov/tao/elnino/el-nino-story.html

Trenberth, Kevin E. "The definition of el nino." Bulletin of the American Meteorological Society 78.12 (1997): 2771-2777.

Walker, G.T. and Bliss, E.W., 1932. World Weather V, Memoirs of the Royal Meteorological Society, 4, (36) , 53-84.

Walker, G.T. and Bliss, E.W., 1930. World Weather IV, Memoirs of the Royal Meteorological Society, 3, (24), 81-95.

Walker, G.T. and Bliss, E.W., 1937. World Weather VI, Memoirs of the Royal Meteorological Society, 4, (39), 119-139.

O Exercício Ilegal da Profissão de Meteorologista

2016-03-09T09:33:00.000-03:00

Nas últimas décadas a profissão de Meteorologista tem chamado cada vez mais a atenção da população brasileira. Infelizmente muitas vezes o profissional de Meteorologia somente é lembrado quando algum desastre natural grave acontece, tais como grandes enchentes, secas severas, eventos de ventos fortes e etc.

Ao longo do tempo os estudos sobre Mudanças Climáticas tem colocado a Meteorologia em destaque no mundo. Grande eventos, grandes campanhas e intensa repercussão em todo tipo de mídias aumentou significativamente, o número de pós-graduandos e pesquisadores estudando sobre assuntos relacionados a Meteorologia (Mudanças Climáticas).

"Este interesse da sociedade em relação a Meteorologia tem atraído mais a atenção dos estudantes de ensino médio para o curso de Meteorologia, claro que a procura ainda é pequena e oscila muito ao longo dos anos. Por se tratar de uma ciência relativamente nova no mundo e ainda mais jovens no Brasil as pessoas não tem muito conhecimento sobre a formação de um Meteorologista. Esta característica é responsável pela desistência de muitos alunos pois possuem uma visão incompleta do curso de graduação em Meteorologia, ou seja, ao se deparar com disciplinas de Cálculo e Física por exemplo alguns desistem. Em um outro post eu irei abordar um pouco sobre as características da carreira de um Meteorologista link."

Acompanhado do aumento do interesse da sociedade, dos alunos de ensino médio e principalmente dos alunos de pós-graduação surgiu um problema que já fora observado em algumas profissões no Brasil e este problema é grave e a legislação tipifica como "Exercício Ilegal da Profissão". Vale uma ressalva aqui: obviamente eu não sou jurista e não tenho conhecimento total da lei e da interpretação dela, portanto agradeço as opiniões e correções. Ultimamente temos observado uma inquietação por parte da comunidade meteorológica em relação a algumas pessoas que estão exercendo ilegalmente a profissão de Meteorologista. A identificação dos casos de exercício ilegal da Profissão está associada principalmente a prática da Previsão de Tempo e Clima.

Segundo a Lei Federal nº 6835 de 14 de Outubro de 1980 Art. 7º - São atribuições do meteorologista:

...
.
.
.

d) executar previsões meteorológicas;
.
.
.

...

Pessoas que alegam serem entusiasta ou apaixonados pela Meteorologia tem repetidamente produzido/reproduzido conteúdos na internet/radio, entre outras mídias, relacionados a práticas atribuídas apenas a Meteorologistas e em alguns relatos estas práticas são remuneradas ($$$). E isto obviamente tem incomodado a comunidade de Meteorologia e seus representantes.

O que sabemos sobre o Exercício Ilegal da Profissão?

Segundo a Seção III do Capítulo 1 da LEI Nº 5.194, DE 24 DEZ 1966, que trata do Exercício Profissional da Engenharia, da Arquitetura e da Agronomia:

Art. 6º - Exerce ilegalmente a profissão de engenheiro, arquiteto ou engenheiro-agrônomo:

a) a pessoa física ou jurídica que realizar atos ou prestar serviços, públicos ou privados, reservados aos profissionais de que trata esta Lei e que não possua registro nos Conselhos Regionais:

b) o profissional que se incumbir de atividades estranhas às atribuições discriminadas em seu registro;

c) o profissional que emprestar seu nome a pessoas, firmas, organizações ou empresas executoras de obras e serviços sem sua real participação nos trabalhos delas;

d) o profissional que, suspenso de seu exercício, continue em atividade;

e) a firma, organização ou sociedade que, na qualidade de pessoa jurídica, exercer atribuições reservadas aos profissionais da Engenharia, da Arquitetura e da Agronomia, com infringência do disposto no parágrafo único do Art. 8º desta Lei.

Portanto, quem faz previsão do tempo e do clima sem registro no CREA exerce ilegalmente a Profissão de Meteorologista!

Qual a função do Conselho Regional de Engenharia e Agronomia (CREA)?

A principal função do CREA é verificar, orientar e fiscalizar os exercícios profissionais com o objetivo de defender a sociedade das práticas ilegais dos ofícios que são abrangidos pelo sistema CONFEA/CREA. Além de promover a valorização profissional e garantir a primazia dos exercícios das atividades profissionais.

Portanto, de maneira simplificada o CONFEA/CREA protege a sociedade das práticas ilegais dos profissionais que são filiados aos conselhos e não o contrário.

Mas e o Exercício Ilegal da Profissão?

Bom, quando a obra não conta com responsável técnico, ou quando o "responsável técnico" identificado é um leigo, o CREA constata a irregularidade, procede à lavratura da Notificação e, quando necessário, do Auto de Infração. Quando este Auto não é atendido dentro das exigências da Lei, o CREA gera um processo administrativo, conforme a tipificação pertinente à atividade e/ou irregularidade encontrada.

Exercer Ilegalmente uma Profissão não é crime? Ninguém vai preso?

Segundo o artigo 282 do Código Penal da legislação brasileira somente médicos, dentistas e farmacêuticos cometem um crime ao exercer ilegalmente a profissão. Para todas as outras profissões, o exercício ilegal é apenas uma contravenção penal. E segundo o artigo 47 do decreto-lei 3.688 de 1941 que trata de contravenções a pena para esta prática é de prisão simples, de 15 dias a três meses, ou multa. Porém, existe um agravante. Se você provocar um dano material ou físico existe neste caso outro delito.

Portanto, o exercício ilegal da profissão é uma contravenção penal e precisa ser denunciada ao CREA e os conselhos regionais precisam colocar em prática uma das suas obrigações que é a de fiscalizar e autuar as pessoas que exercem estas práticas.

A comunidade Meteorológica pode e deve fortalecer e utilizar a Sociedade Brasileira de Meteorologia (SBMET) para promover debates/reuniões que possibilitem aumentar o conhecimento dos profissionais de Meteorologia sobre as legislações que regem a profissão do Meteorologista para que desta forma todos possam cobrar os conselhos regionais que fiscalizem e principalmente que analisem as denuncias feita sobre o tema. Ressaltando que a SBMET não é uma instituição com poder de Polícia ou um orgão regulador, mas uma SBMET bem representada e participativa dará coro as opiniões da comunidade Meteorológica e assim poderá cobrar aos CREAs uma fiscalização e apuração dos casos de Exercício Ilegal da Profissão de Meteorologista.

Este artigo não encerra a discussão é apenas o inicio e eu gostaria de saber oque você acha, por isso convido-o a comentar.

Boa Sorte e Bom Trabalho a Todos!

Interpolando dados do CMIP5 com CDO

2016-03-03T09:45:00.000-03:00

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!!

GrADS Preenchimento

2015-08-24T13:50:00.000-03:00

Olá a todos,

Para usuário do GrADS-2.1 é possível plotar um gráfico em que um determinado intervalo de dados pode ser destacado (hachurado) e isto se tornou possível graças a biblioteca gráfica chamada de Cairo.

A utilização desta biblioteca está bem documentada em alguns sites na internet, mas a pedidos segue um dica simples de como utilizar o comando 'set tile' para destacar um determinado intervalo dos dados. Isto pode ser útil para quem deseja destacar os dados estatisticamente significativo.

A ideia do comando "tile" é desenhar pequenas imagens para cobrir uma determinada área. O GrADS na versão 2.1 suporta a utilização do comando "tile" automaticamente, versões anteriores não possuem suporte.

Segue a descrição do comando:

'set tile nu tipo largura altura llargura fcor bcor'

onde:

nu é o número do tile que pode variar de 0 até 2000;

tipo é o tipo do tile que você quer utilizar; pode variar de 0 a 8 e ao escolher o 0 significa que você irá utilizar um arquivo com os tile (externo);

1 - sólido ;
2 - pontos ;
3 - linhas diagonais esquerda \\\\\\ ;
4 - linhas diagonais direita ////// ;
5 - linhas diagonais esquerda sobrepostas diagonais direita;
6 - linhas verticais ;
7 - linhas horizontais ;
8 - linhas horizontais sobrepostas as linhas verticais.

largura é a largura do "tile" gerada em pixels. O padrão é 9.
altura é a altura do "tile" gerada em pixels. O padrão é 9.
llargura é a espessura da linha utilizada para desenhar a linha especificada na imagem.
fcor é a cor usada para desenhar a linha
bcor é a cor do fundo utilizada para encher o fundo da imagem.

Para os "tile" gerados, uma pequena imagem é criada usando a largura e altura especificada, e, em seguida, ou é preenchido com uma cor sólida, ou é preenchido com a cor de fundo e, em seguida, uma ou duas linhas são desenhadas na imagem usando a cor do primeiro plano e a espessura da linha especificada. Note que as linhas são desenhadas de canto a canto (para linhas diagonais) ou no centro da imagem (diagonal / ou diagonal \ ). Assim, o tamanho da imagem, controlado pelos comando de largura e altura, vai controlar o espaçamento das linhas para o padrão final, e a relação de aspecto da imagem vai controlar o ângulo das linhas. Uma imagem pequena resultará em linhas espaçadas de forma mais estreita no padrão final, e uma imagem maior irá resultar num padrão mais espaçados.

Note que a espessura padrão é 3, a cor de primeiro plano padrão é o número de cores 1, e a cor de fundo padrão é um preto totalmente transparente (vermelho=0, verde=0, azul=0, alfa=0).

O "tile" será descrito como uma combinação do comando 'set tile' e 'set rgb' para definir as cores do primeiro plano e do fundo. Eles podem ser definidos em qualquer ordem, mas deve ser definido antes de qualquer desenho (display) feito na página atual.

Para utilizar um padrão de cores (rgb) ao "tile" que foi criado o comando é o seguinte:

'set rgb cc tile nu'

onde:

cc é o número do RGB criado;
nu é o número do tile.

Ao utilizar um padrão de cores a um "tile" criado não poderá ser utilizado a outro "tile".

Exemplo de utilização:

'reinit'
'sdfopen uwnd.mon.mean.nc'
'set display color white';'c'
'set grads off';'set grid off'

i = 60
n = 60
cmd = 'set rbcols'
while (i<256)
'set rgb 'n' 'i' 'i' 'i
cmd = cmd%' '%n
i = i + 15
n = n + 1
endwhile
cmd

* definindo duas cores
'set rgb 20 255 0 0'
'set rgb 21 0 0 255'

* definindo dois "tile"
'set tile 0 2 6 6 3 20'
'set tile 1 2 6 6 3 21'
'set rgb 22 tile 0'
'set rgb 23 tile 1'

'set gxout shaded'
'd uwnd'

* Desenhando pontilhados azuis e vermelhos para valores de vento superiores a 5 e inferiores a 5.
'set clevs 5'
'set ccols -1 22'
'd uwnd'
'set clevs -5'
'set ccols 23 -1'
'd uwnd'

Exemplo e:

'reinit'
'sdfopen uwnd.mon.mean.nc'
'set display color white';'c'
'set grads off';'set grid off'
'set tile 20 8 5 5 2 0'
'set rgb 30 tile 20'
'set gxout shaded'
'd uwnd'
'cbar'
'set clevs 5'
'set ccols -1 30'
'd uwnd'
'set clevs -5'
'set ccols 30 -1'
'd uwnd'

Abaixo segue um exemplo que mostra uma série de padrões diferentes. As especificações do comando 'tile' que foi usado para criar o padrão é mostrado em cada quadrado (sem os valores de cor). As definições de cor não são mostradas. Os comandos usados para desenhar o último quadrado são mostrados abaixo. Os comandos utilizados para os outros quadrados são semelhantes.

'set rgb 36 255 255 0'
'set rgb 40 255 0 255'
'set tile 21 2 9 9 3 40 36 '
'set rgb 31 tile 21'
'set line 31'
'draw recf 7 1 9 3 '
'draw string 8 0.8 set tile 21 2 9 9 3'

Boa Sorte e Bom Trabalho a Todos!

O ifort voltou a ter versão Grátis/Free

2015-07-20T14:00:00.000-03:00

A intel voltou a disponibilizar uma versão gratuita do seu compilador Fortran ifort. Depois de um ano sem fornecer uma versão gratuita a intel decidiu por voltar a disponibilizar uma versão free para Pesquisador, Estudantes, Professores e Desenvolvedores de Código Aberto.

Algumas mudanças aconteceram como por exemplo:

Não haverá mais uma versão para sistemas 32 bits.
Os Sistemas de derivação Linux suportados são:

Debian* 6.0, 7.0

Fedora* 20, 21

Red Hat Enterprise Linux* 5.0, 6.0, 7.0

SUSE Linux Enterprise Server* 11, 12

Ubuntu* 12.04 LTS, 14.04 LTS

Minimo de 2Gb de memória RAM
100Gb de espaço em HD livre, durante a instalação o software chega a ocupar 12Gb.
...

Para baixar o ifort é só seguir este link. Em breve pretendo criar um tutorial mostrando como instalar e as possíveis novidades do software.

Boa Sorte e Bom Trabalho a Todos!!

Plot TSM, Gelo Marinho e Cobertura de Neve - GrADS

2015-06-22T12:00:00.000-03:00

Você dificilmente vai fazer uma determinada figura se você não precisa dela. Recentemente eu precisei plotar a cobertura de Gelo Marinho a partir de dados de TSM. Vasculhando a internet encontrei um exemplo de um plot bem legal:

Os dados utilizados neste exemplo são do GFS e para ler os dados mais recentes eu irei utilizar o exemplo mostrado aqui no blog no post anterior.

Para que este script funcione perfeitamente é necessário ter as seguinte funções do Kodama em sua maquina: color.gs, xcbar.gs e colormap.gs.

Os primeiros comandos irão pegar a informação de data no sistema e abrir os dados de TSM e Gelo do GFS com resolução 0.25.

'reinit'
'! date +%Y%m%d > date.txt'
date=read('date.txt')
vardate=sublin(date,2)

arqgfs='http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs'vardate'/gfs_0p25_00z'
arqgelo='http://nomads.ncep.noaa.gov:9090/dods/ice/ice'vardate'/ice.00z'

'sdfopen 'arqgfs
'sdfopen 'arqgelo

'!rm date.txt'

Com os dados lidos iremos "setar" algumas funções para plotar os dados utilizando a projeção robinson.

'set gxout shaded'

'set mpdset hires'

'set lon -180 180'

'set lat -90 90'

'set grads off'

'set map 0 1 6'

'set mproj robinson'

'colormaps -l 272 307 0.5 -map jet'

'd tmpsfc'

'xcbar -fs 4'

Após estes comandos o plot exibido é +/- assim:

Agora iremos criar uma mascara para os continentes:

'set rgb 73 80 80 80'

'basemap L 73 0 M'

O que geramos a partir destes comandos é algo assim:

Em seguida iremos plotar o Gelo Marinho, para isto utilizamos os comandos:

'set gxout shaded'

'set dfile 2'

'set z 1'

'set t 1'

'set map 0 1 6'

'color 0 0.6 0.1 -kind dimgray->seashell->white'

'd maskout(icecmsl,icecmsl-0.05)'

O que iremos obter é o seguinte plot:

Por último iremos plotar a Cobertura de Neve e os últimos comandos do script são os seguintes:

'set dfile 1'
'color 5 100 5 -kind dimgray->seashell->white'
'set z 1'
'set t 1'
'set map 0 1 6'
'd maskout(weasdsfc,weasdsfc-5)'

E a última figura exibida tem o seguinte formato:

Espero que seja útil para vocês.

Boa Sorte e Bom trabalho a Todos!

Obtendo a saída do comando Shel - GrADS

2015-06-15T12:00:00.000-03:00

Olá a Todos,

Algumas pessoas utilizam o GrADS para ler os dados do GFS direto no GrADS Data Server. Em alguns casos pode ser necessário executar um script operacionalmente para efetuar esta tarefa. Um problema que pode existir para ler este tipo de dado (GrADS Data Server) operacionalmente é que por ser um dado de modelo os nomes dos arquivos mudam constantemente e torna-se necessário criar uma variável data para atualizar o nome do dado e efetuar a leitura corretamente.

Por exemplo:

Se você deseja ler o dado do GFS 0.25 deg starting from 00Z07jun2015 a url é a seguinte:

http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs20150607/gfs_0p25_00z

Observe que um dos diretórios onde está armazenado o dado contém a data (em negrito) em que o dado foi gerado ou o modelo foi rodado. Para gerar esta data automaticamente você pode criar um script SHELL e utilizar o comando date e resolver este problema, mas você pode utilizar o GrADS e resolver este problema diretamente no script .gs.

Para isso você pode utilizar o comando date no próprio script .gs . O script fica assim:

'reinit'
'! date +%Y%m%d > date.txt'
date=read('date.txt')
vardate=sublin(date,2)

*Neste caso vardate=20150607
*Observe que o comando date produz uma data diferente pois busca a data do sistema

file='http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs'vardate'/gfs_0p25_18z'

'sdfopen 'file

'!rm date.txt'

Eu não conheço uma outra forma de transformar o comando SHELL em variável, a única forma que conheço é escrevendo a saída do comando SHELL em um arquivo e depois ler este arquivo em seguida criando uma variável. Após a leitura do conteúdo do arquivo date.txt eu sugiro remover ('!rm date.txt') este arquivo para não ocupar espaço.

PS: Os dados do GFS são gerados de 6 em 6 horas para criar uma variável para ler os 4 arquivos você pode gerar as variáveis e utilizar um loop (WHILE) para ler cada arquivo.

Boa Sorte e Bom Trabalho a Todos!!

Duas variáveis em um mesmo campo - GrADS

2015-06-08T12:00:00.000-03:00

Olá a todos,

É importante percebermos que ao criar uma figura que será inserida em um artigo para ser publicado precisamos dedicar um pouco de esforço para fazer uma figura com a maior qualidade possível. Muitas vezes seu artigo será mais lembrado pelas figuras do que pelo conteúdo (o conteúdo é mais importante, não é isto que esta em discussão).

PS. Dica para quem tem o GrADS versão 2.

Esta dica é bastante simples, mas pode ser bastante útil quando se deseja criar figuras com alta qualidade. Existem vários comandos no GrADS que te possibilitam gerar uma figura com alta qualidade. No entanto, quando plotamos duas variáveis (sobrepostas) no mesmo campo esquecemos que no GrADS a segunda variável apesar de sobreposta é uma segunda figura (com novos eixos, grades e etc) e por isso é importante inibir a plotagem dos labels dos eixos, as linhas que representam as grades e etc desta segunda variável. Ao inibir o plot do label dos eixos por exemplo você mantêm a qualidade da font das Latitudes e Longitudes por exemplo. Para que isso não ocorra use os comandos abaixo sempre antes do display da próxima variável exibida:

'd press'
'set xlab off'
'set ylab off'
'set grid off'
'set mpdraw off'; *(Este comando inibe o plot do contorno dos continentes)
'd temp'

*Com o comando 'd temp' apenas os dados de temperatura serão plotados em seu campo. Espero que ajude.

Boa Sorte e Bom Trabalho a Todos!!

Imagem de Satelite no Grads

2015-04-30T17:22:00.003-03:00

Pergunta do Arthur Lucas,

Olá Cristiano,

tenho visto em alguns artigos com campos meteorológicos sobrepostos a imagens de Satélites, é possível fazer isso no Grads ?

Obrigado,

É sim possível desde que o dado da imagem de satélite seja disponibilizado no formato que seja possível o Grads ler, por exemplo: Binário, Netcdf e etc.

Alguns Institutos disponibilizam algumas imagens em formato binário. Durante meu doutorado eu utilizei dados de temperatura de brilho e consegui plotar algumas figuras de TB com vento e etc. Segue um exemplo: Eu plotei a temperatura de brilho com o shaded do Grads mas e com um RGB em tons de cinza que fica "quase igual" a uma imagem de satélite.

Abraço e Bom Trabalho a Todos!!

Gerando mascara com o cdo

2015-04-03T10:00:00.000-03:00

Ontem eu mostrei aqui como criar uma mascara (maskout) no Grads para plotar vetor vento significativo.

O Guilherme Martins sempre que possível contribui com alguma sugestão ou script. Desta vez ele compartilha conosco uma forma de resolver este mesmo problema sendo que com o cdo (Climate Data Operator).

PS: Não deixem de olhar a Apostila do próprio Guilherme sobre cdo.

Sem mais delongas, segue:

Dado para fazer este exemplo:

u -> uwnd.nc

v -> vwnd.nc

1) Juntando os arquivos de u e v para calcular a velocidade do vento

cdo merge uwnd.nc vwnd.nc uv.nc

2) Calculando a velocidade do vento

cdo expr,'vel=sqrt(uwnd*uwnd+vwnd*vwnd);' uv.nc vel.nc

3) Criando a máscara. Defino o intervalo da velocidade entre 0 e 4 m/s.

cdo setrtomiss,0,4 vel.nc mascara.nc

4) Aplicando a máscara na componente zonal (u) e meridional (v)

cdo ifthen mascara.nc uwnd.nc vento.u.nc

cdo ifthen mascara.nc vwnd.nc vento.v.nc

Ao plotar o que obtemos é o seguinte:

Boa Sorte e Bom Trabalho a Todos!

Grads Plot do Vento Significante

2015-04-02T17:09:00.000-03:00

Recentemente precisei plotar somente os vetores de vento no Grads significativos, ou seja, dentro de uma área/região eu necessitava que apenas um determinado limiar de magnitude do vetor vento fosse exibido.

Uma solução para este problema é utilizar uma máscara (maskout) para "mascarar" um determinado limiar de magnitude que você não quer exibir. Portanto, o script ou parte dele fica mais ou menos assim:

'open /usr/local/opengrads/Resources/SampleDatasets/model.ctl'

'set display color white'
'c'
'set lev 850'
'set t 1'
'set lon -100 50'
'set lat 0 50'

'wind = mag(ua,va)'

** Escolha do limiar que se deseja "mascarar"

limiar_vento = 8

'sigu = maskout(ua,wind-'limiar_vento')'

'sigv = maskout(va,wind-'limiar_vento')'

'sigwind = mag(sigu,sigv)'

'set gxout shaded'

'd sigwind'

'cbar'

'set gxout vector'

'd sigu;sigv'

Este seria o resultado.

Boa Sorte e Bom Trabalho a Todos!

Alterando os níveis verticais com CDO

2015-02-18T15:02:00.000-02:00

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!

cdo merge levels

2014-10-28T13:15:00.002-02:00

Olá a Todos,

Recentemente estive tentando usar o comando merge do cdo para criar um arquivo com vários níveis. O problema era o seguinte:

Eu possuo 10 arquivos de vento em formato binário sendo cada arquivo contendo 100 tempos e um nível.

A primeira coisa que eu fiz foi criar arquivos netcdf a partir dos arquivos descritores .ctl através do comando cdo import_binary -f nc *.ctl. No entanto desta forma o cdo não definia o nível referente ao arquivo e sempre atribuía nível 0 (zero). Com os arquivos criados com a descrição de nível 0 (zero) ao juntar (merge) os 10 arquivos o resultado era sempre um arquivo com somente um nível.

Eu não consegui resolver este problema de maneira mais simples o que fiz foi o seguinte:

Utilizei o comando sdfwrite do GrADS para gerar os arquivos netcdf que precisava.

Ex:

set x 1 144
set y 1 73
set t 1 100
set z 1
define uwnd = uwnd
set sdfwrite -4d arquivo1.nc
sdfwrite uwnd
quit

OBS1: É necessário inserir a função -4d do comando sdfwrite, desta forma você força que o arquivo de saída tenha 4 dimensões (lat, lon, tempo, nível)

OBS2: O nível no arquivo ctl precisa ser definido de maneira correta, por exemplo:

zdef 1 linear 200 1

Qualquer definição diferente disso deu errado.

Agora que possuía os 10 arquivos .nc utilizei o merge para gerar meu arquivo com os 10 níveis e 100 tempos:

cdo merge arquivo1.nc arquivo2.nc ... arquivo10.nc saida.nc

Boa Sorte e Bom Trabalho a Todos!!

Dado Texto ou ASCII no GrADS

2014-10-07T21:30:00.000-03:00

Olá Pessoal,

Através do GrADS é possível escrever e ler arquivos textos e em algum momento isto pode ser de grande utilidade. Recentemente a Giulia de Salve me escreveu com uma dúvida exatamente sobre este contexto.

O Guilherme Martins em sua Apostila já colocou um exemplo bem explicativo e também existe um guia no site do GrADS que é mais explicativo ainda. O que eu vou colocar aqui é exatamente o que esta na página neste guia do site GrADS.

PS. Eu irei utilizar como exemplo o arquivo model.ctl para quem possui o GrADS 2 (opengrads) instalado ele fica na pasta (grads-2.0.1/Contents/Resources/SampleDatasets/).

O GrADS não possui uma função intrínseca para ler arquivos texto como o sdfopen, por exemplo, por isso é necessário alguns ajustes para isso:

Imagine que você possua um arquivo texto/ascii (my_ascii_file.txt) que tenha uma coluna com valores de 1 a 30:

...

Você pode transformar este arquivo em um binário e criar um ctl para ler no GrADS.

'open model.ctl'

file = 'my_ascii_file.txt'

datafile = 'my_binary_file.dat'
'set gxout fwrite'
'set fwrite -ap 'datafile
'!/bin/rm -f 'datafile
while (1)
         res = read(file)
         line1 = sublin(res,1)
         line2 = sublin(res,2)
         rc1 = subwrd(line1,1)
         if (rc1); break; endif
              val = subwrd(line2,1)
              'd 'val
endwhile
rc = close(file)
'disable fwrite'

O ctl para ler o dado gerado (my_binary_file.dat), imaginando que você queira plotar os dados com relação ao eixo X, ficaria assim:

dset ^my_binary_file.dat
title Sample of ASCII data converted to binary
undef -9.99e8
xdef 100 linear 1 1
ydef 1 linear 1 1
zdef 1 linear 1 1
tdef 1 linear 01jan0001 1dy
vars 1
a 0 99 ascii variable
endvars

Se você abrir o ctl acima e plotar vai gerar uma figura semelhante a esta:

Uma outra opção é criar arquivos textos a partir de qualquer dado binário através da função 'set gxout print'. A formatação desta função pode ser feita através da função 'set prnopts'. A saída é plotada no terminal mas pode ser armazenada em uma variável e posteriormente através do script ser escrita em um arquivo. O exemplo a seguir escreve 72 números por linha que representa os 72 pontos de grade em X.

'open model.ctl'

'set x 1 72'

outfile = my_ascii_file.txt'

'set gxout print'

'set prnopts %7.3f 72 1'

'd ts'

rc = write(outfile,result)

Desta forma o arquivo my_ascii_file.txt irá conter algo como isso:

Printing Grid -- 3312 Values -- Undef = -9.99e+08

258.493 258.493 258.493 258.493 258.493 258.493 258.493 258.493 258.493 258.493 258.493 <etc.>

No caso anterior, observe que a primeira linha do arquivo de saída contem algumas informações sobre o número de valores escritos e valores indefinidos que você pode não querer que seja escrito. O exemplo a seguir omiti essa linha, ou seja, pula a primeira linha de informações de diagnóstico. Além disso neste exemplo a saída será delimitada por vírgulas para importar para uma planilha (excell) por fim a variável apresentada agora é uma expressão em vez de um nome.

'open model.ctl'

'set x 1 72'

outfile='my_ascii_file.txt'

'!/bin/rm -f 'outfile

'set gxout print'

'set prnopts %g, 72 0'

'd ts-273.15'

i=1

while (1)

line = sublin(result,i)

if (line = ''); break; endif

if (i>1)

rc = write(outfile,line,append)

endif

i=i+1

endwhile

Neste caso o arquivo my_ascii_file.txt conterá algo como:

-14.6567,-14.6567,-14.6567,-14.6567,-14.6567,-14.6567,-14.6567,-14.6567,-14.6567,<etc.>

Obs. A função append serve para inserir os dados no final do arquivo caso ele exista. Se você não utilizar a função append e existir um arquivo de mesmo nome ele será sobrescrito. Para mais informações leia "append".

Por fim é possível escrever informações sobre a grade (lat x lon) além dos dados de mais de uma variável, para isso é necessário salvar o resultado da exibição de cada variável/expressão (incluindo longitude e latitude), para cada ponto de grade individualmente, organizar a estrutura que será escrito no arquivo texto de saída e depois escrevê-lo. Neste exemplo será escrito a longitude, latitude, temperatura da superfície, pressão de superfície e precipitação:

'open model.ctl'

'set x 1 72'

outfile='my_ascii_file.txt'

'!/bin/rm -f 'outfile

'set gxout print'

fmt='%8.3f'

numcols=72

'set prnopts 'fmt' 'numcols' 1'

'd lon'

lon_data = result

'd lat'

lat_data = result

'd ts'

v1_data = result

'd ps'

v2_data = result

'd p'

v3_data = result

i=1

while (1)

lons = sublin(lon_data,i)

lats = sublin(lat_data,i)

line1 = sublin(v1_data,i)

line2 = sublin(v2_data,i)

line3 = sublin(v3_data,i)

if (lons='' | lats='' | line1='' | line2='' | line3=''); break; endif

if (i>1)

j=1

while (j<=numcols)

str = subwrd(lons,j); lon = math_format(fmt,str)

str = subwrd(lats,j); lat = math_format(fmt,str)

str = subwrd(line1,j); v1 = math_format(fmt,str)

str = subwrd(line2,j); v2 = math_format(fmt,str)

str = subwrd(line3,j); v3 = math_format(fmt,str)

record = lon' 'lat' 'v1' 'v2' 'v3

rc = write(outfile,record,append)

j=j+1

endwhile

endif

i=i+1

endwhile

Neste caso a saída my_ascii_file.txt apresentará a seguinte estrutura:

0.000 -90.000 258.493 669.911 0.000

5.000 -90.000 258.493 669.911 0.000

10.000 -90.000 258.493 669.911 0.000

15.000 -90.000 258.493 669.911 0.000

20.000 -90.000 258.493 669.911 0.000

25.000 -90.000 258.493 669.911 0.000

30.000 -90.000 258.493 669.911 0.000

35.000 -90.000 258.493 669.911 0.000

40.000 -90.000 258.493 669.911 0.000

A um tempo atrás coloquei a disposição uma função que ajuda a escrever um arquivo texto sem muitas complicações. No entanto o arquivo de saída não possui formatação, veja aqui.

Boa Sorte e Bom Trabalho a Todos!

Cristiano

Terminator: Dividindo seu terminal em várias janelas

2014-10-02T17:32:00.000-03:00

Segue uma excelente sugestão do Mercel José e o Thomas Martin.

Uma "dificuldade" encontrada por algumas pessoas que trabalha com Linux é a necessidade de abrir vários terminais. O terminal atual do Linux possibilita que você abra várias abas e desta forma você precisa migrar para cada aba sempre que precisa efetuar uma atividade e etc. Com um terminal chamado de Terminator você consegue dividir seu terminal em vários terminais e isto pode ser de grande ajuda.

O Terminator é um software leve não ocupa muita memoria, no entanto em telas de computadores pequenas deve ficar um pouco difícil de enxergar. O grande diferencial é ter uma janela dividida em vários terminais.

Para instalar é fácil:

Ubuntu ou Debian:

>sudo apt-get install terminator

Fedora:

>sudo yum install terminator

Para utilizar também é simples:

Ctrl+Shift+X -  Alterna entre mostrar todos os terminais ou apenas um;

Ctrl+Shift+O - Divide a janela em terminais horizontais;

Ctrl+Shift+E - Divide a janela em terminais verticais;

Ctrl+Shift+Right - Redimensiona o terminal para direita;
Ctrl+Shift+Left - Redimensiona o terminal para esquerda;
Ctrl+Shift+Up - Redimensiona o terminal para cima;
Ctrl+Shift+Down - Redimensiona o terminal para baixo;
Ctrl+Shift+W - Fecha o terminal que esta sendo utilizado;
Ctrl+Shift+T -   Abre uma nova aba;
Ctrl+Plus (+) -  Aumenta o tamanho da fonte;
Ctrl+Minus (-) - Diminui o tamanho da fonte;
Ctrl+Zero (0) - Restaura o tamanho original da fonte;


Boa Sorte e Bom Trabalho a Todos!!

Cristiano

O ifort agora é Pago!

2014-09-24T18:43:00.001-03:00

Quem aproveitou, aproveitou!

Recentemente descobrimos que o Compilador fortran e C++ da intel ifort e icc agora são pagos. Seguramente o compilador fortran mais rápido e eficiente já lançado. Diz a lenda que o desenvolvimento do ifort evoluiu junto ao desenvolvimento do Modelo Regional WRF. A lenda conta que ao desenvolver o ifort a intel necessitava de executar um processamento ao qual exigisse o máximo de performance e eles conseguiram ao executar o modelo WRF.

No entanto o desenvolvimento do ifort foi aos poucos migrando para programação orientada ao objeto e muitas funções de programação em C++ foram adicionadas. A utilização do sistema via terminal (modo texto) foi sendo abandonada aos poucos e desde 2011 a intel vinha desenvolvendo o intell Studio XE (modo gráfico) que combinava em um mesmo software o ifort e o icc. O inevitável acabou acontecendo o investimento no software foi aumentando cada vez mais o software foi ficando cada vez mais robusto, consequentemente o intel studio foi adquirindo cada vez mais adepto. A consequência disso tudo é que a partir de 2014 os compiladores ifort e icc se tornarão pagos.

O software pode ser encontrado em 3 níveis de desenvolvimento Cluster Edition, Professional Edition e Composer Edition. O valor varia de acordo com o tipo de comprador (Individual, empresa e etc, Estudante, profissional e etc). A licença individual mais barata e mais simples para fins acadêmico custa 249 dólares a mais cara e mais completa e para fins comerciais 1699 dólares. É possível baixar uma versão trial com duração de 30 dias.

É uma pena, mas já era esperado. Portanto, se você possui uma versão free do ifort e do icc na sua maquina zele por ela :)

Segue o link da página para compra do ifort.

Boa Sorte e Bom Trabalho a Todos!

Cristiano

Barra de cores "pulando" intervalos

2014-09-11T13:13:00.002-03:00

Está é mais uma sugestão de barra de cores (shaded) para usuários Grads. O interessante desta barra é que é possível "saltar" os valores (números) referentes aos contornos da variável.

O autor deste script chama-se Stephen R. McMillan

*
* Script to plot a colorbar
*
* The script will assume a colorbar is wanted even if there is
* not room -- it will plot on the side or the bottom if there is
* room in either place, otherwise it will plot along the bottom and
* overlay labels there if any. This can be dealt with via
* the 'set parea' command. In version 2 the default parea will
* be changed, but we want to guarantee upward compatibility in
* sub-releases.
*
* modifications by mike fiorino 940614
* rvp added by diane stokes 20040729
* cinc added by diane stokes 20050130
*
* - the extreme colors are plotted as triangles
* - the colors are boxed in white
* - input arguments in during a run execution:
*
* run cbarn cinc sf vert xmid ymid rvp
*
* cinc - increment for labels
* sf - scale the whole bar 1.0 = original 0.5 half the size, etc.
* vert - 0 FORCES a horizontal bar = 1 a vertical bar
* xmid - the x position on the virtual page the center the bar
* ymid - the x position on the virtual page the center the bar
* rvp - 1 resets virtual page (one colorbar for several plots)
*
* if cinc is not specified, every color interval is labeled
* if vert,xmid,ymid are not specified, they are selected
* as in the original algorithm
*
* if rvp is not specified, it is 0 (do not reset virtual page)

function colorbar (args)

* Parse command line arguments
* ----------------------------
cinc=subwrd(args,1)
sf=subwrd(args,2)
vert=subwrd(args,3)
xmid=subwrd(args,4)
ymid=subwrd(args,5)
rvp=subwrd(args,6)
if(cinc=''|cinc=0);cinc=1;endif
if(sf='');sf=1.0;endif

* Check shading information
* -------------------------
'query shades'
shdinfo = result
if (subwrd(shdinfo,1)='None')
say 'Cannot plot color bar: No shading information'
return
endif

* Check if reset of virtual page is desired
* -----------------------------------------
if(rvp = 1)
say 'Colorbar based on colors from last plot drawn'
'set vpage off'
endif

* Get plot size info
* ------------------
'query gxinfo'
rec2 = sublin(result,2)
rec3 = sublin(result,3)
rec4 = sublin(result,4)
xsiz = subwrd(rec2,4)
ysiz = subwrd(rec2,6)
ylo = subwrd(rec4,4)
xhi = subwrd(rec3,6)
xd = xsiz - xhi

ylolim=0.6*sf
xdlim1=1.0*sf
xdlim2=1.5*sf
barsf=0.8*sf
yoffset=0.2*sf
stroff=0.05*sf
strxsiz=0.12*sf
strysiz=0.13*sf

* Decide if horizontal or vertical color bar
* and set up constants.
* ------------------------------------------
if (ylo<ylolim & xd<xdlim1)
say "Not enough room in plot for a colorbar"
return
endif
cnum = subwrd(shdinfo,5)

* Logic for setting the bar orientation with user overides
* --------------------------------------------------------
if (ylo<ylolim | xd>xdlim1)
vchk = 1
if(vert = 0) ; vchk = 0 ; endif
else
vchk = 0
if(vert = 1) ; vchk = 1 ; endif
endif

* Vertical bar
* ------------
if (vchk = 1 )
if(xmid = '') ; xmid = xhi+xd/2 ; endif
xwid = 0.2*sf
ywid = 0.5*sf
xl = xmid-xwid/2
xr = xl + xwid
if (ywid*cnum > ysiz*barsf)
ywid = ysiz*barsf/cnum
endif
if(ymid = '') ; ymid = ysiz/2 ; endif
yb = ymid - ywid*cnum/2
'set string 1 l 5'
vert = 1

* Horizontal bar
* --------------
else

ywid = 0.4
xwid = 0.8
if(ymid = '') ; ymid = ylo/2-ywid/2 ; endif
yt = ymid + yoffset
yb = ymid
if(xmid = '') ; xmid = xsiz/2 ; endif
if (xwid*cnum > xsiz*barsf)
xwid = xsiz*barsf/cnum
endif
xl = xmid - xwid*cnum/2
'set string 1 tc 5'
vert = 0
endif

* Plot colorbar
* -------------
'set strsiz 'strxsiz' 'strysiz
num = 0
while (num<cnum)
rec = sublin(shdinfo,num+2)
col = subwrd(rec,1)
hi = subwrd(rec,3)
if (vert)
yt = yb + ywid
else
xr = xl + xwid
endif

* Draw the left/bottom triangle
* -----------------------------
if (num = 0)
if(vert = 1)
xm = (xl+xr)*0.5
'set line 'col
'draw polyf 'xl' 'yt' 'xm' 'yb' 'xr' 'yt' 'xl' 'yt
'set line 1 1 5'
'draw line 'xl' 'yt' 'xm' 'yb
'draw line 'xm' 'yb' 'xr' 'yt
'draw line 'xr' 'yt' 'xl' 'yt
else
ym = (yb+yt)*0.5
'set line 'col
'draw polyf 'xl' 'ym' 'xr' 'yb' 'xr' 'yt' 'xl' 'ym
'set line 1 1 5'
'draw line 'xl' 'ym' 'xr' 'yb
'draw line 'xr' 'yb' 'xr' 'yt
'draw line 'xr' 'yt' 'xl' 'ym
endif
endif

* Draw the middle boxes
* ---------------------
if (num!=0 & num!= cnum-1)
'set line 'col
'draw recf 'xl' 'yb' 'xr' 'yt
'set line 1 1 5'
'draw rec 'xl' 'yb' 'xr' 'yt
endif

* Draw the right/top triangle
* ---------------------------
if (num = cnum-1)
if (vert = 1)
'set line 'col
'draw polyf 'xl' 'yb' 'xm' 'yt' 'xr' 'yb' 'xl' 'yb
'set line 1 1 5'
'draw line 'xl' 'yb' 'xm' 'yt
'draw line 'xm' 'yt' 'xr' 'yb
'draw line 'xr' 'yb' 'xl' 'yb
else
'set line 'col
'draw polyf 'xr' 'ym' 'xl' 'yb' 'xl' 'yt' 'xr' 'ym
'set line 1 1 5'
'draw line 'xr' 'ym' 'xl' 'yb
'draw line 'xl' 'yb' 'xl' 'yt
'draw line 'xl' 'yt' 'xr' 'ym
endif
endif

* Put numbers under each segment of the color key
* but only if this level is to be printed (based on cinc)
* -------------------------------------------------------
if(num < cnum-1)
* base prt on hi if increment is based on data value (less reliable)
* prt = math_fmod(hi,cinc);
* ------------------------------------------------------------------
prt = math_fmod(num,cinc);
if(prt = 0)
if (vert)
xp=xr+stroff
'draw string 'xp' 'yt' 'hi
else
yp=yb-stroff
'draw string 'xr' 'yp' 'hi
endif
endif
endif

* Reset variables for next loop execution
* ---------------------------------------
if (vert)
yb = yt
else
xl = xr
endif
num = num + 1

endwhile
return

Para utilizar é fácil:

Após abrir o arquivo desejado, definir a função gxout shaded e exibir (display) a variável, basta digitar:

ga-> run cbarnskip.gs 2

Eu chamei de cbarnskip.gs você escolhe o nome que achar melhor. O número 2 está relacionado aos valores que serão "removidos" ou alternados.

O resultado pode ser algo do tipo:

Boa Sorte e Bom Trabalho a Todos!

GxYAT um comando para gerar Figuras com alta resolução

2014-02-06T14:57:00.000-02:00

A pedido de Luana Pampuch segue um exemplo para salvar no GrADS figuras em .pdf ou .eps e etc.

No OpenGrADS existe um comando chamado GxYAT que você pode criar figuras em PDF, PNG, Postscript e SVG.

Para quem escreve textos em Latex (meu caso) é excelente pois as figuras ficam com excelente resolução e durante a impressão as figuras não perdem qualidade.

Portanto, para quem esta escrevendo TCC, Dissertação ou Tese fica a dica, aprendam Latex salvem suas figuras em PDF e se livrem de serem cobrados pela má qualidade das figuras durante as defesas :).

como usar no GrADS:

'enable print nome_da_figura.gmf'
.
.
.
'print'
'disable print'

'quit'

No Terminal Linux digite:

> gxyat -o nome_da_figura.pdf nome_da_figura.gmf
É isso!!

Boa Sorte e Bom Trabalho a Todos

Script GrADS Barra em formato de Diamante

2014-02-04T19:39:00.002-02:00

Olha o GrADS surpreendendo novamente. Recentemente encontrei um script que plota uma barra de cores em formato de diamante e irei compartilhar com vocês. Acredito que fazer figuras para apresentações com este tipo de barra pode ser interessante, vamos lá:

Para usar as opções são as seguintes:

'run diamondbar.gs -r -marg 0.1 -font 0.09 -fs 2'

Este comando serve para plotar a figura acima seguindo o exemplo disponibilizado abaixo.

Opções:

-help = Ajuda;
-loc = Escolhe um local para a colorbar (r=direita-vertical, l=esquerda-vertical,
b=base-horizontal, t=topo-horizontal),
você também pode chamar o comando através de (-l, -r, -b, -t);
-size/s = Escolha o tamanho dos diamantes;
-marg = Escolha a margem para o eixo;
-fs/int = Escolha o intervalo para cada número plotado ao longo da colorbar;
-font/fontsize = Escolha o tamanho da fonte dos números;
-scol = Escolha a cor dos números;
-line = Desenha o contorno de cada diamante caso deseje.

Scripts utilizados:

barradiamante.gs
Exemplo.gs

Simples assim, espero que funcione.

Boa Sorte e Bom Trabalho a Todos!!

Instalação PyFerret com Python 2.6

2013-11-22T15:58:00.001-02:00

Mais uma contribuição do Paulo Henrique e Marcelo Rodrigues, encaminhada pelo Dr. Antonio Carlos Santana dos Santos - Prof. Adjunto Vice-Coordenador do Curso de Física da Universidade Estadual do Ceará.

Um Tutorial da instalação do PyFerret para Debian 7 com o Python 2.6.

PS. Segundo os autores esta instalação funcionará para qualquer distribuição Linux desde que tenha o Python 2.6 instalado.

Clique aqui.

Mai uma vez obrigado ao Paulo Henrique e Marcelo Rodrigues pelo desenvolvimento e ao Dr. Antonio Carlos Santana dos Santos pela divulgação e contribuição.

Boa Sorte e Bom Trabalho a Todos!!

Script Download Dados ECMWF/ERA-Interin - Parte 1

2013-11-11T13:43:00.004-02:00

Quando eu atualizei o link do "novo" servidor de download de dados da Reanalise do Era-Interim o Helber comentou de um script que possibilitava o download. Segue como fazer:

PS. Para entender os parâmetros e como modificá-los não esqueça de lê o segundo artigo (Script Download Dados ECMWF/ERA-Interim - Parte 2), link no final deste post.

O ECMWF está disponibilizando uma API (Aplication Programming Interface) nas linguagens (Python, Perl e Java) para acesso/download dos seus dados via script/terminal. É uma ferramenta extremamente útil para quem esta familiarizado com este tipo de sistema e não quer baixar os dados via site. Eu recomendo!

Eu irei tentar ensinar como usar a API desenvolvida em Python, escolha a linguagem que você melhor trabalha. Todas estas informações estão disponíveis no site da API no ECMWF.

Para instalar a biblioteca (ecmwfapi) você precisará primeiro instalar o programa pip.

O programa pip é uma excelente alternativa para instalar bibliotecas desenvolvidas em python de maneira fácil, através do terminal. Para instalar o pip no Ubuntu execute no terminal os seguintes comandos:

> sudo apt-get install python-pip python-dev build-essential
> sudo pip install --upgrade pip

Agora com o pip instalado execute no terminal: (tudo junto separado por espaço)

> sudo pip install https://software.ecmwf.int/wiki/download/attachments/23694554/ecmwf-api-client-python.tgz

Este comando irá instalar a biblioteca, você pode verificar ela instalada no diretório:

/usr/local/lib/python2.7/dist-packages/

Feito isso, agora você irá precisar criar a sua chave de acesso. Quando precisávamos baixar os dados da reanálise do ECMWF assim que escolhíamos os dados o site solicitava os dados pessoais (Nome a Instituição e o email), agora com essa API você gera uma chave, insere no script, e não precisará fazer esta autorização sempre que for baixar, agora ela feita automaticamente via API.

Para criar a chave você precisa se cadastrar no site:

http://apps.ecmwf.int/registration/

Um email de conformação será enviado com seu login e senha, confirme clicando no link. Ao fazer login na pagina abaixo será disponibilizado sua chave de acesso:

https://apps.ecmwf.int/auth/login/

Agora que sua chave já esta disponível é preciso criar um arquivo com o nome (.ecmwfapirc) no seu /home/usuario para o acesso ser automático, Obs: o nome do arquivo é .ecmwfapirc com (.) para ficar oculto e sem extensão.

Neste arquivo você irá inserir:

{
"url" : "https://api.ecmwf.int/v1",
"key" : "XXXXXXXXXXXXXXXXXXXXXX",
"email" : "john.smith@example.com"
}

,onde vocé irá substituir XXXXXXXXXXXXXXXXXXXXXX pela sua chave e john.smith@example.com pelo email que você utilizou para efetuar o cadastro.

Feito isso você já pode baixar os dados. O ecmwf disponibiliza um script exemplo para teste. Portanto, crie um script em PYTHON e dentro dele insira as linhas: eu chamei o meu script de exemplo.py

#!/usr/bin/env python
from ecmwfapi import ECMWFDataServer

server = ECMWFDataServer()

server.retrieve({
'dataset' : "tigge",
'step' : "24/to/120/by/24",
'number' : "all",
'levtype' : "sl",
'date' : "20071001/to/20071003",
'time' : "00/12",
'origin' : "all",
'type' : "pf",
'param' : "tp",
'area' : "70/-130/30/-60",
'grid' : "2/2",
'target' : "data.grib"
})

Salve e dê permissão de execução:

> chmod +x exemplo.py

Ao executar este script:

>./exemplo.py

Aparecerá uma série de informações no terminal sobre o sistema, dado e o andamento do Download.

Script Download Dados ECMWF/ERA-Interim - Parte 2

Boa Sorte e Bom Trabalho a Todos.