<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Laura,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The main issue is that your "data" array is ordered lon x lat, and not lat x lon. Also, the longitudes go from 0 to 360, while on the destination grid they go from -180 to 180.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">To fix both of these problems, I used transpose and lonFlip:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><font face="monospace, monospace">data = <b>transpose</b>(fSrc->precipitationCal)</font></div><div class="gmail_default"><font face="monospace, monospace">data = <b>lonFlip</b>(data)</font></div><div class="gmail_default" style="font-size:small"><font face="monospace, monospace">printVarSummary(data)</font></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The second issue is that you don't need to officially define your source grid, because "data" is a rectilinear grid and already has the lat/lon arrays attached to it. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">This means these five lines can be removed:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace">lonSrc = fSrc->lon</font></div><div class="gmail_default"><font face="monospace, monospace">latSrc = fSrc->lat</font></div><div class="gmail_default"><font face="monospace, monospace">opt@SrcGridType = "0.10x0.10" ; destination grid</font></div><div class="gmail_default"><font face="monospace, monospace">opt@SrcGridLon = lonSrc ; source grid lon input.</font></div><div class="gmail_default"><font face="monospace, monospace">opt@SrcGridLat = latSrc ; source grid lat input.</font></div><div class="gmail_default"><br></div><div class="gmail_default">With the destination grid, you are getting one fewer grid point in both the lat and lon directions because ESMF_regrid is trying to construct the destination lat/lon grid using this information:</div><div class="gmail_default"><div><font face="monospace, monospace"><br></font></div></div></div><div class="gmail_default" style="font-size:small"><font face="monospace, monospace">opt@DstGridType = "0.30x0.30" ; destination grid<br></font></div><div class="gmail_default"><div class="gmail_default" style="font-size:small"><font face="monospace, monospace">opt@DstLLCorner = (/-89.85d, 0.15d/)<br></font></div><div class="gmail_default" style="font-size:small"><font face="monospace, monospace">opt@DstURCorner = (/ 89.85d,359.85d/)</font></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you instead use the "lat" and "lon" arrays that you calculate yourself, then you will get the size you expect:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><div class="gmail_default" style="font-size:small"><font face="monospace, monospace">opt@DstGridLon = lon ; destination grid lon input.</font></div><div class="gmail_default" style="font-size:small"><font face="monospace, monospace">opt@DstGridLat = lat ; destination grid lat input.</font></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">See the attached modified script and two PNG images (from setting MAP_ZOOM to True or False)</div><div class="gmail_default" style="font-size:small"><br></div>Regarding your third question: you should be able to use the <a href="http://destination.GPM.to.0.30deg_by_0.30deg.nc">destination.GPM.to.0.30deg_by_0.30deg.nc</a> file that gets created. Just simply set opt@SkipDstGrid = True before you call ESMF_regrid. You will still need to generate the MPAS source file, and the weights file.</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><div class="gmail_default"><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 26, 2017 at 12:33 PM, Laura Fowler <span dir="ltr"><<a href="mailto:laura@ucar.edu" target="_blank">laura@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi:<br>
<br>
I have two questions regarding ESMF_regrid:<br>
<br>
1. I would like to use ESMF_regrid to regrid the 0.1 deg. GPM grid<br>
(3600 x 1800) to a 0.3 deg. grid (1200 x 600), but I get a core dump<br>
at the very end although the debug file writes "Completed weight<br>
generation successfully. I can also dump the content of some of the<br>
variables too.<br>
Can somebody help me understand why the core file?<br>
<br>
The ncl script is in /glade/p/work/laura/DATA/GPM (see<br>
regrid.GPM.to.0.3deg_0.3deg.<wbr>ncl).<br>
<br>
<br>
2. Looking at the destination file, the grid size is 718201 but I was<br>
expecting 720000 since I decreased the grid resolution by 3. Is there<br>
a way to have the destination grid size equal to 720000 (or 1200 x<br>
600). May that has to do with the fact that I got a core file in the<br>
end?<br>
<br>
<br>
3. Finally, I would like to regrid outputs from my MPAS runs so that<br>
my regridded outputs are on exactly the same grid as the regridded GPM<br>
outputs. Is there a way that I can use information from the GPM<br>
destination file?<br>
<br>
<br>
Thanks for your help.<br>
Laura<br>
<br>
<br>
--<br>
!-----------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------------<br>
Laura D. Fowler<br>
Mesoscale and Microscale Meteorology Division (MMM)<br>
National Center for Atmospheric Research<br>
P.O. Box 3000, Boulder CO 80307-3000<br>
<br>
e-mail: <a href="mailto:laura@ucar.edu">laura@ucar.edu</a><br>
phone: <a href="tel:303-497-1628" value="+13034971628">303-497-1628</a><br>
<br>
!-----------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------------<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>
</blockquote></div><br></div>