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

[Sergey Varlamov]

Dear Douglas,

Please try to change

pdef 67420 1 file ...

to you real grid size

pdef ix jy file ...

where ix and jy are dimensions of your original curvilinear grid,
such that ix*jy=67420.
It must still work in version 1.9 and could work in versions 2 of GrADS.

Sergey Varlamov,

JAMSTEC, Yokohama


Clark, Douglas B. wrote:
> 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
>
>
>
>
>   


-- 

Sergey Varlamov

JAMSTEC, Yokohama