<div dir="ltr"><div>Hi Peirong,<br><br><br></div><div>The problem turns out to be an unfortunate bug in handling 3D shapefiles (your data are perhaps the first 3D shapefile I've encountered). The bug has been fixed and should be available in an upcoming 6.2.1 release of NCL. We'll work with you to try to get you a fixed binary in the meantime.<br>
<br></div><div>Thanks for your help in pointing out the issue and providing an example!<br><br></div><div>Rick<br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Aug 6, 2014 at 10:11 AM, Peirong Lin <span dir="ltr"><<a href="mailto:peironglinlin@gmail.com" target="_blank">peironglinlin@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Thanks anyway! Hope to hear from you soon. <div><br></div><div>Best,</div><div>Peirong</div></div><div class="HOEnZb">
<div class="h5"><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Aug 6, 2014 at 9:50 AM, Rick Brownrigg <span dir="ltr"><<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@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 dir="ltr"><div><div>Hi Peirong,<br><br></div>I am looking into this issue with the files you sent. This appears to be a serious bug somewhere, but I as of yet can not tell where it is originating. I've verified separately, using GDAL utilities, that the COMID field is indeed non-zero. I will post an update on progress later, but this may take much time to track down.<span><font color="#888888"><br>
<br></font></span></div><span><font color="#888888">Rick<br><br></font></span></div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Aug 5, 2014 at 3:23 PM, Rick Brownrigg <span dir="ltr"><<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@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 dir="ltr"><div><div>Hi Peirong,<br><br>Are you able to share those shapefiles, and if so, would you mind sending them to me? If they are too big to be emailed, there's our ftp site:<br>
<br><a href="http://ncl.ucar.edu/report_bug.shtml" target="_blank">http://ncl.ucar.edu/report_bug.shtml</a><br>
<br></div>Thanks,<br></div>Rick<br><div><div><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote"><div>On Tue, Aug 5, 2014 at 3:13 PM, Peirong Lin <span dir="ltr"><<a href="mailto:peironglinlin@gmail.com" target="_blank">peironglinlin@gmail.com</a>></span> wrote:<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr">Dear NCL,<div><br></div><div>I am using the following NCL code to read a shapefile attribute table (which allows further computation), but when I printed out all the "COMIDs", they are all "0". But when I open the shapefile in ArcGIS, the COMIDs can be correctly shown. Can anybody provide some hints to resolve this? </div>
<div><br></div><div><b>begin<br></b></div><div>
<p><b>shp_file = "./NHDFlowline.shp"</b></p>
<p><b>f = addfile(shp_file,"r")</b></p>
<p><b>varnames = getfilevarnames(f)</b></p>
<p><b>print("varnames = "+varnames)</b></p>
<p><b>COMIDs = f->COMID</b></p>
<p><b>print(COMIDs)</b></p>
<p><b>end</b></p></div><div><br></div><div>I've already made sure the COMID is a variable in the shapefile, but the values just cannot be read correctly. Thank you in advance for replying!!</div><div><br></div>
<div>Best,</div><div>Peirong</div><span><font color="#888888"><div><br clear="all"><div><br></div>-- <br><span style="font-family:arial,sans-serif;font-size:13px;border-collapse:collapse"><div style="color:rgb(255,102,0)">
<b>Peirong Lin</b> (Ms.)<br>
Graduate Research Assistant <span style="font-family:arial,sans-serif;font-size:13px;border-collapse:collapse"><span style="font-family:arial,sans-serif;font-size:13px;border-collapse:collapse;color:rgb(204,102,0)">(JGB 5.322)</span></span><br>
Jackson School of Geosciences</div><div><font color="#ff6600">The University of Texas at Austin</font><br><font color="#ff6600">Contact: <a href="tel:512-201-3731" value="+15122013731" target="_blank">512-201-3731</a> </font><br>
<a href="mailto:prlin@utexas.edu" style="color:rgb(255,102,0)" target="_blank">prlin@utexas.edu</a><font color="#ff6600">; </font><a href="mailto:peironglinlin@gmail.com" style="color:rgb(255,102,0)" target="_blank">peironglinlin@gmail.com</a><span style="font-family:arial,sans-serif;font-size:13px;border-collapse:collapse"></span></div>
</span>
</div></font></span></div>
<br></div></div><div>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></div></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><span style="font-family:arial,sans-serif;font-size:13px;border-collapse:collapse"><div style="color:rgb(255,102,0)"><b>Peirong Lin</b> (Ms.)<br>Graduate Research Assistant <span style="font-family:arial,sans-serif;font-size:13px;border-collapse:collapse"><span style="font-family:arial,sans-serif;font-size:13px;border-collapse:collapse;color:rgb(204,102,0)">(JGB 5.322)</span></span><br>
Jackson School of Geosciences</div><div><font color="#ff6600">The University of Texas at Austin</font><br><font color="#ff6600">Contact: <a href="tel:512-201-3731" value="+15122013731" target="_blank">512-201-3731</a> </font><br>
<a href="mailto:prlin@utexas.edu" style="color:rgb(255,102,0)" target="_blank">prlin@utexas.edu</a><font color="#ff6600">; </font><a href="mailto:peironglinlin@gmail.com" style="color:rgb(255,102,0)" target="_blank">peironglinlin@gmail.com</a><span style="font-family:arial,sans-serif;font-size:13px;border-collapse:collapse"></span></div>
</span>
</div>
</div></div></blockquote></div><br></div>