grads-2.0.a1 segmentation fault

[Steven Weiss]

Hi Jennifer,

Thanks, when would the next release be available? I dont think I'm going to be building from source. I will have to ignore this variable for the meantime.

Regards
Steven
  ----- Original Message ----- 
  From: Jennifer Adams 
  To: GRADSUSR at LIST.CINECA.IT 
  Sent: Thursday, April 03, 2008 3:29 PM
  Subject: Re: grads-2.0.a1 segmentation fault


  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/0d1f2c9b/attachment.html