Segmentation fault (core dumped)

From: Vishali P. <vishali.prat_at_nyahnyahspammersnyahnyah>
Date: Tue Jul 05 2011 - 04:27:17 MDT

Hi again.

I'm using 4x Daily NCEP/DOE AMIP-II Reanalysis wind data 1998 to 2004. I
have a script to calculate unweighted seasonal means.
As i was suggested before, I've tried to select region and levels, but the
error stills the same.

"Segmentation fault (core dumped)".

Can anyone help me?

Thank you in advance

Vishali

*
[1] My script*

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

   yrStrt = 1998
   yrLast = 2004
   nyrs = yrLast-yrStrt+1
   klev = 2
   nlat = 8
   mlon = 8
   ntim = 12*nyrs ; number of months

   wMon = new ( (/ntim,klev,nlat,mlon/), "float", 1e20)
   tMon = new ( ntim, "double", "No_FillValue")

   diru = "/home/Vishali/Documents/dadosVentos/SpeedTest/"
   filu = systemfunc("cd "+diru+" ; ls wind*red.nc")
   nfilu = dimsizes(filu)

                    ; calculate monthly means from 6hrly
   nmoStrt = 0
   nmoLast = 11
   do nf=0,nfilu-1

    print("nf="+nf+" "+filu(nf))

      f = addfile(diru+filu(nf), "r")
      w = f->WSpeed
      wMon(nmoStrt:nmoLast,:,:,:) = calculate_monthly_values(w, "avg",0,
False)
       nmoStrt = nmoStrt+12
       nmoLast = nmoLast+12
      delete(w) ; may change size next iteration (leap yr)
   end do
   delete(wMon&time)

   printVarSummary(wMon) ; wMon(84,17,73,144)

   time = yyyymm_time(yrStrt, yrLast, "integer")
   wMon&time = time

   wMon!1 = "level" ; minor bug in "calculate_monthly_values"
   wMon&level = f->level

   printVarSummary(wMon)
                          ; unweighted seasonal means
   wAllSea = wgt_runave_n_Wrap(wMon, 3, 0, 0)
   printVarSummary(wAllSea)

   wSea = wAllSea(::3,:,:,:) ; (28,14,73,144)
   printVarSummary(wSea)

*
[2] Output with "ncl -x" command *

Loading file "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
+
+
+ yrStrt = 1998
+ yrLast = 2004
+ nyrs = yrLast-yrStrt+1
+ klev = 2
+ nlat = 8
+ mlon = 8
+ ntim = 12*nyrs ; number of months
+
+ wMon = new ( (/ntim,klev,nlat,mlon/), "float", 1e20)
+ tMon = new ( ntim, "double", "No_FillValue")
+
+
+ diru = "/home/Vishali/Documents/dadosVentos/SpeedTest/"
+ filu = systemfunc("cd "+diru+" ; ls wind*red.nc")
+ nfilu = dimsizes(filu)
+
+ ; calculate monthly means from 6hrly
+ nmoStrt = 0
+ nmoLast = 11
+ do nf=0,nfilu-1
+
+ print("nf="+nf+" "+filu(nf))
+
+ f = addfile(diru+filu(nf), "r")
+ w = f->WSpeed
+ wMon(nmoStrt:nmoLast,:,:,:) = calculate_monthly_values(w, "avg",0,
False)
+ nmoStrt = nmoStrt+12
+ nmoLast = nmoLast+12
+ delete(w) ; may change size next iteration (leap yr)
+ end do
(0) nf=0 wind.1998.red.nc
*Segmentation fault (core dumped)*

*
[3] Output of 1st data file *

Variable: f (file variable)

filename: wind.1998.red
path: wind.1998.red.nc
   file global attributes:
      CDI : Climate Data Interface version 1.4.3
      Conventions : CF-1.0
      history : Tue Jul 05 10:48:34 2011: cdo sellevel,1000,10
wind.1998.reduzido.nc wind.1998.red.nc
Wed Jun 29 13:32:04 2011: cdo sellonlatbox,25,45,-28,-10 wind.1998.nc
wind.1998.reduzido.nc
Mon Oct 18 11:24:43 2010: cdo mergetime wind1998a.nc wind1998b.nc
wind.1998.nc
      CDO : Climate Data Operators version 1.4.3 (
http://www.mpimet.mpg.de/cdo)
   dimensions:
      lon = 8
      lat = 8
      lev = 2
      time = 1460 // unlimited
   variables:
      float lon ( lon )
         long_name : Longitude
         units : degrees_east
         standard_name : longitude

      float lat ( lat )
         long_name : Latitude
         units : degrees_north
         standard_name : latitude

      double lev ( lev )
         long_name : pressure
         units : Pa
         axis : z

      double time ( time )
         units : hours since 1800-01-01 00:00:00
         calendar : standard

      float WSpeed ( time, lev, lat, lon )
         long_name : Wind Speed
         standard_name : eastward_wind
         units : m/s
         _FillValue : -32767
         valid_range : ( -140, 174.989990234375 )
         parent_stat : Other
         statistic : Individual Obs
         level_desc : Pressure Levels
         dataset : NCEP/DOE AMIP-II Reanalysis (Reanalysis-2)
         var_desc : u-wind
         GRIB_name : UGRD
         GRIB_id : 33
         least_significant_digit : 1
         precision : 2
         actual_range : ( -87.30000305175781, 129.5 )
         unpacked_valid_range : ( -140, 175 )

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jul 5 04:27:29 2011

This archive was generated by hypermail 2.1.8 : Tue Jul 12 2011 - 15:03:53 MDT