Hi,
I was trying to calculate temporal power spectra at a number of longitude
points. When using the 'specx_anal'
function in a do-loop, it causes a Segmentation fault, i.e, could not do the
following:
Rspec = new( (/ntimh,nlon/), "float")
do i = 0, nlon-1
RtD = dtrend(Rt(:,i),False) ; detrend data
spec = specx_anal(RtD,iopt,jave,pct) ; calculate spectrum
Rspec(:,i) = spec@spcx
end do
Where, Rt is Rt(ntim,nlon) and ntimh = ntim/2. (Using specx_anal only once
works fine). I've tried all sorts of things to avoid the above problem, but
with no results (the problem still persists if I delete the last line:
Rspec(:,i) = spec@spcx). Does anyone have any suggestion how to get around
this problem? Thanks.
-- Harun=================================================================== Dr. M. Harun Ar Rashid Bureau of Meteorology Research Centre Email: h.rashid@bom.gov.au PO Box 1289K Tel : +61 3 9669 4550 Melbourne, Victoria 3001 Fax : +61 3 9669 4660 Australia http://www.bom.gov.au/bmrc ===================================================================
_______________________________________________ ncl-talk mailing list ncl-talk@ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Wed May 04 2005 - 08:07:51 MDT