<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 "attached" 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" style><ul style><li style>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 "gsSegments" resource which allows you to add shapefile outlines with a single call. So instead of this code:</li></ul><div><br></div><div>
<p class=""><font face="monospace, monospace"># 10. Add Natura2000 polylines to map </font><span class="" style="font-family:monospace,monospace">for</span><span style="font-family:monospace,monospace"> i </span><span class="" 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 class="" style="font-family:monospace,monospace">for</span><span style="font-family:monospace,monospace"> seg </span><span class="" 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 class="" 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 class=""><br></p><p class="">It would look something like this::</p><p class=""><br></p><p class=""><font face="monospace, monospace">LON1D = numpy.ravel(LON)<br>LAT1D = numpy.ravel(LAT)</font></p></div></div><div>
<p class=""><font face="monospace, monospace">natres.gsSegments = segm<br>id = Ngl.add_polyline(wks, mp, LON1D, LAT1D, natres)</font></p><p class=""><br></p></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 13, 2015 at 4:40 AM, Diego Bruciaferri <span dir="ltr"><<a href="mailto:diego.bruciaferri@ingv.it" target="_blank">diego.bruciaferri@ingv.it</a>></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><div>Dear all,<br><br></div>I'm trying to plot several polylines from a shapefile over a map which contains both scalar (oil concentration) and vector fields (current and stokes' drift fields).<br><br></div><div>As I'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't work (I don't see the lines) for me.<br><br></div><div>I've tried to "play" 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 class="HOEnZb"><font color="#888888"><br><br></font></span></div><span class="HOEnZb"><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>_______________________________________________<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>