[gradsusr] Plotting Station data

[Alexandre Santos]

Hello,

How do for plot lines contour with values?

Thanks


2010/3/19 Jennifer Adams <jma at cola.iges.org>

> Hi, Lukas --
> This question has come up before, but the archives are hard to search
> through at the moment. Attached below is the question and the script
> fragment and image sample that illustrate how it works.
> --Jennifer
>
> On Mar 19, 2010, at 9:53 AM, Lukas Mosimann wrote:
>
> Hello everybody
>
> I have already used Grads for a few months but I have never tried to plot
> station data. And now I tried to do so - but I didn't have much effort ;)
>
> I would like to plot the temperatures of some weather stations in different
> colors depending on the value. I tried the following
>
> ...
>
> set gxout value
>
> set ccolor rainbow
>
> d var
>
>
> But I only got white values - but I would like to get colored values (not
> points. When I use "set gxout stnmark" instead of value, I don't have any
> problems).
> Is there any possibility to plot such a map with station data? And what did
> I do wrong?
>
> Thanks a lot.
>
> Lukas
>
>
> -------------------------------------------------
>
> On Aug 10, 2009, at 3:55 AM, Jennifer Faber wrote:
>
> Hi Jennifer
>
> I was wondering if you and your team are considering adding a feature fto
> GrADS to colour station value according to their value.
> Until so far there was been no good solution to my problem.
>
>
> Here is a script fragment I composed, followed by the image output. I think
> this is what you are trying to do...
>
> 'set lat 22 55'
> 'set lon -130 -60'
> var = 'ts'
> levs = '-15 -10 -5 0 5 10 15'
> cols = '9  4   5  3 7 8  2  6'
> ncols=8
> c=1
> while (c<=ncols)
>   col=subwrd(cols,c)
>   maxlev=subwrd(levs,c)
>   if (c>1)
>     minlev=subwrd(levs,c-1)
>   endif
>   'set ccolor 'col
>   'set gxout value'
>   'set stat on'
>   if (c=1)
>     belowmax='(maskout('var',('maxlev'-0.0001)-'var'))'
>     'd 'belowmax
>   else
>     minandabove='(maskout('var','var'-('minlev')))'
>     if (c=ncols)
>       'd 'minandabove
>     else
>       belowmax='maskout('minandabove',('maxlev'-0.0001)-'minandabove')'
>       'd 'belowmax
>     endif
>   endif
>   c=c+1
> endwhile
>
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100319/eca36380/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/png
Size: 20490 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20100319/eca36380/attachment-0003.png