Re: extract variables from grib: Segmentation fault

This is not enough information.

You are getting data for one point?

rain = gribfile->Rainf_GDS0_SFC(416,1174)

How big is the grid? Are these subscripts within range?

======

==================================
 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

 gribfiles = systemfunc("ls /d1/zhangg/GLDAS/*A"+year+day3+"*.grb")
 setfileoption("grb","SingleElementDimensions","Initial_time")

  f = addfiles(gribfiles, "r")
  rain = f[:]->Rainf_GDS0_SFC(:,416,1174)
  snow = f[:]->Snowf_GDS0_SFC(:,416,1174)

  printVarSummary(rain)
  printVarSUmmary(snow)

  precip = rain + snow
  copy_VarMeta(rain, precip)
  precip@long_name = "precipitation"
  printVarSummary(precip)

> Hi,
>  
> I'm trying extract variables from grib data.   When I just output precip,
> it works pretty well. But when I added ' outfile->rain = rain
> outfile->snow = snow' , I got error: Segmentation fault. How can I fix it?
> Thanks.
>  
> Guo
>  
> The part of my script is like this:  
>  
> precip = new((endday-startday+1)*8,float)
> gribfiles  = systemfunc("ls /d1/zhangg/GLDAS/*A"+year+day3+"*.grb")
>      do ihr = 0,7
>        print("Starting: "+gribfiles(ihr))
>        gribfile =   addfile(gribfiles(ihr), "r")
>        rain = gribfile->Rainf_GDS0_SFC(416,1174)
>        snow = gribfile->Snowf_GDS0_SFC(416,1174)
>        precip((iday-startday)*8+ihr) = rain + snow
>      end do
>      delete(gribfiles)
>    
>    outfile  =
> addfile("output_precip_"+year+"_"+startday3+"_"+endday3+".nc", "c")
>    outfile->precip = precip
>    outfile->rain  = rain
>    outfile->snow  = snow_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Apr 13 16:26:56 2012