[gradsusr] Hello to all users of Grads

[Mohsen Soltani]

Hi Avi,

I'm not quite sure what to say about your script. But you may use the
following commands to draw average trend+ boxes of standard deviation for
your variable.

****************
var='slp/100'
minlat='36'
maxlat='37'
minlon='49'
maxlon='55'
tlast=28
'set x 1'
'set y 1'
'set t 1 'tlast
'define mean =
tloop(aave('var',lon='minlon',lon='maxlon',lat='minlat',lat='maxlat'))'
'define pstd = mean +
tloop(aave(sqrt(pow('var'-mean,2)),lon='minlon',lon='maxlon',lat='minlat',lat='maxlat'))'
'define mstd = mean -
tloop(aave(sqrt(pow('var'-mean,2)),lon='minlon',lon='maxlon',lat='minlat',lat='maxlat'))'
'set vrange 1004 1028'
'set t 0.28 'tlast+0.28
'set gxout bar'
'd mstd;pstd'
'set gxout line'
'd mean'
**************************
The attached is the output of the above commands. As you sea, I used SLP
variable for a period of time in my work. But, you can apply it for any
variable at different levels of the atmosphere.

Hope you find it useful!

Mohsen

--
April showers bring May flowers!
--
*Mohsen Soltani, M.Sc.*
Climatology Grad Student,
Department of Climatology, Faculty of Geography,
University of Tehran, Tehran, Iran
voice: (+98) 9119772934
e-mail: soltani.clima at gmail.com


On Tue, Jul 9, 2013 at 1:44 AM, avi kojokro <kojokroavi at gmail.com> wrote:

> I have a script of the Ensemble Average + standard deviation, I want to
> make this line to be automatically:
>
> 'set vrange 'min1' 'max1''
>
> To do something similar to these lines:
>
>
> 'set lev 500'
> 'set t 1'
> 'define mx=max(hgtprs,t=1,t=65)'
> 'd mx'
> m1=subwrd(result,4)
> 'set lev 500'
> 'define mn=min(hgtprs,t=1,t=65)'
> 'd mn'
> m2=subwrd(result,4)
> 'set lev 500'
> 'set vrange 'm2-1' 'm1+1
>
> But I can not attach these lines to the script.
> Maybe someone know if i can arrange the script properly?
>
> I use Windows, and version of Grads 2.0.a9.
>
> Thanks
> Avi.
>
> _______________________________________________
> 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/20130709/5276cccb/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MeanTrend+SD.jpg
Type: image/jpeg
Size: 316571 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130709/5276cccb/attachment-0003.jpg