segmentation fault

hi, i want to extract two variables from two piles of files, the part of script is presented as follows


do i = 0,dimsizes(yrs)-1
; Rainf --> PRECTmms
dRain = systemfunc("ls "+data_indir+"/"+yrs(i)+"/"+"Rainf"+"/"+"*.grb")
setfileoption("grb","SingleElementDimensions","Initial_time")
fRain = addfiles(dRain,"r")
ListSetType(fRain,"cat")
vRain = fRain[:]->VAR_132_GDS0_SFC
printVarSummary(vRain)

; Tair --> TBOT
dTair = systemfunc("ls "+data_indir+"/"+yrs(i)+"/"+"Tair"+"/"+"*.grb")
setfileoption("grb","SingleElementDimensions","Initial_time")
fTair = addfiles(dTair,"r")
ListSetType(fTair,"cat")
vTair = fTair[:]->TMP_GDS0_SFC
printVarSummary(vTair)
.....

end do

after i execute the code with "ncl -x myfile.ncl", there come to the "segmentation fault (core dumped) " after appear the first printVarSummary.

Variable: vRain
Type: float
Total Size: 10091520000 bytes
            2522880000 values
Number of Dimensions: 3
Dimensions and sizes: [initial_time0_hours | 2920] x [g0_lat_1 | 600] x [g0_lon_2 | 1440]
Coordinates:
            initial_time0_hours: [1761936..1770693]
            g0_lat_1: [-59.875..89.875]
            g0_lon_2: [-179.875..179.875]
Number Of Attributes: 10
  center : Zagreb
  long_name : Unknown Variable Name
  units : unknown
  _FillValue : 1e+20
  level_indicator : 1
  gds_grid_type : 0
  parameter_table_version : 1
  parameter_number : 132
  forecast_time : 0
  forecast_time_units : hours
Segmentation fault (core dumped)


i wonder if the system memory resulted in the segmentation fault, i wrote twice
setfileoption("grb","SingleElementDimensions","Initial_time") just before addfles, if this is relevant to the segmentation fault?

thanks








dyjbean

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jun 4 02:02:13 2013