<div dir="ltr"><div>[1] No</div><div>[2] Yes<br></div><div><br></div><div><pre> nlat = 720<br> latitude = <b>fspan</b>(-49.875, 49.875, nlat)<br> latitude@units = "<b>degrees_north</b>"<br> latitude@long_name = "latitude"<br> latitude!0 = "latitude"<br> latitude&latitude = latitude<br></pre><pre> printVarSummary(latitude)<br> printMinMax(latitude,0)<br><br></pre><pre> mlon = 1440
longitude = <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/lonGlobeFo.shtml" target="_blank"><b>lonGlobeFo</b></a>(mlon, "longitude", "longitude", "degrees_east") ; 0->360<br> longitude = (/ longitude - 180. /) ; subtract 180 from all values
longitude&longitude = longitude ; update coordinates printVarSummary(lon)<br> <br> printVarSummary(longitude)<br> printMinMax(longitude,0)<br><br><pre> f = <b>addfile</b>("NPR_SMOPS_CMAP_D20170<a href="http://410.nc" target="_blank"><wbr>410.nc</a>","r")<br> x = <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml" target="_blank"><b>short2flt</b></a>( f->Blended_SM )<br></pre></pre><pre> x&latitude = latitude<br> x&longitude = longitude <br></pre><pre> printVarSummary(x)<br><br>[3]/[4]<br></pre><pre> xo = .... longitudes of land points...<br></pre><pre> yo = ... latitudes ....<br></pre><pre> xpts = <strong>linint2_points</strong>(x&longitude,x&latitude,x, False, xo,yo, 0)<br></pre><pre> printVarSummary(xpts)<br></pre><pre> printMinMax(xpts,0)<br></pre><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 4, 2018 at 5:50 AM, Ehsan Taghizadeh <span dir="ltr"><<a href="mailto:ehsantaghizadeh@yahoo.com" target="_blank">ehsantaghizadeh@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:bookman old style,new york,times,serif;font-size:16px"><div style="font-size:16px"><div style="font-size:16px"><div id="m_-7665546919046516254ydp2d0caccyiv1638557123"><div style="font-size:16px"><div style="font-size:16px"><div id="m_-7665546919046516254ydp2d0caccyiv1638557123ydp541e2203yiv1133628391"><div><div style="font-family:"bookman old style","new york",times,serif"><span style="font-family:"times new roman","new york",times,serif;font-size:medium">Thank you so much for your nice reply, as always.</span><br></div><div style="font-family:"bookman old style","new york",times,serif"><font size="3" face="times new roman, new york, times, serif">I did your comments, however I encountered some difficulties and I'll be thankful if I could have any help about them, too.</font></div><div style="font-family:"bookman old style","new york",times,serif"><font size="3" face="times new roman, new york, times, serif"><br clear="none"></font></div><div style="font-family:"bookman old style","new york",times,serif"><font size="3" face="times new roman, new york, times, serif">1. Should I change</font></div><span class=""><div style="font-family:"bookman old style","new york",times,serif"><font size="3" face="courier new, courier, monaco, monospace, sans-serif"><span style="background-color:rgb(253,248,105)"></span></font><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><font style="background-color:rgb(253,248,105)" size="3" face="courier new, courier, monaco, monospace, sans-serif">latitude = <b>fspan</b>(-49.875, 49.875, nlat)</font></pre></div></span><div><font size="3" face="times new roman, new york, times, serif">to</font></div><div><font style="font-family:"bookman old style","new york",times,serif;font-size:16px" size="2" face="times new roman, new york, times, serif"><span></span></font><pre style="font-family:"bookman old style","new york",times,serif;white-space:pre-wrap;color:rgb(0,0,0)"><font style="background-color:rgb(253,248,105)" size="3" face="courier new, courier, monaco, monospace, sans-serif">latitude = <b>fspan</b>(-89.875, 89.875, nlat)</font></pre><div><font size="3" face="times new roman, new york, times, serif">or like longitude, use</font></div><div style="font-family:"bookman old style","new york",times,serif"><span></span><pre style="margin-top:0px;margin-bottom:0px;padding:0px;line-height:12pt;font-family:courier;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0)"><span style="background-color:rgb(253,248,105)"><font size="3">lat = <strong style="margin:0px;padding:0px">latGlobeFo</strong>(nlat, "lat", "latitude", "degrees_north")</font></span></pre><br clear="none"></div><font size="3" face="times new roman, new york, times, serif">2. Should I change</font></div><div style="font-size:16px"><font style="font-family:"bookman old style","new york",times,serif" size="2" face="times new roman, new york, times, serif"><span style="background-color:rgb(253,248,105)"></span></font><pre style="font-family:"bookman old style","new york",times,serif;white-space:pre-wrap"><font style="background-color:rgb(253,248,105)" face="courier new, courier, monaco, monospace, sans-serif">latitude@units = "degrees_east"</font></pre><span style="background-color:rgb(255,255,255)"><font face="times new roman, new york, times, serif">to</font></span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><font size="2" face="times new roman, new york, times, serif"><span></span></font><pre style="white-space:pre-wrap"><font style="background-color:rgb(253,248,105)" face="courier new, courier, monaco, monospace, sans-serif">latitude@units = "degrees_north"</font></pre><br clear="none"></div><div><font size="3" face="times new roman, new york, times, serif"><span></span></font><div><font size="3" face="times new roman, new york, times, serif">3. How should I control FillValue (-999s), because after running ncl script, it writes -0.1 for Nan value of the data, both in grided and interpolated output files.</font></div><div><font size="3" face="times new roman, new york, times, serif">4. However the most difficulty is the interpolating part, because all interpolated values for all station points got -0.1, means Nan. There is no value except -0.1 for any station points. Could anybody help me about this?</font></div><div><font size="3" face="times new roman, new york, times, serif">I've attached the output of interpolated file. I'll be thankful for any help.</font></div><div><font size="3" face="times new roman, new york, times, serif"><br></font></div><div><font face="times new roman, new york, times, serif">Sincerely<br></font></div></div><div><font size="3" face="times new roman, new york, times, serif">Ehsan</font></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><br clear="none"></div></div><div class="m_-7665546919046516254ydp2d0caccyiv1638557123ydp541e2203yiv1133628391yahoo_quoted" id="m_-7665546919046516254ydp2d0caccyiv1638557123ydp541e2203yiv1133628391yahoo_quoted_6583562785" style="font-family:"bookman old style","new york",times,serif;font-size:16px">
<div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a">
</div></div></div></div><div class="m_-7665546919046516254ydp2d0caccyiv1638557123yqt3572844941" id="m_-7665546919046516254ydp2d0caccyiv1638557123yqtfd53451" style="font-family:"bookman old style","new york",times,serif"></div></div></div></div></div><div class="m_-7665546919046516254yiv1638557123yqt8965435499" id="m_-7665546919046516254yiv1638557123yqtfd57951"><div><div class="h5"><div>
On Friday, August 31, 2018, 6:14:53 PM GMT+4:30, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:
</div>
<div><br clear="none"></div>
<div><br clear="none"></div>
</div></div><div><div id="m_-7665546919046516254yiv1638557123"><div><div><div class="h5"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Sometimes, you have to read documentation associated with the file(s)<br clear="none"><span class="m_-7665546919046516254yiv1638557123gmail-st"><em>Soil Moisture</em> Operational Products System (<b><em>SMOPS</em></b>)</span><br clear="none"><br clear="none">The file has:<br clear="none"><br clear="none">netcdf NPR_SMOPS_CMAP_D20170410 {<br clear="none">dimensions:<br clear="none"><b> Longitude = 1440 ;<br clear="none"> Latitude = 720 ;</b><br clear="none">variables:<br clear="none"> <b>short</b> Blended_SM(<b>Latitude, Longitude</b>) ;<br clear="none"> Blended_SM:long_name = "Blended Soil Moisture" ;<br clear="none"> Blended_SM:units = "m^3/m^3" ;<br clear="none"> Blended_SM:FillValue = -999s ;<br clear="none"> Blended_SM:valid_range = 0s, 10000s ;<br clear="none"> Blended_SM:<b>scale_factor </b>= 0.0001f ;<br clear="none"> Blended_SM:<b>add_offset</b> = 0.f ;<br clear="none"><br clear="none"></div><div>[SNIP]<br clear="none"> :Product_Resolution = "0.25 degree" ;<br clear="none"> :Date_Start = "20170410" ;<br clear="none"> :Date_End = "20170410" ;<br clear="none"></div><div>[SNIP]<br clear="none"></div><div> ====<br clear="none"></div><div>It does not have Latitude[*] and Longitude[*] coordinates. However, <b>*you*</b> can create them. Then, use ESMF or (say)<a rel="nofollow" shape="rect" href="https://www.ncl.ucar.edu/Document/Functions/Contributed/linint2_Wrap.shtml" target="_blank"><b> linint2_Wrap </b></a>or <a rel="nofollow" shape="rect" href="https://www.ncl.ucar.edu/Document/Functions/Contributed/linint2_points_Wrap.shtml" target="_blank"><b>linint2_points_Wrap</b></a> to perform the interpolation.<br clear="none"><br clear="none">====<br clear="none"><pre> nlat = 720<br clear="none"> latitude = <b>fspan</b>(-49.875, 49.875, nlat)<br clear="none"> latitude@units = "degrees_east"<br clear="none"> latitude@long_name = "latitude"<br clear="none"> latitude!0 = "latitude"<br clear="none"> latitude&latitude = latitude<br clear="none"></pre><pre> printVarSummary(latitude)<br clear="none"> printMinMax(latitude,0)<br clear="none"><br clear="none"></pre><pre> mlon = 1440
longitude = <a rel="nofollow" shape="rect" href="https://www.ncl.ucar.edu/Document/Functions/Contributed/lonGlobeFo.shtml" target="_blank"><b>lonGlobeFo</b></a>(mlon, "longitude", "longitude", "degrees_east") ; 0->360<br clear="none"> longitude = (/ longitude - 180. /) ; subtract 180 from all values
longitude&longitude = longitude ; update coordinates printVarSummary(lon)<br clear="none"> <br clear="none"> printVarSummary(longitude)<br clear="none"> printMinMax(longitude,0)<br clear="none"><br clear="none"></pre><pre> f = <b>addfile</b>("NPR_SMOPS_CMAP_D20170<a rel="nofollow" shape="rect" href="http://410.nc" target="_blank"><wbr>410.nc</a>","r")<br clear="none"> x = <a rel="nofollow" shape="rect" href="https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml" target="_blank"><b>short2flt</b></a>( f->Blended_SM )<br clear="none"></pre><pre> x&latitude = latitude<br clear="none"> x&longitude = longitude <br clear="none"></pre><pre> printVarSummary(x)<br clear="none"></pre><pre> <br clear="none"></pre></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><div class="m_-7665546919046516254yiv1638557123gmail_extra"><br clear="none"><div class="m_-7665546919046516254yiv1638557123gmail_quote"><div><div class="h5">On Fri, Aug 31, 2018 at 5:45 AM, Ehsan Taghizadeh <span dir="ltr"><<a rel="nofollow" shape="rect" href="mailto:ehsantaghizadeh@yahoo.com" target="_blank">ehsantaghizadeh@yahoo.com</a>></span> wrote:<br clear="none"></div></div><blockquote class="m_-7665546919046516254yiv1638557123gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="m_-7665546919046516254yiv1638557123yqt7432983737" id="m_-7665546919046516254yiv1638557123yqt23358"><div><div style="font-family:bookman old style,new york,times,serif;font-size:16px"><div>Hi,</div><div>Until now I've asked some questions related to interpolating some dataset on some unstructured points (synoptic stations). Now I have a new data set which I'm not sure about appropriate method to interpolating. I put one of my data (<span>NPR_SMOPS_CMAP_D20170410.nc</span>) in <span style="font-family:courier;font-size:13.3333px;white-space:pre-wrap"><a rel="nofollow" shape="rect" href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a>. <span><span style="color:rgb(0,0,0)">I'll be thankful if some one could refer me to correct script, like on (</span></span></span><a rel="nofollow" shape="rect" href="https://www.ncl.ucar.edu/Applications/ESMF.shtml" target="_blank">https://www.ncl.ucar.edu/ Applications/ESMF.shtml</a><span style="font-family:courier;font-size:13.3333px;white-space:pre-wrap"><span><span style="color:rgb(0,0,0)">) or (</span></span></span><a rel="nofollow" shape="rect" href="https://www.ncl.ucar.edu/Applications/regrid.shtml" target="_blank">https://www.ncl.ucar.edu/ Applications/regrid.shtml</a><span style="font-family:courier;font-size:13.3333px;white-space:pre-wrap"><span><span style="color:rgb(0,0,0)">), which interpolates my data on some special unstructured points.</span></span></span></div><span class=""><div><span style="font-family:courier;font-size:13.3333px;white-space:pre-wrap"><span><span style="color:rgb(0,0,0)">However the printVarSummary(<span>Blended_SM</span>) result is as below:</span></span></span></div><div><span style="font-family:courier;font-size:13.3333px;white-space:pre-wrap"><span><span style="color:rgb(0,0,0)"><br clear="none"></span></span></span></div><div><span style="font-size:13.3333px;white-space:pre-wrap"><span style="color:rgb(0,0,0);font-size:16px"></span></span><div><font face="times new roman, new york, times, serif">Variable: <span style="color:rgb(0,0,0);font-size:16px;white-space:pre-wrap">Blended_SM</span></font></div><div><font face="times new roman, new york, times, serif">Type: float</font></div><div><font face="times new roman, new york, times, serif">Total Size: 4147200 bytes</font></div><div><font face="times new roman, new york, times, serif"> 1036800 values</font></div><div><font face="times new roman, new york, times, serif">Number of Dimensions: 2</font></div><div><font face="times new roman, new york, times, serif">Dimensions and sizes: [Latitude | 720] x [Longitude | 1440]</font></div><div><font face="times new roman, new york, times, serif">Coordinates:</font></div><div><font face="times new roman, new york, times, serif">Number Of Attributes: 5</font></div><div><font face="times new roman, new york, times, serif"> _FillValue : 1e+20</font></div><div><font face="times new roman, new york, times, serif"> long_name : Blended Soil Moisture</font></div><div><font face="times new roman, new york, times, serif"> units : m^3/m^3</font></div><div><font face="times new roman, new york, times, serif"> FillValue : -999</font></div><div><font face="times new roman, new york, times, serif"> valid_range : ( 0, 1 )</font></div><div><br clear="none"></div><br clear="none"></div><div><span style="white-space:pre-wrap">Sincerely</span><span class="m_-7665546919046516254yiv1638557123HOEnZb"><font color="#888888"><br clear="none"></font></span></div><span class="m_-7665546919046516254yiv1638557123HOEnZb"><font color="#888888"></font></span><div><span style="font-family:courier;font-size:13.3333px;white-space:pre-wrap"><span><span style="color:rgb(0,0,0)">Ehsan</span></span></span></div></span></div></div></div><span class=""><br clear="none">______________________________ _________________<br clear="none">
ncl-talk mailing list<br clear="none">
<a rel="nofollow" shape="rect" href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br clear="none">
List instructions, subscriber options, unsubscribe:<br clear="none">
<a rel="nofollow" shape="rect" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/ mailman/listinfo/ncl-talk</a><br clear="none">
<br clear="none"></span></blockquote></div><br clear="none"></div></div></div></div>
</div><div class="m_-7665546919046516254yiv1638557123yqt8965435499" id="m_-7665546919046516254yiv1638557123yqtfd12555"><div class="m_-7665546919046516254yiv1638557123yqt3572844941" id="m_-7665546919046516254yiv1638557123yqtfd22236">
</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>