grads-2.0.a1 segmentation fault

[Jennifer Adams]

Hi, Steven --
Great detective work! I downloaded the file you pointed to and now I
can see the error -- I fixed the lat/lon/time as you suggest and with
'set gxout contour' got :
     Result value = 4.26098e-311
but with 'set gxout print' I got:
     Printing Grid -- 1 Values -- Undef = 9.999e+20
     9.999e+20

This is clearly a bug -- I will look into it and put a fix in the
next release.
Jennifer


On Apr 3, 2008, at 7:28 AM, Steven Weiss wrote:

> 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

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080403/4bccda2a/attachment.html