<div dir="ltr"><div><div>I have looked at this using NCL 6.3.0 (MAC OS)<br><br>From: PET0.RegridWeightGen.Log<br><br>The error is "<font size="2">Concave Element Detected"<br><br></font></div><font size="2">***Someone else will have to suggest a work around.***<br></font></div><font size="2"></font><div><font size="2">---<br></font><div><font size="2"><br>20170328 173330.778 INFO PET0 Running with ESMF Version 6.3.0rp1<br>20170328 173734.737 INFO PET0 Running with ESMF Version 6.3.0rp1<br>20170328 174808.719 ERROR ~~~~~~~~~~~~~~~~~~ Concave Element Detected ~~~~~~~~~~~~~~~~~~<br>20170328 174808.723 ERROR PET0 concave elem. id=24017<br>20170328 174808.723 ERROR PET0 <br>20170328 174808.723 ERROR PET0 concave elem. coords (lon [-180 to 180], lat [-90 to 90]) (x,y,z)<br>20170328 174808.723 ERROR PET0 -----------------------------------------------------------------<br>20170328 174808.723 ERROR PET0 0 (179.815179, 42.405261) (-0.738390, 0.002382, 0.674370)<br>20170328 174808.723 ERROR PET0 1 (179.627955, 42.541923) (-0.736767, 0.004784, 0.676129)<br>20170328 174808.723 ERROR PET0 2 (89.813369, 42.680222) (0.002395, 0.735145, 0.677906)<br>20170328 174808.723 ERROR PET0 3 (-89.999394, 42.542971) (0.000008, -0.736770, 0.676143)<br>20170328 174808.723 ERROR PET0 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br>20170328 174808.723 ERROR PET0 ESMCI_Regrid_F.C:149 c_esmc_regrid_create() Invalid argument - Src contains a concave cell<br>20170328 174808.727 ERROR PET0 ESMCI_Regrid_F.C:434 c_esmc_regrid_create() Invalid argument Internal subroutine call returned Error<br>20170328 174808.727 ERROR PET0 ESMF_Regrid.F90:321 ESMF_RegridStore Invalid argument - Internal subroutine call returned Error<br>20170328 174808.727 ERROR PET0 ESMF_FieldRegrid.F90:1008 ESMF_FieldRegridStoreNX Invalid argument - Internal subroutine call returned Error<br>20170328 174808.727 ERROR PET0 ESMF_RegridWeightGen.F90:1152 ESMF_RegridWeightGenFile Invalid argument - Internal subroutine call returned <br><br>=======<br>/glade/p/work/shea/ESMF/CDR_NSIDC>dir<br><br>-rw-r--r-- 1 shea ncar 3170894556 Mar 28 19:12 <a href="http://destination_grid_file.nc">destination_grid_file.nc</a><br>-rw-r--r-- 1 shea ncar 4277 Mar 28 17:55 g02202_Regrid_v01.ncl<br>-rw-r--r-- 1 shea ncar 1868 Mar 28 19:12 g02202_Regrid_v01.ncl_orig<br>-rw-r--r-- 1 shea ncar 528585368 Mar 28 19:24 <a href="http://masiemask_ims4km.nc">masiemask_ims4km.nc</a><br>-rw-r--r-- 1 shea ncar 7033 Mar 28 18:44 OUT_g02202.txt<br>-rw-r--r-- 1 shea ncar 1849 Mar 28 17:59 PET0.RegridWeightGen.Log<br>-rw-r--r-- 1 shea ncar 3738850 Mar 28 19:24 <a href="http://seaice_conc_daily_nh_f13_20060101_v02r00.nc">seaice_conc_daily_nh_f13_20060101_v02r00.nc</a><br>-rw-r--r-- 1 shea ncar 11440924 Mar 28 19:12 <a href="http://source_grid_file.nc">source_grid_file.nc</a><br></font></div></div><div class="gmail_extra"><br>=========<br><div class="gmail_quote">On Tue, Mar 28, 2017 at 2:47 PM, Greg Deemer <span dir="ltr"><<a target="_blank" href="mailto:greg.deemer@nsidc.org">greg.deemer@nsidc.org</a>></span> wrote:<br><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote"><div style="word-wrap:break-word;color:rgb(0,0,0);font-size:14px;font-family:calibri,sans-serif"><div>Hello, </div><div><br></div><div>I’m attempting to regrid sea-ice concentrations from a coarse spatial resolution onto a higher-resolution grid. Both grids are NSIDC Sea Ice Polar Stereographic North with a WGS84 reference datum. The coarse product (<a target="_blank" href="https://nsidc.org/data/g02202">https://nsidc.org/data/g02202</a><wbr>) is 448X304, or a 25km nominal resolution. The higher-resolution product (<a target="_blank" href="http://nsidc.org/data/docs/noaa/g02186_masie/index.html">http://nsidc.org/data/docs/<wbr>noaa/g02186_masie/index.html</a>) is 6144x6144, or a 4km nominal resolution.</div><div><br></div><div>The grid for the high-res product is a little funky in that it extends well beyond the equator, from the polar view (for added context, it is the 4km variant of the 24km grid handled here: <a target="_blank" href="http://www.ncl.ucar.edu/Applications/nic_ims.shtml">http://www.ncl.ucar.edu/<wbr>Applications/nic_ims.shtml</a>). I’m wondering if this is causing some error in the ESMF_regrid routine, which is throwing two errors:</div><div><br></div><div>1. nonsensical values to the minimum Extlon2d and minimum ExtGridCenter_Lon, as seen in the processing log</div><div><br></div><div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>get_src_grid_info: source lat dims = (448,304)</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>get_src_grid_info: source lon dims = (448,304)</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>get_src_grid_info: source grid type is 'curvilinear'</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>curvilinear_to_SCRIP: calculating grid corners...</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>curvilinear_to_SCRIP: no lat values are at the poles, so</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> calculating grid corners using</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> calc_SCRIP_corners_<wbr>noboundaries...</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>calc_SCRIP_corners_<wbr>noboundaries</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max original lat: 31.1026717524309/89.<wbr>8368159996152</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max original lon: -180/179.813975395493</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>calc_SCRIP_corners_<wbr>noboundaries</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max Extlat2d: 30.8580632088782/89.<wbr>8368159996152</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max Extlon2d: -539.812759557033/180.<wbr>001207917857</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>calc_SCRIP_corners_<wbr>noboundaries</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max ExtGridCenter_lat: 30.98047461780758/89.<wbr>8368159996152</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max ExtGridCenter_lon: -269.812759557033/179.<wbr>8151793919093</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>=====> CPU Elapsed Time: curvilinear_to_SCRIP: 0.109368 seconds <=====</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>get_dst_grid_info: destination lat dims = (6144,6144)</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>get_dst_grid_info: destination lon dims = (6144,6144)</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>curvilinear_to_SCRIP: calculating grid corners...</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>curvilinear_to_SCRIP: no lat values are at the poles, so</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> calculating grid corners using</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> calc_SCRIP_corners_<wbr>noboundaries...</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>calc_SCRIP_corners_<wbr>noboundaries</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max original lat: -21.4848/89.9729</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max original lon: -180/180</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>calc_SCRIP_corners_<wbr>noboundaries</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max Extlat2d: -21.5023/89.9729</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max Extlon2d: -539.995/180</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span>calc_SCRIP_corners_<wbr>noboundaries</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max ExtGridCenter_lat: -21.4935/89.9729</div><div>(0)<span style="white-space:pre-wrap" class="gmail-m_2410564885733441308Apple-tab-span"> </span> min/max ExtGridCenter_lon: -269.988/179.989</div></div><div><br></div><div>And a second error:</div><div>2. NetCDF: One or more variable sizes violate format constraints</div><div><br></div><div>I tried converting longitudes to monotonic, but the output log still spit out large negative values (e.g. <<0). I’ve also added Opt@srcLargeFile and Opt@dstLargFile lines within the code as I see the resulting ‘<a target="_blank" href="http://destination-grid-file.nc">destination-grid-file.nc</a>’ is on the order of 3.15Gb in size.</div><div><br></div><div>Please find script attached. The data files in use for this script can be found here:</div><div>Coarse: <a target="_blank" href="ftp://sidads.colorado.edu/pub/DATASETS/NOAA/G02202_v2/north/daily/2006/seaice_conc_daily_nh_f13_20060101_v02r00.nc">ftp://sidads.colorado.<wbr>edu/pub/DATASETS/NOAA/G02202_<wbr>v2/north/daily/2006/seaice_<wbr>conc_daily_nh_f13_20060101_<wbr>v02r00.nc</a></div><div>Fine: <a target="_blank" href="ftp://sidads.colorado.edu/DATASETS/NOAA/G02186/ancillary/masiemask_ims4km.nc">ftp://sidads.colorado.<wbr>edu/DATASETS/NOAA/G02186/<wbr>ancillary/masiemask_ims4km.nc</a></div><div><br></div><div>Thank you for the attention to this question and for any help.</div><span class="gmail-HOEnZb"><font color="#888888"><div><br></div><div><div id="gmail-m_2410564885733441308MAC_OUTLOOK_SIGNATURE"><div><div>-- </div><div><br></div></div><div><div style="font-family:-webkit-standard"><b>Greg Deemer</b></div><div style="font-family:-webkit-standard">User Services Office</div><div style="font-family:-webkit-standard">National Snow and Ice Data Center</div><div style="font-family:-webkit-standard"><a target="_blank" value="+13037358068" href="tel:%28303%29%20735-8068">(303) 735-8068</a></div><div style="font-family:-webkit-standard"><a target="_blank" href="http://nsidc.org/">http://nsidc.org</a></div><div style="font-family:-webkit-standard">Twitter: @NSIDC</div></div></div></div></font></span></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 target="_blank" rel="noreferrer" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div>