<div dir="ltr"><div>THX for the grib files.</div><div><br></div><div>I modified one of the ESMF examples. The attached script <b>ESMF_stovern.ncl:</b></div><div>(a) reads the desired variable and the source and destination grids</div><div>(b) interpolates using bilinear interpolation</div><div>(c) plots the source and interpolated variables</div><div>(d) creates a netCDF file containing the interpolated variable on the destination grid.</div><div>===</div><div>The initial interpolation took surprisingly long: 72 seconds on my *old* MAC.</div><div>The ESMF software is a separate/external suite of software. This fact is hidden from users.</div><div>===</div><div>The good news is that once the weight file is generated, the interpolation can be performed in less that a second :-)</div><div><br></div><div><br></div><div>;---ESMF_regrid with weights generated in a previous script<br></div><div> wgtFile = "QPF_weight_bilinear".nc"</div><div> var_regrid = <b>ESMF_regrid_with_weights</b>(var,<b>wgtFile</b>,False) ; less than a second<br></div><div>;---<br></div><div></div><div><br></div><div>Cheers</div><div>D<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Apr 29, 2020 at 8:15 PM Dennis Shea <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Diana,</div><div><br></div><div>Bill Schfetic's response is correct.: ESMF is the way to go.</div><div><br></div><div><a href="http://www.ncl.ucar.edu/Applications/ESMF.shtml" target="_blank"><b>http://www.ncl.ucar.edu/Applications/ESMF.shtml</b></a></div><div>ESMF_{16, 20} are two examples.</div><div>---</div><div>Note that the regridding requires the source and destination lat/lon grid values.</div><div>Map projection information is not used.</div><div>---</div><div>If you can make the source file and the destination grid coordinates available [drop box, ftp, ...] , I could take a look.</div><div><br></div><div>ftp <a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a></div><div>anonymous</div><div>your_email</div><div>cd incoming</div><div>put SOURCE_FILE</div><div>quit</div><div><br></div><div><br></div><div>Stay Healthy!</div><div>D</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Apr 29, 2020 at 7:51 PM William Scheftic via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<p>Hi Diana,</p>
<p>I haven't done your specific regridding. But using the ESMF
regridding routines is the best way to go. You can save the
regridding in a special netCDF file that only you requires the
regridding weights to be calculated once and then reloaded for
future needs. There are some examples that deal with
non-rectilinear grids sending in 2-d lat and lon grids, rather
than 1-d.<br>
</p>
<p><a href="http://www.ncl.ucar.edu/Applications/ESMF.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/ESMF.shtml</a></p>
<p><a href="http://www.ncl.ucar.edu/Document/Functions/ESMF.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/ESMF.shtml</a></p>
<p>Cheers,</p>
<p>Bill</p>
<p><br>
</p>
<div>On 4/29/2020 3:17 PM, Diana Stovern -
NOAA Affiliate via ncl-talk wrote:<br>
</div>
<blockquote type="cite">
<p style="text-align:center"><font color="red"><strong>External
Email</strong><br>
</font></p>
<div dir="ltr">Hello,
<div><br>
</div>
<div>I need to regrid rainfall data from a lambert conformal
grid to a polar sterographic grid. The lat and lon data in
both grids are 2D arrays (see filedump output below). Is
there a regridding routine in NCL that allows for both the
input and output lat and lon to be 2D arrays? I see most of
the routines require lat and lon to be 1D arrays.</div>
<div><br>
</div>
<div>Thanks! </div>
<div><br>
</div>
<div>Diana Stovern </div>
<div><br>
</div>
<div><br>
</div>
<div><b>Grid 1:</b> </div>
<div> float gridlat_0 ( ygrid_0, xgrid_0 )<br>
corners : ( 23.117, 19.80368, 45.61871, 53.51192
)<br>
long_name : latitude<br>
grid_type : Polar Sterographic Projection (North
or South)<br>
units : degrees_north<br>
Dy : 4.763<br>
Dx : 4.763<br>
Lov : 255<br>
Lo1 : 240.977<br>
La1 : 23.117<br>
<br>
float gridlon_0 ( ygrid_0, xgrid_0 )<br>
corners : ( -119.023, -80.74672, -59.95137,
-134.043 )<br>
long_name : longitude<br>
grid_type : Polar Sterographic Projection (North
or South)<br>
units : degrees_east<br>
Dy : 4.763<br>
Dx : 4.763<br>
Lov : 255<br>
Lo1 : 240.977<br>
La1 : 23.117 <br>
</div>
<div><br>
</div>
<div><b>Grid 2:</b> </div>
<div> float gridlat_0 ( ygrid_0, xgrid_0 )<br>
corners : ( 20.192, 20.33182, 54.37271, 54.20572
)<br>
long_name : latitude<br>
grid_type : Lambert Conformal (secant, tangent,
conical or bipolar)<br>
units : degrees_north<br>
Latin2 : 25<br>
Latin1 : 25<br>
Dy : 2.539702<br>
Dx : 2.539703<br>
Lov : 265<br>
Lo1 : 238.446<br>
La1 : 20.192<br>
<br>
float gridlon_0 ( ygrid_0, xgrid_0 )<br>
corners : ( -121.554, -69.20839, -59.04248,
-131.9951 )<br>
long_name : longitude<br>
grid_type : Lambert Conformal (secant, tangent,
conical or bipolar)<br>
units : degrees_east<br>
Latin2 : 25<br>
Latin1 : 25<br>
Dy : 2.539702<br>
Dx : 2.539703<br>
Lov : 265<br>
Lo1 : 238.446<br>
La1 : 20.192<br>
</div>
</div>
<br>
<fieldset></fieldset>
<pre>_______________________________________________
ncl-talk mailing list
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>
List instructions, subscriber options, unsubscribe:
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></pre>
</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></blockquote></div>
</blockquote></div>