<div dir="ltr"><div class="gmail_default" style="font-size:small">In order to properly debug this, I would need to have your script and data. Can you provide this?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You can email me offline about this if you need to.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Thanks,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jan 23, 2017 at 4:57 PM, 林祥 <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"><div>Mary,</div><div><br></div><div> Thank you for the reply! I’ve tried uncommenting the <span style="font-size:small">Opt@SrcGridMask before with the following 2 lines</span></div><span class=""><div><span style="font-size:small"><br></span></div><div><span style="font-size:small"> </span> Opt@SrcGridMask = where(.not.ismissing(y(0,:,:))<wbr>,1,0)</div></span><div> Opt@DstGridMask = where(.not.ismissing(hsic(0,:,<wbr>:)),1,0)</div><div><br></div><div> and the output is as the attached plot file. You can see that the land grid have been masked, but the ocean grids after regridding</div><div><br></div><div> are still inconsistent with the original source grids ( the grids outside 90W-90E are different with the source grids) . </div><div><br></div><div> Besides, I found same problem occurred when I use the “bilinear” and “patch” method. But the combination of using “neareststod”</div><div><br></div><div>method and uncommenting the above 2 lines seem produce the correct results ( see anther attached files). But I still feel puzzled.</div><div><br></div><div><br></div><div>Thanks</div><div><br></div><div><br></div><div>Lin</div><div><br></div><div><br></div><div><br></div><div><img id="m_-7250572314183157786EF3D746F-BB4C-418A-A41F-F7655E1A3455" height="384" width="592" src="cid:1A4F7A64-449F-421B-94F9-F409EDAC944F"></div><div><img id="m_-7250572314183157786F8DEE1DD-F35D-496B-87B2-225D5A8B216C" height="384" width="592" src="cid:D90F13F8-446E-4FB4-9491-BEDC05D558A6"></div><div><div class="h5"><div><br></div><div><br></div><div><br></div><br><div><blockquote type="cite"><div>在 2017年1月24日，02:09，Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> 写道：</div><br class="m_-7250572314183157786Apple-interchange-newline"><div><div dir="ltr" style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><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="m_-7250572314183157786gmail-p1"><span class="m_-7250572314183157786gmail-s1"> </span><span class="m_-7250572314183157786gmail-s2">;Opt@SrcGridMask = where(.not.ismissing(y(0,:,:))<wbr>,1,0) </span></p><p class="m_-7250572314183157786gmail-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="m_-7250572314183157786gmail-p1">--Mary</p><p class="m_-7250572314183157786gmail-p1"><br></p></div></div><div class="gmail_extra" style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><br><div class="gmail_quote">On Mon, Jan 23, 2017 at 12:49 AM, 林祥<span class="m_-7250572314183157786Apple-converted-space"> </span><span dir="ltr"><<a href="mailto:xianglin72@icloud.com" target="_blank">xianglin72@icloud.com</a>></span><span class="m_-7250572314183157786Apple-converted-space"> </span>wro<wbr>te:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style: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<span class="m_-7250572314183157786Apple-converted-space"> </span><a href="ftp://sidads.colorado.edu/pub/DATASETS/NOAA/G02202_v2" target="_blank">ftp://sidads.colorado.<wbr>edu/pub/DATASETS/NOAA/G02202_<wbr>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/h<wbr>adobs/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> <span class="m_-7250572314183157786Apple-converted-space"> </span>or reset the longitudes to 0~360E using the following lines:</div><div> </div><div><div> <span class="m_-7250572314183157786Apple-converted-space"> </span>lon = where(lon.lt.0,lon+360,lon) ; -180 ~ 180 longitude order</div><div> <span class="m_-7250572314183157786Apple-converted-space"> </span>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" target="_blank">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/mailma<wbr>n/listinfo/ncl-talk</a></blockquote></div></div></div></blockquote></div><br></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>