Voce dificilmente vai fazer uma determinada figura se voce nao 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 sao 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 e necessario ter as seguinte funcoes do Kodama em sua maquina: color.gs, xcbar.gs e colormap.gs.Os primeiros comandos irao pegar a informacao de data no sistema e abrir os dados de TSM e Gelo do GFS com resolucao 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 funcoes para plotar os dados utilizando a projecao 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'Apos estes comandos o plot exibido e +/- 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 e 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 e o seguinte plot:Por ultimo iremos plotar a Cobertura de Neve e os ultimos comandos do script sao 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 ultima figura exibida tem o seguinte formato:Espero que seja util para voces. Boa Sorte e Bom trabalho a Todos! Ola a todos,Para usuario do GrADS-2.1 e possivel plotar um grafico em que um determinado intervalo de dados pode ser destacado (hachurado) e isto se tornou possivel gracas a biblioteca grafica chamada de Cairo.A utilizacao desta biblioteca esta 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 util para quem deseja destacar os dados estatisticamente significativo.A ideia do comando "tile" e desenhar pequenas imagens para cobrir uma determinada area. O GrADS na versao 2.1 suporta a utilizacao do comando "tile" automaticamente, versoes anteriores nao possuem suporte.Segue a descricao do comando:'set tile nu tipo largura altura llargura fcor bcor'onde:nu e o numero do tile que pode variar de 0 ate 2000;tipo e o tipo do tile que voce quer utilizar; pode variar de 0 a 8 e ao escolher o 0 significa que voce ira utilizar um arquivo com os tile (externo);1 - solido ;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 e a largura do "tile" gerada em pixels. O padrao e 9.altura e a altura do "tile" gerada em pixels. O padrao e 9.llargura e a espessura da linha utilizada para desenhar