[gradsusr] Plot wind time series
Davide Sacchetti
davide.sacchetti at arpal.gov.it
Thu Jan 18 04:14:42 EST 2018
it's tricky
don't fix x and y; just fix one of them
eg: you want to display wind series of point x=100 y=100:
'set x 100
'set y 99.5 100.5'
'set t 1 last'
'set xyrev on'
'd u;v;mag(u,v)'
actually on Y axis you will have latitude, but it works ...
Bye
Davide
On Wed, 2018-01-17 at 10:01 -0500, MIGUEL FIDENCIO SAAVEDRA JUAN DE
DIOS wrote:
> Hello all:
> I want to get wind graphs as the attached figure (windvecSabine.png).
>
>
> With the grads script that I attached, the results are not as
> expected (vector.png) .
> Please someone have a suggestion?
>
> $ cat vector-time.gs
> 'open total2015.ctl'
> 'set display color white'
> 'set lat -10'
> 'set lon -81'
> 'set t 1 10'
> 'set gxout vector'
> *'set vpage -5 5 -8 8'
> 'd const(u,0);v;mag(u,v)'
> *'d u;v;mag(u,v)d u;v;mag(u,v)'
> *'d u-u;v;mag(u-u,v)'
> 'printim vector.png'
> '!display vector.png'
>
> my ctl is:
> cat total2015.ctl
> DSET ^2015%m2%d2.dat
> OPTIONS TEMPLATE
> UNDEF -9.99e+08
> XDEF 1440 LINEAR 0.125 0.25
> YDEF 641 LINEAR -79.875 0.25
> ZDEF 1 LEVELS 10
> TDEF 360 LINEAR 12Z01JAN2008 1dy
> VARS 2
> u 0 99 eastward_wind
> v 0 99 northward_wind
> ENDVARS
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
More information about the gradsusr
mailing list