<div dir="ltr"><div class="gmail_default" style="font-size:small">[Sorry, I hit send too early]</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I didn&#39;t have a chance to tab the new code correctly, but hopefully you can see from the previous message that you don&#39;t need the double &quot;for&quot; loop. You simply set the special &quot;gsSegments&quot; resources to the &quot;segm&quot; array, and make sure you pass in LON and LAT as 1D arrays (maybe it would work without having to do this).</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Our PyNGL website is down at the moment due to scheduled semi-annual maintenance, but I think that &quot;shapefiles4.py&quot; and/or &quot;shapefiles2.py&quot; on the gallery page use this new resources. I&#39;ve attached shapefiles4.py here in case the website doesn&#39;t have the latest version.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you are interested in trying 1.5.0-beta of PyNGL, we have some binaries at (when the website is back up, hopefully Sunday):</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style><a href="http://www.pyngl.ucar.edu/Training/SEA2015/">http://www.pyngl.ucar.edu/Training/SEA2015/</a><br></div><div class="gmail_default" style><br></div><div class="gmail_default" style>which we created for the tutorial. Note that we also have PyNIO binaries listed here, which also has a number of new features. When we release 1.5.0 of both packages, we hope to do so under &quot;conda&quot;, so you don&#39;t have to download them this way.</div><div class="gmail_default" style><br></div><div class="gmail_default" style>Some of the binaries may have shared library issues, which are usually easy to fix.  Let me know if you try to download one of the binaries and have issues with it.</div><div class="gmail_default" style><br></div><div class="gmail_default" style>--Mary</div><div class="gmail_default" style><br></div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Apr 18, 2015 at 12:10 PM, Mary Haley <span dir="ltr">&lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Diego,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Sorry for the delay in response. Ironically, the delay was caused by us being at a Python conference this week at NCAR.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The problem with using Ngl.polyline, is that the lines are not &quot;attached&quot; to the plot, and hence when you do the maximize_plot call later, the lines are getting ignored and probably discarded.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">There are a couple of things you can try:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><ul><li>Move the #10 part of your code to *after* the #11 code, so that the lines get drawn after you maximize the plot and draw the map.</li></ul><div>or</div><div><br></div><div><ul><li>Updated to PyNGL 1.5.0-beta, and use the special &quot;gsSegments&quot; resource which allows you to add shapefile outlines with a single call. So instead of this code:</li></ul><div><br></div><div>







<p><font face="monospace, monospace"># 10. Add Natura2000 polylines to map                                                          </font><span style="font-family:monospace,monospace">for</span><span style="font-family:monospace,monospace"> i </span><span style="font-family:monospace,monospace">in</span><span style="font-family:monospace,monospace"> range(0,numFeatures):<br></span><span style="font-family:monospace,monospace">        startSegment = geom[i, geom_segIndex]<br></span><span style="font-family:monospace,monospace">        numSegments = geom[i, geom_numSegs]<br></span><span style="font-family:monospace,monospace">        </span><span style="font-family:monospace,monospace">for</span><span style="font-family:monospace,monospace"> seg </span><span style="font-family:monospace,monospace">in</span><span style="font-family:monospace,monospace"> range(startSegment, startSegment+numSegments):<br></span><span style="font-family:monospace,monospace">            startPT = segm[seg, segs_xyzIndex]<br></span><span style="font-family:monospace,monospace">            endPT = startPT + segm[seg, segs_numPnts]</span><span style="font-family:monospace,monospace"># - 1                          <br></span><font face="monospace, monospace">            Ngl.polyline(wks, mp, LON[startPT:endPT], LAT[startPT:endPT], natres</font>)</p><p><br></p><p>It would look something like this::</p><p><br></p><p><font face="monospace, monospace">LON1D = numpy.ravel(LON)<br>LAT1D = numpy.ravel(LAT)</font></p></div></div><div>
<p><font face="monospace, monospace">natres.gsSegments  = segm<br>id = Ngl.add_polyline(wks, mp, LON1D, LAT1D, natres)</font></p><p><br></p></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Mon, Apr 13, 2015 at 4:40 AM, Diego Bruciaferri <span dir="ltr">&lt;<a href="mailto:diego.bruciaferri@ingv.it" target="_blank">diego.bruciaferri@ingv.it</a>&gt;</span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div><div>Dear all,<br><br></div>I&#39;m trying to plot several polylines from a shapefile over a map which contains both scalar (oil concentration) and vector fields (current and stokes&#39; drift fields).<br><br></div><div>As I&#39;ve seen in yours example scripts to plot shapefiles, I want to use the Ngl.polyline function over the Ngl.add_polyline because the first is faster, but it doesn&#39;t work (I don&#39;t see the lines) for me.<br><br></div><div>I&#39;ve tried to &quot;play&quot; with the DrawOrder resources, without any effect.<br></div><div>Can you please suggest me some possible solution?<br><br></div><div>Attached, you can find the code. <br></div><div>Thanks in advance!<span><font color="#888888"><br><br></font></span></div><span><font color="#888888"><div>Diego<br></div><div><div><div><div><br>-- <br><div><div dir="ltr">Diego Bruciaferri
<br>Istituto Nazionale di Geofisica e Vulcanologia (INGV)
<br>National Institute of Geophysics and Volcanology
<br>
<br>Viale Aldo Moro 44
<br>40127 Bologna
<br>Italy
<br>Tel: <a href="tel:%2B39%20051%203782665" value="+390513782665" target="_blank">+39 051 3782665</a> <br></div></div>
</div></div></div></div></font></span></div>
<br></div></div>_______________________________________________<br>
pyngl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/pyngl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/pyngl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>