<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>