grads-2.0.a1 segmentation fault

[Neil Smith]

Any way you can post the corrected gagx.c file for download?
-Neil
--
Neil R. Smith, Comp. Sys. Mngr. neils at tamu.edu
Atmospheric Sciences, Texas A&M Univ. 979/845-6272
--

On Apr 3, 2008, at 9:16 AM, Jennifer Adams wrote:

> Oh, how I wish all bugs could be this quick to find and repair!!!
> Here's the patch in case you are building from source and want to
> fix it right away:
> diff gagx.c97c97,100<         sprintf (pout,"Result value = %g
> \n",pgr->rmin);--->         if (pgr->umin==1)>           sprintf
> (pout,"Result value = %g \n",pgr->rmin);>         else>
> sprintf (pout,"Result value = %g \n",pgr->undef);
>
> Jennifer
>
> On Apr 3, 2008, at 9:29 AM, Jennifer Adams wrote:
>> 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
>>
>>
>>
>
> --
> 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/fbfd16cd/attachment.html