Comandos do LINUX no Grads

duvanzan
(usa Linux Mint)

Enviado em 10/02/2019 - 11:10h

Bom dia,

Estou usando o grads para plotar alguns dados que preciso. Porém, preciso muito usar um comando no linux grads!

Testei o comando no terminal (fora do grads) e funcionou perfeitamente! Quando coloco ele dentro do grads, não funciona de jeito nenhum! Alguém poderia me ajudar?
No grads, os comandos que estão entre '' são comandos usados do próprio grads.


O comando que usei no linux é o seguinte:

pixel=1

cat dados_xy.txt |sed -n $pixel'p'


Quando estou no grads, uso esse comando dentro de um while, dentro do grads está assim:



##### Script para rodar as figuras
## Setando barra de cores
##light yellow to dark red
'set rgb 21 255 250 170'
'set rgb 22 255 232 120'
'set rgb 23 255 192 60'
'set rgb 24 255 160 0'
'set rgb 25 255 96 0'
'set rgb 26 255 50 0'
'set rgb 27 225 20 0'
'set rgb 28 192 0 0'
'set rgb 29 165 0 0'

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

'reinit'
'open saida_irma_v5_final.ctl'


'set lat 12 35'
'set lon 275 310'

##Looping para gerar todas as imagens
count=1
while (count<60)
#O comando abaixo pega os dados da linha. Nesse caso se fosse 1p, pegaria o comando da linha 1, se fosse 2p da linha 2 e assim por diante.
pixel= (cat dados_xy.txt |sed -n $count'p')

'set display color white'
'c'
'set t ' count


#Automatizando a data no titulo de cada figura ;)
'q dims'
dataarq=sublin(result,5)
dataarq2=subwrd(dataarq,6)
dia=substr(dataarq2,4,2)
hora=substr(dataarq2,1,2)


#Setando o mapa e seus contornos
'set map 1 0 6'

#Setando os níveis de tsm
'set clevs 24 26 28 30 32 34'
'set ccols 21 22 24 24 25 26 27'

'set gxout shaded'
'd sst - 273.15'

#Mascarando o continente para usar a TSM. O numero 15 significa a cor cinza
'basemap L 15 M'

#Setando os níveis de pressão
'set clevs 984 986 988 990 992 994 996 998 1000'
'set gxout contour'
'd psfc/100'

#Setando o caminho exato do IRMA na FIGURA
#Os dados de x e y usados no comando abaixo eu peguei atraves do script cospe_valores, onde ele transforma a lat e lon em valores cartesianos e joga os valores no arq dados_xy

'draw mark 6 'pixel' 0.20'


'cbar.gs'

'draw title WRF Furacao Irma'pixel' - 'dia'/SET/17 'hora'Z \ TSM e PNMM '
'printim 'count'_Irma_'hora'Z'dia'SET2017.png'

count=count+1




Espero que alguem possa me ajudar!