<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>The function works. Read the documentation.</div><div><br></div><div>WRF data are ordered ( [time,] <b>bottom_top, </b>south_north x west_east ) ==> (time,'z','lat','lon')<br></div><div><br></div><div>See appended script. For demonstration only 2 points are used. <br></div><div><br></div><div>Variable: x<br>Number of Dimensions: 3<br>Dimensions and sizes: <b>[bottom_top | 34] </b>x [south_north | 160] x [west_east | 180]<br><br>Variable: lat2d<br>Number of Dimensions: 2<br>Dimensions and sizes: [south_north | 160] x [west_east | 180]<br><br>Variable: lon2d<br>Number of Dimensions: 2<br>Dimensions and sizes: [south_north | 160] x [west_east | 180]<br><br></div><div dir="ltr"><br></div><div dir="ltr">Variable: <b>xNew</b><br>Number of Dimensions: 2<br>Dimensions and sizes: <b>[34] x [2] <== levels, 2 random points</b><br><br>(0,0) 10.26271<br>(0,1) -0.3363198<br>(1,0) 10.02834<br>(1,1) -0.3514601<br>(2,0) 9.985236<br>(2,1) -0.3351685</div><div>[snip]</div><div>(32,0) 163.538<br>(32,1) 168.1733<br>(33,0) 182.0521<br>(33,1) 185.6072</div><div><br></div><div>=============</div><div><br></div><div> nt = 0<br> a = addfile("wrfout_d01_2003-07-15_00:00:<a href="http://00.nc">00.nc</a>","r")<br> x = a->T(nt,:,:,:)<br> printVarSummary( x ) ; [bottom_top | 34] x [south_north | 160] x [west_east | 180]<br> print("-----")<br><br> lat2d = a->XLAT(nt,:,:)<br> lon2d = a->XLONG(nt,:,:)<br><br> printVarSummary(lat2d) ; [south_north | 160] x [west_east | 180]<br> printMinMax(lat2d, True)<br> printVarSummary(lon2d) ; [south_north | 160] x [west_east | 180]<br> printMinMax(lon2d, True)<br> print("-----")<br><br> latPts = (/ 22., 29./)<br> lonPts = (/-100., -93./)<br> xNew = rcm2points(lat2d, lon2d, x, latPts, lonPts, 1)<br> print(xNew)<br><br></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Fri, Nov 9, 2018 at 3:04 PM Raj K. Rai <<a href="mailto:rajkrai@gmail.com">rajkrai@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><pre style="white-space:pre-wrap;color:rgb(0,0,0)">Hi Jim, Thanks for your email. Function rcm2points may work for 2d data, but I want to interpolate variables from 3d wrfoutout field (lat, lon, and z) into the given scattered points in space (given lat, lon, and z). </pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)">Thanks,</pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)">Raj </pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><br></pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)">Hi Raj,
I believe the function rcm2points should work for that.
<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/rcm2points.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/rcm2points.shtml</a>
Jim
On 11/1/2018 10:26 AM, Raj K. Rai wrote:
><i> Hi,
</i>><i> I am trying to interpolate variables from WRF output into the
</i>><i> specified set of points, which are scattered in space. I found
</i>><i> routines that interpolates variables into vertical and horizontal
</i>><i> plane, but could not figure out if there is any routine for this task.
</i>><i> Could you please help me finding this?
</i>><i> Thanks,
</i>><i> Raj
</i>><i>
</i>><i>
</i>><i> _______________________________________________
</i>><i> ncl-talk mailing list
</i>><i> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">ncl-talk at ucar.edu</a>
</i>><i> List instructions, subscriber options, unsubscribe:
</i>><i> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a>
</i>
--
Untitled Document
James D. Means, Ph.D.
Senior Atmospheric & Climate Scientist
California Office
Tele: 619-495-1638 | Fax: 775-636-8430
930 Tahoe Blvd., Suite 802-560
Incline Village, Nevada 89451
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">jim at weatherextreme.com</a> <mailto:<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">jim at weatherextreme.com</a>> | vcard
<<a href="http://www.weatherextreme.com/vcards/James%20Means.vcf" target="_blank">http://www.weatherextreme.com/vcards/James%20Means.vcf</a>>
<a href="http://www.weatherextreme.com" target="_blank">www.weatherextreme.com</a> <<a href="http://www.weatherextreme.com/" target="_blank">http://www.weatherextreme.com/</a>></pre></div><br><div class="gmail_quote"><div dir="ltr">On Thu, Nov 1, 2018 at 10:26 AM Raj K. Rai <<a href="mailto:rajkrai@gmail.com" target="_blank">rajkrai@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi, <div>I am trying to interpolate variables from WRF output into the specified set of points, which are scattered in space. I found routines that interpolates variables into vertical and horizontal plane, but could not figure out if there is any routine for this task. Could you please help me finding this?</div><div>Thanks,</div><div>Raj</div></div>
</blockquote></div>
_______________________________________________<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>
</blockquote></div>