Grads Plot do Vento Significante

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!