<div dir="ltr"><div><div><div>Hi Mary & Dennis<br><br></div><div>Thanks for all your info.<br></div>I will try with ESMF routine.<br></div><br></div>Best regards, Dan<br><div><div><div><div class="gmail_extra"><br><div class="gmail_quote">2017-03-25 2:15 GMT+08:00 Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>The following suggests an approach to achieve your objective. <br><br>A user wanted to interpolate curvilinear (NARR) data to specific locations. The user also needed to vertically interpolate to specific pressure levels prior to interpolating to the locations.<br><br></div>Look at: <a href="http://www.ncl.ucar.edu/Applications/narr.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/narr.shtml</a><br></div>Example 7<br><br></div><div>Several steps were used:<br></div><div><br> [1] A previously generated ESMF weight file was used to regrid (interpolate) the NARR curvilinear grid to a high resolution rectilinear grid.<br><br></div><div> [2] The linint2_points function then used the generated rectilinear grid as the input to linint2_points_Wrap<br></div><div><br></div><div>===<br></div><div>There are several ESMF examples the interpolate 'spectral element' (triangular mesh) to rectilinear grids.<br></div><div>ESMF 10, 11, 17-19<br><br></div><div>[1] Using one of the examples as a base, generate a weight file (say): <br> Opt@WgtFileName="Trimesh_to_<wbr>RecLin.nc"<br><br></div>[2] <br>Based upon the function 'narr_to_points' in the above Example 7, you should create a (say) 'trimesh_to_points' function. If not appropriate, delete the vertical interpolation portion of the function. This function could be reused for other interpolations.<br><br></div><div><br></div><div>Good luck<br></div><div><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Mar 24, 2017 at 9:10 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">Dan,</div><div style="font-size:small"><br></div><div style="font-size:small">You might want to check out the ESMF regridding examples page:</div><div style="font-size:small"><br></div><div><a href="http://www.ncl.ucar.edu/Applications/ESMF.shtml" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/ESMF.shtml</a><br></div><div><br></div><div>Example #14 regrids a triangular mesh to a 5 degree grid:</div><div><br></div><div><a href="http://www.ncl.ucar.edu/Applications/ESMF.shtml#ex14" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/ESMF.shtml#ex14</a><br></div><div><br></div><div>There are also a couple of examples associated with the NCL User Guide:</div><div><br></div><div><a href="http://www.ncl.ucar.edu/Applications/NUG.shtml" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/NUG.shtml</a><br></div><div><br></div><div>Search for "triangular". You can actually download the data files for these examples and run them yourselves. The data file link is at the top of the page.</div><div><br></div><div>--Mary</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_-8393679765623293154h5">On Wed, Mar 22, 2017 at 1:41 AM, Yu, Hao-Cheng <span dir="ltr"><<a href="mailto:danishyo@gmail.com" target="_blank">danishyo@gmail.com</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 class="m_-8393679765623293154h5"><div dir="ltr"><div><div><div><div><div>Hi everyone<br><br></div>From website, I know there're lot of routines to do triangular plot and regriding.<br></div>Since plot can be made by predefined FE mesh by specified sfYArray, sfXArray & sfElementNodes.<br>I will assume there must be some interpolation routine for FE mesh, but I can't find this on website. <br><br></div>Is there any 2D interpolation routine we can use directly by specifying above 3 vars with desired output XY points?<br></div>For example, assume this routine called "tri_interp".<br>Pre-defined FE mesh var are X,Y,ele, field to interpolate is ZZ, output points are outX, outY.<br></div><div>Then we can get output points value by out_interp=tri_interp(X,Y,ele,<wbr>ZZ,outX,outY)<br><br></div><div>For regular grid, like linint2_points can provide similar function and efficient.<br></div><div>Just wonder, is there any similar function for uninstructed grid except Natgrid package (which will be quite time consuming if FE mesh is large)? <br></div><div><div><div><div><div><div><br></div><div>Dan<br></div><div><br></div></div></div></div></div></div></div>
<br></div></div>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
<br>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><br><br></div></div></div></div></div>