<html><head></head><body><div style="font-family:bookman old style, new york, times, serif;font-size:16px;"><div></div>
<div>Dear Dennis</div><div>Thank you so much for you help.</div><div>I've changed my script to attached one by adding lines 79-91, however it seems to have some basic problems! Could I ask help me about them, please? I faced following errors:</div><div><br></div><div><div><span style="background-color: rgb(253, 248, 105);">fatal:syntax error: line 87 in file 2_SPL3SMP_E_Interpolating1.ncl before or near \n</span></div><div><span style="background-color: rgb(253, 248, 105);"> curvilinear_to_SCRIP(src_file,var&LAT,var&LON,Opt)</span></div><div><span style="background-color: rgb(253, 248, 105);">-------------------------------------------------------^</span></div><div><span style="background-color: rgb(253, 248, 105);"><br></span></div><div><span style="background-color: rgb(253, 248, 105);">fatal:syntax error: possibly an undefined procedure</span></div><div><span style="background-color: rgb(253, 248, 105);">fatal:syntax error: line 88 in file 2_SPL3SMP_E_Interpolating1.ncl before or near \n</span></div><div><span style="background-color: rgb(253, 248, 105);"> latlon_to_SCRIP(dstGridName,"5x5",Opt)</span></div><div><span style="background-color: rgb(253, 248, 105);">----------------------------------------^</span></div><div><span style="background-color: rgb(253, 248, 105);"><br></span></div><div><span style="background-color: rgb(253, 248, 105);">fatal:syntax error: possibly an undefined procedure</span></div><div><span style="background-color: rgb(253, 248, 105);">fatal:syntax error: line 89 in file 2_SPL3SMP_E_Interpolating1.ncl before or near \n</span></div><div><span style="background-color: rgb(253, 248, 105);"> ESMF_regrid_gen_weights(src_file,dstGridName,wgtFileName,Opt)</span></div><div><span style="background-color: rgb(253, 248, 105);">---------------------------------------------------------------^</span></div><div><span style="background-color: rgb(253, 248, 105);"><br></span></div><div><span style="background-color: rgb(253, 248, 105);">fatal:syntax error: possibly an undefined procedure</span></div><div><span style="background-color: rgb(253, 248, 105);">fatal:Syntax Error in block, block not executed</span></div><div><span style="background-color: rgb(253, 248, 105);">fatal:error at line 126 in file 2_SPL3SMP_E_Interpolating1.ncl</span></div><div><br></div></div><div><br></div><div>Sincerely</div><div>Ehsan</div><div><br></div>
<div id="yahoo_quoted_8638518042" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
On Tuesday, January 30, 2018, 7:34:40 AM GMT+3:30, Dennis Shea <shea@ucar.edu> wrote:
</div>
<div><br></div>
<div><br></div>
<div><div id="yiv3306159953"><div><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 clear="none"><br clear="none"> <a rel="nofollow" shape="rect" target="_blank" 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 clear="none"></div><div><br clear="none"></div>[2] The SMAP data are on a *curvilinear* grid: [DIM_000 | 406] x [DIM_001 | 964] . Hence, 'linint2_points' will not work.<br clear="none"><br clear="none"></div><div>[3] I suggest that you <br clear="none"> (a) use ESMF regridding to interpolate from the curvilinear grid to a rectilinear grid of approximately the same resolution.<br clear="none"></div><div> (b) then use linint2_points on (3a)<br clear="none"><br clear="none">===<br clear="none"></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 clear="none"><br clear="none"> <a rel="nofollow" shape="rect" target="_blank" href="http://www.ncl.ucar.edu/Applications/narr.shtml">http://www.ncl.ucar.edu/Applications/narr.shtml</a><br clear="none"><br clear="none"></div><div>Good luck<br clear="none"></div><div> <br clear="none"></div><div><span style="background-color:rgb(173,215,115);"></span><br clear="none"><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="yiv3306159953gmail_extra"><br clear="none"><div class="yiv3306159953gmail_quote">On Mon, Jan 29, 2018 at 9:23 AM, Ehsan Taghizadeh via ncl-talk <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:ncl-talk@ucar.edu" target="_blank" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>></span> wrote:<br clear="none"><blockquote class="yiv3306159953gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="yiv3306159953yqt5353952860" id="yiv3306159953yqt90247"><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 clear="none"></div><div class="yiv3306159953m_-7490938470158028763ydpdfb568aayahoo_quoted" id="yiv3306159953m_-7490938470158028763ydpdfb568aayahoo_quoted_7713408982"><div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;"><div><div id="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249"><div style="font-family:bookman old style, new york, times, serif;font-size:16px;"><div class="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249ydpbbc61b17yahoo_quoted" id="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249ydpbbc61b17yahoo_quoted_8100818481"><div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;"><div><div id="yiv3306159953m_-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 clear="none"></div><div><div>printVarSummary(lat2d)<br clear="none"></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 clear="none"></div></div><div><div>printVarSummary(lon2d)<br clear="none"></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 clear="none"></div></div><div>printVarSummary(var):<br clear="none"></div><div><div style=""><span style="background-color:rgb(173,215,115);">Variable: var</span></div><div style=""><span style="background-color:rgb(173,215,115);">Type: float</span></div><div style=""><span style="background-color:rgb(173,215,115);">Total Size: 1565536 bytes</span></div><div style=""><span style="background-color:rgb(173,215,115);"> 391384 values</span></div><div style=""><span style="background-color:rgb(173,215,115);">Number of Dimensions: 2</span></div><div style=""><span style="background-color:rgb(173,215,115);">Dimensions and sizes: [DIM_000 | 406] x [DIM_001 | 964]</span></div><div style=""><span style="background-color:rgb(173,215,115);">Coordinates:</span></div><div style=""><span style="background-color:rgb(173,215,115);">Number Of Attributes: 6</span></div><div style=""><span style="background-color:rgb(173,215,115);"> _FillValue : -9999</span></div><div style=""><span style="background-color:rgb(173,215,115);"> coordinates : /Soil_Moisture_Retrieval_Data_ AM/latitude /Soil_Moisture_Retrieval_Data_ AM/longitude</span></div><div style=""><span style="background-color:rgb(173,215,115);"> long_name : Representative soil moisture measurement for the Earth based grid cell.</span></div><div style=""><span style="background-color:rgb(173,215,115);"> units : cm**3/cm**3</span></div><div style=""><span style="background-color:rgb(173,215,115);"> valid_max : 0.5</span></div><div style=""><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 clear="none"></span></div></div><div>However I've faced below error:<br clear="none"></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 clear="none"></div><div>So I used following line:</div><div><div> <span style="background-color:rgb(210,100,170);">SMAP_fo = linint2_points(lon2d,lat2d, 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, 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 clear="none"></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 clear="none"></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 clear="none"></div><div><div class="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd76790">I'll be thankful if I could have your help, again.</div></div><div class="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd02077"><div><br clear="none"></div><div>Sincerely</div><span class="yiv3306159953HOEnZb"><font color="#888888"></font></span><div>Ehsan<br clear="none"></div></div></div></div></div><div class="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd44909">
</div></div><div class="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="yiv3306159953m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd00362">
</div></div></div></div></div>
</div>
</div></div></div></div><br clear="none">______________________________ _________________<br clear="none">
ncl-talk mailing list<br clear="none">
<a rel="nofollow" shape="rect" ymailto="mailto:ncl-talk@ucar.edu" target="_blank" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br clear="none">
List instructions, subscriber options, unsubscribe:<br clear="none">
<a rel="nofollow" shape="rect" target="_blank" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/ mailman/listinfo/ncl-talk</a><br clear="none">
<br clear="none"></blockquote></div><br clear="none"></div></div></div></div>
</div>
</div></div></body></html>