<div dir="ltr"><div class="gmail_default" style="font-size:small">Lin,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I think you need to set your SrcGridMask option. I noticed you had it commented out:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">
<p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">;Opt@SrcGridMask = where(.not.ismissing(y(0,:,:)),1,0) </span></p>
<p class="gmail-p1">If y contains missing values, then you need to set SrcGridMask. If y contains missing values in the same location for every leftmost element of y (assuming it's 3D), then what you have above should work, but you need to uncomment it.</p><p class="gmail-p1">--Mary</p><p class="gmail-p1"><br></p></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jan 23, 2017 at 12:49 AM, 林祥 <span dir="ltr"><<a href="mailto:xianglin72@icloud.com" target="_blank">xianglin72@icloud.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi,all<div><br></div><div><br></div><div> I am regriding some NSIDC sic data (<a href="http://nsidc.org/data/g02202.html" target="_blank">http://nsidc.org/data/g02202.<wbr>html</a>, downloaded from <a href="ftp://sidads.colorado.edu/pub/DATASETS/NOAA/G02202_v2" target="_blank">ftp://sidads.colorado.edu/pub/<wbr>DATASETS/NOAA/G02202_v2</a>).</div><div><br></div><div>The source data is curvilinear grids with 2-d longitude coordinate array ranging between -180E~180E， and the destination grids is HadISST 2.2 (<a href="http://www.metoffice.gov.uk/hadobs/hadisst2/" target="_blank">http://www.metoffice.gov.uk/<wbr>hadobs/hadisst2/</a>)</div><div><br></div><div>which is rectilinear grids ranging from 0E to 360E. </div><div><br></div><div> I try to use the original longitude of source grids with the following :</div><div><div><br></div><div> fi = addfile(data_path + data_filename,"r")</div><div> lon = fi->longitude</div><div> lat = fi->latitude </div></div><div> </div><div> y@lat2d = lat</div><div><div> y@lon2d = lon </div></div><div><br></div><div> or reset the longitudes to 0~360E using the following lines:</div><div> </div><div><div> lon = where(lon.lt.0,lon+360,lon) ; -180 ~ 180 longitude order</div><div> lon@units= "degrees_east"</div></div><div> </div><div>But I got the same result and the regridded data contain missing values outside the range of -90W~90E.Please see the attached plot file</div><div><br></div><div>I don’t know where the problem come from. And I also found no instruction on the NCL webpages explaining the similar case when the curvilinear</div><div><br></div><div>grids share different longitude ranges with the destinate grids. So anyone can show me a hand?</div><div><br></div><div><br></div><div>The script file, source data file , contour plot files with original data and regridded data have been attached.</div><div><br></div><div><br></div><div>Thanks!</div><div><br></div><div><br></div><div><br></div><div>Lin</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div></div></div><br><div style="word-wrap:break-word"><div></div><br><div></div></div><br><div style="word-wrap:break-word"><div></div><div><br></div><div><br></div><div></div></div><br><div style="word-wrap:break-word"><div></div><div> </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>