Transcript sinotica aula 13082013_JoseLeandro

Previsão e Estudo de Caso
Ferramentas Básicas
J.Leandro P.S.Campos
GFS 4 – GLOBAL FORECAST SYSTEM
GFS – Global Forecast System
• Download:
– http://nomads.ncdc.noaa.gov/data/gfs4/
– Previsões ate 192 horas após a rodada (8 dias)
• Dados no formato Grib2
–
–
–
–
Necessário ter o Perl instalado (Unix)
http://www.perl.org/
Necessário ter o Cygwin e Perl (Cygwin/Windows)
http://www.cygwin.com/
• OpenGrads, IDV, McIdas, Matlab, ... (Unix/Windows)
– http://opengrads.org/
– http://www.unidata.ucar.edu/software/idv/
GFS – Global Forecast System
Usuários do OpenGrads
• Usuários do OpenGrads
– Não é possível abrir os dados sem um arquivo
descritor ( .ctl)
– Necessário ter os seguintes scripts:
• g2ctl.pl
• wgrib2
• http://www.cpc.ncep.noaa.gov/products/wesley/g2ctl.html
– Após ter esses dois scripts no mesmo diretório em
que os arquivos ...
• Para criar um template (arquivo descritor) devemos digitar
os seguintes comandos no shel:
GFS – Global Forecast System
Usuários do OpenGrads
• Para previsões:
$ g2ctl.pl gfs_4_20121108_0000_%f3.grb2 >example.ctl (Unix)
$ perl g2ctl.pl gfs_4_20121108_0000_%f3.grb2 >example.ctl
$ gribmap -i example.ctl
ga -> open example.ctl
• Para análises:
$ g2ctl.pl -0 gfs_4_%y4%m2%d2_%h200_000.grb2 >example.ctl
$ perl g2ctl.pl - 0 gfs_4_%y4%m2%d2_%h2_000.grb2 >example.ctl
$ gribmap -0 -i example.ctl
ga -> open example.ctl
• mais detalhes em:
– http://www.cpc.ncep.noaa.gov/products/wesley/g2ctl.html
CONCEITOS MUITO BÁSICOS
SOBRE OPENGRADS
Básico - Grads
• Abrindo arquivos binários
ga-> open example.ctls
• Abrindo arquivos netcdf
ga-> sdfopen example.nc
• Alguns comandos básicos
ga-> q vars
ga-> q ctlinfo
ga-> d variável
ga-> set display color white
ga-> c
ga-> quit
( variáveis contidas no arquivo )
( informações contidas no ctl )
( plotando uma variável )
( tela branca )
( limpa todos os displays )
( sai do grads )
Básico - Grads
• Definindo latitudes e longitude e níveis específicoss
ga-> set lat -90 90
ga-> set lon -180 180
ga-> set lev 500
• Rodando um script
ga-> run script.gs
ga-> script.gs
• Salvando uma imagem
ga-> printim display.png (.jpeg, .gif, ... )
•
Para mais informações
– http://www.cptec.inpe.br/ManualGrADS/
– http://www.iges.org/grads/gadoc/reference_card.pdf
– http://www.iges.org/grads/gadoc/
Construindo um Script
Básico
•
•
•
Pode ser criado por um bloco de notas
Os comandos do grads devem estar entre aspas simples ( ‘d temp’) ou duplas (“d
temp”).
Exemplo:
‘reinit’
‘open example.ctl’
‘set display color white’
*Isto é um comentário
‘c’
‘set gxout shaded’
‘d temp’
‘cbarn’
‘printim temperatura.png’
‘close 1’
•
Salve o script como “temp.gs” e no grads digite:
ga-> run temp.gs
ga-> temp
Construindo um Script
Criando Loops
•
Veja o exemplo abaixo
‘reinit’
‘open example.ctl’
( ... )
tt = 1
tf = 12
* Aqui começa o loop
WHILE( tt <= tf )
‘set t ‘tt
say tt
‘d temp’
‘cbarn’
(...)
‘tt = tt + 1
ENDWHILE
*Aqui termina o loop
‘close 1’
CRIANDO OU LENDO UM ARQUIVO
BINÁRIO NO FORTRAN 95/2003
Lendo ou Escrevendo um Arquivo
Binário do Grads
•
Veja o Exemplo abaixo
PROGRAM teste
IMPLICIT NONE
( ... )
! Lendo um arquivo binário
OPEN(10,FILE=TRIM(iofile),ACTION='read',STATUS='old',ACCESS='direct',RECL=4*xdim*ydim )
DO k = 1, zdim
READ(10,REC=k) var1(:,:,k)
END DO
( ... )
! Escrevendo um arquivo binário
OPEN(20,FILE=TRIM(iofile),ACTION=‘write',STATUS=‘new',ACCESS='direct',RECL=4*xdim*ydim )
DO k = 1, zdim
WRITE(10,REC=k) var2(:,:,k)
END DO
( ... )
END PROGRAM teste
Alguns Links
• Compiladores Fortran
– http://www.silverfrost.com/default.aspx
– http://software.intel.com/en-us/intel-compilers
– http://www.cse.yorku.ca/~roumani/fortran/ftn.ht
m
• Compilador C/ C++
– http://www.codeblocks.org/
• Foruns
– http://forums.silverfrost.com/index.php