<div dir="ltr"><div><div>Thank you,<br></div>Yes, with those modifications the result is satisfying.<br><br></div> Mathieu<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 17, 2014 at 2:15 AM, David Brown <span dir="ltr"><<a href="mailto:dbrown@ucar.edu" target="_blank">dbrown@ucar.edu</a>></span> wrote:<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I agree that something odd is happening when you try to plot this data<br>
using the 'ele' array in the file to define the triangulation.<br>
However, I was able to plot it quite easily using NCL's built-in<br>
triangulation code. The only change I needed to make from your<br>
original code is to comment out the lines:<br>
<br>
; res2@sfElementNodes = f->ele<br>
; res2@sfFirstNodeIndex = 0<br>
<br>
For the plot I am attaching I also changed<br>
<br>
res2@mpDataBaseVersion = "MediumRes"<br>
<br>
because the high res data does not work well for global plots.<br>
<br>
We will need to investigate why the explicit triangulation in your<br>
data set causes a problem. But you can get a better plot faster by not<br>
setting<br>
sfElementNodes and letting NCL compute the triangulation.<br>
<br>
Hope this helps.<br>
-dave<br>
<div><div class="h5"><br>
<br>
On Tue, Dec 16, 2014 at 10:24 AM, Mathieu Dutour<br>
<<a href="mailto:mathieu.dutour@gmail.com">mathieu.dutour@gmail.com</a>> wrote:<br>
> Dear all,<br>
><br>
> I would like to plot some unstructured grid that encompass all<br>
> the oceans of the world.<br>
><br>
> There are two obstacles to that. One is the periodicity in longitude<br>
> and the other is the singularity at the south and north pole. The<br>
> singularity can be avoided by removing points with latitude not<br>
> in acceptable range. The periodicity can be treated by removing<br>
> the points of longitude +-180.<br>
><br>
> A priori, the drawing should then proceed correctly, but it does not.<br>
> The grid I put as input passes all the sanity checks I can think of.<br>
> What makes me think that there is a bug is that if I restrict to<br>
> positive latitudes (0<lat<90) or negative latitude (-90 < lat < 0) then<br>
> the drawings are ok but it is not when I put everything together.<br>
><br>
> All relevant files for debugging are available in<br>
> <a href="http://mathieudutour.altervista.org/Sending/PLOT_globunstruc.tar.gz" target="_blank">http://mathieudutour.altervista.org/Sending/PLOT_globunstruc.tar.gz</a><br>
> The reduction technique from WholeWorld ---> UpperLatitude, LowerLatitude<br>
> is to remove the points of noncorrect latitude and any triangle that<br>
> contain them.<br>
><br>
> The three pictures are<br>
> <a href="http://mathieudutour.altervista.org/Sending/LowerLatitude.png" target="_blank">http://mathieudutour.altervista.org/Sending/LowerLatitude.png</a> (good)<br>
> <a href="http://mathieudutour.altervista.org/Sending/UpperLatitude.png" target="_blank">http://mathieudutour.altervista.org/Sending/UpperLatitude.png</a> (good)<br>
> <a href="http://mathieudutour.altervista.org/Sending/WholeWorld.png" target="_blank">http://mathieudutour.altervista.org/Sending/WholeWorld.png</a> (bad)<br>
><br>
> Thank you in advance for any help in the debugging and please ask<br>
> if I can help you.<br>
><br>
> Best,<br>
><br>
> Mathieu<br>
><br>
</div></div>> _______________________________________________<br>
> ncl-talk mailing list<br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
</blockquote></div></div>