grads-2.0.a1 segmentation fault

[Steven Weiss]

Hi,

I'm getting a segmentation fault which occurs inconsistantly. I've managed
to narrow it down to a particular variable in the grib file - DISRswsfc that
sometimes has garbage in it for certain time steps. I'm working with the
wavewatch 3 data from ncep/noaa.

Below is an example of what is happening:

Grads Output
------------
Grid Analysis and Display System (GrADS) Version 2.0.a1
Copyright (c) 1988-2008 by Brian Doty and the
Institute for Global Environment and Society (IGES)
GrADS comes with ABSOLUTELY NO WARRANTY
See file COPYRIGHT for more information

Config: v2.0.a1 little-endian readline printim grib2 netcdf hdf4-sds
Issue 'q config' command for more information.
Landscape mode? ('n' for portrait):
GX Package Initialization: Size = 11 8.5
Running in Batch mode
ga-> sw.gs
Result value = -1.04874e+80
Segmentation fault

Downloaded data
---------------
ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/wave/prod/wave.20080403/nww3.t00z.grib.grib2

Control file
-------------
dset ^nww3.all.grb
index ^nww3.all.grb.idx
undef 9.999E+20
title nww3.all.grb
*  produced by g2ctl v0.0.3d
* griddef=1:0:(288 x 157):grid_template=0: lat-lon grid:(288 x 157) units 1e-06
input WE:NS output WE:SN res 48 lat 78.000000 to -78.000000 by 1.000000 lon 0.00
0000 to 358.750000 by 1.250000 #points=45216

dtype grib2
ydef 157 linear -78.000000 1
xdef 288 linear 0.000000 1.250000
tdef 61 linear 00Z03apr2008 3hr
zdef 1 linear 1 1
vars 11
DIRPWsfc  0,1,0   10,0,10 ** surface none Primary Wave Direction [deg]
DIRSWsfc  0,1,0   10,0,12 ** surface none Secondary Wave Direction [deg]
HTSGWsfc  0,1,0   10,0,3 ** surface none Significant Height of Combined Wind
 Waves and Swell [m]
PERPWsfc  0,1,0   10,0,11 ** surface none Primary Wave Mean Period [s]
PERSWsfc  0,1,0   10,0,13 ** surface none Secondary Wave Mean Period [s]
UGRDsfc  0,1,0   0,2,2 ** surface none U-Component of Wind [m/s]
VGRDsfc  0,1,0   0,2,3 ** surface none V-Component of Wind [m/s]
WDIRsfc  0,1,0   0,2,0 ** surface none Wind Direction (from which blowing) [deg]
WVPERsfc  0,1,0   10,0,6 ** surface none Mean Period of Wind Waves [s]
ENDVARS

Piece of code for lat/lon -29 16.25
-----------------------------------
   fmt2 = "%-.2f"
   'd dirswsfc'
   val=subwrd(result,4)
   ln=sublin(result,1)
   say ln
   if (val = "9.999e+20" | val = "1e+20")
      val = "N/A"
   else
      rc = math_format(fmt2,val)
*** the segmentation fault occurs on line above due to a bad number in val
   endif