FW: Re: New Win32 Release of GrADS v1.9.0-rc1

Mark Sponsler msponsler at COMCAST.NET
Wed Apr 9 02:21:36 EDT 2008


Thanks Arlindo. Am looking forward to hearing how the ERSI data works, especially in hi-res situations.

I've been playing with shape files available from some of the other sites described in the wiki this evening.  Much good stuff out there.  Alot of work from many divergent sources.

I have noticed that in all cases the 'set line color style thickness'  does not work with regard to style or thickness of lines though.  Not sure if that is in your domain or part of Graziano's shapefile extension.

Just thought I'd mention that as a candidate for a future fix.

Again, many thanks.  Ver 1.9 is a pleasure to work with.
--
Thanks,
Mark

-------------- Original message --------------
From: Arlindo da Silva <dasilva at ALUM.MIT.EDU>



On Wed, Apr 9, 2008 at 12:10 AM, Mark Sponsler <msponsler at comcast.net> wrote:

Many thanks Arlindo.

I'm not quite sure I understand exactly all that Graziano is saying, but it sound like the outlines of areas should work fine, but that the filling of those shapes probably will not.  I played with it a bit more, trying to fill landmasses that were completely bordered (e.g. the big island of Hawaii), but no luck.

Still the hires coastline (and lakes) found in the  '_lines' files are of much value.

Thanks again for all your help.

If any one has other suggestions or ideas it would be appreciated.


I guess what Graziano was saying is that the free files we get from NOAA do not define polygons that can be filled. Commercial shapefile files available from ESRI are more carefully constructed as to permit the polygon fill. I'll be purchasing the DVD from AGI and report back whether it is worth the price tag (> $100).





--
Thanks,
Mark

-------------- Original message --------------
From: Arlindo da Silva <dasilva at ALUM.MIT.EDU>
Mark,


Graziano who is temporarily unable to post to this list sends the following response:


Hi Mark,
gshhs shapefile contains same data as in the binary file, so we have
mostly arcs and not so much polygons.
To fill a shape, it must be closed.
ESRI data are much more quality checked and with lot of useful geometric
characteristic, but they are non-free (being the same very diffused).

Does anybody here know about a source of global public GIS data in a
format supported by gdal?

Graziano.



On Tue, Apr 8, 2008 at 12:54 PM, Mark Sponsler <msponsler at comcast.net> wrote:

Hi Arlindo,
I tried downloading the gshhs lake files too and received the same result.  Here's the segment of the script that overlays the shapefiles. I've attached the output (png file)

'set line 33 1 6'
'shp_lines gshhs_land'
'shp_polyf gshhs_land'

'shp_lines hshhs_lake'
'shp_polyf hshhs_lake'

Also here is the out put from Grads.  As best I can tell it is reading the shapefile fine and appears to be drawing both the lines and polygons, but the polygons do not appear. Any thoughts would be appreciated.


Starting grads under c:\PROGRA~1\PCGRAD~2\win32 ...
Grid Analysis and Display System (GrADS) Version 1.9.0-rc1
Copyright (c) 1988-2007 by Brian Doty and IGES
Center for Ocean-Land-Atmosphere Studies (COLA)
Institute for Global Environment and Society (IGES)
GrADS comes with ABSOLUTELY NO WARRANTY
See file COPYRIGHT for more information
Config: v1.9.0-rc1 32-bit little-endian readline sdf/xdf netcdf dap dap-stn lats
 athena printim
Issue 'q config' command for more information.
Error opening user function definition table
  File name is: c:\Program Files\PCGrADS19\win32\udf
GX Package Initialization: Size = 11 8.5
Running in Batch mode
Grib = test
reg       = sf
map_reg = sf
map = W140N40
 Displaying Beach sf W140N40
 Define beach1 complete
Beach Overlay Complete
 Displaying Relief sf W140N40
Using shapefile /dat/gshhs_land
Shapefile Type: Arc   # of Shapes: 191811
File Bounds: (-180.000, -78.570) to ( 180.000,  83.633)
Done shapefile /dat/gshhs_land
Using shapefile /dat/gshhs_land
Shapefile Type: Arc   # of Shapes: 191811
File Bounds: (-180.000, -78.570) to ( 180.000,  83.633)
Done shapefile /dat/gshhs_land
Using shapefile /dat/hshhs_lake
Shapefile Type: Arc   # of Shapes: 7287
File Bounds: (-180.000, -67.123) to ( 180.000,  83.327)
Done shapefile /dat/hshhs_lake
Using shapefile /dat/h shhs_lake
Shapefile Type: Arc   # of Shapes: 7287
File Bounds: (-180.000, -67.123) to ( 180.000,  83.327)
Done shapefile /dat/hshhs_lake
Start Print
   ** sf Print Complete **
End of file
Region Loop closed
Processing Complete
ga->

--
Thanks,
Mark

-------------- Forwarded Message: --------------
From: msponsler at comcast.net (Mark Sponsler)
To: GRADSUSR at LIST.CINECA.IT
Subject: Re: New Win32 Release of GrADS v1.9.0-rc1
Date: Mon, 07 Apr 2008 20:24:48 +0000


Hi Arlindo,
Thanks so much for providing these new capabilities.  This is wonderful stuff!

I've been playing with shapefiles from gshhs and ersi using Win32 superpack non-self installing version.  Can get the outline of the coast to display fine from both sources but not the polygons from gshhs. Polygons from ersi works fine though. Also cannot seem to get the coast outline to draw 'thicker'.  Not sure if that is an available option, but here's the code I'm using.

'set line 33 1 6'
'shp_lines gshhs_land'
'shp_polyf gshhs_land'

Am I dong something wrong?
--
Thanks,
Mark


-------------- Original message --------------
From: Arlindo da Silva <dasilva at ALUM.MIT.EDU>

All,


  I just uploaded to sf.net a new Win32 GrADS release:


     http://sourceforge.net/project/showfiles.php?group_id=161773&package_id=270750


Notice that Win32 GrADS now has its own packaged called "grads-windows". See the ChangeLog for additional info.


   This version fixes some bugs with DOS-style text files, and includes Release 2 of the OpenGrADS extensions, with Grib-2 and Shapefile support (as well as an improved "fish"). To quickly checkout the new features:


Grib-2:


   ga-> g2_open ../data/model_25.gdf
   ga-> d g2(tmpprs)


or with the % shorthand


   ga-> % d tmpprs


Shapefile:


   ga-> open ../data/model.ctl
   ga-> d ts
   ga-> shp_polyf admin98


Fish:


   ga-> set lev 200
   ga-> d fish_psi(ua,va)
   ga-> clear
   ga-> d fish_chi(ua,va)


   For documentation,


          http://opengrads.org/doc/



   Let me know of any problems.


     Arlindo




--
Arlindo da Silva
dasilva at alum.mit.edu



--
Arlindo da Silva
dasilva at alum.mit.edu



--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080409/93999fb3/attachment.html 


More information about the gradsusr mailing list