<div dir="ltr">Hi Dz,<div>Just a guess: When interpolating make sure that the ranges of the input latitude/longitude arrays match the ranges of the output latitude/longitude arrays. In other words, if the longitudes for the input arrays (lon2d) run from -180:180E and the longitudes for the regridded output (lonnew) run from 0:360 you can get results like what you describe. </div><div><br></div><div>Hope that helps. If you continue to have issues let the ncl-talk email list know. </div><div>Adam</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Sep 18, 2017 at 12:04 AM, dale zuri <span dir="ltr"><<a href="mailto:dalezuri@gmail.com" target="_blank">dalezuri@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hi,<br></div>After interpolating the irregular grid to regular, the output plot reflects the missing values for most of the region. I attached here an output figure.<br></div>I would appreciate any help and suggestions to fix this issue. <br><br></div>Thanks<br></div>Dz<br><br><br>begin<br><br><br><br> iz = addfile("<a href="http://landsea2_wrf.nc" target="_blank">landsea2_wrf.nc</a>","r")<br> lat2d = iz->lat<br> lon2d = iz->lon<br> mas=iz->mask<br><br> c = addfile("<a href="http://era-down_1989_apr-sep_22z.nc" target="_blank">era-down_1989_apr-<wbr>sep_22z.nc</a>","r")<br> t2 = cds->p<br> <br> t2!1="lon2d"<br> t2!2="lat2d"<br> t2@lon2d = lon2d<br> t2@lat2d = lat2d<br>printVarSummary(t2)<br>aavg= t2(time|:,lat2d|:,lon2d|:)<br>printVarSummary(aavg)<br>print(aavg(0,:,0))<br>printMinMax(aavg,0)<br>printMinMax(lat2d,0)<br>printMinMax(lon2d,0)<br>;*****************************<wbr>*************************<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr>;;;;;;;;;;;;;;;;;;;;;;;<br>f1 = addfile("/import/h/<wbr>marchforecast/mar1982-2010/<a href="http://simpleWunit.2005031212.nc" target="_blank">sim<wbr>pleWunit.2005031212.nc</a>", "r") ; note the "s" of addfile<br>T1 = f1->precSum ; read T from all files<br>printVarSummary (T1)<br>latnew=f1->lat_0<br>lonnew=f1->lon_0<br>;;;;;;Interpolation Step;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr>;;;<br>;printVarSummary(T)<br>;printVarSummary(aavg)<br>;print(aavg)<br>printVarSummary(aavg)<br>newsst = rcm2rgrid_Wrap(lat2d,lon2d,<wbr>aavg,latnew,lonnew,1)<br>printMinMax(newsst,0)<br>printVarSummary(newsst)<br>print(newsst)<br><br><br> system("/bin/rm -f <a href="http://simple.nc" target="_blank">simple.nc</a>") ; remove any pre-existing file<br> ncdf = addfile("<a href="http://simple.nc" target="_blank">simple.nc</a>" ,"c") ; open output netCDF file<br><br> ;=============================<wbr>==============================<wbr>========<br> ; create global attributes of the file (optional)<br> ;=============================<wbr>==============================<wbr>========<br> fAtt = True ; assign file attributes<br> fAtt@title = "NCL Simple Approach to netCDF Creation"<br> fAtt@source_file = "<a href="http://original-file.nc" target="_blank">original-file.nc</a>"<br> fAtt@Conventions = "None"<br> fAtt@creation_date = systemfunc ("date")<br> fileattdef( ncdf, fAtt ) ; copy file attributes<br><br> ;=============================<wbr>==============================<wbr>========<br> ; make time an UNLIMITED dimension; recommended for most applications<br> ;=============================<wbr>==============================<wbr>========<br> filedimdef(ncdf,"time",-1,<wbr>True)<br><br> ;=============================<wbr>==============================<wbr>========<br> ; output variables directly; NCL will call appropriate functions<br> ; to write the meta data associated with each variable<br> ;=============================<wbr>==============================<wbr>========<br>ncdf->newsst = newsst<br><br>end<br> <br></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><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>