<div dir="ltr"><div><div>THX for sending a sample file.<br><br></div>See attached.<br><br>The computational part of your original scrip was 'ok':<br><br>   (a) rcm2rgrid  [better: rcm2rgrid_Wrap]<br>   (b) rgrid2rcm  [better: rgrid2rcm_Wrap]<br><br>However, as you found out ... it takes 'forever'. The underlying code is 'brute force' and not efficient.<br><br>---<br>I used ESMF.  I am sure you can follow the code.<br><br>It plots both the (a) rectilinear divergence and the (b)  WRF divergence.<br><br><div>The advantage of the ESMF is that weight files are generated:<br><br><a href="http://WRF_to_Rect.WgtFile_bilinear.700_x_850.nc" target="_blank">WRF_to_Rect.WgtFile_bilinear.<wbr>700_x_850.nc</a><br><br><a href="http://Rect_to_WRF.WgtFile_bilinear.700_x_850.nc" target="_blank">Rect_to_WRF.WgtFile_bilinear.<wbr>700_x_850.nc</a><br><br></div><div>These can be used on any variable that has the original WRF grid structure.<br><br></div><div>Then you can use:<br><br><a href="https://www.ncl.ucar.edu/Document/Functions/ESMF/ESMF_regrid_with_weights.shtml" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/ESMF/ESMF_<wbr>regrid_with_weights.shtml</a><br><br></div>The regridding will be 'instantaneous'<br><br>============<br>d04_div_50_1530.000001.png      regular contour plot<br>d04_div_50_1530.000002.png      One WRF grid with map<br><br></div>Please be aware that interpolating u and  v independently is 'technically' not correct. Because, the components have moved  to a new geographic location, there is rotation <br><br><pre>   Unew = U*cos(ANGLE) - V*sin(ANGLE)
   Vnew = U*sin(ANGLE) + V*cos(ANGLE)<br><br></pre><pre>Still, what was done is 'good enough'  ! <br></pre></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Apr 10, 2018 at 6:24 AM, Ty Buckingham <span dir="ltr"><<a href="mailto:ty.buckingham@manchester.ac.uk" target="_blank">ty.buckingham@manchester.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Hi there,
<div><br>
</div>
<div>Having a little trouble calculating divergence from my WRF output.</div>
<div><br>
</div>
<div>I'm transforming from the WRF curvilinear grid to a linear grid in order to calculate divergence using finite differences, then transforming back for my plots. </div>
<div><br>
</div>
<div>The WRF output I'm using covers a small area over the UK - however, the produced plot is somehow global. Not entirely sure how that's happened.</div>
<div><br>
</div>
<div>I've attached my sample code and an image of the output.</div>
<div><br>
</div>
<div>Thanks!</div><span class="HOEnZb"><font color="#888888">
<div>Ty</div>
</font></span></div>
</div>

<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>