<div dir="ltr">Mike, Jennifer,<div><br></div><div>Thanks for the illustration Mike--I understood what you did--but it does not accomplish what I&#39;m after.  However, it does give me some ideas for trying something else that may be suitable for automation.</div><div><br></div><div>Jennifer, you said user-provided &#39;set clip&#39; commands...do work as expected with &#39;draw&#39;.  As I mentioned earlier, I have not yet seen an example of its use (not in Doc Index either).  Could you please provide a simple sample that illustrates how to use it with &#39;draw&#39;?  I think the example that Mike provided was an alternate to using &#39;clips&#39;.  Perhaps I simply misunderstood its intended use.</div><div><br></div><div>Stephen Mc</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, May 30, 2017 at 3:54 PM, Jennifer M Adams <span dir="ltr">&lt;<a href="mailto:jadams21@gmu.edu" target="_blank">jadams21@gmu.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I looked at the source code, and it’s true that the user-provided ‘set clip’ commands are always overridden with ‘display’, but as Mike F. pointed out, they do work as expected with ‘draw&#39;. I will put this is my list of bugs, but in the meanwhile Mike’s suggestion to draw a few more rectangles on top of your graphic is a good workaround. —Jennifer<br>
<span class=""><br>
On May 30, 2017, at 1:54 PM, Mike Fiorino - NOAA Federal &lt;<a href="mailto:michael.fiorino@noaa.gov">michael.fiorino@noaa.gov</a>&gt; wrote:<br>
<br>
&gt; hi Steve,<br>
&gt;<br>
</span><span class="">&gt; now I understand... you&#39;ll have to do the clipping by hand using:<br>
&gt;<br>
&gt; &#39;draw recf xlo xhi ylo yhi&#39;<br>
&gt;<br>
&gt; which fills in a rectangle set by xlo .. to the current &#39;line&#39; color<br>
&gt;<br>
&gt; here&#39;s an example of clipping to a block from x=2.0 to 6.0 and y = 3.0 to 6.0<br>
&gt;<br>
&gt; first display something:<br>
&gt;<br>
&gt; &#39;d psl*0.01&#39;<br>
&gt;<br>
&gt; then get the dimensions of the plotting area<br>
&gt;<br>
&gt; &#39;q gxinfo&#39;<br>
&gt;<br>
&gt; Last Graphic = Contour<br>
&gt; Page Size = 11 by 8.5<br>
&gt; X Limits = 0.5 to 10.5<br>
&gt; Y Limits = 1.25 to 7.25<br>
&gt; Xaxis = Lon  Yaxis = Lat<br>
&gt; Mproj = 2<br>
&gt;<br>
&gt; then:<br>
&gt;<br>
&gt; &#39;draw recf 0.5 10.5 1.25 3.0&#39;  # xlo for plot is 0.5 and xhi is 10.5<br>
&gt; &#39;draw recf 6.0 3.00 10.5 7.25&#39; # ylo for plot is 1.25 and xhi is 7.25 -- this is the default parea<br>
&gt; &#39;draw recf 0.5 3.00 2.00 7.25&#39;<br>
&gt; &#39;draw recf 2.0 6.00 6.00 7.25&#39;<br>
&gt;<br>
&gt; # -- draw a box<br>
&gt; &#39;set line 1 0<br>
&gt; &#39;draw rec 2.0 3.0 6.0 6.0&#39;<br>
&gt;<br>
&gt; here&#39;s the plot (ecmwf sea level pressure for 2017052800 analysis -- note the tropical cyclone in the Bay of Bengal (02B)):<br>
&gt;<br>
</span>&gt; &lt;bbknhkmnngojjgok.png&gt;<br>
<span class="">&gt;<br>
&gt; hope this helps...<br>
&gt;<br>
&gt; best /R Mike<br>
&gt;<br>
&gt; On 5/25/17 3:32 PM, Stephen McMillan wrote:<br>
&gt;&gt; Mike,<br>
&gt;&gt;<br>
&gt;&gt; Thanks for the response, but I don&#39;t think you understood what I&#39;m trying to accomplish.<br>
&gt;&gt;<br>
&gt;&gt; Imagine you&#39;ve displayed a variable, using the default plot area, and it displays across most of North America.  Now, let&#39;s say you want to display the same variable, but only including part of the  western US, indicated by the dashed rectangle in attached image.  What I want to do is display the zoomed-in area but fill the screen to the default plot area, not just that part of the screen.<br>
&gt;&gt;<br>
&gt;&gt; If I use &#39;set parea...&#39; as you suggest, e.g., &#39;set parea 2.07 3.44 5.35 6.11&#39; I then get the third image shown.  It displays the entirety of the original plot, not the zoomed-in area.<br>
&gt;&gt;<br>
&gt;&gt; I&#39;ve tried using &#39;set clip...&#39; with a draw command, but no luck.  Could you send a simple code example to illustrate?  I have yet to see how it&#39;s used--the Doc Index doesn&#39;t help in that regard.<br>
&gt;&gt;<br>
&gt;&gt; Steve Mc<br>
&gt;&gt;<br>
</span>&gt;&gt; &lt;image.png&gt;<br>
&gt;&gt;<br>
&gt;&gt; &lt;image.png&gt;<br>
&gt;&gt;<br>
&gt;&gt; &lt;image.png&gt;<br>
<div><div class="h5">&gt;&gt;<br>
&gt;&gt; On Wed, May 24, 2017 at 5:52 PM, Mike Fiorino - NOAA Federal &lt;<a href="mailto:michael.fiorino@noaa.gov">michael.fiorino@noaa.gov</a>&gt; wrote:<br>
&gt;&gt; hi Steve,<br>
&gt;&gt;<br>
&gt;&gt; I only use &#39;set clip&#39; for &#39;draw&#39; commands and it does work, but apparently not for &#39;display&#39; (took a quick look at the source code and it wasn&#39;t obvious why it wouldn&#39;t for &#39;display&#39; except that &#39;set clip&#39; works at the very low level plot commands)<br>
&gt;&gt;<br>
&gt;&gt; have you tried the &#39;parea&#39; command?  this sets the plot area, e.g.,<br>
&gt;&gt;<br>
&gt;&gt; &#39;set parea 1 6 3 7&#39;<br>
&gt;&gt; &#39;d air&#39;<br>
&gt;&gt;<br>
&gt;&gt; note that the labels are outside the parea...<br>
&gt;&gt;<br>
&gt;&gt; best /R Mike<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; On 5/24/17 7:36 PM, Stephen McMillan wrote:<br>
&gt;&gt;&gt; Is the &#39;set clip&#39; command disabled in GrADS?  I can&#39;t get it to work per the Documentation Index instructions.  A search of the User Forum yields no help, but one entry does suggest that it may indeed be disabled without modifying the code.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; For example (file open with variable &quot;air&quot;):<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; &#39;set clip 1 6 3 7&#39;<br>
&gt;&gt;&gt; &#39;d air&#39;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Result is a map displayed to the default plot area, not to the &quot;clip&quot; area.  I would normally use &#39;set mpvals&#39; but cannot get that to work the way I want it, either.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;  I am using v2,1,a2,iga,1 on a Windows 10 PC.  I also cannot get it to work using v2.0.2.oga.2 on a Linux.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; I don&#39;t need suggestions on how to crop images outside GrADS.  I know how to do that.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Stephen McMillan<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; ______________________________<br>
&gt;&gt;&gt; _________________<br>
&gt;&gt;&gt; gradsusr mailing list<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt;&gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;<br>
&gt; ______________________________<wbr>_________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
<br>
</div></div>--<br>
Jennifer Miletta Adams<br>
Center for Ocean-Land-Atmosphere Studies (COLA)<br>
George Mason University<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
<br>
______________________________<wbr>_________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
</div></div></blockquote></div><br></div>