<div dir="ltr"><div class="gmail_default" style="font-size:small">[I had to reject my own ncl-talk post because my previous attachment was too large!]</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default">Adam,</div><div class="gmail_default"><br></div><div class="gmail_default">As I said offline, I have a script that draws your grid, but I'm not sure this is what you want.</div><div class="gmail_default"><br></div><div class="gmail_default">You are drawing the grid from a SCRIP file, and it looks like each boundary is a rectangle and not a pentagon. I recommend drawing the grid from the original source, rather than the SCRIP file.</div><div class="gmail_default"><br></div><div class="gmail_default">Anyway, you should be able to use a modified version of the attached code to draw the original grid. I actually showed three different ways to draw the grid, mainly to illustrate the difference in speed of each method. The timings are as follows:</div><div class="gmail_default"><br></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace">Elapsed time for method 1 using gsn_add_polyline = 15.9526</font></div><div class="gmail_default"><font face="monospace, monospace">Elapsed time for method 2 using gsn_polyline = 19.662</font></div><div class="gmail_default"><font face="monospace, monospace">Elapsed time for method 3 using gsn_add_polyline/gsSegments = 0.850868</font></div></div><div class="gmail_default"><br></div><div class="gmail_default">gsSegments can really speed things up if you have a lot of lines to add.</div><div class="gmail_default"><br></div><div class="gmail_default">I've attached the script and only one of the PNG images (since the other two should be identical).</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 7, 2016 at 11:18 AM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</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 style="font-size:small"><div>Adam,</div><div><br></div><div>As I said offline, I have a script that draws your grid, but I'm not sure this is what you want.</div><div><br></div><div>You are drawing the grid from a SCRIP file, and it looks like each boundary is a rectangle and not a pentagon. I recommend drawing the grid from the original source, rather than the SCRIP file.</div><div><br></div><div>Anyway, you should be able to use a modified version of the attached code to draw the original grid. I actually showed three different ways to draw the grid, mainly to illustrate the difference in speed of each method. The timings are as follows:</div><div><br></div><div><div><font face="monospace, monospace">Elapsed time for method 1 using gsn_add_polyline = 15.9526</font></div><div><font face="monospace, monospace">Elapsed time for method 2 using gsn_polyline = 19.662</font></div><div><font face="monospace, monospace">Elapsed time for method 3 using gsn_add_polyline/gsSegments = 0.850868</font></div></div><div><br></div><div>gsSegments can really speed things up if you have a lot of lines to add.</div><div><br></div><div>I've attached the script and only one of the PNG images (since the other two should be identical).</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>--Mary</div><div><br></div><div><br></div><div><br></div></font></span></div></div><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Thu, Apr 7, 2016 at 9:37 AM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></span> wrote:<br></span><div><div class="h5"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div style="font-size:small">I'll try to take a look soon. I had a tutorial to give yesterday and was not checking email for a bit.</div><div style="font-size:small"><br></div><div style="font-size:small">Usually, tools want the cells be oriented in counter clockwise order. This becomes critically important to have this consistency in order for things like regridding to be applied correctly.</div><span><font color="#888888"><div style="font-size:small"><br></div><div style="font-size:small">--Mary</div></font></span></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 4, 2016 at 11:43 AM, adamrhster . <span dir="ltr"><<a href="mailto:adamrhster@gmail.com" target="_blank">adamrhster@gmail.com</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">Hi Mary,<div><br></div><div>I'm trying to plot the cubed-sphere grid using cam-se data. This doesn't appear to be addressed in the "seam" examples. In the attached plot, the nodes are the red points that I want to connect with ~roughly orthogonal grid lines. The blue lines are a result of the cnFillMode = "CellFill" using vertices that aren't really of interest to me (these are for the interpolating fields to the coupler from file "<a href="http://ne4np4-pentagons.nc" target="_blank">ne4np4-pentagons.nc</a>")...Alternatively, the ESMF re-grid generated a sourcefile with seemingly opposite vertex and center coord data (3 vertices for each center coord, and the vertices appear to be on the actual computational nodes)</div><div><br></div><div>Maybe gsSegments will work...but I don't know a direction to orient the line segments since the nodes are not in any particular order. My script (and data) are on glade:</div><div><br></div><div>/glade/p/work/aherring/grids/seam.ncl</div><div><br></div><div><br></div></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 4, 2016 at 12:35 PM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</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 style="font-size:small">Hi Adam,</div><div style="font-size:small"><br></div><div style="font-size:small">We have a function called "gsn_coordinates" which draws grid points or grid lines for the given data array, but it doesn't handle the special case of an unstructured grid with N cells and NE edges.</div><div style="font-size:small"><br></div><div style="font-size:small">However, if your grid is defined such that you have an array of indexes that represent vertex edges for another array that contains the nodes, then you might be able to use the gsn_add_polyline routine with the special "gsSegments" resource.</div><div style="font-size:small"><br></div><div style="font-size:small">See example "mpas_faster_2.ncl" which shows the gsSegments method, and the mpas_2.ncl example which draws the same thing, but in a slower fashion using basic line draws. The point of this example is to show how much faster gsSegments can make your code.</div><div style="font-size:small"><br></div><div><a href="https://www.ncl.ucar.edu/Applications/mpas.shtml" target="_blank">https://www.ncl.ucar.edu/Applications/mpas.shtml</a><br></div><div style="font-size:small"><br></div><div style="font-size:small">The ICON page that Guido pointed you to also makes use of the gsSegments resource, in example icon_faster_2.ncl.</div><div style="font-size:small"><br></div><div><a href="https://www.ncl.ucar.edu/Applications/icon.shtml" target="_blank">https://www.ncl.ucar.edu/Applications/icon.shtml</a><br></div><div><br></div><div style="font-size:small">If you have an unstructured grid that is different from the ICON and the MPAS grids, then I would love to add this as yet another example on our website. I could help you create the script if necessary, but would need access to your data file. I could perhaps use this as a way to enhance gsn_coordinates to better handle unstructured grids.</div><span><font color="#888888"><div style="font-size:small"><br></div><div style="font-size:small">--Mary</div><div style="font-size:small"> </div><div style="font-size:small"><br></div><div style="font-size:small"><br></div></font></span></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Sun, Apr 3, 2016 at 3:01 PM, Adam Herrington <span dir="ltr"><<a href="mailto:adam.herrington@stonybrook.edu" target="_blank">adam.herrington@stonybrook.edu</a>></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><div dir="ltr">How would one plot grid lines for a data array on an unstructured grid? <div><br></div><div>I am able to plot grid boxes around nodes by setting cnFillMode = "CellFill" and defining the vertices of the grid box. But what if I want grid lines connecting the nodes? And what if I want to plot grid lines on every other node (not that the data arrays have one dimension for the horizontal dimension, containing all the tiles of the sphere --this is a global model)</div><div><br></div><div>Any help would be greatly appreciated</div><span><font color="#888888"><div><br></div><div>Adam</div></font></span></div>
<br></div></div><span>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></span></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div></div></div><br></div>
</blockquote></div><br></div></div>