Segmentation fault with dtrend_msg

From: Barry Baker <bbaker_at_nyahnyahspammersnyahnyah>
Date: Wed, 11 Oct 2006 14:45:04 -0600

Hello,

I'm new to NCL and having a problem with dtrend_msg. When I set return_info
to True I get a segmentation fault. I can't figure out why and was
wondering if anyone else may have run into this problem? I'm running NCL on
Powerbook G4 with Mac OS X 10.4.7 I've included my script is below which was
borrowed from Michael Notaro's web page. Thanks in advance, Barry

;******************************
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/contributed.ncl"

begin
diri = "/Users/bbaker/data/Vulnerability/"
fili = "co_tmax1895-1997.nc"

a = addfile(diri+fili, "r")
yrs = fspan(1.,103.,103)
 lat = a->lat ; get lat
 lon = a->lon ; get lon
 ta1 = a->tmax(:,:,:)/100.0 ; get data
 print(yrs)
 ta1!0 = "time"
 ta1!1 = "lat"
 ta1!2 = "lon"
 ta1&time=yrs
 ta1&lat=lat
 ta1&lon=lon

 ta2 = ta1(lat|:,lon|:,time|:)
 times = ispan(1,103,1)
 yDtrend = dtrend_msg(times,ta2,False,True) ; Calculate the trend

end
;*********************************
_________________________
Barry Baker, Ph.D.
Global Climate Change Initiative
The Nature Conservancy
2424 Spruce St., Suite 100
Boulder, CO 80302

Tel: (303)-541-0322
Fax: (303)-449-4328

http://nature.org/tncscience/scientists/misc/baker.html

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Oct 11 2006 - 14:45:04 MDT

This archive was generated by hypermail 2.2.0 : Thu Oct 12 2006 - 10:22:33 MDT