130 likes | 235 Views
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)
E N D
Previsão e Estudo de CasoFerramentas Básicas J.LeandroP.S.Campos
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 SystemUsuá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 SystemUsuários do OpenGrads • Para previsões: $ g2ctl.pl gfs_4_20121108_0000_%f3.grb2 >example.ctl (Unix) $ perlg2ctl.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
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 ( variáveis contidas no arquivo ) ga-> q ctlinfo ( informações contidas no ctl ) ga-> d variável ( plotando uma variável ) ga-> set display color white ( tela branca ) ga-> c ( limpa todos os displays ) ga-> quit ( 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 ScriptBá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 gxoutshaded’ ‘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 ScriptCriando Loops • Veja o exemplo abaixo ‘reinit’ ‘open example.ctl’ ( ... ) tt = 1 tf = 12 * Aqui começa o loop WHILE( tt <= tf ) ‘set t ‘tt saytt ‘d temp’ ‘cbarn’ (...) ‘tt = tt + 1 ENDWHILE *Aqui termina o loop ‘close 1’
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) DOk = 1, zdim WRITE(10,REC=k)var2(:,:,k) END DO ( ... ) END PROGRAMteste
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.htm • Compilador C/ C++ • http://www.codeblocks.org/ • Foruns • http://forums.silverfrost.com/index.php