grads-2.0.a1 segmentation fault

Jennifer Adams jma at COLA.IGES.ORG
Thu Apr 3 10:16:07 EDT 2008


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.c
97c97,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/76993235/attachment.html 


More information about the gradsusr mailing list