<div dir="ltr"><div>I have a shapefile containing elevation contours I obtained from USGS. They are very high in resolution (20-ft contour intervals). If I run the standard draw command 'draw shp [shapefile]' then <i>most, but not all</i>, of the shapefiles appear to be drawn. See the attached final image as proof. The geographical extent of this image is about 0.1 deg lat/lon.</div><div><br></div><div>The shapefile itself is over 350 MB in size, and I suspect that if I tell GrADS to draw the whole thing, it runs into memory issues, hence the incomplete drawing. However, I cannot confirm that the shapefile itself does not contain the cutouts of the data. It is possible that the shapefile itself is bad somehow. I just don't know how to find out.</div><div><br></div><div>In light of the memory issue, I wrote a script to query the database file and pull out the record numbers corresponding only to a subset of elevations I want (namely, every 100 ft instead). However, when I run my plotting script to individually plot each shapefile using the syntax 'draw shp [shapefile] [record]', no contours are drawn, even though the record numbers being entered appear to be accurate. I am issuing this command in a loop over all of the record numbers I found. For example, to draw all shapefile contours corresponding to the 13,000-ft elevation, the record numbers follow in this table:</div><div><br></div><div>12800 53<br>13000 35<br>13000 60<br>13000 97<br>13000 116<br>13000 129<br>13000 3075<br>13000 3091<br>13000 3167<br>13000 3555<br>13000 3567<br>13000 7884<br>13000 12789<br>13000 12818<br>13000 12850<br>13000 13370<br>13000 13375<br>13000 14231<br>13000 14233<br>13000 14556<br>13000 14645<br>13000 21232<br>13000 21248<br>13000 22460<br>13000 22478<br>13000 22517<br>13000 22552<br>13000 22583<br>13000 22654<br>13000 23316<br>13000 23530<br>13000 24222<br>13000 26865<br>13000 27949<br>13000 28224<br>13000 28300<br>13000 34185<br>13000 36604<br>13000 36621<br>13000 36632<br>13000 36682<br>13000 36695<br>13000 36721<br>13000 36739<br>13000 36802<br>13000 36822<br>13000 37771<br>13000 38627<br>13000 38708</div><div>(48 'draw shp' commands are issued for this subset of the list)</div><div><br></div><div>There are about 45,000 records in the shapefile, and I have used Microsoft Excel to open the dbf file and directly examine the contour values corresponding to some records to check that these numbers do appear correct.</div><div><br></div><div>So why won't GrADS plot these shapefiles when the calls to draw them are specified individually? No error message comes up when I do this.</div><div><br></div><div>Also, I should note that when I examined the contents of the shapefile directly using 'q shp', the records have entries like this:</div><div><br></div><div>4465:  ArcZ  parts=1  vertices=15  XBounds=-105.994:-105.993  YBounds=38.1644:38.1648  ZBounds=0:0  MBounds=-1.79769e+308:-1.79769e+308</div><div><br></div><div>The part that is unexpected is that first entry after the record number: "ArcZ". I am expecting to see "point"/"polygon"/"line" etc. Is that the cause of the failed plotting?</div><div><br></div><div>Thanks.</div><div><br></div><div>Jeff Duda<br></div><div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font size="2"><span>Jeff Duda, Ph.D.</span></font></div><div><font size="2"><span></span>Research scientist with CU/CIRES and NOAA/GSL<br></font></div></div></div></div></div></div></div></div></div></div></div></div>