<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>Regarding my problem drawing counties on an Alaska map... I've tried all the suggestions and still am getting no results.</div><div>Interesting result when I did Jennifer's query on the shape file (which others are using with success), it says:</div><div><br></div><div>UnknownShapeType parts=0 verticies=0 XBounds=0:0 YBounds=0:0 ZBounds=0:0 MBounds= 0:0 </div><div><br></div><div>for every shape.</div><div><br></div><div>I don't have this problem with Admin98. I've tried many different shape files and can't seem to plot anything unless it's from Admin98.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
------------------------------<br>
<br>Date: Sat, 20 Aug 2016 14:01:13 +0000<br>
From: Jennifer M Adams <<a href="mailto:jadams21@gmu.edu">jadams21@gmu.edu</a>><br>
Subject: Re: [gradsusr] Drawing Shape Files<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID: <<a href="mailto:3807099F-678F-45D3-AC16-713AE9A20CFC@gmu.edu">3807099F-678F-45D3-AC16-<wbr>713AE9A20CFC@gmu.edu</a>><br>
Content-Type: text/plain; charset="windows-1252"<br>
<br>
Hi, Andrew ?<br>
You?ve gotten a lot of good help from other users on this question. The only thing I can add is that maybe it?s a longitude wrapping issue. Try testing with the output from ?q shp? ? this returns a list of all the polygons in the file and some metadata about them including their lat/lon boundaries. E.g.<br>
<br>
ga-> q shp admin98<br>
Shapefile Type=Polygon #Shapes=2604 XBounds=-180:180 YBounds=-90:83.6236<br>
0: Polygon parts=1 vertices=19 XBounds=-70.0633:-69.8733 YBounds=12.4111:12.6311 ZBounds=0:0 MBounds=0:0<br>
1: Polygon parts=2 vertices=48 XBounds=-61.8911:-61.6664 YBounds=16.9897:17.725 ZBounds=0:0 MBounds=0:0<br>
2: Polygon parts=1 vertices=462 XBounds=69.9325:74.9157 YBounds=35.4396:38.4721 ZBounds=0:0 MBounds=0:0<br>
<br>
If the XBounds are defined from -180 to 180 and you are drawing 0 to 360, you may not see all the shapes.<br>
?Jennifer<br>
<br>
<br>
On Aug 17, 2016, at 4:10 PM, Andrew Revering <<a href="mailto:andy@f5data.com">andy@f5data.com</a><mailto:<a href="mailto:andy@f5data.com">andy@<wbr>f5data.com</a>>> wrote:<br>
<br>
Thanks for the link Stephen.<br>
I tried the file but still doesn't draw. On the same map I used Admin98 and it drew the Yukon Territory no problem. I also tried putting these shape files into the same path as admin98 but no luck.<br>
<br>
open gfs003.ctl<br>
Scanning description file: gfs003.ctl<br>
Data file gfs003.grib is open as file 1<br>
LON set to 0 360<br>
LAT set to -90 90<br>
LEV set to 1 1<br>
Time values set: 2016:8:17:9 2016:8:17:9<br>
E set to 1 1<br>
<br>
ga-> set lon -179.8 -127.2<br>
LON set to -179.8 -127.2<br>
ga-> set lat 50.9 73.9<br>
LAT set to 50.9 73.9<br>
<br>
ga-> set shp cb_2015_us_county_500k.shp<br>
Shapefile output file name root: cb_2015_us_county_500k.shp<br>
Shapefile output type: line<br>
Shapefile format string is %12.6f<br>
ga-> d tmp2m<br>
Contouring: 264 to 291 interval 3<br>
<br>
ga-> draw shp admin98<br>
ga-> draw shp cb_2015_us_county_500k.shp<br>
ga-> printim test.png<br>
Created PNG file test.png<br>
<br>
result: <a href="http://f5data.net/test.png" rel="noreferrer" target="_blank">http://f5data.net/test.png</a><br>
<br></blockquote></div></div></div>