<div dir="ltr"><div><div><div>Regarding display map outside shp file <br><br></div>I got the same problem ... but i have done it in a different way which includes some manual work . I have taken the lat,lon points of my region and i created a matrix with ones at the required region with required resolution and remaining i made it zero using grads and i made this into grd file inorder to have equal size as that of our original data and multiplied the data with data in grd .. so that it can display only regions which u r interested .<br>
<br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, May 31, 2013 at 9:30 PM, <span dir="ltr"><<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send gradsusr mailing list submissions to<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of gradsusr digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Issue with creating cross section of vertical wind vectors<br>
and reflectivities. (Nimrod Micael)<br>
2. Clean regions outside a shapefile (Roger Rodrigues Torres)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Fri, 31 May 2013 11:38:19 -0400<br>
From: Nimrod Micael <<a href="mailto:nmicael@aggies.ncat.edu">nmicael@aggies.ncat.edu</a>><br>
Subject: [gradsusr] Issue with creating cross section of vertical wind<br>
vectors and reflectivities.<br>
To: <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
Message-ID:<br>
<<a href="mailto:CA%2B_TaRNTxmc9dE9WBwJ4sDwSVBG83R7Q_CC2xd8sRgzojG6V8Q@mail.gmail.com">CA+_TaRNTxmc9dE9WBwJ4sDwSVBG83R7Q_CC2xd8sRgzojG6V8Q@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hey Grad users, I am trying to figure out what went wrong with my script. I<br>
am trying to plot a vertical cross section of mag of u, v wind with max_dbz<br>
at the locations 80.5 W (set lat 276.5)-83.5 W (set lat 279.5) and 34.5<br>
N-37.2 N. I am getting an error below:<br>
<br>
<br>
Errors received:<br>
<br>
Operation error: Invalid dimension environment<br>
Looping dimension does not vary<br>
Operation error: Invalid dimension environment<br>
Looping dimension does not vary<br>
Cannot plot color bar: No shading information<br>
<br>
'reinit'<br>
'clear'<br>
'open 20110416_d03.ctl'<br>
'set display color white'<br>
'set map 1 1 5'<br>
'set annot 1 6'<br>
'set xlopts 1 6 .17'<br>
'set ylopts 1 6 .16'<br>
'set mpdset hires'<br>
'set clab on'<br>
'set grid off'<br>
'set timelab off'<br>
'set gxout grfill'<br>
<br>
<br>
********************************<br>
* Winds Speed At SFC Plot (shaded) in kts<br>
********************************<br>
'set t 169'<br>
'set gxout shaded'<br>
'set lev 1000 100'<br>
'd max_dbz'<br>
<br>
<br>
<br>
**************************************<br>
*<br>
* Settings for vector plot<br>
***************************************<br>
'set lev 1000 100'<br>
'set gxout vector'<br>
'set ccolor 49'<br>
'd skip(u,4,4)*1.94;v*1.94'<br>
<br>
<br>
<br>
<br>
*********************************************<br>
* Draw Title, Run Cbar, Draw Map<br>
**********************************************<br>
'draw map'<br>
'run cbarn'<br>
<br>
<br>
<br>
<br>
--<br>
Nimrod Micael<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20130531/b326d4bc/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20130531/b326d4bc/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Fri, 31 May 2013 12:48:41 -0300<br>
From: Roger Rodrigues Torres <<a href="mailto:roger.torres@cptec.inpe.br">roger.torres@cptec.inpe.br</a>><br>
Subject: [gradsusr] Clean regions outside a shapefile<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<CAEOS777EkjH9rN+e+d9hrGk1FR2TYHfHxvz7E=<a href="mailto:RO%2BcFogmbE%2BA@mail.gmail.com">RO+cFogmbE+A@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Dear Grads Users,<br>
I am displaying a map in Opengrads, and trying to clean the results outside<br>
a specific region defined in a Shapefile,<br>
in a way similar when using "basemap" to clean ocean or land.<br>
I tried "shp_polyf" and "draw shp", but what I could obtain was only the<br>
countour of the region.<br>
Could someone help me please?<br>
Thanks in advance,<br>
Roger Torres<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20130531/931603d6/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20130531/931603d6/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
End of gradsusr Digest, Vol 39, Issue 87<br>
****************************************<br>
</blockquote></div><br><br clear="all"><br>-- <br><div>Umakanth U </div>
<div>Jr. Research Fellow</div>
<div>National Atmospheric research laboratory.<br></div>
<div>Tirupati.</div>
</div>