<div dir="ltr"><div>Hi everyone</div><div><br></div><div>I'm using latest precompile v6.6.2 ncl to plot with unstructured data.</div><div>sfXarray, sfYarray and sfElementNodes were used to do the contour.</div><div>Things are going well with FE mesh under 600~700k nodes (about 1.2 million elements).</div><div>However, if mesh increase to certain amount, say 1.8 million nodes (5 million elements), I have to use multi-threads with rasterfill to get reasonable efficiency, but the figure output is totally not right. <br></div><div>I have no chance to test areafill because after run few hours, no any contour plot can be 
generated.</div><div><br></div><div>I also try with Matlab, and get the correct contour plot. <br></div><div>The attach figures show the difference.</div><div> Matlab took about 90s to produce the plot with 1 core.</div><div> ncl took 130s under 6 threads,  20 mins with 1 thread.</div><div><br></div><div>From ncl website, I saw MPAS example with 2.6 million cell and get quite reasonable efficiency with rasterfill.</div><div>I notice 
sfElementNodes is not used in script.</div><div>When I comment out this resource, ncl can produce the plot about 60s with 6 threads, but without element table, the plot is still wrong.</div><div></div><div></div><div><br></div><div></div><div>Just curious, will it be any chance to get better efficiency and produce correct figure on PyNGL?</div><div><br></div><div>Dan<br></div><div><br></div></div>