<div dir="ltr"><div><div><div><div><br><br></div>Unfortunately (for you) we are teaching the NCL WKshop this week. Hence, we are quite busy.<br><br></div>I did look ...<br><br></div>The rgrid2rcm uses a very simple brute force algorithm. *All* coordinates are required to be monotonically increasing. <br></div><div>However, after the 'where' statement, if you print the 'lon2d' you will see the issue. <br></div><div><br><br></div>(row,col)<br><div>[snip]<br></div><div>(58,40) 352.3439<br>(58,41) 355.4963<br>(58,42) 358.3603<br> ---------------------> not monotonic<br>(58,43) 0.9611969<br>(58,44) 3.324268<br>(58,45) 5.473709<br>(58,46) 7.431938<br><br></div><div>The reason for the gap in your plot is that in the "seam" occurs where the longitudes go from 3?? to 0 fails (not monotonic). It then 'works' again where the longitudes are increasing.<br></div><div><br>---<br></div><div>That said .... I initially thought there is something unusual about the source file. There still may be but I just do not have the time to look. Sometime ... yes but not now. Sorry.<br><br>---<br><br></div><div>I randomly took the ESMF example 9 and more or less force fed your file info info. It is not working but you can use it as a template. Again, just do not have the time.<br><br></div><div>Good luck<br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 9, 2015 at 3:16 PM, Michael Stone <span dir="ltr"><<a href="mailto:mrstone@colorado.edu" target="_blank">mrstone@colorado.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">Hi again,<br><div dir="ltr"><div><div><div><div><br></div>I've continued to investigate the issue and decided to try out using the ESMF regridding package instead. I wrote a script that does a simple interpolation using the same data and grids but using the ESMF regridding/interpolation and ran into the same issue with missing values. This script is attached below and requires the same "geo_em" file in my prior email as well as the data file. I plotted the missing values and they are in the same location. <br><br></div>To me, this indicates some issue with how I'm creating the geo_em file but I am having a hard time finding out why that might be. I'm going to continue to work on this and try a different grid to see if I run into the same issue. If anyone has any ideas I would greatly appreciate it. <br><br></div>Thanks,<br></div>Michael<span class=""><br><div><div><div><div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 8, 2015 at 2:26 PM, Michael Stone <span dir="ltr"><<a href="mailto:mrstone.colorado@gmail.com" target="_blank">mrstone.colorado@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"><div><div><div>Hi All,<br><br></div>I'm having issues interpolating data using the rgrid2rcm function from NCEP to WRF coordinates. After interpolation there is a strip of missing data and I'm having a very difficult time figuring out exactly why that is. <br><br>Attached you'll find the script I'm running to interpolate the data, a plot of where the missing values occur, and the location of the lats/lons that are missing. The NCEP data is just air temperature (although I'm trying to do many variables) and I put a single year on the NCAR ftp server under "<a href="http://air.1950.nc" target="_blank">air.1950.nc</a>". There is also a file "<a href="http://geo_em.d01.nc" target="_blank">geo_em.d01.nc</a>" that represents the WRF grid to interpolate to that is necessary for the script to run located on the ftp. <br><br></div><div>ncl -V<br><br>6.2.0<br><br></div><div>uname -a <br><br>Linux <a href="http://vardar.colorado.edu" target="_blank">vardar.colorado.edu</a> 2.6.32-504.12.2.el6.i686 #1 SMP Sun Feb 1 12:14:25 EST 2015 i686 i686 i386 GNU/Linux<br></div><div><br></div>Thanks for the help.<span><font color="#888888"><br><br></font></span></div><span><font color="#888888">Michael<br></font></span></div>
</blockquote></div><br></div></div></div></div></div></div></div></span></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><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><br></div>