[gradsusr] pdef file works at v1.9b4 but not at v2.0.a9

[Clark, Douglas B.]

Dear All,

I have a dataset that uses a pdef file and that I can use in GrADS v1.9b4 (in both SunOS and RH executables downloaded from the GrADS website), but which results in a segmentation fault at v2.0.a9 (and, from memory, at v2.0.a8) - again in various downloads. Looks like it might be a bug in GrADS...but of course most things turn out to be user error! Any ideas?

At v2.0.a9:

ga-> q config
Config: v2.0.a9 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids geotiff shapefile
Grid Analysis and Display System (GrADS) Version 2.0.a9
Copyright (c) 1988-2010 by Brian Doty and the
Institute for Global Environment and Society (IGES) 
This program is distributed WITHOUT ANY WARRANTY 
See file COPYRIGHT for more information. 

Built Fri Sep  3 17:21:21 UTC 2010 for x86_64-redhat-linux-gnu

This version of GrADS has been configured with the following options:
  o Built on a LITTLE ENDIAN machine
  o Athena Widget GUI DISABLED
  o Command line editing ENABLED 
      http://tiswww.case.edu/php/chet/readline/rltop.html 
  o printim command for image output ENABLED 
      http://www.zlib.net 
      http://www.libpng.org/pub/png/libpng.html 
      http://www.libgd.org/Main_Page 
  o GRIB2 interface ENABLED 
      http://www.ijg.org 
      http://www.ece.uvic.ca/~mdadams/jasper 
      http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 
      g2clib-1.2.0  
  o NetCDF interface ENABLED 
      http://www.unidata.ucar.edu/software/netcdf  
      netcdf 4.1.1 of Dec  6 2010 08:20:12 $  
  o OPeNDAP gridded data interface ENABLED
  o OPeNDAP station data interface DISABLED
  o HDF4 and HDF5 interfaces ENABLED 
      http://hdfgroup.org 
      HDF 4.2r5 
      HDF5 1.8.5 
  o GeoTIFF and KML/TIFF output ENABLED
      http://www.libtiff.org 
      http://geotiff.osgeo.org 
  o KML contour output ENABLED
  o Shapefile interface ENABLED
      http://shapelib.maptools.org 

For additional information please consult http://iges.org/grads

ga-> open Rainf_WFD_GPCC.ctl 
Scanning description file:  Rainf_WFD_GPCC.ctl
PDEF FILE Error: The offsets in the pdef file 
  must be 1-based (i.e., > 0 and <= isize*jsize). 
  The data file was not opened. 
Segmentation fault (core dumped)

-------------------------------------------------------------

The first few lines of the ctl file are:

dset ^../data/Rainf_WFD_GPCC/Rainf_WFD_GPCC_%y4%m2.nc
options template
dtype netcdf
undef 1.0e20
pdef 67420  1 file 1 stream binary-big ^../../../ancil/data/WFD_0p5deg_pdefData.gra
xdef    720 linear -179.75 0.5
ydef    280 linear  -55.75 0.5


-------------------------------------------------------------

Works fine at v1.9.b4:

ga-> open Rainf_WFD_GPCC.ctl 
Scanning description file:  Rainf_WFD_GPCC.ctl
Data file ../data/Rainf_WFD_GPCC/Rainf_WFD_GPCC_%y4%m2.nc is open as file 1
LON set to 0 360 
LAT set to -55.75 83.75 
LEV set to 1 1 
Time values set: 2001:1:1:0 2001:1:1:0 
Notice: Implied interpolation for file Rainf_WFD_GPCC.ctl
 Interpolation will be performed on any data displayed from this file
ga-> d rainf
Notice:  Automatic Grid Interpolation Taking Place
Contouring: 0 to 0.0033 interval 0.0003


-------------------------------------------------------------

Vague notion: The pdef data file is big-endian (as declared in the pdef line of the ctl file) but
I'm using a little-endian machine. I wonder if v2.0.a9 can't cope with this mix?


Cheers,
Doug




-- 
This message (and any attachments) is for the recipient only. NERC
is subject to the Freedom of Information Act 2000 and the contents
of this email and any reply you make may be disclosed by NERC unless
it is exempt from release under the Act. Any material supplied to
NERC may be stored in an electronic records management system.