<div dir="ltr"><div><div>[1] As noted in the documentation, <span style="font-weight:normal"><font size="2">linint2_points</font></span> Interpolates from a *rectilinear grid* to a user specified set of points using
bilinear interpolation.
<br><br>     <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/linint2_points.shtml">http://www.ncl.ucar.edu/Document/Functions/Built-in/linint2_points.shtml</a><br></div><div><br></div>[2] The SMAP data are on a *curvilinear* grid: [DIM_000 | 406] x [DIM_001 | 964] . Hence, 'linint2_points'  will not work.<br><br></div><div>[3] I suggest that you <br>     (a) use ESMF regridding to interpolate from the curvilinear grid to a rectilinear grid of approximately the same resolution.<br></div><div>     (b) then use linint2_points on (3a)<br><br>===<br></div><div>See the NARR Example 5. NARR uses a curvilinear grid. Create an SMAP ESMF weight file [eg: ESMF eamaple 30) ; Use the weight file to generate a rectilinear grid; use linint2_points.<br><br> <a href="http://www.ncl.ucar.edu/Applications/narr.shtml">http://www.ncl.ucar.edu/Applications/narr.shtml</a><br><br></div><div>Good luck<br></div><div>       <br></div><div><span style="background-color:rgb(173,215,115)"></span><br><span style="background-color:rgb(173,215,115)"></span></div><div><span style="background-color:rgb(173,215,115)"></span><span style="font-weight:normal"></span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jan 29, 2018 at 9:23 AM, Ehsan Taghizadeh via ncl-talk <span dir="ltr"><<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</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></div>
            <div><span style="color:rgb(38,40,42)">Hi,</span><br></div><div id="m_-7490938470158028763ydpdfb568aayahoo_quoted_7713408982" class="m_-7490938470158028763ydpdfb568aayahoo_quoted"><div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a"><div><div id="m_-7490938470158028763ydpdfb568aayiv8920612249"><div style="font-family:bookman old style,new york,times,serif;font-size:16px"><div class="m_-7490938470158028763ydpdfb568aayiv8920612249ydpbbc61b17yahoo_quoted" id="m_-7490938470158028763ydpdfb568aayiv8920612249ydpbbc61b17yahoo_quoted_8100818481"><div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a"><div><div id="m_-7490938470158028763ydpdfb568aayiv8920612249ydpbbc61b17yiv0430163723"><div style="font-family:bookman old style,new york,times,serif;font-size:16px"><div>I want to interpolate soil moisture data from SMAP (Soil Moisture Ative/Passive) to some specified locations (in stations_NW.csv attached file). Some information of data are as belows.</div><div><br></div><div><div>printVarSummary(lat2d)<br></div><div><div><span style="background-color:rgb(173,215,115)">Variable: lat2d</span></div><div><span style="background-color:rgb(173,215,115)">Type: float</span></div><div><span style="background-color:rgb(173,215,115)">Total Size: 1565536 bytes</span></div><div><span style="background-color:rgb(173,215,115)">            391384 values</span></div><div><span style="background-color:rgb(173,215,115)">Number of Dimensions: 2</span></div><div><span style="background-color:rgb(173,215,115)">Dimensions and sizes:   [DIM_000 | 406] x [DIM_001 | 964]</span></div><div><span style="background-color:rgb(173,215,115)">Coordinates:</span></div><div><span style="background-color:rgb(173,215,115)">Number Of Attributes: 2</span></div><div><span style="background-color:rgb(173,215,115)">  long_name :   Latitude of the center of the Earth based grid cell.</span></div><div><span style="background-color:rgb(173,215,115)">  units :       degrees_north</span></div><div><br></div></div><div><div>printVarSummary(lon2d)<br></div><div><span style="background-color:rgb(173,215,115)">Variable: lon2d</span></div><div><span style="background-color:rgb(173,215,115)">Type: float</span></div><div><span style="background-color:rgb(173,215,115)">Total Size: 1565536 bytes</span></div><div><span style="background-color:rgb(173,215,115)">            391384 values</span></div><div><span style="background-color:rgb(173,215,115)">Number of Dimensions: 2</span></div><div><span style="background-color:rgb(173,215,115)">Dimensions and sizes:   [DIM_000 | 406] x [DIM_001 | 964]</span></div><div><span style="background-color:rgb(173,215,115)">Coordinates:</span></div><div><span style="background-color:rgb(173,215,115)">Number Of Attributes: 2</span></div><div><span style="background-color:rgb(173,215,115)">  long_name :   Longitude of the center of the Earth based grid cell.</span></div><div><span style="background-color:rgb(173,215,115)">  units :       degrees_east</span></div><div><br></div></div><div>printVarSummary(var):<br></div><div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">Variable: var</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">Type: float</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">Total Size: 1565536 bytes</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">            391384 values</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">Number of Dimensions: 2</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">Dimensions and sizes:   [DIM_000 | 406] x [DIM_001 | 964]</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">Coordinates:</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">Number Of Attributes: 6</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">  _FillValue :  -9999</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">  coordinates : /Soil_Moisture_Retrieval_Data_<wbr>AM/latitude /Soil_Moisture_Retrieval_Data_<wbr>AM/longitude</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">  long_name :   Representative soil moisture measurement for the Earth based grid cell.</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">  units :       cm**3/cm**3</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">  valid_max :   0.5</span></div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><span style="background-color:rgb(173,215,115)">  valid_min :   0.02</span></div></div><div><span style="background-color:rgb(173,215,115)"><br></span></div></div><div>However I've faced below error:<br></div><div><span style="background-color:rgb(228,172,100)">fatal:(lon) is not a dimension name in variable (var), could not determine dimension number</span><br clear="none"></div><div><br clear="none"></div><div>I searched in google and I found may be <span style="background-color:rgb(125,190,241)">"<span style="color:rgb(0,0,0);font-size:medium">There is no 'lon' dimension name. In fact, there are no latitude/longitude </span><span style="color:rgb(0,0,0);font-size:medium">coordinates on the file at all.</span>"</span><br></div><div>So I used following line:</div><div><div>    <span style="background-color:rgb(210,100,170)">SMAP_fo = linint2_points(lon2d,lat2d,<wbr>var(</span><span style="background-color:rgb(132,90,167)">DIM_000 | :, DIM_001 | :)</span><span style="background-color:rgb(210,100,170)">, False , lon_NW, lat_NW, 0)</span></div><div><br clear="none"></div>instead of:</div><div><div>    <span style="background-color:rgb(210,100,170)">SMAP_fo = linint2_points(lon2d,lat2d,<wbr>var(</span><span style="background-color:rgb(132,90,167)">lat | :, lon | :</span><span style="background-color:rgb(210,100,170)">), False , lon_NW, lat_NW, 0)</span></div><div><br clear="none"></div>and after that I faced following error:</div><div><span style="background-color:rgb(228,172,100)">fatal:linint2_points: If xi is not one-dimensional, then it must have one less dimension than fi<br clear="none"></span></div><div><br clear="none"></div><div>After that I added following lines to the script:</div><div><div><span style="background-color:rgb(210,100,170)">    lat1d = ndtooned(lat2d)</span></div><div><span style="background-color:rgb(210,100,170)">    lon1d = ndtooned(lon2d)</span></div><div><br></div>Again I had following error:</div><div><span style="background-color:rgb(228,172,100)">fatal:linint2_points: The rightmost dimensions of fi must be nyi x nxi, where nyi and nxi are the lengths of yi and xi respectively</span><br></div><div><br clear="none"></div><div>And here I'm not sure about what I've done.</div><div>I've attached the script, stations_NW.csv (input for unstructured points) and links.txt (which contains links for downloading SMAP data. Sorry about that but I couldn't attach SMAP files, because of their sizes. However to download them it is necessary to have earthdata user name!).<br></div><div><div class="m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd76790">I'll be thankful if I could have your help, again.</div></div><div class="m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd02077"><div><br clear="none"></div><div>Sincerely</div><span class="HOEnZb"><font color="#888888"><div>Ehsan<br></div></font></span></div></div></div></div><div class="m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd44909">
                </div></div><div class="m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd00362">
            </div></div></div></div></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>