<div dir="ltr"><div class="gmail_default" style="font-size:small">You are treating "lat2d" and "lon2d" like coordinate arrays with the following call:</div><div class="gmail_default" style="font-size:small"><br></div> curvilinear_to_SCRIP(srcGridName,var&lat2d,var&lon2d,Opt)<div><br></div><div><div class="gmail_default" style="font-size:small">These two variables are NOT attached to var, and they couldn't be anyway, because they violate the rule that coordinate arrays be one-dimensional.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">lat2d and lon2d are standalone variables that you read off a file earlier in your script, so you simply want:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"> curvilinear_to_SCRIP(srcGridName,lat2d,lon2d,Opt)</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 class="gmail_default" style="font-size:small"><div class="gmail_default"><br></div><div class="gmail_default"><br></div></div><div><br></div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 5, 2018 at 5:05 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></div>
<div>Dear NCL users,</div><div>I've modified the script and I have some problems with it, yet.</div><div>New outputs and errors are:</div><div><br></div><div><span class=""><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><span style="background-color:rgb(173,215,115)"><br></span></div></span><span class=""><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><span style="background-color:rgb(173,215,115)"><br></span></div></span><span class=""><div><span style="background-color:rgb(173,215,115)">Variable: var</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: 6</span></div><div><span style="background-color:rgb(173,215,115)"> _FillValue : -9999</span></div><div><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><span style="background-color:rgb(173,215,115)"> long_name : Representative soil moisture measurement for the Earth based grid cell.</span></div><div><span style="background-color:rgb(173,215,115)"> units : cm**3/cm**3</span></div><div><span style="background-color:rgb(173,215,115)"> valid_max : 0.5</span></div><div><span style="background-color:rgb(173,215,115)"> valid_min : 0.02</span></div></span><div><span style="background-color:rgb(253,248,105)">fatal:(lat2d) is not a named dimension in variable (var).</span></div><div><span style="background-color:rgb(253,248,105)">fatal:["Execute.c":8640]:<wbr>Execute: Error occurred at or near line 101 in file 2_SPL3SMP_E_Interpolating1.ncl</span></div><div><br></div>I'll be thankful if I could have your help.</div><div><br></div><div><i>Sincerely</i></div><span class="HOEnZb"><font color="#888888"><div><i>Ehsan Taghizadeh</i></div></font></span><div><div class="h5">
<div id="m_3543603993019868521yahoo_quoted_8403328213" class="m_3543603993019868521yahoo_quoted">
<div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a">
<div>
On Sunday, February 4, 2018, 3:32:14 PM GMT+3:30, Ehsan Taghizadeh <<a href="mailto:ehsantaghizadeh@yahoo.com" target="_blank">ehsantaghizadeh@yahoo.com</a>> wrote:
</div>
<div><br></div>
<div><br></div>
<div><div id="m_3543603993019868521yiv4071997559"><div><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 clear="none"></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,<wbr>var&LAT,var&LON,Opt)</span></div><div><span style="background-color:rgb(253,248,105)">------------------------------<wbr>-------------------------^</span></div><div><span style="background-color:rgb(253,248,105)"><br clear="none"></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,"<wbr>5x5",Opt)</span></div><div><span style="background-color:rgb(253,248,105)">------------------------------<wbr>----------^</span></div><div><span style="background-color:rgb(253,248,105)"><br clear="none"></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_<wbr>file,dstGridName,wgtFileName,<wbr>Opt)</span></div><div><span style="background-color:rgb(253,248,105)">------------------------------<wbr>------------------------------<wbr>---^</span></div><div><span style="background-color:rgb(253,248,105)"><br clear="none"></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 clear="none"></div></div><div><br clear="none"></div><div>Sincerely</div><div>Ehsan</div><div><br clear="none"></div>
<div class="m_3543603993019868521yiv4071997559yahoo_quoted" id="m_3543603993019868521yiv4071997559yahoo_quoted_8638518042">
<div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a">
<div class="m_3543603993019868521yiv4071997559yqt1462782316" id="m_3543603993019868521yiv4071997559yqtfd57874"><div>
On Tuesday, January 30, 2018, 7:34:40 AM GMT+3: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 id="m_3543603993019868521yiv4071997559"><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" href="http://www.ncl.ucar.edu/Document/Functions/Built-in/linint2_points.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Document/Functions/Built-in/<wbr>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" href="http://www.ncl.ucar.edu/Applications/narr.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>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="m_3543603993019868521yiv4071997559gmail_extra"><br clear="none"><div class="m_3543603993019868521yiv4071997559gmail_quote">On Mon, Jan 29, 2018 at 9:23 AM, Ehsan Taghizadeh via ncl-talk <span dir="ltr"><<a rel="nofollow" shape="rect" href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>></span> wrote:<br clear="none"><blockquote class="m_3543603993019868521yiv4071997559gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="m_3543603993019868521yiv4071997559yqt5353952860" id="m_3543603993019868521yiv4071997559yqt90247"><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="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayahoo_quoted" id="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayahoo_quoted_7713408982"><div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a"><div><div id="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayiv8920612249"><div style="font-family:bookman old style,new york,times,serif;font-size:16px"><div class="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayiv8920612249ydpbbc61b17yahoo_quoted" id="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayiv8920612249ydpbbc61b17yahoo_quoted_8100818481"><div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a"><div><div id="m_3543603993019868521yiv4071997559m_-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><span style="background-color:rgb(173,215,115)">Variable: var</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: 6</span></div><div><span style="background-color:rgb(173,215,115)"> _FillValue : -9999</span></div><div><span style="background-color:rgb(173,215,115)"> coordinates : /Soil_Moisture_Retrieval_Data_ AM/latitude /Soil_Moisture_Retrieval_Data_ AM/longitude</span></div><div><span style="background-color:rgb(173,215,115)"> long_name : Representative soil moisture measurement for the Earth based grid cell.</span></div><div><span style="background-color:rgb(173,215,115)"> units : cm**3/cm**3</span></div><div><span style="background-color:rgb(173,215,115)"> valid_max : 0.5</span></div><div><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="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd76790">I'll be thankful if I could have your help, again.</div></div><div class="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd02077"><div><br clear="none"></div><div>Sincerely</div><span class="m_3543603993019868521yiv4071997559HOEnZb"><font color="#888888"></font></span><div>Ehsan<br clear="none"></div></div></div></div></div><div class="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayiv8920612249yqtfd44909">
</div></div><div class="m_3543603993019868521yiv4071997559m_-7490938470158028763ydpdfb568aayiv8920612249yqt7366927260" id="m_3543603993019868521yiv4071997559m_-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" 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"></blockquote></div><br clear="none"></div></div></div></div>
</div></div><div class="m_3543603993019868521yiv4071997559yqt1462782316" id="m_3543603993019868521yiv4071997559yqtfd43205">
</div></div></div></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>