Segmentation fault (core dumped)

From: Vishali P. <vishali.prat_at_nyahnyahspammersnyahnyah>
Date: Mon Jun 20 2011 - 02:40:40 MDT

Dear NCL users

I'm using 4x Daily NCEP/DOE AMIP-II Reanalysis wind data 1998 to 2004. I've
calculated wind speed for each year...and created outputs.
I was suggested a script to calculate unweighted seasonal means. And then,
average of each season from 1998 to 2004, ie, average if season 1 from 1998
to 2004,average of season 2 from 1998 to 2004, then 3rd and 4th.

But when I try to run,appears error "Segmentation fault (core dumped)".
Can anyone help me ?
Current version of my NCL is 5.2.1.

Thank you in advance

Vishali

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

   yrStrt = 1998
   yrLast = 2004
   nyrs = yrLast-yrStrt+1
   klev = 17
   nlat = 73
   mlon = 144
   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/Speed/"
   filu = systemfunc("cd "+diru+" ; ls wind*nc")
   nfilu = dimsizes(filu)

                    ; calculate monthly means from 6hrly
   nmoStrt = 0
   nmoLast = 11
   do nf=0,nfilu-1
      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)

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 20 02:40:51 2011

This archive was generated by hypermail 2.1.8 : Mon Jun 20 2011 - 12:30:20 MDT